Parallel Auxiliary Space AMG Solver for $H(div)$ Problems
Abstract
In this paper we present a family of scalable preconditioners for matrices arising in the discretization of $H(div)$ problems using the lowest order Raviart--Thomas finite elements. Our approach belongs to the class of “auxiliary space''--based methods and requires only the finite element stiffness matrix plus some minimal additional discretization information about the topology and orientation of mesh entities. We provide a detailed algebraic description of the theory, parallel implementation, and different variants of this parallel auxiliary space divergence solver (ADS) and discuss its relations to the Hiptmair--Xu (HX) auxiliary space decomposition of $H(div)$ [SIAM J. Numer. Anal., 45 (2007), pp. 2483--2509] and to the auxiliary space Maxwell solver AMS [J. Comput. Math., 27 (2009), pp. 604--623]. An extensive set of numerical experiments demonstrates the robustness and scalability of our implementation on large-scale $H(div)$ problems with large jumps in the material coefficients.
Keywords
MSC codes
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
R. A. Adams and J. F. Fournier, Sobolev Spaces, Pure Appl. Math. 140, Academic Press, Boston, 2003.
2.
C. Amrouche, C. Bernardi, M. Dauge, and V. Girault, Vector potentials in three-dimensional nonsmooth domains, Math. Methods Appl. Sci., 21 (1998), pp. 823--864.
3.
D. N. Arnold, R. S. Falk, and J. Gopalakrishnan, Mixed Finite Element Approximation of the Vector Laplacian with Dirichlet Boundary Conditions, preprint, arXiv:1109.3668, 2011.
4.
D. N. Arnold, R. S. Falk, and R. Winther, Preconditioning in H(div) and applications, Math. Comp., 66 (1998), pp. 957--984.
5.
D. N. Arnold, R. S. Falk, and R. Winther, Multigrid in H(div) and H(curl), Numer. Math., 85 (2000), pp. 197--217.
6.
A. H. Baker, R. D. Falgout, Tz. V. Kolev, and U. M. Yang, Multigrid smoothers for ultraparallel computing, SIAM J. Sci. Comput., 33 (2011), pp. 2864--2887.
7.
N. Bell and L. N. Olson, Algebraic multigrid for k-form Laplacians, Numer. Linear Algebra Appl., 15 (2008), pp. 165--185.
8.
P. Bochev, C. Siefert, R. Tuminaro, J. Xu, and Y. Zhu, Compatible Gauge Approaches for H(div) Equations, Tech. Report SAND 2007-5384P, Sandia National Laboratories, Albuquerque, NM, 2007.
9.
A. Bossavit, Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements, Academic Press, San Diego, 1998.
10.
J. H. Bramble, A proof of the inf-sup condition for the Stokes equations on Lipschitz domains, Math. Model. Methods Appl. Sci., 13 (2003), pp. 361--371.
11.
F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer Ser. in Comput. Math. 15, Springer-Verlag, New York, 1991.
12.
T. A. Brunner and T. V. Kolev, Algebraic multigrid for linear systems obtained by explicit element reduction, SIAM J. Sci. Comput., 33 (2011), pp. 2706--2731.
13.
Z. Cai, R. Lazarov, T. A. Manteuffel, and S. F. McCormick, First-order system least squares for second-order partial differential equations: part I, SIAM J. Numer. Anal., 31 (1994), pp. 1785--1799.
14.
Z. Cai, C. Tong, P. S. Vassilevski, and C. Wang, Mixed finite element methods for incompressible flow: Stationary Stokes equations, Numer. Methods Partial Differential Equations, 26 (2010), pp. 957--978.
15.
M. Christie and M. Blunt, Tenth SPE comparative solution project: A comparison of upscaling techniques, SPE Reservoir Engrg. Eval., 4 (2001), pp. 308--317.
16.
P. Clément, Approximation by finite element functions using local regularization, Rev. Française Automat. Informat. Recherche Opérationnelle Sér. RAIRO Anal. Numér., 9 (1975), pp. 77--84.
17.
L. A. Dávalos and D. Zanette, Fundamentals of Electromagnetism: Vacuum Electrodynamics, Media and Relativity, Springer-Verlag, New York, 1999.
18.
R. D. Falgout and P. S. Vassilevski, On generalizing the algebraic multigrid framework, SIAM J. Numer. Anal., 42 (2004), pp. 1669--1693.
19.
M. Gee, C. Siefert, J. Hu, R. Tuminaro, and M. Sala, ML 5.0 Smoothed Aggregation User's Guide, Tech. Report SAND2006-2649, Sandia National Laboratories, Albuquerque, NM, 2006.
20.
N. A. Gentile, Implicit Monte Carlo diffusion---an acceleration method for Monte Carlo time-dependent radiative transfer simulations, J. Comput. Phys., 172 (2001), pp. 543--571.
21.
F. Graziani and J. LeBlanc, The Crooked Pipe Test Problem, Tech. Report UCRL-MI-143393, Lawrence Livermore National Laboratory, Livermore, CA, 2000.
22.
V. E. Henson and U. M. Yang, BoomerAMG: A parallel algebraic multigrid solver and preconditioner, Appl. Numer. Math., 41 (2002), pp. 155--177.
23.
R. Hiptmair, Multigrid method for H(div) in three dimensions, Electron. Trans. Numer. Anal., 6 (1997), pp. 133-152.
24.
R. Hiptmair, Finite elements in computational electromagnetism, Acta Numer., 11 (2002), pp. 237--339.
25.
R. Hiptmair and J. Xu, Nodal auxiliary space preconditioning in H(curl) and H(div) spaces, SIAM J. Numer. Anal., 45 (2007), pp. 2483--2509.
26.
hypre: High Performance Preconditioners. \newblock http://www.llnl.gov/CASC/hypre/, 2011.
27.
T. V. Kolev and P. S. Vassilevski, Parallel auxiliary space AMG for H(curl) problems, J. Comput. Math., 27 (2009), pp. 604--623.
28.
T. V. Kolev and P. S. Vassilevski, Regular Decompositions for H(div) Spaces, Comput. Methods Appl. Math., 12 (2012), pp. 437--447.
29.
I. Lashuk and P. S. Vassilevski, Element agglomeration coarse Raviart-Thomas spaces with improved approximation properties, Numer. Linear Algebra Appl., 19 (2012), \newblock pp. 414--426.
30.
P. Lin, A sequential regularization method for time-dependent incompressible Navier--Stokes equations, SIAM J. Numer. Anal., 34 (1997), pp. 1051--1071.
31.
P.-H. Maire, A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes, J. Comput. Phys., 228 (2009), pp. 2391--2425.
32.
K.-A. Mardal, J. Schöberl, and R. Winther, A uniformly stable Fortin operator for the Taylor-Hood element, Numer. Math.
33.
K.-A. Mardal and R. Winther, Preconditioning discretizations of systems of partial differential equations, Numer. Linear Algebra Appl., 18 (2011), pp. 1--40.
34.
MFEM: Modular Finite Element Methods Library. \newblock http://mfem.googlecode.com.
35.
P. Monk, Finite Element Methods for Maxwell's Equations, Numer. Math. Sci. Comput., Oxford University Press, New York, 2003.
36.
J. E. Morel, T.-Y. Brian Yang, and J. S. Warsa, Linear multifrequency-grey acceleration recast for preconditioned krylov iterations, J. Comput. Phys., 227 (2007), pp. 244--263.
37.
J.-C. Nédélec, Mixed finite elements in ${\mathbb R}^3$, Numer. Math., 35 (1980), pp. 315--341.
38.
M. A. Olshanskii, J. Peters, and A. Reusken, Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface equations, Numer. Math., 105 (2006), pp. 159--191.
39.
P. Oswald, Multilevel Finite Element Approximation: Theory Applications, Teubner Skr. Numer., Teubner, Stuttgart, 1994.
40.
J. E. Pasciak and P. S. Vassilevski, Exact de Rham sequences of spaces defined on macro-elements in two and three spatial dimensions, SIAM J. Sci. Comput., 30 (2008), pp. 2427--2446.
41.
P. Raviart and J. Thomas, Primal hybrid finite element method for 2nd order elliptic problems, Math. Comp., 31 (1977), pp. 391--413.
42.
SPE Comparative Solution Project, \newblock http://www.spe.org/csp (2000).
43.
H. D. Sterck, R. D. Falgout, J. W. Nolting, and U. M. Yang, Distance-two interpolation for parallel algebraic multigrid, Numer. Linear Algebra Appl., 15 (2008), pp. 115--139.
44.
H. D. Sterck, U. M. Yang, and J. J. Heys, Reducing complexity in parallel algebraic multigrid preconditioners, SIAM J. Matrix Anal. Appl., 27 (2006), pp. 1019--1039.
45.
F. Tampieri, Newell's method for computing the plane equation of a polygon, Academic Press Professional, San Diego, 1992, pp. 231--232.
46.
A. T. Till, T. A. Brunner, and T. S. Bailey, A higher-order face finite element radiation diffusion method for unstructured curvilinear meshes, LLNL-PRES-494711, \newblock Talk presented at Lawrence Livermore National Laboratory, Livermore, CA, 2011.
47.
R. Tuminaro, J. Xu, and Y. Zhu, Auxiliary space preconditioners for mixed finite element methods, in Domain Decomposition Methods in Science and Engrg XVIII, Lecture Notes in Comput. Sci. Engrg. 70, Springer, Berlin, 2009, pp. 99--109.
48.
P. S. Vassilevski, Multilevel Block Factorization Preconditioners: Matrix-based Analysis and Algorithms for Solving Finite Element Equations, Springer, New York, 2008.
49.
P. S. Vassilevski and U. Villa, A Block-Diagonal Algebraic Multigrid Preconditioner for the Brinkman Problem, LLNL-CONF-522863, Lawrence Livermore National Laboratory, Livermore, CA, SIAM J. Sci. Comput., to appear.
50.
J. Xu, The auxiliary space method and optimal preconditioning techniques for unstructured grids, Computing, 56 (1996), pp. 215--235.
51.
J. Xu and Y. Zhu, Uniform convergent multigrid methods for elliptic problems with strongly discontinuous coefficients, Math. Models Methods Appl. Sci., 18 (2008), pp. 77--105.
Information & Authors
Information
Published In
SIAM Journal on Scientific Computing
Pages: A3079 - A3098
DOI: 10.1137/110859361
ISSN (online): 1095-7197
Copyright
© 2012, Society for Industrial and Applied Mathematics.
History
Submitted: 19 December 2011
Accepted: 15 October 2012
Published online: 18 December 2012
Keywords
MSC codes
Authors
Metrics & Citations
Metrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited By
- End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizationsThe International Journal of High Performance Computing Applications, Vol. 37, No. 5 | 1 June 2023
- Parameter-robust preconditioners for Biot’s modelSeMA Journal, Vol. 12 | 5 August 2023
- An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned MeshesSIAM Journal on Scientific Computing, Vol. 0, No. 0 | 19 July 2023
- Parallel Element-Based Algebraic Multigrid for \(\boldsymbol{H}\) (curl) and \(\boldsymbol{H}\) (div) Problems Using the ParELAG LibrarySIAM Journal on Scientific Computing, Vol. 45, No. 3 | 26 June 2023
- Low-Order Preconditioning for the High-Order Finite Element de Rham ComplexSIAM Journal on Scientific Computing, Vol. 45, No. 2 | 26 April 2023
- Fast solution of three‐dimensional elliptic equations with randomly generated jumping coefficients by using tensor‐structured preconditionersNumerical Linear Algebra with Applications, Vol. 151 | 17 November 2022
- Patch‐smoother and multigrid for the dual formulation for linear elasticityInternational Journal for Numerical Methods in Engineering, Vol. 122, No. 24 | 20 October 2021
- Block preconditioners for mixed-dimensional discretization of flow in fractured porous mediaComputational Geosciences, Vol. 25, No. 2 | 4 July 2020
- A New Class of AMG Interpolation Methods Based on Matrix-Matrix MultiplicationsSIAM Journal on Scientific Computing, Vol. 43, No. 5 | 19 July 2021
- Robust Preconditioners for Perturbed Saddle-Point Problems and Conservative Discretizations of Biot's Equations Utilizing Total PressureSIAM Journal on Scientific Computing, Vol. 43, No. 4 | 4 August 2021
- An observation on the uniform preconditioners for the mixed Darcy problemNumerical Methods for Partial Differential Equations, Vol. 36, No. 6 | 10 July 2020
- Robust Preconditioners for a New Stabilized Discretization of the Poroelastic EquationsSIAM Journal on Scientific Computing, Vol. 42, No. 3 | 4 June 2020
- Mixed-Dimensional Auxiliary Space PreconditionersSIAM Journal on Scientific Computing, Vol. 42, No. 5 | 22 October 2020
- A Review of Regular Decompositions of Vector Fields: Continuous, Discrete, and Structure-PreservingSpectral and High Order Methods for Partial Differential Equations ICOSAHOM 2018 | 12 August 2020
- Algebraic Hybridization and Static Condensation with Application to Scalable $H$(div) PreconditioningSIAM Journal on Scientific Computing, Vol. 41, No. 3 | 7 May 2019
- Modifying AMG Coarse Spaces with Weak Approximation Property to Exhibit Approximation in Energy NormSIAM Journal on Matrix Analysis and Applications, Vol. 40, No. 3 | 24 September 2019
- Multigrid methods for H(div) in three dimensions with nonoverlapping domain decomposition smoothersNumerical Linear Algebra with Applications, Vol. 25, No. 5 | 17 June 2018
- Lumping Techniques for Mixed Finite Element Diffusion DiscretizationsJournal of Computational and Theoretical Transport, Vol. 47, No. 4-6 | 8 May 2018
- A Scalable Preconditioner for a Primal Discontinuous Petrov--Galerkin MethodSIAM Journal on Scientific Computing, Vol. 40, No. 2 | 24 April 2018
- The Auxiliary Space Preconditioner for the de Rham ComplexSIAM Journal on Numerical Analysis, Vol. 56, No. 6 | 6 November 2018
- A Smoother Based on Nonoverlapping Domain Decomposition Methods for H(div) Problems: A Numerical StudyDomain Decomposition Methods in Science and Engineering XXIV | 5 January 2019
- A robust optimal preconditioner for the mixed finite element discretization of elliptic optimal control problemsNumerical Linear Algebra with Applications, Vol. 25, No. 1 | 10 November 2017
- Algebraic interface-based coarsening AMG preconditioner for multi-scale sparse matrices with applications to radiation hydrodynamics computationNumerical Linear Algebra with Applications, Vol. 24, No. 2 | 21 December 2016
- A Fluidity-Based First-Order System Least-Squares Method for Ice SheetsSIAM Journal on Scientific Computing, Vol. 39, No. 2 | 23 March 2017
- A Comparison of Finite Element Spaces for $H$(div) Conforming First-Order System Least SquaresSIAM Journal on Scientific Computing, Vol. 39, No. 5 | 26 October 2017
- A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random FieldsSIAM Journal on Scientific Computing, Vol. 39, No. 5 | 26 October 2017
- Space-Time CFOSLS Methods with AMGe UpscalingDomain Decomposition Methods in Science and Engineering XXIII | 18 March 2017
- Parallel Solver for $$\boldsymbol{H}$$ (div) Problems Using Hybridization and AMGDomain Decomposition Methods in Science and Engineering XXIII | 18 March 2017
- On the Robustness and Prospects of Adaptive BDDC Methods for Finite Element Discretizations of Elliptic PDEs with High-Contrast CoefficientsProceedings of the Platform for Advanced Scientific Computing Conference | 8 June 2016
- PCBDDC: A Class of Robust Dual-Primal Methods in PETScSIAM Journal on Scientific Computing, Vol. 38, No. 5 | 27 October 2016
- A BDDC Preconditioner for Problems Posed in H(div) with Deluxe ScalingDomain Decomposition Methods in Science and Engineering XXII | 1 Jan 2016
- A parallel Robin–Robin domain decomposition method for H(div)-elliptic problemsInternational Journal of Computer Mathematics, Vol. 92, No. 2 | 27 March 2014
- A Mixed Formulation for the Brinkman ProblemSIAM Journal on Numerical Analysis, Vol. 52, No. 1 | 4 February 2014
- Analysis of Block Preconditioners for Models of Coupled Magma/Mantle DynamicsSIAM Journal on Scientific Computing, Vol. 36, No. 4 | 19 August 2014
- An Overlapping Schwarz Algorithm for Raviart--Thomas Vector Fields with Discontinuous CoefficientsSIAM Journal on Numerical Analysis, Vol. 51, No. 1 | 31 January 2013
- A Block-Diagonal Algebraic Multigrid Preconditioner for the Brinkman ProblemSIAM Journal on Scientific Computing, Vol. 35, No. 5 | 28 October 2013