Local Fourier Analysis of p-Multigrid for High-Order Finite Element Operators

DOI10.1137/21M1431199MaRDI QIDQ6108158zbMATH OpenFDO

Authors Jed Brown, Yunhui He

Publication date 29 June 2023

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2108.01751

zbMATH Keywords

finite elementshigh-orderlocal Fourier analysisp-multigrid

Mathematics Subject Classification ID

Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Numerical methods in Fourier analysis (65Txx)

Abstract: Multigrid methods are popular for solving linear systems derived from discretizing PDEs. Local Fourier Analysis (LFA) is a technique for investigating and tuning multigrid methods. P-multigrid is popular for high-order or spectral finite element methods, especially on unstructured meshes. In this paper, we introduce LFAToolkit.jl, a new Julia package for LFA of high-order finite element methods. LFAToolkit.jl analyzes preconditioning techniques for arbitrary systems of second order PDEs and supports mixed finite element methods. Specifically, we develop LFA of p-multigrid with arbitrary second-order PDEs using high-order finite element discretizations and examine the performance of Jacobi and Chebyshev smoothing for two-grid schemes with aggressive p-coarsening. A natural extension of this LFA framework is the analysis of h-multigrid for finite element discretizations or finite difference discretizations that can be represented in the language of finite elements. With this extension, we can replicate previous work on the LFA of h-multigrid for arbitrary order discretizations using a convenient and extensible abstraction. Examples in one, two, and three dimensions are presented to validate our LFA of p-multigrid for the Laplacian and linear elasticity.

Recommendations

• Two‐level Fourier analysis of multigrid for higher‐order finite‐element discretizations of the Laplacian
• Local Fourier analysis for multigrid with overlapping smoothers applied to systems of PDEs.
• On local Fourier analysis of multigrid methods for PDEs with jumping and random coefficients
• Multicolor Fourier analysis of the multigrid method for quadratic FEM discretizations
• Local Fourier Analysis of Multigrid for Hybridized and Embedded Discontinuous Galerkin Methods

Cites work

• A Multigrid Tutorial, Second Edition
• Automated local Fourier analysis (aLFA)
• Efficient nonlinear solvers for nodal high-order finite elements in 3D
• High-Order Methods for Incompressible Fluid Flow
• scientific article; zbMATH DE number 2130683 (Why is no real title available?)
• scientific article; zbMATH DE number 3896238 (Why is no real title available?)
• scientific article; zbMATH DE number 1561761 (Why is no real title available?)
• scientific article; zbMATH DE number 3797119 (Why is no real title available?)
• Jacobian-free Newton-Krylov methods: a survey of approaches and applications.
• Local Fourier analysis for multigrid with overlapping smoothers applied to systems of PDEs.
• Local Fourier analysis of balancing domain decomposition by constraints algorithms
• Local Fourier analysis of multigrid methods with polynomial smoothers and aggressive coarsening
• Multi-Level Adaptive Solutions to Boundary-Value Problems
• On local Fourier analysis of multigrid methods for PDEs with jumping and random coefficients
• On the validity of the local Fourier analysis
• Parallel multigrid smoothing: Polynomial versus Gauss--Seidel
• Spectral element multigrid. I: Formulation and numerical results
• The Chebyshev iteration revisited
• Toward a universal h-p adaptive finite element strategy. I: Constrained approximation and data structure
• Toward a universal h-p adaptive finite element strategy. II: A posteriori error estimation
• Toward a universal h-p adaptive finite element strategy. III: Design of h-p meshes
• Tuning Multigrid Methods with Robust Optimization and Local Fourier Analysis
• Two‐level Fourier analysis of multigrid for higher‐order finite‐element discretizations of the Laplacian

Cited in

(4)

• An algebraic generalization of local Fourier analysis for grid transfer operators in multigrid based on Toeplitz matrices
• Two‐level Fourier analysis of multigrid for higher‐order finite‐element discretizations of the Laplacian
• Optimized sparse matrix operations for reverse mode automatic differentiation
• Smoothers with localized residual computations for geometric multigrid methods for higher-order finite elements

This page was built for publication: Local Fourier Analysis of p-Multigrid for High-Order Finite Element Operators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6108158)