Commit 5492ef8f authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

Aplanck au carre poru HLP

parent cdb2e051
......@@ -87,8 +87,8 @@ application testKlass -group=test -s=../src/camel test/testKlass.cc
application testCLHEP -group=test -s=../src/camel/test testCLHEP.cc
application testBAO -group=test -s=../src/camel/test testBAO.cc
application testHiLLiPOP -group=test -s=../src/camel/test testHiLLiPOP.cc
application test_clik -group=test -s=../src/camel test/test_clik.cc
application test_lensing -group=test -s=../src/camel test/test_lensing.cc
application dump_clik -group=test -s=../src/camel test/test_clik.cc
application dump_lensing -group=test -s=../src/camel test/test_lensing.cc
application test_jla -group=test -s=../src/camel test/test_jla.cc
application testPico -group=test -s=../src/camel/pico/exec testPico.cc
......
......@@ -6,7 +6,8 @@ use class v2.4.3
#macro cpp "icpc"
#macro cppflags " -O2 -ip -ansi_alias -align -Wbrief "
macro cppflags " -O2 -pipe -ansi -Woverloaded-virtual "
macro cppflags " -O2 -pipe -ansi -Woverloaded-virtual -Wpedantic -Wunused -Wmaybe-uninitialized -Wno-deprecated -Wno-div-by-zero -Wno-variadic-macros"
#macro cppflags " -O2 -pipe -ansi"
macro cflags " -O2 -fopenmp "
......@@ -86,8 +87,8 @@ application testKlass -group=test -s=../src/camel test/testKlass.cc
application testCLHEP -group=test -s=../src/camel/test testCLHEP.cc
application testBAO -group=test -s=../src/camel/test testBAO.cc
application testHiLLiPOP -group=test -s=../src/camel/test testHiLLiPOP.cc
application test_clik -group=test -s=../src/camel test/test_clik.cc
application test_lensing -group=test -s=../src/camel test/test_lensing.cc
application dump_clik -group=test -s=../src/camel test/test_clik.cc
application dump_lensing -group=test -s=../src/camel test/test_lensing.cc
application test_jla -group=test -s=../src/camel test/test_jla.cc
application testPico -group=test -s=../src/camel/pico/exec testPico.cc
......
......@@ -770,8 +770,8 @@ double HiLLiPOP_3d::computeLikelihood(const vector<unsigned int>& l,vector<doubl
//SP: tmpBeamEigenmodes = pow(1+beta[_XFreq2XSpectra[f][c]]*_beamEigenmodes[_XFreq2XSpectra[f][c]][l],2);
tmpBeamEigenmodes = 1.;
tmpCal = 1+cal[_XSpectra2Maps[_XFreq2XSpectra[f][c]][0]]+cal[_XSpectra2Maps[_XFreq2XSpectra[f][c]][1]];
//SP
tmpCal*=Aplanck;
//SP 15/09/15: au carre
tmpCal*=(Aplanck*Aplanck);
#ifdef PolarizationEfficiency
if( m == 0 ) residual[m][f][l] += tmpWeight*(_ClData[_XFreq2XSpectra[f][c]][modeOffset+l]-tmpCal*tmpBeamEigenmodes*(ClCMB[modeOffset+l]+foregrounds));
if( m == 1 || m == 2 ) {
......
......@@ -500,9 +500,15 @@ ClassEngine::getLensing(const std::vector<unsigned>& lvec, //input
for (size_t i=0;i<lvec.size();i++){
try{
/*
clpp[i]=getCl(ClassEngine::PP,lvec[i])/tr.lcmb_rescale/tr.lcmb_rescale;
cltp[i]=getCl(ClassEngine::TP,lvec[i])/tr.lcmb_rescale;
clep[i]=getCl(ClassEngine::EP,lvec[i])/tr.lcmb_rescale;
*/
clpp[i]=getCl(ClassEngine::PP,lvec[i]);
cltp[i]=getCl(ClassEngine::TP,lvec[i]);
clep[i]=getCl(ClassEngine::EP,lvec[i]);
}
catch(exception &e){
cout << "plantage!" << endl;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment