Changeset 3604 in Sophya for trunk/SophyaLib/BaseTools/randinterf.cc

Timestamp:

Apr 29, 2009, 12:14:04 PM (16 years ago)

Author:

ansari

Message:

Modifs,petites extensions + numero de version/initialisation de l'instance globale de RandomGenerator, Reza 29/04/2009

File:

• trunk/SophyaLib/BaseTools/randinterf.cc (modified) (5 diffs)

trunk/SophyaLib/BaseTools/randinterf.cc

@@ -14 +14 @@
 //-------------------------------------------------------------------------------
 // ------ Definition d'interface des classes de generateurs de nombres aleatoires
+/*!
+   \class RandomGeneratorInterface
+   \ingroup BaseTools
+   \brief Base class for random number generators
+
+   This class defines the interface for random number generator classes and 
+   implements the generation of some specific distributions (Gaussian, Poisson ...) 
+   through generation of random number with a flat distribution in the range [0,1[. 
+
+   The sub classes inheriting from this class should implement the Next() method.
+
+   This base class manages also a global instance of a default generator.
+
+   \sa frand01 drand01 frandpm1 drandpm1 
+   \sa GauRnd PoissRand
+
+*/
+
+
+RandomGeneratorInterface* RandomGeneratorInterface::gl_rndgen_p = NULL;
+
+/*! 
+   \brief: static method to set or change the intance of the global Random Generator object
+
+   This method should be called during initialization, before any call to global 
+   functions for random number generation. The rgp object should be created using new.
+*/
+void RandomGeneratorInterface::SetGlobalRandGenP(RandomGeneratorInterface* rgp)
+{
+  if (rgp == NULL) return;
+  if (gl_rndgen_p) delete gl_rndgen_p;
+  gl_rndgen_p = rgp;
+  return;
+}
 
 RandomGeneratorInterface::RandomGeneratorInterface()
-  : usegaussian_(0), usepoisson_(0), useexpo_(0)
-{
+{
+  SelectGaussianAlgo();
+  SelectPoissonAlgo();
+  SelectExponentialAlgo();
 }
  
@@ -26 +62 @@
 }
 
-void RandomGeneratorInterface::UseGaussian(uint_2 type)
-{
-  usegaussian_ = type;
-}
-
-void RandomGeneratorInterface::UseExpo(uint_2 type)
-{
-  useexpo_ = type;
-}
-
-void RandomGeneratorInterface::UsePoisson(uint_2 type)
-{
-  usepoisson_ = type;
-}
 
 /////////////////////////////////////////////////////////////////////////
@@ -150 +172 @@
 }
 
-
-/////////////////////////////////////////////////////////////////////////
-/////////////////////////////////////////////////////////////////////////
-/////////////////////////////////////////////////////////////////////////
+/////////////////////////////////////////////////////////////////////////
+/////////////////////////////////////////////////////////////////////////
+/////////////////////////////////////////////////////////////////////////
+
+r_8 RandomGeneratorInterface::Gaussian() 
+{
+  switch (usegaussian_) {
+    case C_Gaussian_BoxMuller : 
+      return GaussianBoxMuller();
+      break;
+    case C_Gaussian_RandLibSNorm :
+      return GaussianSNorm();
+      break;
+    case C_Gaussian_PolarBoxMuller :
+      return GaussianPolarBoxMuller();
+      break;
+    case C_Gaussian_RatioUnif :
+      return GaussianRatioUnif();
+      break;
+    case C_Gaussian_LevaRatioUnif :
+      return GaussianLevaRatioUnif();
+      break;
+    default:
+      return GaussianBoxMuller();
+      break;
+  }
+}
 
 //--- Generation de nombre aleatoires suivant une distribution gaussienne 
@@ -374 +419 @@
 /////////////////////////////////////////////////////////////////////////
 
+uint_8 RandomGeneratorInterface::Poisson(double mu, double mumax)
+{
+  switch (usepoisson_) {
+    case C_Poisson_Simple :
+      return PoissonSimple(mu,mumax);
+      break;
+    case C_Poisson_Ahrens :
+      return PoissonAhrens(mu);
+      break;
+    default:
+      return PoissonSimple(mu,mumax);
+      break;
+  }
+}
+
+
 //--- Generation de nombre aleatoires suivant une distribution de Poisson 
 uint_8 RandomGeneratorInterface::PoissonSimple(double mu,double mumax)
@@ -634 +695 @@
 /////////////////////////////////////////////////////////////////////////
 /////////////////////////////////////////////////////////////////////////
+
+r_8 RandomGeneratorInterface::Exponential() 
+{
+  switch (useexpo_) {
+    case C_Exponential_Simple :
+      return ExpoSimple();
+      break;
+    case C_Exponential_Ahrens :
+      return ExpoAhrens();
+      break;
+    default:
+      return ExpoSimple();
+      break;
+  }
+}
 
 r_8 RandomGeneratorInterface::ExpoSimple(void)