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-dimensional arrays (matrices). We build some iterative solvers on available multi-dimensional structures completely avoiding the use of a sparse matrix.