Abstract for the talk at 23.03.2011 (10:00 h)Oberseminar NUMERIK / WISSENSCHAFTLICHES RECHNEN
Jan Valdman (MPI MIS, Leipzig + University of Iceland, Reykjavik)
Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements
We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. The major loops in the code have been vectorized using the so called array operation in MATLAB, and no low level languages like the C or Fortran has been used for the purpose. The implementation is based on having the vectorization part separated, in other words hidden, from the original code thereby preserving its original structure, and its flexibility as a finite element code. The code is fast and scalable. This is a joint work with T. Rahman (Bergen).
Report available: www.mis.mpg.de/publications/other-series/tr/report-1111.html