Daniel Kressner
ETH Zurich
Low rank tensor methods for solving high-dimensional linear systems
Linear systems arising from the discretization of high-dimensional problems
are typically too large to be handled by standard algorithms. On the other
hand, high-dimensional problems often feature a rather simple geometry and
favorable smoothness properties. In such cases, it can be expected that the
solution of the discretized system admits low rank tensor approximations. This
talk presents numerical algorithms that operate with such low rank tensor ap-
proximations and ideally scale only linearly (instead of exponentially) as the
number of dimensions grows. It is shown how a combination of Krylov subspace
and low rank techniques can be used to derive such algorithms. A number of
examples, including parametrized and stochastic partial diff erential equations,
illustrate the use of the newly developed algorithms in applications. This is
joint work with Christine Tobler, ETH Zurich.