Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements. (English) Zbl 1288.65169

Summary: We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. The major loops in the code have been vectorized using the so called array operation in MATLAB, and no low level languages like the C or Fortran has been used for the purpose. The implementation is based on having the vectorization part separated, in other words hidden, from the original code thereby preserving its original structure, and its flexibility as a finite element code. The code is fast and scalable with respect to time.


65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
Full Text: DOI


[1] Alberty, J.; Carstensen, C.; Funken, S. A., Remarks around 50 lines of MATLAB: short finite element implementation, Numer. Algorithms, 20, 117-137, (1999) · Zbl 0938.65129
[2] Alberty, J.; Carstensen, C.; Funken, S. A.; Klose, R., MATLAB implementation of the finite element method in elasticity, Computing, 69, 236-263, (2002) · Zbl 1239.74092
[3] S. Funken, D. Praetorius, P. Wissgott, Efficient implementation of adaptive P1-FEM in MATLAB, ASC Report 19/2008, Institute for Analysis and Scientific Computing, Vienna University of Technology, Wien, 2008. · Zbl 1284.65197
[4] Gockenbach, M. S., Understanding and implementing the finite element method, (2006), SIAM · Zbl 1105.65112
[5] Koko, J., vectorized MATLAB codes for linear two-dimensional elasticity, Sci. Program., 15, 3, 157-172, (2007)
[6] Persson, P.-O.; Strang, G., A simple mesh generation in MATLAB, SIAM Rev., 42, 329-345, (2004) · Zbl 1061.65134
[7] Rahman, T., SERF2D-MATLAB (ver. 1.1) - documentation, (2003), University of Augsburg
[8] Smith, I. M.; Griffiths, D. V., Programming the finite element method, (2004), John Wiley & Sons · Zbl 1091.74002
[9] J. Valdman, Minimization of Functional Majorant in A Posteriori Error Analysis based on H(div) Multigrid-Preconditioned CG Method, Advances in Numerical Analysis, vol. 2009, Article ID 164519, 2009. · Zbl 1200.65095
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.