Added wrappers for CGLS and LSMR least squares solvers.
Overloaded the sum function to allow spot operators as input.