Fast implementation of mixed RT0 finite elements in MATLAB
No Thumbnail Available
Permanent Link
Author/Creator
Author/Creator ORCID
Date
Type of Work
Department
Program
Citation of Original Publication
Theodore Weinberg; Fast implementation of mixed RT0 finite elements in MATLAB; SIAM Undergraduate Research Online (SIURO), vol. 12; https://www.siam.org/Portals/0/Publications/SIURO/Vol12/S01743.pdf;
Rights
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.
Abstract
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.