Changeset 3617 in Sophya

Timestamp:

May 2, 2009, 11:43:17 AM (16 years ago)

Author:

cmv

Message:

chgt nom Gaussian -> GaussianRand etc.. pour ambiguite, cmv 02/05/2009

Location:

trunk

Files:

• Cosmo/SimLSS/cmvtstagn.cc (modified) (1 diff)
• Cosmo/SimLSS/genefluct3d.cc (modified) (2 diffs)
• SophyaLib/SysTools/cexpre.cc (modified) (1 diff)
• SophyaLib/SysTools/rpneval.cc (modified) (1 diff)
• SophyaProg/Tests/lpk.cc (modified) (2 diffs)
• SophyaProg/Tests/tfft.cc (modified) (1 diff)
• SophyaProg/Tests/tmtrnd.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvtstagn.cc

@@ -103 +103 @@
    for(long i=0;i<npix;i++) {
      if(i%l==0) cout<<"...trial: "<<i<<endl;
-     uint_8 n = Poisson(mu,10.);
+     uint_8 n = PoissonRand(mu,10.);
      hnobj.Add((double)n);
      sumn += n; sumn2 += n*n;

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -694 +694 @@
        // Explication de la division par 2: voir perandom.cc
        // ou egalement Coles & Lucchin, Cosmology formula 13.7.2 p279
-       T_(l,j,i) = ComplexGaussian(sqrt(pk/2.));
+       T_(l,j,i) = ComplexGaussianRand(sqrt(pk/2.));
      }
    }
@@ -1347 +1347 @@
    double v = data_[ip];
    if(v>0.) {
-     uint_8 dn = Poisson(v,10.);
+     uint_8 dn = PoissonRand(v,10.);
      data_[ip] = (double)dn;
      sum += (double)dn;

trunk/SophyaLib/SysTools/cexpre.cc

@@ -220 +220 @@
 static double _CE_rand01() { return drand01(); }
 static double _CE_randpm1() { return drandpm1(); }
-static double _CE_gaurand() { return Gaussian(1.,0.); }
+static double _CE_gaurand() { return GaussianRand(1.,0.); }
 //---------------------------------------------------------
 CE_FuncExp::CE_FuncExp(string const & func)

trunk/SophyaLib/SysTools/rpneval.cc

@@ -238 +238 @@
     }
     else if (args[k] == "gaurand") { 
-      double rnd = Gaussian(1.,0.);
+      double rnd = GaussianRand(1.,0.);
       rpnstack_.push(rnd);
     }

trunk/SophyaProg/Tests/lpk.cc

@@ -92 +92 @@
   Matrix a(n,n);
   for(i=0; i<n; i++)
-    for(j=0; j<n; j++)  a(j,i) = Gaussian(1.,0.);
+    for(j=0; j<n; j++)  a(j,i) = GaussianRand(1.,0.);
   
   Vector x(n), b;
@@ -98 +98 @@
   if (prtlev > 0) 
     cout << " ------------ Vector X = \n " << x << "\n" << endl; 
-  for(i=0; i<n; i++) x(i) = Gaussian(1.5,2.);
+  for(i=0; i<n; i++) x(i) = GaussianRand(1.5,2.);
   b = a*x;
 

trunk/SophyaProg/Tests/tfft.cc

@@ -70 +70 @@
   if (inp_typ_random) 
     for (i=0; i<num ; i++){
-      ino[i] = in[i] = Gaussian(1.,0.);
+      ino[i] = in[i] = GaussianRand(1.,0.);
       inc[i] = complex<T> (in[i], 0.);
     }

trunk/SophyaProg/Tests/tmtrnd.cc

@@ -83 +83 @@
     Timer tm("f0_tmtrnd-Gaussian/drand01");
     if (fgg) 
-      for(sa_size_t k=0; k<VSZ; k++)  DATA(k) = Gaussian(1.,0.);
+      for(sa_size_t k=0; k<VSZ; k++)  DATA(k) = GaussianRand(1.,0.);
     else 
       for(sa_size_t k=0; k<VSZ; k++)  DATA(k) = drand01();