Abstract:We consider the problem of approximating the solution to $A(\mu) x(\mu) = b$ for many different values of the parameter $\mu$. Here we assume $A(\mu)$ is large, sparse, and nonsingular with a nonlinear dependence on $\mu$. Our method is based on a companion linearization derived from an accurate Chebyshev interpolation of $A(\mu)$ on the interval $[-a,a]$, $a \in \mathbb{R}$. 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(\mu)$ for many different values of the parameter $\mu \in [-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 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.

Restrictively Preconditioned Chebyshev Method for Solving Systems of Linear Equations

Restrictive Preconditioners for Conjugate Gradient Methods for Symmetric Positive Definite Linear Systems

Restrictive Preconditioning for Convection-Diffusion Distributed Control Problems

Chebyshev-Type Methods and Preconditioning Techniques

Practical Restrictively Preconditioned Conjugate Gradient Methods for a Class of Block Two-by-Two Linear Systems

Preconditioned Gauss-Seidel Type Iterative Method for Solving Linear Systems

On a Chebyshev Accelerated Splitting Iteration Method with Application to Two‐by‐two Block Linear Systems

Constraint Preconditioners for Symmetric Indefinite Matrices

On Hybrid Preconditioning Methods for Large Sparse Saddle-Point Problems

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

The Upper Jacobi and Upper Gauss–Seidel Type Iterative Methods for Preconditioned Linear Systems

A Preconditioned MINRES Method for Block Lower Triangular Toeplitz Systems

BCCB preconditioners for solving linear systems from delay differential equations

A Preconditioned Iteration Method for Solving Sylvester Equations.

A CONSTRAINT PRECONDITIONER FOR SOLVING SYMMETRIC POSITIVE DEFINITE SYSTEMS AND APPLICATION TO THE HELMHOLTZ EQUATIONS AND POISSON EQUATIONS

Generalized Tridiagonal Preconditioners for Solving Linear Systems

Preconditioned Chebyshev BiCG for parameterized linear systems

RESTRICTIVELY PRECONDITIONED CONJUGATE GRADIENT METHOD FOR A SERIES OF CONSTANTLY AUGMENTED LEAST SQUARES PROBLEMS

A preconditioned iteration method for solving saddle point problems

Preconditioned iterative method for nonsymmetric saddle point linear systems