I am trying to run the generic sampler on the latest version of cosmomc.
The last version I ran used 'num_hard' in settings.f90 to set the number of parameters. However, I see this is now redundant.
When I try to run the program I get the following error:
Code: Select all
file_root:FR/chains
Random seeds: 9367, 2262 rand_inst: 1
Change the Parameterization type implementation to use Parameter names
Computing tensors: F
Doing CMB lensing: F
Doing non-linear Pk: 0
lmax = 6500
lmax_computed_cl = 10
max_eta_k = 20.0000000000000
transfer kmax = 0.800000011920929
Number of C_ls = 4
parameter ranges not found: 1
Code: Select all
param1 = 0. -10 10. 0.01 0.01
param2 = 0. -10 10. 0.01 0.01
param3 = 0. -10 10. 0.01 0.01
param4 = 0. -10 10. 0.01 0.01
param5 = 0. -10 10. 0.01 0.01
use_CMB = F
use_HST = F
use_mpk = F
use_clusters = F
use_BBN = F
use_Age_Tophat_Prior = F
use_SN = F
use_min_zre = 0
Code: Select all
function GenericLikelihoodFunction(Params)
type(ParamSet) Params
real(mcp) :: GenericLikelihoodFunction
real(mcp), allocatable, save :: covInv(:,:)
real(mcp) X(num_params_used)
alpha1 = Params%P(1)
alpha2 = Params%P(2)
beta1 = Params%P(3)
beta2 = Params%P(4)
beta3 = Params%P(5)
GenericLikelihoodFunction = sfr_like()
call MpiStop('GenericLikelihoodFunction: need to write this function!')
end function GenericLikelihoodFunction
So basically, its not clear to me how you tell cosmomc how many paramaters you are using in generic mode and thus why I get error message:
parameter ranges not found: 1
I presume its to do with .paranames file but not sure how you link it. Any help would be appreciated!
cheers