Block filtering decomposition

Riadh Fezzani 1 Laura Grigori 2, 3 Frédéric Nataf 2, 3 Ke Wang 4
2 ALPINES - Algorithms and parallel tools for integrated numerical simulations
LJLL - Laboratoire Jacques-Louis Lions, Inria Paris-Rocquencourt, Institut National des Sciences Mathématiques et de leurs Interactions
Abstract : This paper introduces a new preconditioning technique that is suitable for matrices arising from the discretization of a system of PDEs on unstructured grids. The preconditioner satisfies a so-called filtering property, which ensures that the input matrix is identical with the preconditioner on a given filtering vector. This vector is chosen to alleviate the effect of low-frequency modes on convergence and so decrease or eliminate the plateau that is often observed in the convergence of iterative methods. In particular, the paper presents a general approach that allows to ensure that the filtering condition is satisfied in a matrix decomposition. The input matrix can have an arbitrary sparse structure. Hence, it can be reordered using nested dissection, to allow a parallel computation of the preconditioner and of the iterative process. We show the efficiency of our preconditioner through a set of numerical experiments on symmetric and nonsymmetric matrices.
Type de document :
Article dans une revue
Numerical Linear Algebra with Applications, Wiley, 2014, 21 (6), pp.18. 〈10.1002/nla.1921〉
Liste complète des métadonnées
Contributeur : Laura Grigori <>
Soumis le : vendredi 5 décembre 2014 - 04:18:21
Dernière modification le : vendredi 25 mai 2018 - 12:02:06

Lien texte intégral




Riadh Fezzani, Laura Grigori, Frédéric Nataf, Ke Wang. Block filtering decomposition. Numerical Linear Algebra with Applications, Wiley, 2014, 21 (6), pp.18. 〈10.1002/nla.1921〉. 〈hal-01091252〉



Consultations de la notice