Changeset 3354 in Sophya

Timestamp:

Oct 22, 2007, 10:29:41 AM (18 years ago)

Author:

cmv

Message:

ajout options -8 sigmar,R , cmv 22/10/2007

Location:

trunk/Cosmo/SimLSS

Files:

• cmvobserv3d.cc (modified) (3 diffs)
• genefluct3d.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -47 +47 @@
      <<" -R schmassdist.ppf : read mass distribution for trials from file"<<endl
      <<"               instead of computing it (ONLY if \"-Q\" option is activated)"<<endl
+     <<" -8 sigmaR,R : normalisation du spectre de puissance, R en Mpc"<<endl
+     <<"               (default sigmaR=1, R=8/h100 Mpc)"<<endl
      <<" -W : write cube in FITS format (complex cube is coded as real cube)"<<endl
      <<" -P : write cube in PPF format"<<endl
@@ -131 +133 @@
 
  char c;
- while((c = getopt(narg,arg,"ha0PWSV2UG:F:x:y:z:s:Z:M:A:T:N:Q:R:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWSV2UG:F:x:y:z:s:Z:M:A:T:N:Q:R:8:")) != -1) {
   int nth = 0;
   switch (c) {
@@ -185 +187 @@
     ////sscanf(optarg,"%lf,%lf,%lf",&lfjy_agn,&lsigma_agn,&powlaw_agn);
     ////break;
+  case '8' :
+    sscanf(optarg,"%lf,%lf",&sigmaR,&R);
+    break;
   case 'V' :
     compvarreal = true;

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -1023 +1023 @@
 // Par definition: SigmaR^2 = <(M-<M>)^2>/<M>^2
 //                 ou M = masse dans sphere de rayon R
+// --- ATTENTION: la variance calculee a une tres grande dispersion
+//  (surtout si le volume du cube est petit). Pour verifier
+//  que le sigmaR calcule par cette methode est en accord avec
+//  le sigmaR en input, il faut faire plusieurs simulations (~100)
+//  et regarder la moyenne des sigmaR reconstruits
 {
   if(lp_>0) cout<<"--- VarianceFrReal R="<<R<<endl;