Fast FEM assembly: edge elements

Vectorized assembly of FEM matrices.
2.8K Downloads
Updated 3 Aug 2018

View License

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.

Cite As

Jan Valdman (2026). Fast FEM assembly: edge elements (https://nl.mathworks.com/matlabcentral/fileexchange/46635-fast-fem-assembly-edge-elements), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2012b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Verification, Validation, and Test in Help Center and MATLAB Answers
Version Published Release Notes
1.2.0.0

Title a summary updated.

1.1.0.0

New code version, description updated with a link to arxiv version of the explaining paper.
New feature added: piecewise constant coefficients in FEM matrices.
New feature added: piecewise constant coefficients in FEM matrices.
New feature added: piecewise constant coefficients in FEM matrices.
New feature added: piecewise constant coefficients in FEM matrices.

1.0.0.0