Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements
Immanuel Anjam, Jan Valdman
TL;DR
This paper introduces a fast, scalable MATLAB method for assembling finite element matrices specifically for edge elements like Raviart-Thomas and Nedelec, improving efficiency in 2D and 3D discretizations.
Contribution
It presents a novel vectorization approach for assembling FEM matrices in MATLAB tailored for edge elements, enhancing speed and flexibility.
Findings
Significant reduction in assembly time for FEM matrices.
Method is scalable for large 2D and 3D problems.
Provides a freely available MATLAB implementation.
Abstract
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.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.