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