Changeset 4018 in Sophya for trunk/SophyaLib/BaseTools/randinterf.h

Timestamp:

Sep 21, 2011, 6:21:17 PM (14 years ago)

Author:

cmv

Message:

• introduction du tirage gaussien par la methode de la ziggurat
• ajout argument slim par defaut dans GaussianTail pour proteger d'une reccursion infinie genere par le tirage gaussien ziggurat dans le cas ou GaussianTail applique la methode par rejection (sdev<slim)

cmv, 21/09/2011

File:

• trunk/SophyaLib/BaseTools/randinterf.h (modified) (3 diffs)

trunk/SophyaLib/BaseTools/randinterf.h

@@ -29 +29 @@
   C_Gaussian_PolarBoxMuller = 2,
   C_Gaussian_RatioUnif = 3,
-  C_Gaussian_LevaRatioUnif = 4
+  C_Gaussian_LevaRatioUnif = 4,
+  C_Gaussian_Ziggurat128 = 5
 };
 
@@ -88 +89 @@
   virtual r_8 GaussianRatioUnif(); 
   virtual r_8 GaussianLevaRatioUnif();
+  virtual r_8 GaussianZiggurat128(); 
   /*! \brief Return a random number following a gaussian distribution "sigma", (mean=0)*/ 
   inline r_8 Gaussian(double sigma) {return sigma*Gaussian();}
@@ -94 +96 @@
 
   /*! \brief Return a random number following a gaussian tail distribution for x>sdev */ 
-  virtual r_8 GaussianTail(double sdev);
+  virtual r_8 GaussianTail(double sdev,double slim=1.);
 
   // --- Le tirage sur une distribution de poisson