Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements
It provides an efficient implementation for researchers using edge finite elements in MATLAB, though the approach is incremental.
The paper presents a fast and flexible MATLAB assembly method for finite element stiffness and mass matrices, specifically for edge elements (Raviart-Thomas and Nédélec), achieving scalable performance.
We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time.