CosmoCoffee Forum Index CosmoCoffee

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

Tutorial in CosmoMC
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: March 03 2017  Reply with quote

Hi everyone,

I am learning CosmoMc to get best parameter fit to my model of CMB tempertaure Data(Planck Data).My model includes 3 user defined parameters and 6 cosmological parameters as in Lambda-CDM model.So I need to make a best fit of 9 parameters.As I am new how to use cosmoMC, Can u suggest a basic tutorial of CosmoMC.I have installed it and compiled it and it compiled without any errors.What to do next?

Thanks in Advance!!

Purnalingam
Back to top
View user's profile [ Hidden ] Visit poster's website
Erick Almaraz



Joined: 10 Jun 2012
Posts: 24
Affiliation: University of Mexico

PostPosted: March 03 2017  Reply with quote

Check out Antonio Cuesta's slides

Just google Antonio Cuesta CosmoMC
Back to top
View user's profile [ Hidden ]
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: March 04 2017  Reply with quote

I have seen all the three power points.I have some doubts.

1. What is the difference in parameter extraction with (action=0) and (action=2).

2.There is no params.ini file in CosmoMC main directory.Only test.ini is present.

3.I am using planck likelihood and planck data.After building planck likelihood I get a file named test_planck.ini. When I try compiling cosmomc using this file I get errors.

4.What does building plank likelihood actually do? Does it just take plank data or is the likelihood function is itself different?

5.For adding my own parameters,I have edited my camb code.Where should i edit my CosmoMC code?

Thanks in advance !!
Back to top
View user's profile [ Hidden ] Visit poster's website
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: March 04 2017  Reply with quote

One more thing,
what is the difference between disttest.ini and distparams.ini
Back to top
View user's profile [ Hidden ] Visit poster's website
Antony Lewis



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

PostPosted: March 04 2017  Reply with quote

What the different action options do is described in the readme:

http://cosmologist.info/cosmomc/readme.html
Back to top
View user's profile [ Hidden ] Visit poster's website
Bob Zou



Joined: 21 Jul 2016
Posts: 1
Affiliation: Beijing Institute of Technology

PostPosted: March 15 2017  Reply with quote

Hi, everyone
I have followed Antonio Cuesta's tutorial "Cosmology with CosmoMC
Solution to Exercises" , succeed in running a parameterization=background case. The tutorial is so good, thanks @Antonio Cuesta.
My question is about function "rombint" used in Planckprior.f90, chapter "Extra exercise: Background parameterization" of Antonio Cuesta's tutorial , which is the Romberg method for integration. In the module calling "rombint", I have not found the defination of function “rombint”, neither in the modules Planckprior.f90 used. But the integration works well, giving right values every time when called. I do not understand How and where CosmoMC call function "rombint".

Thanks for replying.
Good Luck.
Back to top
View user's profile  
Antony Lewis



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

PostPosted: March 15 2017  Reply with quote

It's in CAMB (subroutines.f90).
Back to top
View user's profile [ Hidden ] Visit poster's website
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: April 03 2017  Reply with quote

I am using Nov 2016 CosmoMc. Can you please give me the procedure to add a new theoritical parameter in COSMOMC.

Thanks
Back to top
View user's profile [ Hidden ] Visit poster's website
Sofie Marie Koksbang



Joined: 27 Feb 2017
Posts: 9
Affiliation: Aarhus University

PostPosted: April 03 2017  Reply with quote

Hi,

A description of how to add theory parameters is given e.g. at http://cosmocoffee.info/viewtopic.php?t=1149&postdays=0&postorder=asc&start=15
Back to top
View user's profile  
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: April 05 2017  Reply with quote

I have set MPI_Converge_Stop = 0.1 but the code is still continuing after R−1 has reached 0.1.Why?

updating proposal density
Chain:4 drag accpt: 8.1330597E−02 fast/slow 17.71744 slow: 50432
Chain:5 drag accpt: 8.2306042E−02 fast/slow 17.70989 slow: 50378
Chain:1 drag accpt: 8.1103392E−02 fast/slow 17.70531 slow: 50382
Chain:7 drag accpt: 8.0609366E−02 fast/slow 17.71154 slow: 50007
Chain 6 MPI communicating
Chain 2 MPI communicating
Chain 8 MPI communicating
Chain 5 MPI communicating
Chain 7 MPI communicating
Chain 3 MPI communicating
Chain 4 MPI communicating
Chain 1 MPI communicating
Current convergence R−1 = 8.7707318E−02 chain steps = 82017
slow changes 50723 power changes 0
omegabh2 lim err 0.092
omegabh2 lim err 1.042
omegach2 lim err 0.216
omegach2 lim err 0.396
theta lim err 0.401
theta lim err 0.112
logA lim err 0.077
logA lim err 0.274
ns lim err 0.055
ns lim err 0.051
calPlanck lim err 0.092
calPlanck lim err 0.121
BBdust lim err 0.097
BBdust lim err 0.100
BBbetadust lim err 0.095
BBbetadust lim err 0.094
Current limit err = 1.0421 for omegabh2; samps = 81931
Back to top
View user's profile [ Hidden ] Visit poster's website
Antony Lewis



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

PostPosted: April 05 2017  Reply with quote

Probably also have MPI_Check_Limit_Converge =T?
Back to top
View user's profile [ Hidden ] Visit poster's website
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: April 05 2017  Reply with quote

yep MPI_Check_Limit_Converge =T but what does this do?
Back to top
View user's profile [ Hidden ] Visit poster's website
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: April 05 2017  Reply with quote

I have added 3 new parameters w0,g0,a0 and did the necessary modifications in camb and cosmomc but I get error as the following(action = 2)
Number of MPI processes: 2
file_root:test
Random seeds: 1656, 1087 rand_inst: 1
Random seeds: 1756, 1087 rand_inst: 2
Doing non-linear Pk: F
Doing CMB lensing: T
Doing non-linear lensing: T
TT lmax = 2500
EE lmax = 2500
ET lmax = 2500
BB lmax = 2500
PP lmax = 2500
lmax_computed_cl = 2500
Computing tensors: F
max_eta_k = 14000.00
transfer kmax = 5.000000
adding parameters for: smica_g30_ftl_full_pp
adding parameters for: BKPlanck_detset_comb_dust
Fast divided into 2 blocks
Block breaks at: 10
11 parameters ( 8 slow ( 0 semi-slow), 3 fast ( 0 semi-fast))
skipped unused params: acib217 xi asz143 aps100 aps143 aps143217 aps217 aksz kgal100 kgal143 kgal143217 kgal217 cal0 cal2
finding best fit point...
minmizing fast parameters
minimizing 3 with rhobeg, rhoend = 3.000000 5.0000001E−02
minimization points: 7 chi2 converge tol: 5.0000001E−02
Compile with CosmoRec to use fdm
MpiStop: 0
————————————————————————–
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 538976288.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
————————————————————————–
skipped unused params: acib217 xi asz143 aps100 aps143 aps143217 aps217 aksz kgal100 kgal143 kgal143217 kgal217 cal0 cal2
minmizing fast parameters
minimizing 3 with rhobeg, rhoend = 3.000000 5.0000001E−02
minimization points: 7 chi2 converge tol: 5.0000001E−02
forrtl: error (78): process killed (SIGTERM)
Image PC Routine Line Source
cosmomc 0000000000793DCD Unknown Unknown Unknown
cosmomc 0000000000791C67 Unknown Unknown Unknown
cosmomc 0000000000720864 Unknown Unknown Unknown
cosmomc 0000000000720676 Unknown Unknown Unknown
cosmomc 00000000006A7149 Unknown Unknown Unknown
cosmomc 00000000006AF94E Unknown Unknown Unknown
Unknown 00007F9D8208E330 Unknown Unknown Unknown
cosmomc 0000000000703D65 Unknown Unknown Unknown
cosmomc 000000000042087E Unknown Unknown Unknown
cosmomc 000000000042B08C Unknown Unknown Unknown
cosmomc 0000000000441119 Unknown Unknown Unknown
cosmomc 0000000000502146 Unknown Unknown Unknown
cosmomc 000000000043BD32 Unknown Unknown Unknown
cosmomc 000000000043BCDF Unknown Unknown Unknown
cosmomc 00000000005A3D43 Unknown Unknown Unknown
cosmomc 00000000005A2AA6 Unknown Unknown Unknown
cosmomc 00000000005A88A0 Unknown Unknown Unknown
cosmomc 00000000005A8BBB Unknown Unknown Unknown
cosmomc 00000000005A65F9 Unknown Unknown Unknown
cosmomc 00000000004DF53E Unknown Unknown Unknown
cosmomc 00000000004DE42F Unknown Unknown Unknown
cosmomc 00000000004D3C33 Unknown Unknown Unknown
cosmomc 00000000004B8803 Unknown Unknown Unknown
cosmomc 00000000004B8637 Unknown Unknown Unknown
cosmomc 00000000004E20D2 Unknown Unknown Unknown
cosmomc 00000000004E3CA1 Unknown Unknown Unknown
cosmomc 00000000005B09DC Unknown Unknown Unknown
cosmomc 000000000041009E Unknown Unknown Unknown
libc.so.6 00007F9D81CD7F45 Unknown Unknown Unknown
cosmomc 000000000040FFA9 Unknown Unknown Unknown
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.