CosmoMC: estimate_propose_matrix with new parameters
Posted: September 16 2005
Hi,
I am trying to use the above option set to true in CosmoMC to find a new proposal distribution for a model in which I have added several fast parameters. In one case I have a pretty good idea of the marginalized standard deviation of the parameters, in the others, not so much. But in all cases, the code crashes after a while with errors that look like what is pasted below.
I can always try the MPI methods to learn the propose matrix that are recommended in previous posts in this forum. However, I would really like to get the code working as it should in a serial environment before migrating to a cluster. The error looks like there is some weird memory leak happening, so I would really like to know its cause, even if it ends up that I can't use this option for some reason. I would be really grateful if anyone can provide any suggestions.
Thanks a lot!
Hiranya
For initial lower limit on step length (lambda1), 0.0000000000E+00 , gradient.dirn is -483.36
40747
For initial upper estimate (lambda2), 0.8407531306E-02 , gradient.dirn is 57083.14844
Locating the minimum in this direction...
Finding the gradient numerically...
gradient wrt 1 th element is: 58.58262634
gradient wrt 2 th element is: -70.12072754
gradient wrt 3 th element is: -6.720314026
gradient wrt 4 th element is: -149.6531372
gradient wrt 5 th element is: -698.9159546
gradient wrt 6 th element is: 1242.675415
gradient wrt 7 th element is: 380.3808289
after 1 bisections, lambda2,df2: 0.4203765653E-02 30061.18359
Finding the gradient numerically...
gradient wrt 1 th element is: 28.08916283
gradient wrt 2 th element is: -33.91596985
gradient wrt 3 th element is: -4.793701172
gradient wrt 4 th element is: -75.75560760
gradient wrt 5 th element is: -350.4357605
gradient wrt 6 th element is: 631.8175049
gradient wrt 7 th element is: 185.3844452
after 2 bisections, lambda2,df2: 0.2101882827E-02 15189.93262
Finding the gradient numerically...
cosmomc(15553) malloc: *** vm_allocate(size=897024) failed (error code=3)
cosmomc(15553) malloc: *** error: can't allocate region
cosmomc(15553) malloc: *** set a breakpoint in szone_error to debug
1525-108 Error encountered while attempting to allocate a data object. The program will stop.
I am trying to use the above option set to true in CosmoMC to find a new proposal distribution for a model in which I have added several fast parameters. In one case I have a pretty good idea of the marginalized standard deviation of the parameters, in the others, not so much. But in all cases, the code crashes after a while with errors that look like what is pasted below.
I can always try the MPI methods to learn the propose matrix that are recommended in previous posts in this forum. However, I would really like to get the code working as it should in a serial environment before migrating to a cluster. The error looks like there is some weird memory leak happening, so I would really like to know its cause, even if it ends up that I can't use this option for some reason. I would be really grateful if anyone can provide any suggestions.
Thanks a lot!
Hiranya
For initial lower limit on step length (lambda1), 0.0000000000E+00 , gradient.dirn is -483.36
40747
For initial upper estimate (lambda2), 0.8407531306E-02 , gradient.dirn is 57083.14844
Locating the minimum in this direction...
Finding the gradient numerically...
gradient wrt 1 th element is: 58.58262634
gradient wrt 2 th element is: -70.12072754
gradient wrt 3 th element is: -6.720314026
gradient wrt 4 th element is: -149.6531372
gradient wrt 5 th element is: -698.9159546
gradient wrt 6 th element is: 1242.675415
gradient wrt 7 th element is: 380.3808289
after 1 bisections, lambda2,df2: 0.4203765653E-02 30061.18359
Finding the gradient numerically...
gradient wrt 1 th element is: 28.08916283
gradient wrt 2 th element is: -33.91596985
gradient wrt 3 th element is: -4.793701172
gradient wrt 4 th element is: -75.75560760
gradient wrt 5 th element is: -350.4357605
gradient wrt 6 th element is: 631.8175049
gradient wrt 7 th element is: 185.3844452
after 2 bisections, lambda2,df2: 0.2101882827E-02 15189.93262
Finding the gradient numerically...
cosmomc(15553) malloc: *** vm_allocate(size=897024) failed (error code=3)
cosmomc(15553) malloc: *** error: can't allocate region
cosmomc(15553) malloc: *** set a breakpoint in szone_error to debug
1525-108 Error encountered while attempting to allocate a data object. The program will stop.