Changeset 4061 in Sophya for trunk/Cosmo/SimLSS/cmvtstpk.cc

Timestamp:

Apr 17, 2012, 11:21:23 AM (13 years ago)

Author:

cmv

Message:

cmv, 17/04/2012

File:

• trunk/Cosmo/SimLSS/cmvtstpk.cc (modified) (3 diffs)

trunk/Cosmo/SimLSS/cmvtstpk.cc

@@ -80 +80 @@
  if(optind<narg) zval = atof(arg[optind]);
 
- cout<<"h100="<<h100<<" Om0="<<Om0<<" Ob0="<<Ob0<<" Ol0="<<Ol0<<" w0="<<w0<<endl;
+ cout<<"h100="<<h100<<" Om0="<<Om0<<" Ob0="<<Ob0<<" Ocdm="<<Om0-Ob0<<" Ol0="<<Ol0<<" w0="<<w0<<endl;
  cout<<"lkmin="<<lkmin<<" lkmax="<<lkmax<<" npt="<<npt<<" dlk="<<dlk<<endl;
  cout<<"scale="<<scale<<endl;
@@ -91 +91 @@
 
  TransfertEisenstein tf(h100,Om0-Ob0,Ob0,T_CMB_Par,false);
- cout<<"kpeak="<<tf.KPeak()<<endl;
+   double tf0 = tf(0);
+   cout<<"kpeak="<<tf.KPeak()<<endl;
  TransfertEisenstein tfnosc2(tf); tfnosc2.SetNoOscEnv(2);
+   double tfnosc20 = tfnosc2(0);
  TransfertEisenstein tfnosc1(tf); tfnosc1.SetNoOscEnv(1);
+   double tfnosc10 = tfnosc1(0);
  TransfertEisenstein tfnob(h100,Om0,0.,T_CMB_Par,true);
+   double tfnob0 = tfnob(0);
+ cout<<"tf(0)="<<tf0<<" tfnosc2(0)="<<tfnosc20<<" tfnosc1(0)="<<tfnosc10<<" tfnob(0)="<<tfnob0<<endl;
 
  GrowthEisenstein d1(Om0,Ol0);
- cout<<"GrowthFactor: "<<d1(zval)<<endl;
+ cout<<"GrowthFactor: D1(0)="<<d1(0.)<<" D1("<<zval<<") = "<<d1(zval)<<endl;
 
  PkSpecCalc pk0(pkini,tf,d1,0.);
@@ -173 +178 @@
    fprintf(fdata,"# z= %g : k(Mpc^-1) pkini, tf pk0 pkz, tfnosc2 pk0nosc2 pkznosc2, ",zval);
    fprintf(fdata,"tfnosc1 pk0nosc1 pkznosc1, tfnob pk0nob pkznob\n");
+   int n = 0;
    for(double lk=lkmin;lk<lkmax+dlk/2.;lk+=dlk) {
      double k = pow(10.,lk);
-     fprintf(fdata,"%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n",
-             k,pkini(k),
+     fprintf(fdata,"%d %e %e %e %e %e %e %e %e %e %e %e %e %e %e\n",
+             n,k,pkini(k),
              tf(k),pk0(k),pkz(k,zval),
              tfnosc2(k),pk0nosc2(k),pkznosc2(k,zval),
              tfnosc1(k),pk0nosc1(k),pkznosc1(k,zval),
              tfnob(k),pk0nob(k),pkznob(k,zval));
+     n++;
    }
    fclose(fdata);