2000 character limit reached
Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements (1409.4618v2)
Published 16 Sep 2014 in cs.MS and math.NA
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.