Computation of the Response Surface in the Tensor Train data format

Abstract

We apply the Tensor Train (TT) approximation to construct the Polynomial Chaos Expansion (PCE) of a random field, and solve the stochastic elliptic diffusion PDE with the stochastic Galerkin discretization. We compare two strategies of the polynomial chaos expansion: sparse and full polynomial (multi-index) sets. In the full set, the polynomial orders are chosen independently in each variable, which provides higher flexibility and accuracy. However, the total amount of degrees of freedom grows exponentially with the number of stochastic coordinates. To cope with this curse of dimensionality, the data is kept compressed in the TT decomposition, a recurrent low-rank factorization. PCE computations on sparse grids sets are extensively studied, but the TT representation for PCE is a novel approach that is investigated in this paper. We outline how to deduce the PCE from the covariance matrix, assemble the Galerkin operator, and evaluate some post-processing (mean, variance, Sobol indices), staying within the low-rank framework. The most demanding are two stages. First, we interpolate PCE coefficients in the TT format using a few number of samples, which is performed via the block cross approximation method. Second, we solve the discretized equation (large linear system) via the alternating minimal energy algorithm. In the numerical experiments we demonstrate that the full expansion set encapsulated in the TT format is indeed preferable in cases when high accuracy and high polynomial orders are required.

28 pages. Submitted to SIAM J. of Uncertainty Quantification

References in corpus (1)

• Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation

Cited by in corpus (10)

• Spectral tensor-train decomposition
• Enabling High-Dimensional Hierarchical Uncertainty Quantification by ANOVA and Tensor-Train Decomposition
• To be or not to be intrusive? The solution of parametric and stochastic equations --- Proper Generalized Decomposition
• Tensor Numerical Methods for High-dimensional PDEs: Basic Theory and Initial Applications
• Prospects of tensor-based numerical modeling of the collective electrostatic potential in many-particle systems
• Adaptive non-intrusive reconstruction of solutions to high-dimensional parametric PDEs
• Partial inversion of the elliptic operator to speed up computation of likelihood in Bayesian inference
• Tucker tensor method for fast grid-based summation of long-range potentials on 3D lattices with defects
• Quantification of airfoil geometry-induced aerodynamic uncertainties - comparison of approaches
• Tensor Approximation of Advanced Metrics for Sensitivity Analysis