Page 1 of 1

Problem with python vs fortran getdist

Posted: August 24 2017
by Sunny Vagnozzi
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

Code: Select all

RuntimeWarning: divide by zero encountered in divide
  corrected = (bins2D[ix] * A + xP * Ax + yP * Ay) / denom
whereas the fortran version does not.

Any ideas on what is causing this and how to solve it?

Cheers,
Sunny

Re: Problem with python vs fortran getdist

Posted: August 24 2017
by Antony Lewis
Are you using the latest getdist python version?

Re: Problem with python vs fortran getdist

Posted: August 24 2017
by Sunny Vagnozzi
Antony Lewis wrote:Are you using the latest getdist python version?
Yep, downloaded the latest version. Not sure what's going wrong though.

Re: Problem with python vs fortran getdist

Posted: August 24 2017
by Antony Lewis
If you can email me a link to the chains I can have a look at what's going on with the bottom-left plot.

Problem with python vs fortran getdist

Posted: August 24 2017
by Sunny Vagnozzi
Thanks, I've just sent you an email but if can be useful for anyone else, here is a public link to the chains and distparams.ini used: https://uploadfiles.io/nrg8f