Mainly bugfixes. One was that the partial derivatives of the observation operator cost function were not updated if the parameter was set as CONSTANT. There was a summation over all possible sites to be performed that was missing. Other stuff includes some generalisation of the hessians.
Solved some sparse definitions, added a quick sparse hessian for the prior
Added a per band fit for the time series object. This is the same kind of thing as the image GP function.
Ooops, cost and derivative were undefined!!!
So the der_cost method is now able to use per band emulators, which should result in faster performance.
A couple of small fixes. It turns out that for the time series stuff, the per band is overkill most of the time
Weird happenings, examples not quite working. commiting and going back to see what's different