Preconditioned Chebyshev BiCG for parameterized linear systems

12/08/2022
by   Siobhán Correnty, et al.
0

The biconjugate gradient method (BiCG) is one of the most popular short-term recurrence methods for solving non-symmetric linear systems of equations. The objective of this paper is to provide an efficient adaption of BiCG to parameterized linear systems. More precisely, we consider the problem of approximating the solution to A(μ) x(μ) = b for many different values of the parameter μ. Here we assume A(μ) is large, sparse, and nonsingular with a nonlinear dependence on μ. Our method is based on a companion linearization derived from an accurate Chebyshev interpolation of A(μ) on the interval [-a,a], a ∈ℝ. The solution to the linearization is approximated in a preconditioned BiCG setting for shifted systems, where the Krylov basis matrix is formed once. This process leads to a short-term recurrence method, where one execution of the algorithm produces the approximation to x(μ) for many different values of the parameter μ∈ [-a,a] simultaneously. In particular, this work proposes one algorithm which applies a shift-and-invert preconditioner exactly as well as an algorithm which applies the same preconditioner inexactly. The competitiveness of the algorithms are illustrated with large-scale problems arising from a finite element discretization of a Helmholtz equation with parameterized material coefficient. The software used in the simulations is publicly available online, and thus all our experiments are reproducible.

READ FULL TEXT

page 14

page 15

research
06/10/2022

Preconditioned infinite GMRES for parameterized linear systems

We are interested in obtaining approximate solutions to parameterized li...
research
02/08/2021

Infinite GMRES for parameterized linear systems

We consider linear parameter-dependent systems A(μ) x(μ) = b for many di...
research
09/24/2021

Compound Krylov subspace methods for parametric linear systems

In this work, we propose a reduced basis method for efficient solution o...
research
10/30/2020

Numerical solution of the wave propagation problem in a plate

In this work, the propagation of an ultrasonic pulse in a thin plate is ...
research
06/15/2020

Alternating Energy Minimization Methods for Multi-term Matrix Equations

We develop computational methods for approximating the solution of a lin...
research
08/01/2023

Anderson Accelerated PMHSS for Complex-Symmetric Linear Systems

This paper presents the design and development of an Anderson Accelerate...
research
04/21/2023

A Deep Learning algorithm to accelerate Algebraic Multigrid methods in Finite Element solvers of 3D elliptic PDEs

Algebraic multigrid (AMG) methods are among the most efficient solvers f...

Please sign up or login with your details

Forgot password? Click here to reset