Changeset 3768 in Sophya for trunk/Cosmo/SimLSS/cmvtgrowth.cc

Timestamp:

May 3, 2010, 4:08:34 PM (15 years ago)

Author:

cmv

Message:

• refonte du code pour creer uniquement des conditions initiales
• introduction du tirage des vitesse LOS pour les redshift-distortion

cmv 03/05/2010

File:

• trunk/Cosmo/SimLSS/cmvtgrowth.cc (modified) (4 diffs)

trunk/Cosmo/SimLSS/cmvtgrowth.cc

@@ -14 +14 @@
 
 void usage(void);
-void usage(void) {cout<<"cmvtgrowth z1,z2,dz [Omatter0,Lambda0]"<<endl;}
-
-void tstprint(CosmoCalc& univ,double z1,double z2,double dz);
-void tstspeed(CosmoCalc& univ,double z1,double z2,double dz);
-void tstntuple(CosmoCalc& univ,double z1,double z2,double dz);
+void usage(void) {cout<<"cmvtgrowth z1,z2,nz [Omatter0,Lambda0]"<<endl;}
 
 int main(int narg,char *arg[])
 {
 
- double z1=0., z2=10., dz=1.;
- if(narg>1) sscanf(arg[1],"%lf,%lf,%lf",&z1,&z2,&dz);
- cout<<"z1="<<z1<<"  z2="<<z2<<"  dz="<<dz<<endl;
+ double z1=0., z2=10.;
+ int nz = 100;;
+ if(narg>1) sscanf(arg[1],"%lf,%lf,%d",&z1,&z2,&nz);
+ if(nz<=0) nz = 100;
+ double dz = (z2-z1)/nz;
+ cout<<"z1="<<z1<<"  z2="<<z2<<"  nz="<<nz<<"  dz="<<dz<<endl;
 
  double om0=0.267804, ol0=0.73;
@@ -34 +33 @@
  cout<<"D1(z=0) = "<<growth(0.)<<endl;
 
- const int n = 2;
- const char *vname[n] = {"z","d1"};
+ const int n = 3;
+ const char *vname[n] = {"z","d1","d1dz"};
  NTuple nt(n,vname);
  double xnt[n];
@@ -41 +40 @@
    xnt[0]  = z;
    xnt[1]  = growth(z);
+   xnt[2]  = growth.DsDz(z,dz/2.);
    nt.Fill(xnt);
  }
@@ -59 +59 @@
 set cut z<5
 
+# --- growth
 n/plot nt.d1%z $cut ! "nsta connectpoints"
 n/plot nt.1./(1.+z)%z $cut ! "nsta connectpoints red same"
+
+# --- growth'/growth
+n/plot nt.d1dz/d1%z $cut ! "nsta connectpoints"
+n/plot nt.-1./(1.+z)%z $cut ! "nsta connectpoints red same"
+
+# --- d(growth)/dz
+zmin = 0.
+zmax = 10.
+set npt ${dv.size}
+dd = ($zmax-$zmin)/(${npt}-1.)
+exptovec dv nt d1
+c++exec for(int i=1;i<dv.Size(); i++) dv(i-1) = dv(i)-dv(i-1);
+
+n/plot nt.-1./pow(1.+z,2.)%z $cut ! "nsta connectpoints red"
+n/plot nt.d1dz%z $cut ! "nsta connectpoints same"
+n/plot dv.val/${dd}%${zmin}+(n+1)*$dd ! ! "nsta plusmarker5 green same"
+
 */