Added a linear operator base class for LHS and Preconditioner operations with the conj grad solver, identity operator for testing purposes, mass matrix operator without optimisation
This MR has been partially split into: