An efficient way to assemble finite element matrices in vector languages
From MaRDI portal
Publication:329010
DOI10.1007/S10543-015-0587-4zbMATH Open1351.65088arXiv1401.3301OpenAlexW2963397658MaRDI QIDQ329010FDOQ329010
Authors: François Cuvelier, Caroline Japhet, Gilles Scarella 
Publication date: 21 October 2016
Published in: BIT (Search for Journal in Brave)
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.
Full work available at URL: https://arxiv.org/abs/1401.3301
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
performancealgorithmnumerical resultfinite elementsvectorizationMatlabPythonmatrix assemblyOctavevector languages
Complexity and performance of numerical algorithms (65Y20) Packaged methods for numerical algorithms (65Y15) Boundary value problems for second-order elliptic equations (35J25) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)
Cites Work
- Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities
- New development in freefem++
- Direct Methods for Sparse Linear Systems
- Title not available (Why is that?)
- The finite element methods for elliptic problems.
- Title not available (Why is that?)
- Galerkin Finite Element Methods for Parabolic Problems
- Title not available (Why is that?)
- Finite Element Methods and Their Applications
- Efficient implementation of adaptive P1-FEM in Matlab
- Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
- Finite element method.
- An efficient way to assemble finite element matrices in vector languages
- Matlab implementation of the finite element method in elasticity
- Title not available (Why is that?)
- Scientific Computing with MATLAB and Octave
- Numerical modeling for differential problems
- Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements
Cited In (16)
- Efficient P1-FEM for any space dimension in \textsc{Matlab}
- An unstructured mesh finite difference/finite element method for the three-dimensional time-space fractional Bloch-Torrey equations on irregular domains
- The generalized finite difference method with third- and fourth-order approximations and treatment of ill-conditioned stars
- Fast assembly of Galerkin matrices for 3D solid laminated composites using finite element and isogeometric discretizations
- 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
- Optimizing the Evaluation of Finite Element Matrices
- An efficient way to assemble finite element matrices in vector languages
- FEM and BEM simulations with the Gypsilab framework
- A posteriori stopping criteria for space-time domain decomposition for the heat equation in mixed formulations
- Vectorizing matrix operations arising from PDE discretization on 9-point stencils
- Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements
- Fast native-MATLAB stiffness assembly for SIPG linear elasticity
- On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems
- A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries
- FEM-BEM coupling methods for tokamak plasma axisymmetric free-boundary equilibrium computations in unbounded domains
- Efficient multi-frequency solutions of FE-be coupled structural-acoustic problems using Arnoldi-based dimension reduction approach
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)
