Fast assembly of stiffness and matrices in finite element method

This code demonstrates vectorization concepts from a paper

Talal Rahman and Jan Valdman:

Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements,

Applied Mathematics and Computation 219, 7151–7158 (2013)

Original technical report from 2011 is available at

http://www.mis.mpg.de/publications/other-series/tr/report-1111.html

To run the code, call

"start_assembly_P1_3D" or

"start_assembly_P1_2D"

in the Matlab window (for 3D or 2D cases)

Considering a given 3D or 2D triangulation, mass and stiffness matrices corresponding to linear nodal elements are generated in an efficient vectorized and fast way.

Recently, functionality was extended to generate mass and stiffness matrices for a (vector) linear elasticity problem.

To run the code, call

"start_assembly_P1_3D_elasticity" or

"start_assembly_P1_2D_elasticity"