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

Abdellah Chkifa, Albert Cohen, Giovanni Migliorati, Fabio Nobile, Raul Tempone

Research output: Contribution to journalArticlepeer-review

65 Scopus citations


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.
Original languageEnglish (US)
Pages (from-to)815-837
Number of pages23
JournalESAIM: Mathematical Modelling and Numerical Analysis
Issue number3
StatePublished - Apr 8 2015


Dive into the research topics of 'Discrete least squares polynomial approximation with random evaluations − application to parametric and stochastic elliptic PDEs'. Together they form a unique fingerprint.

Cite this