CosmoMC: out of bounds in MatterPowerAt
-
- Posts: 23
- Joined: February 21 2006
- Affiliation: McGill University
CosmoMC: out of bounds in MatterPowerAt
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
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
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: CosmoMC: out of bounds in MatterPowerAt
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).
-
- Posts: 23
- Joined: February 21 2006
- Affiliation: McGill University
CosmoMC: out of bounds in MatterPowerAt
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
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
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: CosmoMC: out of bounds in MatterPowerAt
You should be able to increase num_matter_power.