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.
翻译:在这项工作中,对高维函数的指数 $\ exp(h) 美元进行压缩表示的一般计算方法为高维函数 $h 美元。这种指数函数在不确定性量化的若干问题中起着重要作用,例如对正对随机字段的近似或对巴耶西亚后天测量的评估。通常,这些高维天体在数字上是难以操作的,只能通过抽样方法获得。相反,拟议方法通过利用其性质作为普通差分方程的解决方案来构建指数的功能代表。对这个方程采用Petrov-Galerkin 方案提供了一种高压阵列的表达方法,我们由此产生了一个高效和可靠的后天误标数。用日志随机字段的数值实验和Bayesian 的可能性说明了该方法相对于最近对各应用程序的低级表述的绩效。虽然目前的工作只考虑具体的差异方程式,但所提出的方法可以在更笼统的设置中应用。我们所介绍的方法可以显示,在普通方程式中,普通方程式的公式的构成可以与高方程式的公式函数相适应。