Optimizing the Evaluation of Finite Element Matrices
From MaRDI portal
Abstract: Assembling stiffness matrices represents a significant cost in many finite element computations. We address the question of optimizing the evaluation of these matrices. By finding redundant computations, we are able to significantly reduce the cost of building local stiffness matrices for the Laplace operator and for the trilinear form for Navier-Stokes. For the Laplace operator in two space dimensions, we have developed a heuristic graph algorithm that searches for such redundancies and generates code for computing the local stiffness matrices. Up to cubics, we are able to build the stiffness matrix on any triangle in less than one multiply-add pair per entry. Up to sixth degree, we can do it in less than about two. Preliminary low-degree results for Poisson and Navier-Stokes operators in three dimensions are also promising.
Recommendations
- Optimizing finite element matrix calculations using the general technique of element vectorization
- Geometric Optimization of the Evaluation of Finite Element Matrices
- scientific article; zbMATH DE number 4133438
- Topological Optimization of the Evaluation of Finite Element Matrices
- Optimal matrix approximation in finite element model updating
- Combinatorial optimization of matrix-vector multiplication in finite element assembly
- scientific article; zbMATH DE number 4155773
- scientific article; zbMATH DE number 4074338
- scientific article; zbMATH DE number 495827
- An efficient way to assemble finite element matrices in vector languages
Cited in
(17)- Optimizing finite element matrix calculations using the general technique of element vectorization
- Geometric Optimization of the Evaluation of Finite Element Matrices
- A stencil scaling approach for accelerating matrix-free finite element implementations
- On Memory Traffic and Optimisations for Low-order Finite Element Assembly Algorithms on Multi-core CPUs
- Numerical integration on GPUs for higher order finite elements
- A discrete-mechanical approach for computation of three-dimensional flows
- Combinatorial optimization of matrix-vector multiplication in finite element assembly
- Automating the finite element method
- Automatic implementation of finite strain anisotropic hyperelastic models using hyper-dual numbers
- High-performance evaluation of finite element variational forms via commuting diagrams and duality
- Finite element numerical integration for first order approximations on multi- and many-core architectures
- Locally assembled stiffness matrix: a novel method to obtain global stiffness matrix
- Automation of primal and sensitivity analysis of transient coupled problems
- The Target-Matrix Optimization Paradigm for High-Order Meshes
- ConvPDE-UQ: convolutional neural networks with quantified uncertainty for heterogeneous elliptic partial differential equations on varied domains
- Topological Optimization of the Evaluation of Finite Element Matrices
- An algorithm for the optimization of finite element integration loops
This page was built for publication: Optimizing the Evaluation of Finite Element Matrices
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5470339)