Research reports

Years: 2024 2023 2022 2021 2020 2019 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2004 2003 2002 2001 2000 1999 1998 1997 1996 1995 1994 1993 1992 1991

Higher order QMC Petrov-Galerkin discretization for affine parametric operator equations with random field inputs

by J. Dick and F.Y. Kuo and Q.T. Le Gia and D. Nuyens and Ch. Schwab

(Report number 2013-29)

Abstract
We construct quasi-Monte Carlo methods to approximate the expected values of linear functionals of Petrov-Galerkin discretizations of parametric operator equations which depend on a possibly infinite sequence of parameters. Such problems arise in the numerical solution of differential and integral equations with random field inputs. We analyze the regularity of the solutions with respect to the parameters in terms of the rate of decay of the fluctuations of the input field. If $p\in (0,1]$ denotes the ``summability exponent'' corresponding to the fluctuations in affine-parametric families of operators, then we prove that deterministic ``interlaced polynomial lattice rules'' of order $\alpha = \lfloor 1/p \rfloor+1$ in $s$ dimensions with $N$ points can be constructed using a fast component-by-component algorithm, in $\cal O(\alpha\,s\, N\log N + \alpha^2\,s^2 N)$ operations, to achieve a convergence rate of $\cal O(N^{-1/p})$, with the implied constant independent of $s$. This dimension-independent convergence rate is superior to the rate $\cal O(N^{-1/p+1/2})$ for $2/3\leq p\leq 1$, which was recently established for randomly shifted lattice rules under comparable assumptions. In our analysis we use a non-standard Banach space setting and introduce ``smoothness-driven product and order dependent (SPOD)'' weights for which we develop a new fast CBC construction.

Keywords: Quasi Monte-Carlo methods, interlaced polynomial lattice rules, higher order digital nets, parametric operator equations, infinite dimensional quadrature, Galerkin discretization

@Techreport{DKLNS13_526,
  author = {J. Dick and F.Y. Kuo and Q.T. Le Gia and D. Nuyens and Ch. Schwab},
  title = {Higher order QMC Petrov-Galerkin discretization for affine parametric operator equations with random field inputs},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2013-29},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2013/2013-29.pdf },
  year = {2013}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).