Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree
<p>The Riesz maps of the 𝐿<sup>2</sup> de Rham complex frequently arise as subproblems in the construction of fast preconditioners for more complicated problems. In this work, we present multigrid solvers for high-order finite-element discretizations of these Riesz maps with the sa...
Main Authors: | , |
---|---|
Format: | Journal article |
Language: | English |
Published: |
Society for Industrial and Applied Mathematics
2024
|
_version_ | 1811140371276627968 |
---|---|
author | Brubeck, PD Farrell, PE |
author_facet | Brubeck, PD Farrell, PE |
author_sort | Brubeck, PD |
collection | OXFORD |
description | <p>The Riesz maps of the 𝐿<sup>2</sup> de Rham complex frequently arise as subproblems in the construction of fast preconditioners for more complicated problems. In this work, we present multigrid solvers for high-order finite-element discretizations of these Riesz maps with the same time and space complexity as sum-factorized operator application, i.e., with optimal complexity in polynomial degree in the context of Krylov methods. The key idea of our approach is to build new finite elements for each space in the de Rham complex with orthogonality properties in both the 𝐿<sup>2</sup>- and 𝐻(d)-inner products (d∈{grad,curl,div}) on the reference hexahedron. The resulting sparsity enables the fast solution of the patch problems arising in the Pavarino, Arnold–Falk–Winther, and Hiptmair space decompositions in the separable case. In the nonseparable case, the method can be applied to an auxiliary operator that is sparse by construction. With exact Cholesky factorizations of the sparse patch problems, the application complexity is optimal, but the setup costs and storage are not. We overcome this with the finer Hiptmair space decomposition and the use of incomplete Cholesky factorizations imposing the sparsity pattern arising from static condensation, which applies whether static condensation is used for the solver or not. This yields multigrid relaxations with time and space complexity that are both optimal in the polynomial degree.</p> |
first_indexed | 2024-03-07T08:13:00Z |
format | Journal article |
id | oxford-uuid:ae5bbfb4-38c4-4169-81bb-1d416ed87175 |
institution | University of Oxford |
language | English |
last_indexed | 2024-09-25T04:20:55Z |
publishDate | 2024 |
publisher | Society for Industrial and Applied Mathematics |
record_format | dspace |
spelling | oxford-uuid:ae5bbfb4-38c4-4169-81bb-1d416ed871752024-08-07T17:09:18ZMultigrid solvers for the de Rham complex with optimal complexity in polynomial degreeJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:ae5bbfb4-38c4-4169-81bb-1d416ed87175EnglishSymplectic ElementsSociety for Industrial and Applied Mathematics2024Brubeck, PDFarrell, PE<p>The Riesz maps of the 𝐿<sup>2</sup> de Rham complex frequently arise as subproblems in the construction of fast preconditioners for more complicated problems. In this work, we present multigrid solvers for high-order finite-element discretizations of these Riesz maps with the same time and space complexity as sum-factorized operator application, i.e., with optimal complexity in polynomial degree in the context of Krylov methods. The key idea of our approach is to build new finite elements for each space in the de Rham complex with orthogonality properties in both the 𝐿<sup>2</sup>- and 𝐻(d)-inner products (d∈{grad,curl,div}) on the reference hexahedron. The resulting sparsity enables the fast solution of the patch problems arising in the Pavarino, Arnold–Falk–Winther, and Hiptmair space decompositions in the separable case. In the nonseparable case, the method can be applied to an auxiliary operator that is sparse by construction. With exact Cholesky factorizations of the sparse patch problems, the application complexity is optimal, but the setup costs and storage are not. We overcome this with the finer Hiptmair space decomposition and the use of incomplete Cholesky factorizations imposing the sparsity pattern arising from static condensation, which applies whether static condensation is used for the solver or not. This yields multigrid relaxations with time and space complexity that are both optimal in the polynomial degree.</p> |
spellingShingle | Brubeck, PD Farrell, PE Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title | Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title_full | Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title_fullStr | Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title_full_unstemmed | Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title_short | Multigrid solvers for the de Rham complex with optimal complexity in polynomial degree |
title_sort | multigrid solvers for the de rham complex with optimal complexity in polynomial degree |
work_keys_str_mv | AT brubeckpd multigridsolversforthederhamcomplexwithoptimalcomplexityinpolynomialdegree AT farrellpe multigridsolversforthederhamcomplexwithoptimalcomplexityinpolynomialdegree |