Year: 2010
Numerical Mathematics: Theory, Methods and Applications, Vol. 3 (2010), Iss. 3 : pp. 276–294
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.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/nmtma.2010.33.2
Numerical Mathematics: Theory, Methods and Applications, Vol. 3 (2010), Iss. 3 : pp. 276–294
Published online: 2010-01
AMS Subject Headings:
Copyright: COPYRIGHT: © Global Science Press
Pages: 19
Keywords: Hybrid direct/iterative solver domain decomposition incomplete/partial factorization Schur approximation scalable preconditioner convection-diffusion large $3D$ problems parallel scientific computing High Performance Computing.