Abstract:In this work we study convergence properties of sparse polynomial approximations for a class of affine parametric saddle point problems. Such problems can be found in many computational science and engineering fields, including the Stokes equations for viscous incompressible flow, mixed formulation of diffusion equations for heat conduction or groundwater flow, time-harmonic Maxwell equations for electromagnetics, etc. Due to the lack of knowledge or intrinsic randomness, the coefficients of such problems are uncertain and can often be represented or approximated by high- or countably infinite-dimensional random parameters equipped with suitable probability distributions, and the coefficients affinely depend on a series of either globally or locally supported basis functions, e.g., Karhunen--Loève expansion, piecewise polynomials, or adaptive wavelet approximations. Consequently, we are faced with solving affine parametric saddle point problems. Here we study sparse polynomial approximations of the parametric solutions, in particular sparse Taylor approximations, and their convergence properties for these parametric problems. With suitable sparsity assumptions on the parametrization, we obtain the algebraic convergence rates $O(N^{-r})$ for the sparse polynomial approximations of the parametric solutions, in cases of both globally and locally supported basis functions. We prove that $r$ depends only on a sparsity parameter in the parametrization of the random input, and in particular does not depend on the number of active parameter dimensions or the number of polynomial terms $N$. These results imply that sparse polynomial approximations can effectively break the curse of dimensionality, thereby establishing a theoretical foundation for the development and application of such practical algorithms as adaptive, least-squares, and compressive sensing constructions.

Dynamical low-rank tensor approximations to high-dimensional parabolic problems: existence and convergence of spatial discretizations

On low-rank approximability of solutions to high-dimensional operator equations and eigenvalue problems

A space-time adaptive low-rank method for high-dimensional parabolic partial differential equations

Adaptive Low-Rank Methods: Problems on Sobolev Spaces

Preconditioned low-rank Riemannian optimization for linear systems with tensor product structure

A tensor approximation method based on ideal minimal residual formulations for the solution of high-dimensional problems

Collocation methods for nonlinear differential equations on low-rank manifolds

Solving differential Riccati equations: A nonlinear space-time method using tensor trains

Parametric PDEs: Sparse or Low-Rank Approximations?

A low-rank isogeometric solver based on Tucker tensors

Improved ParaDiag via low-rank updates and interpolation

On the Global Convergence of the Alternating Least Squares Method for Rank-One Approximation to Generic Tensors

Cross Interpolation for Solving High-Dimensional Dynamical Systems on Low-Rank Tucker and Tensor Train Manifolds

A literature survey of low-rank tensor approximation techniques

Sparse polynomial approximations for affine parametric saddle point problems

Low-rank plus diagonal approximations for Riccati-like matrix differential equations

Time integration of rank-constrained Tucker tensors

Interpolatory dynamical low-rank approximation for the 3+3d Boltzmann-BGK equation

Unconditional convergence for discretizations of dynamical optimal transport

Tensor-Sparsity of Solutions to High-Dimensional Elliptic Partial Differential Equations

Weighted sparsity and sparse tensor networks for least squares approximation