Changeset 3599 in Sophya for trunk/SophyaLib

Timestamp:

Apr 22, 2009, 12:56:23 AM (16 years ago)

Author:

cmv

Message:

redimensionnement du buffer RandomGenerator, cmv 22/04/2009

Location:

trunk/SophyaLib/BaseTools

Files:

• stsrand.cc (modified) (1 diff)
• stsrand.h (modified) (3 diffs)

trunk/SophyaLib/BaseTools/stsrand.cc

@@ -57 +57 @@
   }
 }
- 
+
 
 RandomGenerator::~RandomGenerator(void) 
 {
   // rien a faire 
+}
+
+void RandomGenerator::SetBuffSize(size_t n)
+// redimensionnement du buffer
+{
+  if(fg_nothrsafe) return;
+  if (n < 1) n = 1024;
+  rseq_.ReSize(n, false);
+  idx_ = n;
 }
 

trunk/SophyaLib/BaseTools/stsrand.h

@@ -22 +22 @@
   RandomGenerator(RandomGenerator const & rg);
 
+  void SetBuffSize(size_t n);
+
   virtual ~RandomGenerator();
 
@@ -39 +41 @@
   inline r_8 Flatpm1() {return 2.*Next()-1.;}
 
-  /*! \brief Return a random number following a gaussian distribution (sigma=1, mean=0*/ 
+  /*! \brief Return a random number following a gaussian distribution (sigma=1, mean=0)*/ 
   virtual r_8 Gaussian(); 
 
   /*! \brief Return a random number following a gaussian distribution "sigma", (mean=0)*/ 
-  virtual r_8 Gaussian(double sigma) 
-    { return sigma*Gaussian(); }
+  virtual r_8 Gaussian(double sigma)  {return sigma*Gaussian();}
  /*! \brief Return a random number following a gaussian distribution with mean=mu, and sigma */ 
-  virtual r_8 Gaussian(double sigma,double mu) 
-    { return sigma*Gaussian()+mu; }
+  virtual r_8 Gaussian(double sigma,double mu) {return sigma*Gaussian()+mu;}
 
   /*! \brief alias for Gaussian() */ 
-  inline r_8 NorRand()
-    { return Gaussian(); }
+  inline r_8 NorRand() {return Gaussian();}
 
   /*! \brief Return a random number following a poisson distribution with mean mu */ 
@@ -57 +56 @@
   
   /*! \brief Return a random number following a poisson distribution with mean mu */ 
-  inline uint_8 Poisson(double mu)
-    { return Poisson(mu, -1.); }
+  inline uint_8 Poisson(double mu) {return Poisson(mu, -1.);}
 
    /*! \brief Return a random number following a poisson distribution with mean mu */ 
-  inline uint_8  PoissRandLimit(double mu,double mumax=10.)
-    { return Poisson(mu, mumax); }
+  inline uint_8  PoissRandLimit(double mu,double mumax=10.) {return Poisson(mu, mumax);}
+
+  /*! \brief Returns a random complex number such that real and imaginary parts are gaussians with variance sig^2 */
+  inline complex< r_8 > ComplexGaussRan(void)
+         {return complex< r_8 >(NorRand(),NorRand());}
+  inline complex< r_8 > ComplexGaussRan(double sig)
+         {return complex< r_8 >(sig*NorRand(),sig*NorRand());}
+  /*! \brief Returns the module of a random complex number generated by ComplexGaussRan */
+  inline double ModComplexGaussRan(double sig=1.)
+         {double r=-log(1.-Flat01()); return sig*sqrt(2.*r);}
+
  
   //  Pour la gestion de persistance PPF