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.
DOI : 10.1051/m2an/2014050
Mots clés : Approximation theory, polynomial approximation, least squares, parametric and stochastic PDEs, high-dimensional approximation
@article{M2AN_2015__49_3_815_0, author = {Chkifa, Abdellah and Cohen, Albert and Migliorati, Giovanni and Nobile, Fabio and Tempone, Raul}, title = {Discrete least squares polynomial approximation with random evaluations \ensuremath{-} application to parametric and stochastic elliptic {PDEs}}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {815--837}, publisher = {EDP-Sciences}, volume = {49}, number = {3}, year = {2015}, doi = {10.1051/m2an/2014050}, zbl = {1318.41004}, mrnumber = {3342229}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2014050/} }
TY - JOUR AU - Chkifa, Abdellah AU - Cohen, Albert AU - Migliorati, Giovanni AU - Nobile, Fabio AU - Tempone, Raul TI - Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2015 SP - 815 EP - 837 VL - 49 IS - 3 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2014050/ DO - 10.1051/m2an/2014050 LA - en ID - M2AN_2015__49_3_815_0 ER -
%0 Journal Article %A Chkifa, Abdellah %A Cohen, Albert %A Migliorati, Giovanni %A Nobile, Fabio %A Tempone, Raul %T Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2015 %P 815-837 %V 49 %N 3 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2014050/ %R 10.1051/m2an/2014050 %G en %F M2AN_2015__49_3_815_0
Chkifa, Abdellah; Cohen, Albert; Migliorati, Giovanni; Nobile, Fabio; Tempone, Raul. Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 49 (2015) no. 3, pp. 815-837. doi : 10.1051/m2an/2014050. http://www.numdam.org/articles/10.1051/m2an/2014050/
A stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Numer. Anal. 45 (2007) 1005–1034 . | DOI | MR | Zbl
, and ,Galerkin finite element approximations of stochastic elliptic partial differential equations. SIAM J. Numer. Anal. 42 (2004) 800–825. | DOI | MR | Zbl
, and ,On the optimal polynomial approximation of stochastic PDEs by Galerkin and Collocation methods. Math. Model. Methods Appl. Sci. 22 (2012) 1250023. | DOI | MR | Zbl
, , and ,Convergence of quasi-optimal Stochastic Galerkin Methods for a class of PDES with random coefficients. Comput. Math. Appl. 67 (2014) 732–751. | DOI | MR | Zbl
, , , ,Sparse adaptive Taylor approximation algorithms for parametric and stochastic elliptic PDEs. Math. Model. Numer. Anal. 47 (2013) 253–280. | DOI | Numdam | MR | Zbl
, , and ,High-dimensional adaptive sparse polynomial interpolation and applications to parametric PDEs. Found. Comput. Math. 14 (2014) 601–633. | DOI | MR | Zbl
, and ,A. Chkifa, A. Cohen and C. Schwab, Breaking the curse of dimensionality in sparse polynomial approximation of parametric PDEs. To appear in J. Math. Pures Appl. (2015). | MR | Zbl
On the stability and accuracy of least squares approximations. Found. Comput. Math. 13 (2013) 819–834. | DOI | MR | Zbl
, and ,Convergence rates of best -term Galerkin approximations for a class of elliptic sPDEs, Found. Comput. Math. 10 (2010) 615–646. | DOI | MR | Zbl
, and ,Analytic regularity and polynomial approximation of parametric and stochastic PDE’s. Anal. Appl. 9 (2011) 1–37. | DOI | MR | Zbl
, and ,The growth of polynomials bounded at equally spaced points. SIAM J. Math. Anal. 23 (1992) 970–983. | DOI | MR | Zbl
and ,Ph.J. Davis, Interpolation and Approximation. Blaisdell Publishing Company (1963). | MR | Zbl
Computational aspects of polynomial interpolation in several variables. Math. Comput. 58 (1992) 705–727. | DOI | MR | Zbl
and ,An adaptive stochastic Galerkin method for random elliptic operators. Math. Comput. 82 (2013) 1515–1541. | DOI | MR | Zbl
,M. Kleiber and T. D. Hien, The stochastic finite element methods. John Wiley & Sons, Chichester (1992). | MR | Zbl
J. Kuntzman, Méthodes numériques - Interpolation, dérivées. Dunod, Paris (1959). | MR | Zbl
Solvability problems of bivariate interpolation I. Constructive Approximation 2 (1986) 153–169. | DOI | MR | Zbl
and ,G.Migliorati, Polynomial approximation by means of the random discrete projection and application to inverse problems for PDEs with stochastic data. Ph.D. thesis, Dipartimento di Matematica “Francesco Brioschi” Politecnico di Milano, Milano, Italy, and Centre de Mathématiques Appliquées, École Polytechnique, Palaiseau, France (2013).
Multivariate Markov-type and Nikolskii-type inequalities for polynomials associated with downward closed multi-index sets. J. Approx. Theory 189 (2015) 137–159. | DOI | MR | Zbl
,Analysis of discrete projection on polynomial spaces with random evaluations. Found. Comput. Math. 14 (2014) 419–456. | MR | Zbl
, , and ,Approximation of Quantities of Interest in stochastic PDEs by the random discrete projection on polynomial spaces. SIAM J. Sci. Comput. 35 (2013) A1440–A1460. | DOI | MR | Zbl
, , and ,High order Galerkin approximations for parametric, second order elliptic partial differential equations, Report 2012-22, Seminar for Applied Mathematics, ETH Zürich. Math. Models Methods Appl. Sci. 23 (2013). | MR | Zbl
and ,A sparse grid stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Num. Anal. 46 (2008) 2309–2345. | DOI | MR | Zbl
, and ,An anisotropic sparse grid stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Num. Anal. 46 (2008) 2411–2442. | DOI | MR | Zbl
, and ,Cité par Sources :