This solver assembles and solves the FEM for the 2D scalar Helmholtz equation, using P1 triangular elements. The assembly is fully vectorized for efficiency.
For tutoring, the script of a "pedagogic" naive assembly is also provided in comments.
The resolution is performed using the Matlab \ operator (here leads to LU factorization).
Five different geometries are provided using built-in Matlab tools, but the solver is also compatible with arbitrary geometries. In particular, external mesh files (for instance from GMSH) can be used, while satisfying the PET format.
david gasperini (2021). FEM solver for 2D Helmholtz equation (https://www.mathworks.com/matlabcentral/fileexchange/91695-fem-solver-for-2d-helmholtz-equation), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!