A stencil scaling approach for accelerating matrix-free finite element implementations

DOI10.1137/17M1148384MaRDI QIDQ4560707zbMATH OpenOpenAlexWikidataFDO

Authors

Publication date 7 December 2018

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

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

zbMATH Keywords

matrix-free finite elementsoptimal order a priori estimatesstencil scalingvariable coefficientsvariational crime analysis

Mathematics Subject Classification ID

A priori estimates in context of PDEs (35B45) PDEs in connection with fluid mechanics (35Q35) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Error bounds for boundary value problems involving PDEs (65N15) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Complexity and performance of numerical algorithms (65Y20) Flows in porous media; filtration; seepage (76S05)

Abstract: We present a novel approach to fast on-the-fly low order finite element assembly for scalar elliptic partial differential equations of Darcy type with variable coefficients optimized for matrix-free implementations. Our approach introduces a new operator that is obtained by appropriately scaling the reference stiffness matrix from the constant coefficient case. Assuming sufficient regularity, an a priori analysis shows that solutions obtained by this approach are unique and have asymptotically optimal order convergence in the $H^1$- and the $L^2$-norm on hierarchical hybrid grids. For the pre-asymptotic regime, we present a local modification that guarantees uniform ellipticity of the operator. Cost considerations show that our novel approach requires roughly one third of the floating-point operations compared to a classical finite element assembly scheme employing nodal integration. Our theoretical considerations are illustrated by numerical tests that confirm the expectations with respect to accuracy and run-time. A large scale application with more than a hundred billion ($1.6cdot10^{11}$) degrees of freedom executed on 14,310 compute cores demonstrates the efficiency of the new scaling approach.

Recommendations

• Stencil scaling for vector-valued PDEs on hybrid grids with applications to generalized Newtonian fluids
• A two-scale approach for efficient on-the-fly operator assembly in massively parallel high performance multigrid codes
• Optimizing the Evaluation of Finite Element Matrices
• High-performance evaluation of finite element variational forms via commuting diagrams and duality
• Fast matrix-free evaluation of discontinuous Galerkin finite element operators

Cites work

• 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 two-scale approach for efficient on-the-fly operator assembly in massively parallel high performance multigrid codes
• Cubature of integrands containing derivatives
• Efficient nonlinear solvers for nodal high-order finite elements in 3D
• Element-by-element linear and nonlinear solution schemes
• Hierarchical hybrid grids: data structures and core algorithms for multigrid
• scientific article; zbMATH DE number 410743 (Why is no real title available?)
• scientific article; zbMATH DE number 5602968 (Why is no real title available?)
• scientific article; zbMATH DE number 3596197 (Why is no real title available?)
• scientific article; zbMATH DE number 1036632 (Why is no real title available?)
• scientific article; zbMATH DE number 762608 (Why is no real title available?)
• scientific article; zbMATH DE number 770271 (Why is no real title available?)
• Maximum principle and uniform convergence for the finite element method
• Optimal a priori estimates for higher order finite elements for elliptic interface problems
• Parallel octree-based finite element method for large-scale earthquake ground motion simulation
• Tetrahedral grid refinement
• Towards Textbook Efficiency for Parallel Multigrid
• Transfinite element methods: Blending-function interpolation over arbitrary curved element domains
• Weakened acute type condition for tetrahedral triangulations and the discrete maximum principle
• Weighted Poincare inequalities

Cited in

(9)

• Double-grid quadrature with interpolation-projection (DoGIP) as a novel discretisation approach: an application to FEM on simplexes
• Matrix-free multigrid solvers for phase-field fracture problems
• GPU-based matrix-free finite element solver exploiting symmetry of elemental matrices
• A two-scale approach for efficient on-the-fly operator assembly in massively parallel high performance multigrid codes
• Textbook efficiency: massively parallel matrix-free multigrid for the Stokes system
• A Massively Parallel Eulerian-Lagrangian Method for Advection-Dominated Transport in Viscous Fluids
• Stencil scaling for vector-valued PDEs on hybrid grids with applications to generalized Newtonian fluids
• The Surrogate Matrix Methodology: A Priori Error Estimation
• A Matrix-Free ILU Realization Based on Surrogates

Describes a project that uses

Uses Software

• likwid

This page was built for publication: A stencil scaling approach for accelerating matrix-free finite element implementations

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