Tensor-Sparsity of Solutions to High-Dimensional Elliptic Partial Differential Equations

A recurring theme in attempts to break the curse of dimensionality in the numerical approximation of solutions to high-dimensional partial differential equations (PDEs) is to employ some form of sparse tensor approximation. Unfortunately, there are only a few results that quantify the possible advantages of such an approach. This paper introduces a class of functions, which can be written as a sum of rank-one tensors using a total of at most parameters, and then uses this notion of sparsity to prove a regularity theorem for certain high-dimensional elliptic PDEs. It is shown, among other results, that whenever the right-hand side of the elliptic PDE can be approximated with a certain rate in the norm of by elements of , then the solution can be approximated in from to accuracy for any . Since these results require knowledge of the eigenbasis of the elliptic operator considered, we propose a second "basis-free" model of tensor-sparsity and prove a regularity theorem for this second sparsity model as well. We then proceed to address the important question of the extent to which such regularity theorems translate into results on computational complexity. It is shown how this second model can be used to derive computational algorithms with performance that breaks the curse of dimensionality on certain model high-dimensional elliptic PDEs with tensor-sparse data.