Efficient approximation of high-dimensional exponentials by tensornetworks

05/19/2021
by   Martin Eigel, et al.
0

In this work a general approach to compute a compressed representation of the exponential exp(h) of a high-dimensional function h is presented. Such exponential functions play an important role in several problems in Uncertainty Quantification, e.g. the approximation of log-normal random fields or the evaluation of Bayesian posterior measures. Usually, these high-dimensional objects are intractable numerically and can only be accessed pointwise in sampling methods. In contrast, the proposed method constructs a functional representation of the exponential by exploiting its nature as a solution of an ordinary differential equation. The application of a Petrov–Galerkin scheme to this equation provides a tensor train representation of the solution for which we derive an efficient and reliable a posteriori error estimator. Numerical experiments with a log-normal random field and a Bayesian likelihood illustrate the performance of the approach in comparison to other recent low-rank representations for the respective applications. Although the present work considers only a specific differential equation, the presented method can be applied in a more general setting. We show that the composition of a generic holonomic function and a high-dimensional function corresponds to a differential equation that can be used in our method. Moreover, the differential equation can be modified to adapt the norm in the a posteriori error estimates to the problem at hand.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/19/2018

A New Bound on Hrushovski's Algorithm for Computing the Galois Group of a Linear Differential Equation

The complexity of computing the Galois group of a linear differential eq...
research
01/21/2022

Robust a posteriori estimates for the stochastic Cahn-Hilliard equation

We derive a posteriori error estimates for a fully discrete finite eleme...
research
07/03/2018

Improving the approximation of the first and second order statistics of the response process to the random Legendre differential equation

In this paper, we deal with uncertainty quantification for the random Le...
research
10/29/2019

Invariant Galerkin Ansatz Spaces and Davison-Maki Methods for the Numerical Solution of Differential Riccati Equations

The differential Riccati equation appears in different fields of applied...
research
03/08/2023

Determining the Rolle function in Hermite interpolatory approximation by solving an appropriate differential equation

We determine the pointwise error in Hermite interpolation by numerically...
research
12/06/2020

Estimating Vector Fields from Noisy Time Series

While there has been a surge of recent interest in learning differential...
research
06/08/2022

A Lanczos-type procedure for tensors

The solution of linear non-autonomous ordinary differential equation sys...

Please sign up or login with your details

Forgot password? Click here to reset