An efficient way to perform the assembly of finite element matrices in Matlab and Octave

Abstract

We describe different optimization techniques to perform the assembly of finite element matrices in Matlab and Octave, from the standard approach to recent vectorized ones, without any low level language used. We finally obtain a simple and efficient vectorized algorithm able to compete in performance with dedicated software such as FreeFEM++. The principle of this assembly algorithm is general, we present it for different matrices in the P1 finite elements case and in linear elasticity. We present numerical results which illustrate the computational costs of the different approaches

Inria: No: RR-8305 (2013)

Cited by in corpus (4)

• Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements
• Fast Matlab compatible sparse assembly on multicore computers
• The optimal multilevel Monte-Carlo approximation of the stochastic drift-diffusion-Poisson system
• On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems