recently I switched from using the fortran version of getdist to using the python one. However, I noticed differences between the results produced by the two. Specifically, the triangle plots I obtain after analysing chains with the Python version are sometimes quite different. Please see the files fortran.pdf and python.pdf where I compare the results obtained after using the fortran getdist and the python getdist respectively (same distparams.ini, same chains)!
By the way, the python getdist gives me the following warning
Code: Select all
RuntimeWarning: divide by zero encountered in divide corrected = (bins2D[ix] * A + xP * Ax + yP * Ay) / denom
Any ideas on what is causing this and how to solve it?