Seg fault in May version of cosmomc
Posted: May 14 2008
The May release of Cosmomc (Union SN sample, etc.) seg faults for me
out of the box (i.e., running ./cosmomc params.ini with the supplied versions
does not work).
The issue seems to be one of variable initialization. The seg fault
is when it tries to write chain results in WriteParams (in params_CMB.f90),
as nuisance_params_used has a very large random value.
Tracing this through the code, I note that
ReadDataSet in cmbdata.F90 does -not- initialize aset%nuisance_parameters
for the WMAP case. Adding a line doing so (i.e.,
aset%nuisance_parameters=0
inside the IF (aname == 'WMAP' ...
statement) seems to fix things, although I'm a little leery
of this solution since I don't understand exactly how the code
is trying to handle the WMAP data.
out of the box (i.e., running ./cosmomc params.ini with the supplied versions
does not work).
The issue seems to be one of variable initialization. The seg fault
is when it tries to write chain results in WriteParams (in params_CMB.f90),
as nuisance_params_used has a very large random value.
Tracing this through the code, I note that
ReadDataSet in cmbdata.F90 does -not- initialize aset%nuisance_parameters
for the WMAP case. Adding a line doing so (i.e.,
aset%nuisance_parameters=0
inside the IF (aname == 'WMAP' ...
statement) seems to fix things, although I'm a little leery
of this solution since I don't understand exactly how the code
is trying to handle the WMAP data.