CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

Problem with python vs fortran getdist
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Sunny Vagnozzi



Joined: 15 Aug 2016
Posts: 41
Affiliation: Oskar Klein Centre, Stockholm University

PostPosted: August 24 2017  Reply with quote

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:
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
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1300
Affiliation: University of Sussex

PostPosted: August 24 2017  Reply with quote

Are you using the latest getdist python version?
Back to top
View user's profile [ Hidden ] Visit poster's website
Sunny Vagnozzi



Joined: 15 Aug 2016
Posts: 41
Affiliation: Oskar Klein Centre, Stockholm University

PostPosted: August 24 2017  Reply with quote

Antony Lewis wrote:
Are you using the latest getdist python version?


Yep, downloaded the latest version. Not sure what's going wrong though.
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1300
Affiliation: University of Sussex

PostPosted: August 24 2017  Reply with quote

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.
Back to top
View user's profile [ Hidden ] Visit poster's website
Sunny Vagnozzi



Joined: 15 Aug 2016
Posts: 41
Affiliation: Oskar Klein Centre, Stockholm University

PostPosted: August 24 2017  Reply with quote

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
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.