Abstract:The solution of large linear systems with multiple right-hand sides given simultaneously is required in many large-scale scientific and engineering applications modelled by either partial differential or boundary integral equations. Block Krylov subspace methods are attractive to use for this problem class as they can overcome the memory bottleneck of direct methods and they perform block matrix-vector products, achieving high computational efficiency on modern cache-based computer architectures. In this paper we introduce variants of the block GMRES method that combine initial deflation and eigenvalue recycling strategies to remedy some typical convergence problems of block Krylov solvers. The new class of block iterative solvers has the ability to handle the approximate linear dependence of the block of right-hand sides and exploits approximate invariant subspaces recycled over the iterations to mitigate the bad effects that small eigenvalues can have on the convergence, by adapting an existing preconditioner. We illustrate the numerical behavior of the spectrally updated and initially deflated block GMRES method on a set of linear systems arising from the discretization of the Dirac equation and of boundary integral equations in electromagnetics scattering.

Block Compressed Sparse Row of Multiple Dof and Large Sparse Linear System Real Domain Solution

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

A Scalable HPC-Based Domain Decomposition Method for Multiphysics Modeling of RF Devices

An Efficient Domain Decomposition Method in Multi-physics Simulation of Package Systems

Block Variants of the Cocg and Cocr Methods for Solving Complex Symmetric Linear Systems with Multiple Right-Hand Sides

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

A Method for Solving Sparse Linear Equations of Power Systems Based on GPU

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

Scalable Linear Solvers for Sparse Linear Systems from Large-Scale Numerical Simulations.

Fast Finite Element Electrostatic Analysis with Domain Decomposition Method

Limited Memory Block Krylov Subspace Optimization for Computing Dominant Singular Value Decompositions.

A NEW ALGORITHM OF LARGE-SCALE QUADRATIC PROGRAMMING SOLUTION FOR FINITE ELEMENT METHOD

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

Combining the CORS and BiCORSTAB Iterative Methods with MLFMA and SAI Preconditioning for Solving Large Linear Systems in Electromagnetics

Relaxed Incomplete Cholesky Factorization Preconditioned Symmetric Lbcg Algorithm for Solving Electromagnetic Problems

A compression scheme for domain decomposition method in solving electromagnetic problems

A spectrally preconditioned and initially deflated variant of the restarted block GMRES method for solving multiple right-hand sides linear systems

Implementation of a Parallel Sparse Direct Solver on Vector Architecture

Nonlinear Asynchronous Block Iterative Method

An Improved ICCG Method for Large Scale Sparse Linear Equations

A parallel sparse linear system solver for large-scale circuit simulation based on Schur Complement