Fast MATLAB assembling functions for 2D/3D FEM Matrices
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)
引用格式
Jonas Koko (2024). Fast MATLAB assembling functions for 2D/3D FEM Matrices (https://www.mathworks.com/matlabcentral/fileexchange/59616-fast-matlab-assembling-functions-for-2d-3d-fem-matrices), MATLAB Central File Exchange. 检索时间: .
MATLAB 版本兼容性
平台兼容性
Windows macOS Linux标签
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!KPDE/
版本 | 已发布 | 发行说明 | |
---|---|---|---|
1.0.0.0 |