An efficient way to assemble finite element matrices in vector languages
From MaRDI portal
Publication:329010
Abstract: Efficient Matlab codes in 2D and 3D have been proposed recently to assemble finite element matrices. In this paper we present simple, compact and efficient vectorized algorithms, which are variants of these codes, in arbitrary dimension, without the use of any lower level language. They can be easily implemented in many vector languages (e.g. Matlab, Octave, Python, Scilab, R, Julia, C++ with STL,...). The principle of these techniques is general, we present it for the assembly of several finite element matrices in arbitrary dimension, in the P1 finite element case. We also provide an extension of the algorithms to the case of a system of PDE's. Then we give an extension to piecewise polynomials of higher order. We compare numerically the performance of these algorithms in Matlab, Octave and Python, with that in FreeFEM++ and in a compiled language such as C. Examples show that, unlike what is commonly believed, the performance is not radically worse than that of C : in the best/worst cases, selected vector languages are respectively 2.3/3.5 and 2.9/4.1 times slower than C in the scalar and vector cases. We also present numerical results which illustrate the computational costs of these algorithms compared to standard algorithms and to other recent ones.
Recommendations
- Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
- Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements
- \(o\) FEM: an object oriented finite element package for Matlab
- Fast native-MATLAB stiffness assembly for SIPG linear elasticity
- scientific article; zbMATH DE number 4197325
Cites work
- scientific article; zbMATH DE number 5719436 (Why is no real title available?)
- scientific article; zbMATH DE number 635667 (Why is no real title available?)
- scientific article; zbMATH DE number 1131501 (Why is no real title available?)
- scientific article; zbMATH DE number 5179816 (Why is no real title available?)
- An efficient way to assemble finite element matrices in vector languages
- Direct Methods for Sparse Linear Systems
- Efficient implementation of adaptive P1-FEM in Matlab
- Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements
- Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
- Finite Element Methods and Their Applications
- Finite element method.
- Galerkin Finite Element Methods for Parabolic Problems
- Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities
- Matlab implementation of the finite element method in elasticity
- New development in freefem++
- Numerical modeling for differential problems
- Scientific computing with MATLAB and Octave
- The finite element methods for elliptic problems.
Cited in
(17)- FEM and BEM simulations with the Gypsilab framework
- Efficient P1-FEM for any space dimension in \textsc{Matlab}
- FEM-BEM coupling methods for tokamak plasma axisymmetric free-boundary equilibrium computations in unbounded domains
- Vectorizing matrix operations arising from PDE discretization on 9-point stencils
- Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements
- An unstructured mesh finite difference/finite element method for the three-dimensional time-space fractional Bloch-Torrey equations on irregular domains
- Fast assembly of Galerkin matrices for 3D solid laminated composites using finite element and isogeometric discretizations
- A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries
- Fast native-MATLAB stiffness assembly for SIPG linear elasticity
- Vectfem: a generalized MATLAB-based vectorized algorithm for the computation of global matrix/force for finite elements of any type and approximation order in linear elasticity
- An efficient way to assemble finite element matrices in vector languages
- On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems
- Efficient multi-frequency solutions of FE-be coupled structural-acoustic problems using Arnoldi-based dimension reduction approach
- A posteriori stopping criteria for space-time domain decomposition for the heat equation in mixed formulations
- Optimizing the Evaluation of Finite Element Matrices
- The generalized finite difference method with third- and fourth-order approximations and treatment of ill-conditioned stars
- Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
Describes a project that uses
Uses Software
This page was built for publication: An efficient way to assemble finite element matrices in vector languages
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q329010)
