Rahman and Valdman (2013) introduced a vectorized way to assemble finite element stiffness and mass matrices in MATLAB. Local element matrices are computed all at once by array operations and stored in multi-dimentional arrays (matrices). We build some iterative solvers on available multi-dimentional structures completely avoiding the use of a sparse matrix. The implementation of those solvers is directly based on MATLAB codes from Rahman and Valdman (2013).
Type 'start_tensor_solver_Dirichlet' to run the code.
MATLAB Implementation of Element-based Solvers by Leszek Marcinkowski and Jan Valdman.
Inspired by: Fast FEM assembly: nodal elements
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!