source: Sophya/trunk/SophyaLib/BaseTools/randr48.h@ 3615

#ifndef RANDR48_H_SEEN
#define RANDR48_H_SEEN

// Classes DR48RandGen
// Generateur aleatoire utilisant l'aleatoire drand48()
//
// R. Ansari          (C) UPS+LAL IN2P3/CNRS  2009
// C. Magneville      (C) DAPNIA/SPP  CEA     2009


#include "machdefs.h"
#include "objfio.h"
#include "ndatablock.h"
#include "randinterf.h"


/* ------- SPECIFIQUE Mac Darwin OSX  <= 10.2 (gcc < 3.x) -------- */
#if defined(Darwin) &&  defined(__GNUC__) && (__GNUC__ < 3)
#include "osx_values.h"
/* Declaration de drand48 et srand48 etc , mis ds SysSpec, Reza 11/02/2003 */
#ifdef __cplusplus
extern "C" {
#endif
unsigned short int * seed48(unsigned short int seed16v[3]);
void srand48(long seedval);
#ifdef __cplusplus
}
#endif
#define drand48() ((double)(random())/LONG_MAX)
#endif
/* ------- FIN SPECIFIQUE Mac OSX / Darwin V <= 10.2 --------- */


namespace SOPHYA {

//! Implementation of the RandomGeneratorInterface class using drand48() functions 
class DR48RandGen : public RandomGeneratorInterface {

 public:
  DR48RandGen(long int seed=12345);
  virtual ~DR48RandGen();

  virtual void SetSeed(long int seed);
  virtual void SetSeed(uint_2 seed[3]);
  virtual void GetSeed(uint_2 seed[3]);
  virtual void AutoInit(int lp=0);
  virtual void ShowRandom();

  friend class ObjFileIO<DR48RandGen> ;

 protected:
  virtual r_8 Next();

};


//--------------------------------------------------------------------------------
//--------------------------------------------------------------------------------
//--------------------------------------------------------------------------------

//! Version thread-safe de RandomGeneratorInterface avec drand48() functions 
class ThSDR48RandGen : public DR48RandGen {

 public:
  ThSDR48RandGen(size_t n=1024,bool tsafe=true);
  ThSDR48RandGen(ThSDR48RandGen const & rg);
  virtual ~ThSDR48RandGen();

  void SetBuffSize(size_t n);

  virtual void SetSeed(long int seed);
  virtual void SetSeed(uint_2 seed[3]);
  virtual void GetSeed(uint_2 seed[3]);
  virtual void AutoInit(int lp=0);
  virtual void ShowRandom();

  friend class ObjFileIO<ThSDR48RandGen> ;

 protected:
  // ---- protected data members
  NDataBlock<r_8>  rseq_;  
  size_t idx_;
  bool fg_nothrsafe;  // if true --> NOT thread-safe
  // ---- protected methods 
  void GenSeq();
  virtual r_8 Next();

  // Non thread-safe version of Init() and GetSeed()
  void SetSeed_P(uint_2 seed[3]);
  void GetSeed_P(uint_2 seed[3]);

};

// Classe pour la gestion de persistance PPF :  ObjFileIO<ThSDR48RandGen>

/*! Writes the random generator object state in the POutPersist stream \b os */
inline POutPersist& operator << (POutPersist& os, ThSDR48RandGen & obj)
{ ObjFileIO<ThSDR48RandGen> fio(&obj);  fio.Write(os);  return(os); }
/*! Reads the random generator object state from the PInPersist stream \b is */
inline PInPersist& operator >> (PInPersist& is, ThSDR48RandGen & obj)
{ ObjFileIO<ThSDR48RandGen> fio(&obj); is.SkipToNextObject(); fio.Read(is); return(is); }

// Classe pour la gestion de persistance PPF :  ObjFileIO<DR48RandGen>

/*! Writes the random generator object state in the POutPersist stream \b os */
inline POutPersist& operator << (POutPersist& os, DR48RandGen & obj)
{ ObjFileIO<DR48RandGen> fio(&obj);  fio.Write(os);  return(os); }
/*! Reads the random generator object state from the PInPersist stream \b is */
inline PInPersist& operator >> (PInPersist& is, DR48RandGen & obj)
{ ObjFileIO<DR48RandGen> fio(&obj); is.SkipToNextObject(); fio.Read(is); return(is); }

} /* namespace SOPHYA */

#endif