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" |