Log In Sign Up

A parameter-dependent smoother for the multigrid method

by   Lars Grasedyck, et al.

The solution of parameter-dependent linear systems, by classical methods, leads to an arithmetic effort that grows exponentially in the number of parameters. This renders the multigrid method, which has a well understood convergence theory, infeasible. A parameter-dependent representation, e.g., a low-rank tensor format, can avoid this exponential dependence, but in these it is unknown how to calculate the inverse directly within the representation. The combination of these representations with the multigrid method requires a parameter-dependent version of the classical multigrid theory and a parameter-dependent representation of the linear system, the smoother, the prolongation and the restriction. A derived parameter-dependent version of the smoothing property, fulfilled by parameter-dependent versions of the Richardson and Jacobi methods, together with the approximation property prove the convergence of the multigrid method for arbitrary parameter-dependent representations. For a model problem low-rank tensor formats represent the parameter-dependent linear system, prolongation and restriction. The smoother, a damped Jacobi method, is directly approximated in the low-rank tensor format by using exponential sums. Proving the smoothing property for this approximation guarantees the convergence of the parameter-dependent method. Numerical experiments for the parameter-dependent model problem, with bounded parameter value range, indicate a grid size independent convergence rate.


page 1

page 2

page 3

page 4


Bayesian inversion for electromyography using low-rank tensor formats

The reconstruction of the structure of biological tissue using electromy...

Low-rank tensor methods for Markov chains with applications to tumor progression models

Continuous-time Markov chains describing interacting processes exhibit a...

Tensor train based isogeometric analysis for PDE approximation on parameter dependent geometries

This work develops a numerical solver based on the combination of isogeo...

Local convergence of alternating low-rank optimization methods with overrelaxation

The local convergence of alternating optimization methods with overrelax...

A Low-rank Method for Parameter-dependent Fluid-structure Interaction Discretizations With Hyperelasticity

In aerospace engineering and boat building, fluid-structure interaction ...

A Low-rank Approach for Nonlinear Parameter-dependent Fluid-structure Interaction Problems

Parameter-dependent discretizations of linear fluid-structure interactio...

Tensor-Krylov method for computing eigenvalues of parameter-dependent matrices

In this paper we extend the Residual Arnoldi method for calculating an e...