CosmoMC: out of bounds in MatterPowerAt

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Loison Hoi
Posts: 23
Joined: February 21 2006
Affiliation: McGill University

CosmoMC: out of bounds in MatterPowerAt

Post by Loison Hoi » November 03 2006

Hi All,

I have a problem when I run chains with Oct (20th) 2006 CosmoMC. I used WMAP3 + (SDSS or LUAQS Lyman-alpha). Sometimes the chains were stopped and said:

k/h out of bounds in MatterPowerAt ( 3.842072 )

And I found the corresponding the statement in cmbtypes.f90:

x = log(kh/matter_power_minkh) / matter_power_dlnkh
if (x < 0 .or. x >= num_matter_power-1) then
write (*,*) ' k/h out of bounds in MatterPowerAt (',kh,')'
stop
end if

Since kh = 3.842072 in my case, and matter_power_minkh = 0.999e-4 and matter_power_dlnkh = 0.143911568 by default, so x = 73.36 >= num_matter_power - 1 = 74 -1 = 73. Therefore, the code runs into the stop statement. Does anyone know how to avoid this problem? Thanks a lot!

Loison Hoi
2 Nov 2006

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

Re: CosmoMC: out of bounds in MatterPowerAt

Post by Antony Lewis » November 05 2006

I would guess this is a problem only with ly-alpha data since this is a rather small scale? I think the default ranges should be OK for the the defauly lya data (SDSS).

Loison Hoi
Posts: 23
Joined: February 21 2006
Affiliation: McGill University

CosmoMC: out of bounds in MatterPowerAt

Post by Loison Hoi » November 07 2006

Hi Antony,

Can I avoid this problem by just modifying one number? For example, let num_matter_power be 100 or more, but keeping others unchanged. Or are there relations that the following numbers must satisfy?

num_matter_power = 74 !number of points computed in matter power spectrum
matter_power_minkh = 1e-4 !minimum value of k/h to store
matter_power_dlnkh = 0.143911568 !log spacing in k/h
matter_power_maxz = 0. !6.0
matter_power_lnzsteps = 1 !20

Thanks a lot!

Loison Hoi
6 Nov 2006

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

Re: CosmoMC: out of bounds in MatterPowerAt

Post by Antony Lewis » November 07 2006

You should be able to increase num_matter_power.

Post Reply