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 by: Fast FEM assembly: nodal elements, inttet, calc_meshdata
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!Create scripts with code, output, and formatted text in a single executable document.