CAMB Aug06: Scalar + Tensor?
Posted: September 13 2006
Hi,
I may be missing something, but with the last CAMB version, both the standalone and the one in cosmomc, I cannot get simultaneously scalar and tensor Cls output files.
For instance, if in params.ini I set:
get_scalar_cls = T
get_vector_cls = F
get_tensor_cls = T
I get only this output file: test_tensCls.dat.
No sure if this is the correct things to do, but in camb.f90, everything works fine again by uncommenting the following lines and moving them out of all the "if endif" scalar/tensor/vector constructs [in CAMB_GetResult()] and pasting them just before the end of the routine.
!Was there a good reason for doing this??
!CP%Transfer%high_precision = Params%Transfer%high_precision
!CP%WantTransfer = Params%WantTransfer
!CP%WantTensors = Params%WantTensors
!CP%WantVectors = Params%WantTensors
!Params = CP
Anthony, could you confirm that I am not screwing up something else by doing this?
Cheers,
Chris.
PS: these lines are uncommented in camb within cosmomc, but still inside some if endif construct. Therefore, the default Params values seem to not be restored before calling the output cls function in module.f90
I may be missing something, but with the last CAMB version, both the standalone and the one in cosmomc, I cannot get simultaneously scalar and tensor Cls output files.
For instance, if in params.ini I set:
get_scalar_cls = T
get_vector_cls = F
get_tensor_cls = T
I get only this output file: test_tensCls.dat.
No sure if this is the correct things to do, but in camb.f90, everything works fine again by uncommenting the following lines and moving them out of all the "if endif" scalar/tensor/vector constructs [in CAMB_GetResult()] and pasting them just before the end of the routine.
!Was there a good reason for doing this??
!CP%Transfer%high_precision = Params%Transfer%high_precision
!CP%WantTransfer = Params%WantTransfer
!CP%WantTensors = Params%WantTensors
!CP%WantVectors = Params%WantTensors
!Params = CP
Anthony, could you confirm that I am not screwing up something else by doing this?
Cheers,
Chris.
PS: these lines are uncommented in camb within cosmomc, but still inside some if endif construct. Therefore, the default Params values seem to not be restored before calling the output cls function in module.f90