$\ell$FEM: An efficient loop-free Matlab implementation of isoparametric bulk and surface finite elements
Numerical Analysis
2026-05-15 v1 Numerical Analysis
Abstract
The FEM MATLAB package provides a simple, efficient, and flexible implementation of isoparametric finite elements in bulk domains and on surfaces. The finite element matrix assemblies are based on MATLAB's paged operators and therefore completely loop-free. We give a short and conscious description of high-order isoparametric surface finite elements, which is then used to describe the assembly process and the implementation. We report on relevant numerical experiments (runtime comparisons, modifications for non-linear problems, etc.), and on additional functions, examples, and a testing unit which are all part of the FEM package.
Keywords
Cite
@article{arxiv.2605.14035,
title = {$\ell$FEM: An efficient loop-free Matlab implementation of isoparametric bulk and surface finite elements},
author = {Balázs Kovács and Michael Lantelme},
journal= {arXiv preprint arXiv:2605.14035},
year = {2026}
}