This code demonstrates vectorization concepts from a report by
Talal Rahman and Jan Valdman
on
Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
Technical report 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"
|