Sparse approximations of the Schur complement for parallel algebraic hybrid solvers in 3D

L. Giraud, A. Haidar, Y. Saad

Research output: Contribution to journalArticlepeer-review

16 Scopus citations

Abstract

In this paper we study the computational performance of variants of an algebraic additive Schwarz preconditioner for the Schur complement for the solution of large sparse linear systems. In earlier works, the local Schur complements were computed exactly using a sparse direct solver. The robustness of the preconditioner comes at the price of this memory and time intensive computation that is the main bottleneck of the approach for tackling huge problems. In this work we investigate the use of sparse approximation of the dense local Schur complements. These approximations are computed using a partial incomplete LU factorization. Such a numerical calculation is the core of the multi-level incomplete factorization such as the one implemented in pARMS. The numerical and computing performance of the new numerical scheme is illustrated on a set of large 3D convection-diffusion problems; preliminary experiments on linear systems arising from structural mechanics are also reported.

Original languageEnglish (US)
Pages (from-to)276-294
Number of pages19
JournalNumerical Mathematics
Volume3
Issue number3
DOIs
StatePublished - Aug 2010

Keywords

  • Convection-diffusion
  • Domain decomposition
  • High performance computing
  • Hybrid direct/iterative solver
  • Incomplete/partial factorization
  • Large 3D problems
  • Parallel scientific computing
  • Scalable preconditioner
  • Schur approximation

Fingerprint Dive into the research topics of 'Sparse approximations of the Schur complement for parallel algebraic hybrid solvers in 3D'. Together they form a unique fingerprint.

Cite this