def getClsAxionCamb(out_pre,spec,templateIni,params,AccuracyBoost=False,seed=0):
option = 0 # for regular CAMB call
CAMB = cambInterface(out_pre,templateIni,cambRoot=os.environ['AXIONCAMB_DIR'],option=option,seed=seed)
for key in params:
CAMB.setParam(key,params[key])
CAMB.call(suppress=False)
if spec == 'lensed':
Cls = CAMB.getCls(lensed=True)
else:
Cls = CAMB.getCls(lensed=False)
CAMB.done()
return Cls
def getClsCambRed(out_pre, templateIni, params, AccuracyBoost=False, savefid=False, seed=0):
option = 1
cRed = cambInterface(out_pre, templateIni, cambRoot=os.environ["CAMBREDNN_DIR"], option=option, seed=seed)
for key in params:
cRed.setParam(key, params[key])
# sys.exit()
cRed.call(suppress=False)
# print 'call cambRed done!'
Cls = cRed.getCls()
if savefid:
os.system("cp $CAMBREDNN_DIR/" + out_pre + "_params.ini $FISHER_DIR/output/" + out_pre + "_galaxy_fid.csv")
cRed.done()
return Cls
def getClsCamb(out_pre,
spec,
templateIni,
params,
cambRoot=os.environ['AXIONCAMB_DIR'],
AccuracyBoost=False,
seed=0):
option = 0 # for regular CAMB call
CAMB = cambInterface(out_pre,
templateIni,
cambRoot=cambRoot,
option=option,
seed=seed)
for key in params:
CAMB.setParam(key, params[key])
CAMB.call(suppress=False)
if spec == 'lensed':
Cls = CAMB.getCls(lensed=True)
else:
Cls = CAMB.getCls(lensed=False)
CAMB.done()
return Cls
def getClsCambRed(out_pre,
templateIni,
params,
AccuracyBoost=False,
savefid=False,
seed=0):
option = 1
cRed = cambInterface(out_pre,
templateIni,
cambRoot=os.environ['CAMBREDNN_DIR'],
option=option,
seed=seed)
for key in params:
cRed.setParam(key, params[key])
#sys.exit()
cRed.call(suppress=False)
#print 'call cambRed done!'
Cls = cRed.getCls()
if savefid:
os.system('cp $CAMBREDNN_DIR/' + out_pre +
'_params.ini $FISHER_DIR/output/' + out_pre +
'_galaxy_fid.csv')
cRed.done()
return Cls
