Quasi-kernel polynomials and convergence results for quasi-minimal residual iterations

Recently, Freund and Nachtigal have proposed a novel polynominal-based iteration, the quasi-minimal residual algorithm (QMR), for solving general nonsingular non-Hermitian linear systems. Motivated by the QMR method, we have introduced the general concept of quasi-kernel polynomials, and we have shown that the QMR algorithm is based on a particular instance of quasi-kernel polynomials. In this paper, we continue our study of quasi-kernel polynomials. In particular, we derive bounds for the norms of quasi-kernel polynomials. These results are then applied to obtain convergence theorems both for the QMR method and for a transpose-free variant of QMR, the TFQMR algorithm

The Bramble-Pasciak preconditioner for saddle point problems

The Bramble-Pasciak Conjugate Gradient method is a well known tool to solve linear systems in saddle point form. A drawback of this method in order to ensure applicability of Conjugate Gradients is the need for scaling the preconditioner which typically involves the solution of an eigenvalue problem. Here, we introduce a modified preconditioner and inner product which without scaling enable the use of a MINRES variant and can be used for the simplified Lanczos process. Furthermore, the modified preconditioner and inner product can be combined with the original Bramble-Pasciak setup to give new preconditioners and inner products. We undermine the new methods by showing numerical experiments for Stokes problems

QMR: A Quasi-Minimal Residual method for non-Hermitian linear systems

The biconjugate gradient (BCG) method is the natural generalization of the classical conjugate gradient algorithm for Hermitian positive definite matrices to general non-Hermitian linear systems. Unfortunately, the original BCG algorithm is susceptible to possible breakdowns and numerical instabilities. A novel BCG like approach is presented called the quasi-minimal residual (QMR) method, which overcomes the problems of BCG. An implementation of QMR based on a look-ahead version of the nonsymmetric Lanczos algorithm is proposed. It is shown how BCG iterates can be recovered stably from the QMR process. Some further properties of the QMR approach are given and an error bound is presented. Finally, numerical experiments are reported

Recent advances in Lanczos-based iterative methods for nonsymmetric linear systems

In recent years, there has been a true revival of the nonsymmetric Lanczos method. On the one hand, the possible breakdowns in the classical algorithm are now better understood, and so-called look-ahead variants of the Lanczos process have been developed, which remedy this problem. On the other hand, various new Lanczos-based iterative schemes for solving nonsymmetric linear systems have been proposed. This paper gives a survey of some of these recent developments

Image reconstruction for diffuse optical tomography using bi-conjugate gradient and transpose-free quasi minimal residual algorithms and comparison of them

Diffuse optical tomography (DOT) is a new emerging modality in the diagnosis of soft tissue abnormalities. DOT image quality substantially depends on the recon struction stage. In the literature, there are many reconstruction algorithms used in DOT systems. However, some algorithms were improved for solving specific cases but still need to be improved. The bi-conjugate gradient (BiCG) enhanced is one of the conjugate gradient (CG)-based reconstruction techniques for non-Hermitian systems. The BiCG provides a solution to a non-Hermitian system. However, it has erratic convergence in some cases. Therefore, DOT images reconstructed by BiCG can be at the wrong location and is inaccurate in some cases. In this study, we used continuous-wave diffuse optical tomography (CW-DOT) to acquire mea surements from breast tissue phantoms with single or double inclusion at different depths and center-to-center separations and we have used the transpose free quasi minimal residual (TFQMR) reconstruction algorithm, improved as an alternative to BiCG for the first time in the CW-DOT system. Moreover, we have experimen tally proved that TFQMR is superior to BiCG in some specific cases for the first time in CW-DOT. Therefore, we concluded that TFQMR has the potential to be able to be used in the reconstruction stage in CW-DOT

BiCGCR2: A new extension of conjugate residual method for solving non-Hermitian linear systems

In the present paper, we introduce a new extension of the conjugate residual (CR) for solving non-Hermitian linear systems with the aim of developing an alternative basic solver to the established biconjugate gradient (BiCG), biconjugate residual (BiCR) and biconjugate A-orthogonal residual (BiCOR) methods. The proposed Krylov subspace method, referred to as the BiCGCR2 method, is based on short-term vector recurrences and is mathematically equivalent to both BiCR and BiCOR. We demonstrate by extensive numerical experiments that the proposed iterative solver has often better convergence performance than BiCG, BiCR and BiCOR. Hence, it may be exploited for the development of new variants of non-optimal Krylov subspace methods