Fast FEM assembly: edge elements

Vectorized assembly of FEM matrices.

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.