BUGFIX You can now pass a vector/array of gamma/regularisation parameters for the spatial regularisation, and it will get picked up
Provide timing information, as BFGS is a bit shit in this respect
use scipy.optimize.minimize rather than the actual solver. This opens the possibility of a dictionary of configuration parameters etc. but not yet implemented