Yes. I figured this. A bit too late though 
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.
Yes. I figured this. A bit too late though 
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.