Guaranteed a posteriori error bounds for low rank tensor approximate solutions

Abstract

We propose a guaranteed and fully computable upper bound on the energy norm of the error in low-rank Tensor Train (TT) approximate solutions of (possibly) high dimensional reaction-diffusion problems. The error bound is obtained from Euler-Lagrange equations for a complementary flux reconstruction problem, which are solved in the low-rank TT representation using the block Alternating Linear Scheme. This bound is guaranteed to be above the energy norm of the total error, including the discretization error, the tensor approximation error, and the error in the solver of linear algebraic equations, although quadrature errors, in general, can pollute its evaluation. Numerical examples with the Poisson equation and the Schroedinger equation with the Henon-Heiles potential in up to 40 dimensions are presented to illustrate the efficiency of this approach.

References in corpus (7)

• Alternating minimal energy methods for linear systems in higher dimensions. Part I: SPD systems
• Alternating minimal energy methods for linear systems in higher dimensions. Part II: Faster algorithm and application to nonsymmetric systems
• Computation of extreme eigenvalues in higher dimensions using block tensor train format
• Quasioptimality of maximum-volume cross interpolation of tensors
• Robust error bounds for finite element approximation of reaction-diffusion problems with non-constant reaction coefficient in arbitrary space dimension
• Simple and robust equilibrated flux a posteriori estimates for singularly perturbed reaction-diffusion problems
• A Simple Approach to Reliable and Robust A Posteriori Error Estimation for Singularly Perturbed Problems