The Surrogate Matrix Methodology: A Priori Error Estimation
From MaRDI portal
(Redirected from The Surrogate Matrix Methodology: A Priori Error Estimation)
a priori analysisfinite element methodsgeometric multigridlow ordermatrix-freesurrogate numerical methods
Numerical interpolation (65D05) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Parallel numerical computation (65Y05) Complexity and performance of numerical algorithms (65Y20)
Abstract: We give the first mathematically rigorous analysis of an emerging approach to finite element analysis (see, e.g., Bauer et al. [Appl. Numer. Math., 2017]), which we hereby refer to as the surrogate matrix methodology. This methodology is based on the piece-wise smooth approximation of the matrices involved in a standard finite element discretization. In particular, it relies on the projection of smooth so-called stencil functions onto high-order polynomial subspaces. The performance advantage of the surrogate matrix methodology is seen in constructions where each stencil function uniquely determines the values of a significant collection of matrix entries. Such constructions are shown to be widely achievable through the use of locally-structured meshes. Therefore, this methodology can be applied to a wide variety of physically meaningful problems, including nonlinear problems and problems with curvilinear geometries. Rigorous a priori error analysis certifies the convergence of a novel surrogate method for the variable coefficient Poisson equation. The flexibility of the methodology is also demonstrated through the construction of novel methods for linear elasticity and nonlinear diffusion problems. In numerous numerical experiments, we demonstrate the efficacy of these new methods in a matrix-free environment with geometric multigrid solvers. In our experiments, up to a twenty-fold decrease in computation time is witnessed over the classical method with an otherwise identical implementation.
Recommendations
- An error estimate for matrix equations
- Estimating the error in matrix function approximations
- Matrix strategies for computing the least trimmed squares estimation of the general linear and SUR models
- \textit{A priori} error estimation for the stochastic perturbation method
- Surrogate-based methods
- scientific article; zbMATH DE number 468800
- A posteriori error estimates of Krylov subspace approximations to matrix functions
- Surrogate methods for linear inequalities
- Precision matrix estimation using penalized generalized Sylvester matrix equation
- The least-squares method for matrices that depend on a parameter
Cites work
- A fully asynchronous multifrontal solver using distributed dynamic scheduling
- A generic interface for parallel cell-based finite element operator application
- A scalable multi-level preconditioner for matrix-free µ-finite element analysis of human bone structures
- A scalable, matrix-free multigrid preconditioner for finite element discretizations of heterogeneous Stokes flow
- A stencil scaling approach for accelerating matrix-free finite element implementations
- A two-scale approach for efficient on-the-fly operator assembly in massively parallel high performance multigrid codes
- Coupled variational formulations of linear elasticity and the DPG methodology
- Curved Elements in the Finite Element Method. I
- Efficient nonlinear solvers for nodal high-order finite elements in 3D
- Element-by-element linear and nonlinear solution schemes
- Finite Element Approximation of the p-Laplacian
- Finite Element Interpolation of Nonsmooth Functions Satisfying Boundary Conditions
- Hierarchical hybrid grids: data structures and core algorithms for multigrid
- scientific article; zbMATH DE number 4138299 (Why is no real title available?)
- scientific article; zbMATH DE number 192847 (Why is no real title available?)
- scientific article; zbMATH DE number 1219577 (Why is no real title available?)
- scientific article; zbMATH DE number 1036632 (Why is no real title available?)
- scientific article; zbMATH DE number 1069512 (Why is no real title available?)
- scientific article; zbMATH DE number 3416931 (Why is no real title available?)
- scientific article; zbMATH DE number 3081828 (Why is no real title available?)
- Parallel octree-based finite element method for large-scale earthquake ground motion simulation
- Stencil scaling for vector-valued PDEs on hybrid grids with applications to generalized Newtonian fluids
- Tetrahedral grid refinement
- The surrogate matrix methodology: low-cost assembly for isogeometric analysis
- Towards Textbook Efficiency for Parallel Multigrid
Cited in
(6)- The surrogate matrix methodology: accelerating isogeometric analysis of waves
- Fast and multiscale formation of isogeometric matrices of microstructured geometric models
- The surrogate matrix methodology: low-cost assembly for isogeometric analysis
- Automatic stabilization of finite-element simulations using neural networks and hierarchical matrices
- A Matrix-Free ILU Realization Based on Surrogates
- Mixed-Precision Computing for Scientific Discovery: Formats, Co-Design, and Responsible Approximation
This page was built for publication: The Surrogate Matrix Methodology: A Priori Error Estimation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5204013)
