We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.
MiS Preprint
22/2010
A Projection Method to Solve Linear Systems in Tensor Format
Jonas Ballani and Lars Grasedyck
Abstract
In this paper we propose a method for the numerical solution of linear systems of equations in low rank tensor format. Such systems may arise from the discretisation of PDEs in high dimensions but our method is not limited to this type of application. We present an iterative scheme which is based on the projection of the residual to a low dimensional subspace. The subspace is spanned by vectors in low rank tensor format which --- similarly to Krylov subspace methods --- stem from the subsequent (approximate) application of the given matrix to the residual. All calculations are performed in hierarchical Tucker format which allows for applications in high dimensions. The mode size dependency is treated by a multilevel method. We present numerical examples that include high-dimensional convection-diffusion equations and shift-invert eigenvalue solvers.
