Problem with python vs fortran getdist
Posted: August 24 2017
Hi Antony et al.,
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
whereas the fortran version does not.
Any ideas on what is causing this and how to solve it?
Cheers,
Sunny
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?
Cheers,
Sunny