• Media type: Text; Report; E-Book
  • Title: Efficient approximation of high-dimensional exponentials by tensor networks
  • Contributor: Eigel, Martin [Author]; Farchmin, Nando [Author]; Heidenreich, Sebastian [Author]; Trunschke, Philipp [Author]
  • imprint: Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021
  • Issue: published Version
  • Language: English
  • DOI: https://doi.org/10.34657/8600; https://doi.org/10.20347/WIAS.PREPRINT.2844
  • ISSN: 2198-5855
  • Keywords: Uncertainty quantification ; dynamical system approximation ; a posteriori error bounds ; tensor train format ; Petrov--Galerkin ; tensor product methods ; holonomic functions ; log-normal random field ; Bayesian likelihoods
  • Origination:
  • Footnote: Diese Datenquelle enthält auch Bestandsnachweise, die nicht zu einem Volltext führen.
  • Description: 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.
  • Access State: Open Access