Tutorial in CosmoMC

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » March 03 2017

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

Erick Almaraz
Posts: 25
Joined: June 10 2012
Affiliation: University of Mexico

Tutorial in CosmoMC

Post by Erick Almaraz » March 03 2017

Check out Antonio Cuesta's slides

Just google Antonio Cuesta CosmoMC

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » March 04 2017

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 !!

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » March 04 2017

One more thing,
what is the difference between disttest.ini and distparams.ini

Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Tutorial in CosmoMC

Post by Antony Lewis » March 04 2017

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

http://cosmologist.info/cosmomc/readme.html

Bob Zou
Posts: 5
Joined: July 21 2016
Affiliation: Beijing Institute of Technology

Tutorial in CosmoMC

Post by Bob Zou » March 15 2017

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.

Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Tutorial in CosmoMC

Post by Antony Lewis » March 15 2017

It's in CAMB (subroutines.f90).

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » April 03 2017

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

Thanks

Sofie Marie Koksbang
Posts: 9
Joined: February 27 2017
Affiliation: Aarhus University

Tutorial in CosmoMC

Post by Sofie Marie Koksbang » April 03 2017

Hi,

A description of how to add theory parameters is given e.g. at http://cosmocoffee.info/viewtopic.php?t ... c&start=15

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » April 05 2017

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

Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Tutorial in CosmoMC

Post by Antony Lewis » April 05 2017

Probably also have MPI_Check_Limit_Converge =T?

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » April 05 2017

yep MPI_Check_Limit_Converge =T but what does this do?

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

Tutorial in CosmoMC

Post by R Purnalingam » April 05 2017

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

Post Reply