Fast FEM assembly: edge elements

Version 1.2.0.0 (80.4 KB) by Jan Valdman

Vectorized assembly of FEM matrices.

5.0

(7)

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://www.mathworks.com/matlabcentral/fileexchange/46635-fast-fem-assembly-edge-elements), MATLAB Central File Exchange. Retrieved January 9, 2026.

MATLAB Release Compatibility

Created with R2012b

Compatible with any release

Platform Compatibility

Windows macOS Linux

Tags Add Tags

Acknowledgements

Inspired by: Fast FEM assembly: nodal elements, inttet, calc_meshdata

Inspired: Continuum undergoing combined elasto-plasto-damage transformation., Fast FEM evaluation of nonlinear energies: hp elements, Efficient vectorization techniques for FEM computations, p-Laplace equation solver using 1D, 2D FEM, Fast-Implementation-Mixed-FEM, Fast FEM evaluation of nonlinear energies: nodal elements, Hyper elasticity with a non-penetration condition, Implementation of C1 FEM, Simulation of von Kármán viscoelastic plates, Hierarchical hp FEM on rectangles in 2D

software/debug

software/example_comparison

software/example_eddycurrent

software/example_eddycurrent/test2D_1

software/example_eddycurrent/test2D_2

software/example_eddycurrent/test2D_3

software/example_eddycurrent/test3D_1

software/example_eddycurrent/test3D_2

software/example_majorant

software/example_majorant/test2D

software/example_majorant/test3D

software/path

software/path/library_fem

software/path/library_integration

software/path/library_vectorization

Version	Published	Release Notes
1.2.0.0	3 Aug 2018	Title a summary updated.	Download
1.1.0.0	30 Oct 2015	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.	Zip Toolbox
1.0.0.0	16 May 2014		Download