Row scaling as a preconditioner for some nonsymmetric linear systems with discontinuous coefficients

Dan Gordon, Rachel Gordon

Research output: Contribution to journalArticlepeer-review

Abstract

Linear systems with large differences between the coefficients, called "discontinuous coefficients", often arise when physical phenomena in heterogeneous media are modeled by partial differential equations (PDEs). Such problems are usually solved by domain decomposition techniques, but these can be difficult to implement when subdomain boundaries are complicated or the grid is unstructured. It is known that for such systems, diagonal scaling can sometimes improve the eigenvalue distribution and the convergence properties of some algorithm/preconditioner combinations. However, there seems to be no study outlining both the usefulness and limitations of this approach. It is shown that L2-scaling of the equations is a generally useful preconditioner for such problems when the system matrices are nonsymmetric, but only when the off-diagonal elements are small to moderate. Tests were carried out on several nonsymmetric linear systems with discontinuous coefficients derived from convectiondiffusion elliptic PDEs with small to moderate convection terms. It is shown that L2-scaling improved the eigenvalue distribution of the system matrix by reducing their concentration around the origin very significantly. Furthermore, such scaling improved the convergence properties of restarted GMRES and Bi-CGSTAB, with and without the ILU(0) preconditioner. Since ILU(0) is theoretically oblivious to diagonal scaling, these results indicate that L2-scaling also improves the runtime numerical stability.

Original languageEnglish
Pages (from-to)3480-3495
Number of pages16
JournalJournal of Computational and Applied Mathematics
Volume234
Issue number12
DOIs
StatePublished - 15 Oct 2010

Keywords

  • Bi-CGSTAB
  • Diagonal scaling
  • Discontinuous coefficients
  • Domain decomposition
  • Geometric scaling
  • GMRES
  • GRS
  • L-norm
  • Linear equations
  • Nonsymmetric equations
  • Partial differential equations

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Row scaling as a preconditioner for some nonsymmetric linear systems with discontinuous coefficients'. Together they form a unique fingerprint.

Cite this