Element-based Solvers

2D Poisson's solver using only local FEM matrices without the global matrix assembly.

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.