Added linear operator, mass and identity, and ConjGrad operators, and implementations

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:

• Identity matrix implementation (see !36 (merged))
• Mass matrix implementation (see !35 (merged))