Rebecca-Rumney

I've been looking into how to do this and this is my idea: It requires changing the `__call__` function of each ProblemLogLikelihood and adding 2 new functions so that: def...

@DavAug That's a good point. What I've written there only saves the last step's log-likelihoods (so of size N) rather than the whole N x M x K matrix and...

I have been using scipy optimisation, and `fmin` is a legacy function. The current `minimize` function uses fatol and xatol for absolute differences and uses ftol and xtol for relative...