File Exchange

image thumbnail

Fast MATLAB assembling functions for 2D/3D FEM Matrices

version (19.7 KB) by Jonas Koko
MATLAB assembling functions for 2D/3D Poisson and linear elasticity problems


Updated 09 Oct 2016

View License

We propose a MATLAB implementation of the P 1 finite element method for the numerical solutions of the Poisson problem and the linear elasticity problem in two-dimensional (2D) and three-dimensional (3D). The code consists of vectorized (and short) assembling functions for the matrices (mass and stiffness) and the right-hand sides. Since for the P 1 finite element, the element mass matrix and right-hand side are simple, the implementation uses only the MATLAB function sparse on the elements volume. For the stiffness matrix, to obtain a MATLAB implementation close to the standard form, cell-arrays are used to store the gradients of the element basis functions. The assembling procedure can then use matrix/vector products on small size cell-arrays. Numerical experiments show that our implementation is fast and scalable with respect to time.
To cite
Koko J., Fast MATLAB assembly of FEM matrices in 2D and 3D using cell-array approach, International Journal of Modeling, Simulation, and Scientific Computing Vol. 7, No. 3 (2016)

Cite As

Jonas Koko (2020). Fast MATLAB assembling functions for 2D/3D FEM Matrices (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (0)

MATLAB Release Compatibility
Created with R2014b
Compatible with any release
Platform Compatibility
Windows macOS Linux