Geodesic
Adaptive reduction-based multigrid for nearly singular and highly disordered physical systems
Electronic transactions on numerical analysis, Tome 37 (2010), pp. 276-295.

Voir la notice de l'article provenant de la source Electronic Library of Mathematics

Summary: Classical multigrid solution of linear systems with matrices that have highly variable entries and are nearly singular is made difficult by the compounding difficulties introduced by these two model features. Efficient multigrid solution of nearly singular matrices is known to be possible, provided the so-called Brandt-McCormick (or eigenvector approximation) criterion is satisfied, which requires building interpolation to fit the near-null-space modes with high accuracy. When these modes are known, traditional multigrid approaches may be very effective. In this paper, we consider the case of matrices describing highly disordered systems, such as those that arise in lattice quantum chromodynamics (QCD), where the near-null modes cannot be easily expressed in closed form. We develop a variational adaptive reduction-based algebraic multigrid preconditioner for such systems and present a two-level convergence theory for the approach for Hermitian and positive-definite systems. The proposed method is applied to a two-dimensional model known as the Gauge Laplacian, a common test problem for development of solvers in quantum dynamics applications, showing promising numerical results. The proposed reduction-based setup uses compatible relaxation coarsening together with a sparse approximation to the so-called ideal interpolation operator to recursively construct the coarse spaces.
Classification : 65N55, 65N22, 65F10
Keywords: algebraic multigrid, lattice QCD, gauge Laplacian, iterative methods, adaptive multigrid 
@article{ETNA_2010__37__a8,
     author = {Brannick, J. and Frommer, A. and Kahl, K. and Maclachlan, S. and Zikatanov, L.},
     title = {Adaptive reduction-based multigrid for nearly singular and highly disordered physical systems},
     journal = {Electronic transactions on numerical analysis},
     pages = {276--295},
     publisher = {mathdoc},
     volume = {37},
     year = {2010},
     language = {en},
     url = {http://geodesic.mathdoc.fr/item/ETNA_2010__37__a8/}
}
TY  - JOUR
AU  - Brannick, J.
AU  - Frommer, A.
AU  - Kahl, K.
AU  - Maclachlan, S.
AU  - Zikatanov, L.
TI  - Adaptive reduction-based multigrid for nearly singular and highly disordered physical systems
JO  - Electronic transactions on numerical analysis
PY  - 2010
SP  - 276
EP  - 295
VL  - 37
PB  - mathdoc
UR  - http://geodesic.mathdoc.fr/item/ETNA_2010__37__a8/
LA  - en
ID  - ETNA_2010__37__a8
ER  - 
%0 Journal Article
%A Brannick, J.
%A Frommer, A.
%A Kahl, K.
%A Maclachlan, S.
%A Zikatanov, L.
%T Adaptive reduction-based multigrid for nearly singular and highly disordered physical systems
%J Electronic transactions on numerical analysis
%D 2010
%P 276-295
%V 37
%I mathdoc
%U http://geodesic.mathdoc.fr/item/ETNA_2010__37__a8/
%G en
%F ETNA_2010__37__a8
Brannick, J.; Frommer, A.; Kahl, K.; Maclachlan, S.; Zikatanov, L. Adaptive reduction-based multigrid for nearly singular and highly disordered physical systems. Electronic transactions on numerical analysis, Tome 37 (2010), pp. 276-295. http://geodesic.mathdoc.fr/item/ETNA_2010__37__a8/