MATLAB Implementation of Element-based Solvers

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.