Changeset 3344 in Sophya for trunk/Cosmo/SimLSS/cmvdefsurv.cc

Timestamp:

Oct 9, 2007, 10:43:47 AM (18 years ago)

Author:

cmv

Message:

set schechter mmin=1e8 -> 1e7 cmv 09/10/2007

File:

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

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -299 +299 @@
  cout<<"mstar*sch(mstar) = "<<sch(mstar)<<" Msol/Mpc^3/Msol"<<endl;
  int npt = 10000;
- double lnx1=log10(1e+6), lnx2=log10(1e+14), dlnx=(lnx2-lnx1)/npt;
+ double lnx1=log10(1.e+6), lnx2=log10(1.e+13), dlnx=(lnx2-lnx1)/npt;
  double masshimpc3 = IntegrateFuncLog(sch,lnx1,lnx2,0.001,dlnx,10.*dlnx,6);
  cout<<"Mass density: "<<masshimpc3<<" Msol/Mpc^3"<<endl;
@@ -312 +312 @@
  cout<<"\nsch(mstar) = "<<sch(mstar)<<" /Mpc^3/Msol"<<endl;
  cout<<"Galaxy number density:"<<endl;
- for(double x=lnx1; x<10.5; x+=1.) {
+ for(double x=lnx1; x<lnx2-0.5; x+=1.) {
    double n = IntegrateFuncLog(sch,x,lnx2,0.001,dlnx,10.*dlnx,6);
    cout<<"  m>"<<x<<" Msol: "<<n<<" /Mpc^3, "<<n*dvol<<" /pixel, "