I am writing a likelihood where I need to call H(z) at different redshifts, which from my understanding can be easily done using the function dtauda provided by CAMB.
However, I am getting some problems. Whatever argument I give to dtauda, it always returns the same value, 463506.174001024 . For instance, dtauda(1.0) = 463506.174001024, dtauda(0.5) = 463506.174001024, and so on.
In my likelihood (called CC), I am using the following modules
Code: Select all
module CC
use settings
use CosmologyTypes
use CosmoTheory
use Calculator_CAMB
use Calculator_Cosmology
use Likelihood_Cosmology
use CAMB
Code: Select all
real(mcp) function CC_LnLike(this, CMB)
Class(CCLikelihood) :: this
Class(CMBParams) CMB
real(mcp):: chisq
real(dl) :: dtauda
external dtauda
Code: Select all
DATAMODULES += $(OUTPUT_DIR)/CC.o
[...]
$(OUTPUT_DIR)/CC.o: $(OUTPUT_DIR)/Calculator_CAMB.o $(OUTPUT_DIR)/Calculator_Cosmology.o $(OUTPUT_DIR)/CosmoTheory.o $(OUTPUT_DIR)/CosmologyTypes.o $(OUTPUT_DIR)/Likelihood_Cosmology.o
Cheers,
Sunny