Unconsistent logp, Evaluating the results of find_MAP on Custom Likelihood

Yes. I figured this. A bit too late though :slight_smile:
On a side note, I got distracted by the fact that the logp was not deterministic, I am getting different results at each run. Going through the source code, find_MAP seems to use scipy’s optimizer.