Study of Aggressive Coarsening and Multipass Interpolation in Algebraic Multigrid

Hans De Sterck

Department of Applied Mathematics, University of Waterloo, Ontario, Canada

Ulrike Meier Yang, Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, California, USA

Abstract

Algebraic multigrid is a very efficent algorithm for solving large linear systems on unstructured grids. Use of coarsening schemes such as parallel variants of the standard coarsening algorithm by Ruge and Stueben [1] or CLJP coarsening [2], a method based on parallel maximal independent set algorithms, can lead to high complexities with regard to memory usage as well as computation time, which adversely affect scalability.
In recent work [3], we have proposed two new parallel AMG coarsening schemes, that are based on solely enforcing a maximum independent set property, resulting in sparser coarse grids. The new coarsening techniques remedy memory and execution time complexity growth for various large three-dimensional (3D) problems. If used within AMG as a preconditioner for Krylov subspace methods, the resulting iterative methods tend to converge fast. For some difficult problems, however, these methods still produce complexities that are too high, or don't converge well enough, and further remedies in terms of coarsening and interpolation need to be found in order to obtain scalable methods.
In this paper we describe the combination of these various coarsening methods with "aggressive coarsening" techniques, which require long range "multipass" interpolation [4], and empirically study complexity and convergence properties of the resulting iterative methods. The resulting AMG methods, implemented in the hypre solver library [5], are applied to first-order system least-squares (FOSLS) discretizations of elliptic PDE systems. Parallel scalability of the combined FOSLS-AMG method is investigated for large-scale three-dimensional applications.

*This work was performed under the auspices of the U.S. Department of Energy by University of California Lawrence Livermore National Laboratory under contract number W-7405-Eng-48 and subcontract number B545391.

[1] J. Ruge, K. Stueben, Algebraic multigrid (AMG), in: S. McCormick, ed., Multigrid Methods, vol. 3 of Frontires in Applied Mathematics (SIAM, 1987) 73-130.
[2] V. E. Henson, U. M. Yang, BoomerAMG: a Parallel Algebraic Multigrid Solver and Preconditioner, Applied Numerical Mathematics, 41 (2002) 155-177.
[3] H. De Sterck, U. M. Yang, and J. J. Heys, Reducing Complexity in Parallel Algebraic Multigrid Preconditioners, submitted to SIAM Journal on Matrix Analysis and Applications, 2004.
[4] K. Stueben, Algebraic multigrid (AMG): an introduction with applications, in: U. Trottenberg, C. Osterlee and A. Schueller, eds., Multigrid (Academic Press, 2000) 413-532.
[5] R. Falgout, U. M. Yang, hypre: a Library of High Performance Preconditioners, in Computational Science - ICCS 2002 Part III, P. Sloot, C. Tan, J. Dongarra and A. Hoekstra, eds., Lecture Notes in Computer Science, 2331 (Springer, 2002) 632-641.