The following Matlab project contains the source code and Matlab examples used for fast assembly of stiffness and matrices in finite element method using nodal elements.
This code demonstrates vectorization concepts from 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).
The source code and files included in this project are listed in the project files section, please make sure whether the listed source code meet your needs there.
Project Files:
| File Name | Size |
|---|---|
| deleterepeatedrows.m | 1668 |
| elements.mat | 228 |
| getEdges.m | 1784 |
| getFaces.m | 1545 |
| hexmesh.mat | 898280 |
| integration_point_transformation.m | 564 |
| aminv.m | 1731 |
| amsm.m | 271 |
| amsv.m | 235 |
| amt.m | 382 |
| amtam.m | 280 |
| astam.m | 231 |
| avtam.m | 213 |
| phider.m | 1710 |
| shapeder.m | 8191 |
| shapefun.m | 4725 |
| smamt.m | 273 |
| svamt.m | 260 |
| license.txt | 1501 |
| mass_matrixP1_2D.m | 214 |
| mass_matrixP1_2D_elasticity.m | 475 |
| mass_matrixP1_3D.m | 1001 |
| mass_matrixP1_3D_elasticity.m | 1480 |
| mesh3D.mat | 816 |
| mesh_ring.mat | 2862 |
| refinement_uniform.m | 1632 |
| refinement_uniform3D.m | 425 |
| show_mesh.m | 1272 |
| start_assembly_P1_2D.m | 1812 |
| start_assembly_P1_2D_elasticity.m | 2144 |
| start_assembly_P1_3D.m | 1072 |
| start_assembly_P1_3D_elasticity.m | 1028 |
| start_assembly_Q1_3D_elasticity.m | 1567 |
| start_compute_areas.m | 1991 |
| stifness_matrixP1_2D.m | 735 |
| stifness_matrixP1_2D_elasticity.m | 985 |
| stifness_matrixP1_3D.m | 846 |
| stifness_matrixP1_3D_elasticity.m | 1182 |
| stifness_matrixQ1_3D_elasticity.m | 1562 |
| tetrarefine3.m | 2619 |
