Sparse matrix and vector classes, solvers (Krylov methods, AMG, ...)