This code demonstrates vectorization concepts from the the paper
Immanuel Anjam, Jan Valdman: Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements. Applied Mathematics and Computation 267, 252–263 (2015)
We extended techniques from the 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)
to a fast assembly of FEM matrices using edge elements - Raviart-Thomas elements for Hdiv problems and Nedelec elements for Hcurl problems. In addition, vectorized higher order quadratures were added.
A link to the paper can be found at the author web page located at http://sites.google.com/site/janvaldman/publications
Please cite the paper if you find the code useful.
To compare the assembly times, call
"start_2D" or "start_3D" in the "example_comparison" directory.
You can also call
"start_2D" or "start_3D" in "example_majorant" and "example_eddycurrect" directories
to obtain solution of the functional majorant minimization in Hdiv space and a solution of a eddy current problem in Hcurl space.
Jan Valdman (2021). Fast FEM assembly: edge elements (https://www.mathworks.com/matlabcentral/fileexchange/46635-fast-fem-assembly-edge-elements), MATLAB Central File Exchange. Retrieved .
Inspired: Continuum undergoing combined elasto-plasto-damage transformation., Fast-Implementation-Mixed-FEM, Hyper elasticity with a non-penetration condition , Implementation of C1 FEM, Simulation of von Kármán viscoelastic plates
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!