### adding parameters, empty chains

Posted:

**January 16 2009**Dear all,

I have added three new slow parameters following the information in the Readme file of the COSMOMC and in several similar topics here. I summarize the changes I have made at the end of this message. I tried running COSMOMC with both options; asking for estimating a new covariance matrix and also with the old matrix filled with zeros for the added entries. In any case I got the following output on the screen (this is an example using 6 processors)

===

Number of MPI processes: 6

Random seeds: 2354, 15909 rand_inst: 1

Random seeds: 12426, 30014 rand_inst: 5

Random seeds: 30414, 24716 rand_inst: 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Random seeds: 17023, 772 rand_inst: 6

Random seeds: 5732, 22170 rand_inst: 4

Random seeds: 8613, 796 rand_inst: 2

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

===

and nothing else, even after a couple of days. I have checked the chains and all of they are empty as well as the log files.

I would appreciate any help.

Best,

César.

====

Summary of modifications to COSMOMC

====

The new parameters are defined in modules.f90 of CAMB as

type CAMBparams

real(dl) :: qVpar1, qVpar2, qVpar3

---

In COSMOMC they are added with

In settings.f90

integer, parameter :: num_hard =10

---

In cmbtypes.f90

Type CMBParams

real Q1, Q2, Q3

---

In params_CMB.f90

subroutine SetForH(Params,CMB,H0)

CMB%Q1 = Params(8)

CMB%Q2 = Params(9)

CMB%Q3 = Params(10)

subroutine CMBParamsToParams(CMB, Params)

Params(8) = CMB%Q1

Params(9) = CMB%Q2

Params(10) = CMB%Q3

---

In CMB_Cls_simple.f90

subroutine CMBToCAMB(CMB,P)

P%qVpar1 = CMB%Q1

P%qVpar2 = CMB%Q2

P%qVpar3 = CMB%Q3

---

In params.ini

param8 = 10 100 1000 10 10

param9 = 0.01 0.001 0.1 0.01 0.01

param10 = 0 0 0 0 0

#n_s

param11 = 0.95 0.5 1.5 0.02 0.01

#n_t

param12 = 0 0 0 0 0

#n_run

param13 = 0 0 0 0 0

#log[10^10 A_s]

param14 = 3 2.7 4 0.01 0.01

#amp_ratio

param15 = 0 0 0 0 0

#SZ amplitude, as in WMAP analysis

param16 = 1 0 2 0.4 0.4

I have added three new slow parameters following the information in the Readme file of the COSMOMC and in several similar topics here. I summarize the changes I have made at the end of this message. I tried running COSMOMC with both options; asking for estimating a new covariance matrix and also with the old matrix filled with zeros for the added entries. In any case I got the following output on the screen (this is an example using 6 processors)

===

Number of MPI processes: 6

Random seeds: 2354, 15909 rand_inst: 1

Random seeds: 12426, 30014 rand_inst: 5

Random seeds: 30414, 24716 rand_inst: 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Random seeds: 17023, 772 rand_inst: 6

Random seeds: 5732, 22170 rand_inst: 4

Random seeds: 8613, 796 rand_inst: 2

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

Computing tensors: F

Doing CMB lensing: T

lmax = 2100

Number of C_ls = 3

===

and nothing else, even after a couple of days. I have checked the chains and all of they are empty as well as the log files.

I would appreciate any help.

Best,

César.

====

Summary of modifications to COSMOMC

====

The new parameters are defined in modules.f90 of CAMB as

type CAMBparams

real(dl) :: qVpar1, qVpar2, qVpar3

---

In COSMOMC they are added with

In settings.f90

integer, parameter :: num_hard =10

---

In cmbtypes.f90

Type CMBParams

real Q1, Q2, Q3

---

In params_CMB.f90

subroutine SetForH(Params,CMB,H0)

CMB%Q1 = Params(8)

CMB%Q2 = Params(9)

CMB%Q3 = Params(10)

subroutine CMBParamsToParams(CMB, Params)

Params(8) = CMB%Q1

Params(9) = CMB%Q2

Params(10) = CMB%Q3

---

In CMB_Cls_simple.f90

subroutine CMBToCAMB(CMB,P)

P%qVpar1 = CMB%Q1

P%qVpar2 = CMB%Q2

P%qVpar3 = CMB%Q3

---

In params.ini

param8 = 10 100 1000 10 10

param9 = 0.01 0.001 0.1 0.01 0.01

param10 = 0 0 0 0 0

#n_s

param11 = 0.95 0.5 1.5 0.02 0.01

#n_t

param12 = 0 0 0 0 0

#n_run

param13 = 0 0 0 0 0

#log[10^10 A_s]

param14 = 3 2.7 4 0.01 0.01

#amp_ratio

param15 = 0 0 0 0 0

#SZ amplitude, as in WMAP analysis

param16 = 1 0 2 0.4 0.4