An Odyssey into Local Refinement and Multilevel Preconditioning III: Implementation and Numerical Experiments

DOI10.1137/S1064827502407676MaRDI QIDQ4442480zbMATH OpenFDO

Authors Burak Aksoylu, Stephen D. Bond, Michael Holst

Publication date 20 January 2004

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

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

zbMATH Keywords

2D and 3DalgorithmsBPXcomputational complexitydata structuresfinite element methodsiterative methodslocal mesh refinementmultilevel preconditioningnumerical experimentsred and red-green refinementstabilizations

Mathematics Subject Classification ID

Boundary value problems for second-order elliptic equations (35J25) Iterative numerical methods for linear systems (65F10) Numerical computation of matrix norms, conditioning, scaling (65F35) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50) Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55)

Abstract: In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on BPX-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis method (HB), and their use in the local mesh refinement setting. In this article, we describe in detail the implementation of these types of algorithms, including detailed discussions of the datastructures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard datatypes available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. Our implementations are performed in both C and MATLAB using the Finite Element ToolKit (FETK), an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement.

Recommendations

• Optimality of Multilevel Preconditioners for Local Mesh Refinement in Three Dimensions
• Algebraic multilevel preconditioning methods. I
• scientific article; zbMATH DE number 1167599
• Performance Enhancements and Parallel Algorithms for Two Multilevel Preconditioners
• Multilevel finite element preconditioning for \(\sqrt{3}\) refinement

Cited in

(9)

• Goal-oriented adaptivity and multilevel preconditioning for the Poisson-Boltzmann equation
• Arbitrary-level hanging nodes and automatic adaptivity in the hp-FEM
• Robust multigrid preconditioners for the high-contrast biharmonic plate equation
• A finite element method with mesh-separation-based approximation technique and its application in modeling crack propagation with adaptive mesh refinement
• Multilevel finite element preconditioning for \(\sqrt{3}\) refinement
• Optimality of Multilevel Preconditioners for Local Mesh Refinement in Three Dimensions
• Multiscale modeling of beam and plates using customized second-generation wavelets
• Robust multigrid preconditioners for cell-centered finite volume discretization of the high-contrast diffusion equation
• Optimal multilevel methods for graded bisection grids

Describes a project that uses

Uses Software

• FEtk
• Matlab

This page was built for publication: An Odyssey into Local Refinement and Multilevel Preconditioning III: Implementation and Numerical Experiments

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