z-logo
Premium
An iterated tangential filtering decomposition
Author(s) -
Achdou Y.,
Nataf F.
Publication year - 2003
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.326
Subject(s) - preconditioner , mathematics , discretization , iterated function , factorization , matrix (chemical analysis) , incomplete lu factorization , multigrid method , matrix decomposition , diagonal , dimension (graph theory) , linear system , algorithm , partial differential equation , eigenvalues and eigenvectors , mathematical analysis , pure mathematics , geometry , materials science , physics , quantum mechanics , composite material
Large linear systems arising from the discretization of partial differential equations with finite differences or finite elements on structured grids in dimension d ( d ⩾ 3) require efficient preconditioners. For a symmetric and positive definite (SPD) matrix, we propose a SPD block LDL T preconditioner whose factorized form requires a smaller amount of memory than the original matrix. Moreover, the computing time for the preconditioner solves is linear with respect to the number of unknowns. The preconditioner is built in d stages: in a first stage, we use the tangential filtering decomposition of Wittum et al . and obtain a preconditioner which remains rather difficult to factorize. Then, in a second stage, we apply tangential filtering decomposition recursively to the diagonal blocks of this first preconditioner. The final stage consists of factorizing exactly the blocks corresponding to one dimensional problems. Such preconditioners can also be computed adaptively and combined in a multiplicative way. A generic programming implementation is discussed and numerical tests are presented, in particular for problems with highly heterogeneous media. Copyright © 2003 John Wiley & Sons, Ltd.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here