  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.
