Changeset 3249 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

May 16, 2007, 4:39:35 PM (18 years ago)

Author:

cmv

Message:

print cmv 16/05/2007

File:

• trunk/Cosmo/SimLSS/cmvtstsch.cc (modified) (2 diffs)

trunk/Cosmo/SimLSS/cmvtstsch.cc

@@ -91 +91 @@
  cout<<"sch(mstar) = "<<sch(mstar)<<" /Mpc^3/Msol"<<endl;
  sch.SetOutValue(0);
- for(int i=0;i<=5;i++) {
-   double num = IntegrateFuncLog(sch,lnx1-i,lnx2+i,perc,dlxinc,dlxmax,glorder);
-   cout<<"["<<lnx1-i<<","<<lnx2+i<<"] integrated number = "<<num<<" Mpc^-3"<<endl;
+ for(double lnx=lnx1-5.;lnx<=lnx2-1.;lnx+=1.) {
+   double num = IntegrateFuncLog(sch,lnx,lnx2,perc,dlxinc,dlxmax,glorder);
+   cout<<"["<<lnx<<","<<lnx2<<"] integrated number = "<<num<<" Mpc^-3"<<endl;
  }
  Histo hdndm(lnx1,lnx2,npt); hdndm.ReCenterBin();
@@ -102 +102 @@
  sch.SetOutValue(1);
  cout<<"mstar*sch(mstar) = "<<sch(mstar)<<" Msol/Mpc^3/Msol"<<endl;
- for(int i=0;i<=5;i++) {
-   double sum = IntegrateFuncLog(sch,lnx1-i,lnx2+i,perc,dlxinc,dlxmax,glorder);
-   cout<<"["<<lnx1-i<<","<<lnx2+i<<"] integrated mass = "<<sum<<" Msol.Mpc^-3"<<endl;
+ for(double lnx=lnx1-5.;lnx<=lnx2-1.;lnx+=1.) {
+   double sum = IntegrateFuncLog(sch,lnx,lnx2,perc,dlxinc,dlxmax,glorder);
+   cout<<"["<<lnx<<","<<lnx2<<"] integrated mass = "<<sum<<" Msol.Mpc^-3"<<endl;
  }
  Histo hmdndm(lnx1,lnx2,npt); hmdndm.ReCenterBin();