Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
Talal Rahman and Jan Valdman
Contact the author: Please use for correspondence this email.
Submission date: 02. Mar. 2011 (revised version: August 2011)
published in: Applied mathematics and computation, 219 (2013) 13, p. 7151-7158
DOI number (of the published article): 10.1016/j.amc.2011.08.043
MSC-Numbers: 68N30, 65M60
Keywords and phrases: MATLAB vectorization, finite elements, elasticity
Download full preprint: PDF (307 kB)
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.