Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs

Motivated by the numerical treatment of parametric and stochastic PDEs, we analyze the least-squares method for polynomial approximation of multivariate functions based on random sampling according to a given probability measure. Recent work has shown that in the univariate case, the least-squares method is quasi-optimal in expectation in [A. Cohen, M A. Davenport and D. Leviatan. Found. Comput. Math. 13 (2013) 819–834] and in probability in [G. Migliorati, F. Nobile, E. von Schwerin, R. Tempone, Found. Comput. Math. 14 (2014) 419–456], under suitable conditions that relate the number of samples with respect to the dimension of the polynomial space. Here “quasi-optimal” means that the accuracy of the least-squares approximation is comparable with that of the best approximation in the given polynomial space. In this paper, we discuss the quasi-optimality of the polynomial least-squares method in arbitrary dimension. Our analysis applies to any arbitrary multivariate polynomial space (including tensor product, total degree or hyperbolic crosses), under the minimal requirement that its associated index set is downward closed. The optimality criterion only involves the relation between the number of samples and the dimension of the polynomial space, independently of the anisotropic shape and of the number of variables. We extend our results to the approximation of Hilbert space-valued functions in order to apply them to the approximation of parametric and stochastic elliptic PDEs. As a particular case, we discuss “inclusion type” elliptic PDE models, and derive an exponential convergence estimate for the least-squares method. Numerical results confirm our estimate, yet pointing out a gap between the condition necessary to achieve optimality in the theory, and the condition that in practice yields the optimal convergence rate.

[1]  Fabio Nobile,et al.  Analysis of Discrete $$L^2$$L2 Projection on Polynomial Spaces with Random Evaluations , 2014, Found. Comput. Math..

[2]  P. Revesz.  Interpolation and Approximation , 2010 .

[3]  Albert Cohen,et al.  Breaking the curse of dimensionality in sparse polynomial approximation of parametric PDEs , 2015 .

[4]  Simone Vantini,et al.  Analysis of AneuRisk65 data: $k$-mean alignment , 2014 .

[5]  Raul Tempone,et al.  MATHICSE Technical Report : Analysis of the discrete $L^2$ projection on polynomial spaces with random evaluations , 2011 .

[6]  Albert Cohen,et al.  Sparse adaptive Taylor approximation algorithms for parametric and stochastic elliptic PDEs , 2011 .

[7]  T. J. Rivlin,et al.  The growth of polynomials bounded at equally spaced points , 1992 .

[8]  Humberto Contreras,et al.  The stochastic finite-element method , 1980 .

[9]  Fabio Nobile,et al.  Computers and Mathematics with Applications Convergence of Quasi-optimal Stochastic Galerkin Methods for a Class of Pdes with Random Coefficients , 2022 .

[10]  R. Tempone,et al.  ON THE OPTIMAL POLYNOMIAL APPROXIMATION OF STOCHASTIC PDES BY GALERKIN AND COLLOCATION METHODS , 2012 .

[11]  C. D. Boor,et al.  Computational aspects of polynomial interpolation in several variables , 1992 .

[12]  Raul Tempone,et al.  An anisotropic sparse grid stochastic collocation method for elliptic partial differential equations with random input data , 2007 .

[13]  Albert Cohen,et al.  Convergence Rates of Best N-term Galerkin Approximations for a Class of Elliptic sPDEs , 2010, Found. Comput. Math..

[14]  Victor Nistor,et al.  High order Galerkin appoximations for parametric second order elliptic partial differential equations , 2012 .

[15]  Fabio Nobile,et al.  A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[16]  BabuskaIvo,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007 .

[17]  Fabio Nobile,et al.  An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[18]  Claude Jeffrey Gittelson,et al.  An adaptive stochastic Galerkin method , 2011 .

[19]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[20]  Albert Cohen,et al.  Breaking the curse of dimensionality in sparse polynomial approximation of parametric PDEs , 2013 .

[21]  Albert Cohen,et al.  High-Dimensional Adaptive Sparse Polynomial Interpolation and Applications to Parametric PDEs , 2013, Foundations of Computational Mathematics.

[22]  George G. Lorentz,et al.  Solvability problems of bivariate interpolation I , 1986 .

[23]  Albert Cohen,et al.  On the Stability and Accuracy of Least Squares Approximations , 2011, Foundations of Computational Mathematics.

[24]  Claude Jeffrey Gittelson,et al.  An adaptive stochastic Galerkin method for random elliptic operators , 2013, Math. Comput..

[25]  Raúl Tempone,et al.  Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations , 2004, SIAM J. Numer. Anal..

[26]  Giovanni Migliorati,et al.  Polynomial approximation by means of the random discrete L2 projection and application to inverse problems for PDEs with stochastic data , 2013 .

[27]  Fabio Nobile,et al.  Approximation of Quantities of Interest in Stochastic PDEs by the Random Discrete L2 Projection on Polynomial Spaces , 2013, SIAM J. Sci. Comput..

[28]  G. Migliorati,et al.  Multivariate Markov-type and Nikolskii-type inequalities for polynomials associated with downward closed multi-index sets , 2015, J. Approx. Theory.

[29]  M. Newman,et al.  Interpolation and approximation , 1965 .

[30]  Victor Nistor,et al.  HIGH-ORDER GALERKIN APPROXIMATIONS FOR PARAMETRIC SECOND-ORDER ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS , 2013 .

[31]  R. DeVore,et al.  ANALYTIC REGULARITY AND POLYNOMIAL APPROXIMATION OF PARAMETRIC AND STOCHASTIC ELLIPTIC PDE'S , 2011 .

[32]  Steen Krenk,et al.  Stochastic Finite Element Methods , 2004 .

[33]  R. Tempone,et al.  A sparse grid stochastic collocation method for elliptic partial differential equations with random input data , 2006 .