SPARSE HIERARCHICAL PRECONDITIONERS USING PIECEWISE SMOOTH APPROXIMATIONS OF EIGENVECTORS.

Saved in:
Bibliographic Details
Title: SPARSE HIERARCHICAL PRECONDITIONERS USING PIECEWISE SMOOTH APPROXIMATIONS OF EIGENVECTORS.
Authors: KLOCKIEWICZ, BAZYLI, DARVE, ERIC
Source: SIAM Journal on Scientific Computing; 2020, Vol. 42 Issue 6, pA3907-A3931, 25p
Subject Terms: SMOOTHNESS of functions, ELLIPTIC equations, FLOW simulations, LINEAR equations, RIGID bodies
Abstract: When solving linear systems arising from PDE discretizations, iterative methods (such as conjugate gradient (CG), GMRES, or MINRES) are often the only practical choice. To converge in a small number of iterations, however, they have to be coupled with an efficient preconditioner. One approach to preconditioning is provided by the hierarchical approximate factorization methods. However, to guarantee sufficient accuracy on the eigenvectors corresponding to the smallest eigenvalues, these methods typically have to be performed at very stringent accuracies, making the preconditioner expensive to apply. On the other hand, for a large class of problems, including many elliptic equations, the eigenvectors corresponding to small eigenvalues are smooth functions of the PDE grid. In this paper, we describe a hierarchical approximate factorization approach which focuses on improving accuracy on the smooth eigenvectors. The improved accuracy is achieved by preserving the action of the factorized matrix on piecewise polynomial functions of the grid. Based on the factorization, we propose a family of sparse preconditioners with Ϭ(n) or Ϭ (n log n) construction complexities. Our methods exhibit rapid convergence of CG in benchmarks run on large elliptic problems, arising for example in flow or mechanical simulations. In the case of the linear elasticity equation the preconditioners are exact on the near-kernel rigid body modes. [ABSTRACT FROM AUTHOR]
Copyright of SIAM Journal on Scientific Computing is the property of Society for Industrial & Applied Mathematics and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Complementary Index
Description
Abstract:When solving linear systems arising from PDE discretizations, iterative methods (such as conjugate gradient (CG), GMRES, or MINRES) are often the only practical choice. To converge in a small number of iterations, however, they have to be coupled with an efficient preconditioner. One approach to preconditioning is provided by the hierarchical approximate factorization methods. However, to guarantee sufficient accuracy on the eigenvectors corresponding to the smallest eigenvalues, these methods typically have to be performed at very stringent accuracies, making the preconditioner expensive to apply. On the other hand, for a large class of problems, including many elliptic equations, the eigenvectors corresponding to small eigenvalues are smooth functions of the PDE grid. In this paper, we describe a hierarchical approximate factorization approach which focuses on improving accuracy on the smooth eigenvectors. The improved accuracy is achieved by preserving the action of the factorized matrix on piecewise polynomial functions of the grid. Based on the factorization, we propose a family of sparse preconditioners with Ϭ(n) or Ϭ (n log n) construction complexities. Our methods exhibit rapid convergence of CG in benchmarks run on large elliptic problems, arising for example in flow or mechanical simulations. In the case of the linear elasticity equation the preconditioners are exact on the near-kernel rigid body modes. [ABSTRACT FROM AUTHOR]
ISSN:10648275
DOI:10.1137/20M1315683