CosmoCoffee Forum Index CosmoCoffee

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

CosmoMC: Updated PPF Module for w_0, w_a constraints
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Ben Brown



Joined: 12 Sep 2011
Posts: 4
Affiliation: University of Pittsburgh

PostPosted: September 20 2011  Reply with quote

I've tried to apply Wenjuan's PPF module found here: http://camb.info/ppf/. It no longer compiles, due to updates to the vanilla Camb implementation (esp equations.f90).

How should I modify the new Camb? Should I try to re-implement the modifications in the source to allow w to vary? Or is there an implementation of this already available?

Edit: "new Camb" is the June 2011 version of equations.f90, packaged with the July 2011 version of camb.

Thanks,
Ben
Back to top
View user's profile  
Ben Brown



Joined: 12 Sep 2011
Posts: 4
Affiliation: University of Pittsburgh

PostPosted: September 28 2011  Reply with quote

I'm sure that this is a common interest. I assumed that the old, modified camb source files found here http://camb.info/ppf/ should not be expected to compile, and therefore didn't post compilation errors, etc.

Details:

After following the directions to replace (as opposed to modifying the makefile found in camb/) equations.f90, etc:

Code:
$ make camb
cd ../camb && make --file=Makefile_main libcamb.a
make[1]: Entering directory `/home/ben/packages/cosmomc/camb'
ifort -O2 -W0 -WB -fpp -Vaxlib -c equations.f90
equations.f90(1352): error #6404: This name does not have a type, and must have an explicit type.   [DLFDLQ]
           y(EV%iq0+off_ix+i-1)=-0.25_dl*dlfdlq(i)*InitVec(i_clxr)
-----------------------------------------^
equations.f90(1797): error #6404: This name does not have a type, and must have an explicit type.   [DQ]
                q=(i-0.5_dl)*dq
-----------------------------^
equations.f90(1806): error #6404: This name does not have a type, and must have an explicit type.   [DLFDLQ]
                     -k*2._dl/15._dl*sigma*dlfdlq(i)
-------------------------------------------^
equations.f90(2232): error #6404: This name does not have a type, and must have an explicit type.   [DLFDLQ]


However, getting the old ppf module to work in this manner with the otherwise updated camb is probably not optimal to the community, since the improvements to CAMB are welcome ones.
Back to top
View user's profile  
Jannis Bielefeld



Joined: 03 Oct 2011
Posts: 1
Affiliation: Dartmouth

PostPosted: October 04 2011  Reply with quote

Hi,

I am experiencing the same problem and think that this should be addressed soon, as this package has a broad common interest.

I am also using the July 2011 version of CAMB.
Back to top
View user's profile  
Savvas Nesseris



Joined: 05 Apr 2005
Posts: 50
Affiliation: UAM/IFT

PostPosted: October 04 2011  Reply with quote

I agree that it is very useful.

@Antony: Perhaps the PPF module could be updated and become part of the official code, so that it does not become obsolete?
Back to top
View user's profile [ Hidden ] Visit poster's website
Wenjuan Fang



Joined: 21 Dec 2007
Posts: 3
Affiliation: University of Michigan

PostPosted: October 04 2011  Reply with quote

Thank you all for your interests in the code. I'll take a look at the problem and try to update it as soon as I can.
Back to top
View user's profile [ Hidden ]
Ben Brown



Joined: 12 Sep 2011
Posts: 4
Affiliation: University of Pittsburgh

PostPosted: October 04 2011  Reply with quote

Wenjuan-
Thank you so much for your attention and contribution!
-Ben
Back to top
View user's profile  
Ben Brown



Joined: 12 Sep 2011
Posts: 4
Affiliation: University of Pittsburgh

PostPosted: October 26 2011  Reply with quote

Wenjuan-
Have you had a chance to look at this yet?

I might simply use an older version of cosmomc to get the likelihood constraints I'm calculating using the original ppf module.
-Ben
Back to top
View user's profile  
Wenjuan Fang



Joined: 21 Dec 2007
Posts: 3
Affiliation: University of Michigan

PostPosted: April 27 2012  Reply with quote

Hi guys,

The PPF module has been updated finally (check the PPF webpage). Sorry for such a delay, I was held up by various things related to job application. Right now, the updated code does not support massive neutrinos. If you want to use the PPF module to study massive neutrinos, you may want to use the old code, or check the PPF webpage at a later date. In the updated code, setcgammappf is called in init_background, instead of in inidriver.F90 as in the old code. This ensures that c_Gamma_ppf is set to the correct value while CosmoMC and the modified CAMB are run together. Thanks Michael Mortonson, Antonio Cuesta, and Kushal Mehta for making this suggestion! Let me know if you find any problem with the updated code.

-Wenjuan
Back to top
View user's profile [ Hidden ]
Savvas Nesseris



Joined: 05 Apr 2005
Posts: 50
Affiliation: UAM/IFT

PostPosted: May 08 2012  Reply with quote

Many thanks! If you don't mind, I will check if I can integrate it with my CAMB mathematica interface, which is long overdue for an update anyway :D

Cheers
Back to top
View user's profile [ Hidden ] Visit poster's website
Hamzeh Alavi



Joined: 27 Oct 2011
Posts: 4
Affiliation: KIT

PostPosted: November 19 2012  Reply with quote

Wenjuan Fang wrote:
Hi guys,

The PPF module has been updated finally (check the PPF webpage). Sorry for such a delay, I was held up by various things related to job application. Right now, the updated code does not support massive neutrinos. If you want to use the PPF module to study massive neutrinos, you may want to use the old code, or check the PPF webpage at a later date. In the updated code, setcgammappf is called in init_background, instead of in inidriver.F90 as in the old code. This ensures that c_Gamma_ppf is set to the correct value while CosmoMC and the modified CAMB are run together. Thanks Michael Mortonson, Antonio Cuesta, and Kushal Mehta for making this suggestion! Let me know if you find any problem with the updated code.

-Wenjuan


Hi Wenjuan

I'm trying to add w0_ppf and wa_ppf to cosmomc, but I got this error as mentioned in http://cosmocoffee.info/viewtopic.php?p=5589#5589

Random seeds: 1717, 15182 rand_inst: 0
Computing tensors: F
Doing CMB lensing: T
lmax = 2100
Number of Cls = 4
Varying 9 parameters ( 0 fast)
starting Monte-Carlo
Initialising BBN Helium data...
bad xa input in splint

how should I fix it?
Back to top
View user's profile  
Antony Lewis



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

PostPosted: November 19 2012  Reply with quote

Does it work if you run LCDM model with your modified code? I have been running cosmomc chains uing equations_ppf.f90 included with the latest CAMB version no problem.
Back to top
View user's profile [ Hidden ] Visit poster's website
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.