Fast implementation of mixed RT0 finite elements in MATLAB


Author/Creator ORCID




Citation of Original Publication

Theodore Weinberg; Fast implementation of mixed RT0 finite elements in MATLAB; SIAM Undergraduate Research Online (SIURO), vol. 12;;


This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.


We develop a fast implementation of the mixed finite element method for the Darcy’s problem discretized by lowest-order Raviart-Thomas finite elements using Matlab. The implementation is based on the so-called vectorized approach applied to the computation of the finite element matrices and assembly of the global finite element matrix. The code supports both 2D and 3D domains, and the finite elements can be triangular, rectangular, tetrahedral or hexahedral. The code can also be easily modified to import user-provided meshes. We comment on our freely available code and present a performance comparison with the standard approach.