High-Dimensional Adaptive Sparse Polynomial Interpolation and Applications to Parametric PDEs

Abdellah Chkifa (1) , Albert Cohen (2) , Christoph Schwab (3)

1 ORNL - Oak Ridge National Laboratory [Oak Ridge]
2 LJLL - Laboratoire Jacques-Louis Lions
3 ETH-SAM - Seminar for Applied Mathematics

We consider the problem of Lagrange polynomial interpolation in high or countably infinite dimension, motivated by the fast computation of solutions to partial differential equations (PDEs) depending on a possibly large number of parameters which result from the application of generalized polynomial chaos discretizations to random and stochastic PDEs. In such applications there is a substantial advantage in considering polynomial spaces that are sparse and anisotropic with respect to the different parametric variables. In an adaptive context, the polynomial space is enriched at different stages of the computation. In this paper, we study an interpolation technique in which the sample set is incremented as the polynomial dimension increases, leading therefore to a minimal amount of PDE solving. This construction is based on standard principle of tensorization of a one dimensional interpolation scheme and sparsifi-cation. We derive bounds on the Lebesgue constants for this interpolation process in terms of their univariate counterpart. For a class of model elliptic parametric PDE's, we have shown in [11] that certain polynomial approximations based on Taylor expansions converge in terms of the polynomial dimension with an algebraic rate that is robust with respect to the parametric dimension. We show that this rate is preserved when using our interpolation algorithm. We also propose a greedy algorithm for the adaptive selection of the polynomial spaces based on our interpolation scheme, and illustrate its performance both on scalar valued functions and on parametric elliptic PDE's.

Domaines: Analyse numérique [math.NA]

Abdellah Chkifa, Albert Cohen, Christoph Schwab. High-Dimensional Adaptive Sparse Polynomial Interpolation and Applications to Parametric PDEs. Foundations of Computational Mathematics, 2014, 14 (4), pp.601-633. ⟨10.1007/s10208-013-9154-z⟩. ⟨hal-01350602⟩