The current code calculates the Hessian contribution of the spatial smoother, and this seems to work fine. Need to do a better job with the Hessian calculations of the ImageGP part, but that's work for when i grow up...
We also return the forward modelled observations, which will be returned to the user for inspection etc.
The GP-based methods store the forward modelled observations, so that at the end of the processing, these are available for each operator (I think!)
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
Move the fwd model calculations to the init method
Fixed a couple of typos
Now we save the intermediate Hessians as a pickle, as converting them to numpy arrays breaks all hell. Theres a new option to pass an output name (otherwise, it gives a default filename)
ImageOp Hessian calculations appear to work but take forever :(