Conjugate gradient-type methods for linear systems with complex symmetric coefficient matrices. (English) Zbl 0761.65018

The work employs a generalization of the non-symmetric Lanczos method based on the minimization of the residual norm. Quasiminimization implies that the minimum residual norm is sought not over the entire Krylov subspace but in some specially determined unidimensional subspace.
In the considered algorithm similarly to the standard Lanczos method the normal process of computation can break down due to the vanishing of some scalar product. Suggestions about the continuation of the iterative process in case of breakdown are given. The paper contains a considerable number of numerical results that show the advantages of the suggested method.


65F10 Iterative numerical methods for linear systems
65F50 Computational methods for sparse matrices
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs


Full Text: DOI Link