The following Matlab project contains the source code and Matlab examples used for fast assembly of stiffness and matrices in finite element method using edge elements.
This code demonstrates vectorization concepts from the prepared paper
Immanuel Anjam and Jan Valdman: Fast MATLAB assembly of FEM matrices in 2D and 3D: edge elements.
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 |
|---|---|
| license.txt | 1533 |
| add_paths.m | 244 |
| mesh3D.mat | 819 |
| start_2D.m | 5875 |
| start_3D.m | 5907 |
| start_assembly_Nedelec0_2D.m | 1482 |
| start_assembly_Nedelec0_3D.m | 1256 |
| start_assembly_P1_2D.m | 1711 |
| start_assembly_P1_3D.m | 1484 |
| start_assembly_RT0_2D.m | 1473 |
| start_assembly_RT0_3D.m | 1290 |
| add_paths.m | 244 |
| num_error.m | 1548 |
| solver_eddycurrent.m | 2190 |
| start_2D.m | 3344 |
| start_3D.m | 2246 |
| exact.m | 376 |
| loadfun.m | 301 |
| mesh.mat | 356 |
| exact.m | 474 |
| loadfun.m | 308 |
| mesh.mat | 819 |
| add_paths.m | 244 |
| majorant.m | 2137 |
| norm_L2.m | 1000 |
| num_error.m | 1338 |
| saveTightFigure.m | 766 |
| solver_majorant.m | 2848 |
| solver_poisson.m | 1728 |
| solver_poisson_orig.m | 2356 |
| start_2D.m | 5876 |
| start_3D.m | 3576 |
| const.mat | 179 |
| exact.m | 427 |
| loadfun.m | 283 |
| mesh.mat | 356 |
| const.mat | 179 |
| exact.m | 557 |
| loadfun.m | 375 |
| mesh.mat | 819 |
| affine_transformations.m | 1919 |
| copy_triu.m | 600 |
| deleterepeatedrows.m | 1668 |
| getBoundaryFaces.m | 1429 |
| getDirichletEdges.m | 1450 |
| getEdges.m | 1828 |
| getFaces.m | 1046 |
| basis_Nedelec0.m | 3132 |
| basis_P1.m | 2304 |
| basis_RT0.m | 2347 |
| mass_matrix_Nedelec0.m | 1830 |
| mass_matrix_P1.m | 1568 |
| mass_matrix_P1_2D.m | 644 |
| mass_matrix_P1_3D.m | 1583 |
| mass_matrix_RT0.m | 1706 |
| solution_Nedelec0.m | 3529 |
| solution_P1.m | 2452 |
| solution_RT0.m | 2804 |
| stifness_matrix_Nedelec0.m | 2438 |
| stifness_matrix_P1.m | 1874 |
| stifness_matrix_P1_2D.m | 1243 |
| stifness_matrix_P1_3D.m | 1281 |
| stifness_matrix_RT0.m | 1569 |
| intquad.m | 815 |
| inttet.m | 24937 |
| inttri.m | 28635 |
| 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 | 7355 |
| shapefun.m | 4144 |
| smamt.m | 273 |
| svamt.m | 260 |
| refinement_uniform.m | 1634 |
| refinement_uniform3D.m | 439 |
| show_constant_scalar.m | 595 |
| show_mesh.m | 569 |
| show_nodal_scalar.m | 460 |
| signs_edges.m | 1512 |
| signs_faces.m | 2768 |
| tetrarefine3.m | 2619 |
