#include <RanshiEngine.h>

Inheritance diagram for CLHEP::RanshiEngine:

Public Member Functions
	RanshiEngine ()

	RanshiEngine (std::istream &is)

	RanshiEngine (long seed)

	RanshiEngine (int rowIndex, int colIndex)

virtual	~RanshiEngine ()

double	flat ()

void	flatArray (const int size, double *vect)

void	setSeed (long seed, int)

void	setSeeds (const long *seeds, int)

void	saveStatus (const char filename[]="RanshiEngine.conf") const

void	restoreStatus (const char filename[]="RanshiEngine.conf")

void	showStatus () const

	operator double ()

	operator float ()

	operator unsigned int ()

virtual std::ostream &	put (std::ostream &os) const

virtual std::istream &	get (std::istream &is)

virtual std::istream &	getState (std::istream &is)

std::string	name () const

std::vector< unsigned long >	put () const

bool	get (const std::vector< unsigned long > &v)

bool	getState (const std::vector< unsigned long > &v)

Public Member Functions inherited from CLHEP::HepRandomEngine
	HepRandomEngine ()

virtual	~HepRandomEngine ()

bool	operator== (const HepRandomEngine &engine)

bool	operator!= (const HepRandomEngine &engine)

virtual double	flat ()=0

virtual void	flatArray (const int size, double *vect)=0

virtual void	setSeed (long seed, int)=0

virtual void	setSeeds (const long *seeds, int)=0

virtual void	saveStatus (const char filename[]="Config.conf") const =0

virtual void	restoreStatus (const char filename[]="Config.conf")=0

virtual void	showStatus () const =0

virtual std::string	name () const =0

virtual std::ostream &	put (std::ostream &os) const

virtual std::istream &	get (std::istream &is)

virtual std::istream &	getState (std::istream &is)

virtual std::vector< unsigned long >	put () const

virtual bool	get (const std::vector< unsigned long > &v)

virtual bool	getState (const std::vector< unsigned long > &v)

long	getSeed () const

const long *	getSeeds () const

virtual	operator double ()

virtual	operator float ()

virtual	operator unsigned int ()

Static Public Member Functions
static std::string	beginTag ()

static std::string	engineName ()

Static Public Member Functions inherited from CLHEP::HepRandomEngine
static std::string	beginTag ()

static HepRandomEngine *	newEngine (std::istream &is)

static HepRandomEngine *	newEngine (const std::vector< unsigned long > &v)

Additional Inherited Members
Static Protected Member Functions inherited from CLHEP::HepRandomEngine
static double	exponent_bit_32 ()

static double	mantissa_bit_12 ()

static double	mantissa_bit_24 ()

static double	mantissa_bit_32 ()

static double	twoToMinus_32 ()

static double	twoToMinus_48 ()

static double	twoToMinus_49 ()

static double	twoToMinus_53 ()

static double	nearlyTwoToMinus_54 ()

static bool	checkFile (std::istream &file, const std::string &filename, const std::string &classname, const std::string &methodname)

Protected Attributes inherited from CLHEP::HepRandomEngine
long	theSeed

const long *	theSeeds

Detailed Description

Author

Definition at line 53 of file RanshiEngine.h.

Constructor & Destructor Documentation

◆ RanshiEngine() [1/4]

CLHEP::RanshiEngine::RanshiEngine ( )

Definition at line 67 of file RanshiEngine.cc.

: HepRandomEngine(),
  halfBuff(0), numFlats(0) 
{
  int numEngines = numberOfEngines++;
  int i = 0;
  while (i < numBuff) {    
    buffer[i] = (unsigned int)((numEngines+19780503L*(i+1))& 0xffffffff);
    ++i;
  }
  theSeed = numEngines+19780503L*++i;
  redSpin = (unsigned int)(theSeed & 0xffffffff);
 
  for( i = 0; i < 10000; ++i) flat();  // Warm-up by running thorugh 10000 nums
}

Referenced by setSeed().

◆ RanshiEngine() [2/4]

CLHEP::RanshiEngine::RanshiEngine ( std::istream & is )

Definition at line 83 of file RanshiEngine.cc.

: HepRandomEngine(),
  halfBuff(0), numFlats(0) 
{
  is >> *this;
}

◆ RanshiEngine() [3/4]

CLHEP::RanshiEngine::RanshiEngine ( long seed )

Definition at line 90 of file RanshiEngine.cc.

: HepRandomEngine(),
  halfBuff(0), numFlats(0) 
{
  for (int i = 0; i < numBuff; ++i) {
    buffer[i] = (unsigned int)seed&0xffffffff;
  }
  theSeed = seed;
  redSpin = (unsigned int)(theSeed & 0xffffffff);
  int j;
  for (j = 0; j < numBuff*20; ++j) {      // "warm-up" for engine to hit
    flat();                               //  every ball on average 20X.
  }
}

◆ RanshiEngine() [4/4]

CLHEP::RanshiEngine::RanshiEngine	(	int	rowIndex,
		int	colIndex
	)

Definition at line 105 of file RanshiEngine.cc.

: HepRandomEngine(),
  halfBuff(0), numFlats(0) 
{
  int i = 0;
  while( i < numBuff ) {
    buffer[i] = (unsigned int)((rowIndex + (i+1)*(colIndex+8))&0xffffffff);
    ++i;
  }
  theSeed = rowIndex;
  redSpin = colIndex & 0xffffffff;
  for( i = 0; i < 100; ++i) flat();    // Warm-up by running thorugh 100 nums
}

◆ ~RanshiEngine()

CLHEP::RanshiEngine::~RanshiEngine ( )

virtual

Definition at line 119 of file RanshiEngine.cc.

119{ }

Member Function Documentation

◆ beginTag()

std::string CLHEP::RanshiEngine::beginTag ( )

static

Definition at line 329 of file RanshiEngine.cc.

                                    { 
  return "RanshiEngine-begin"; 
}

◆ engineName()

static std::string CLHEP::RanshiEngine::engineName ( )

inlinestatic

Definition at line 97 of file RanshiEngine.h.

97{return "RanshiEngine";}

Referenced by restoreStatus().

◆ flat()

double CLHEP::RanshiEngine::flat ( )

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 121 of file RanshiEngine.cc.

                          {
  unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
  unsigned int blkSpin     = buffer[redAngle] & 0xffffffff;
  unsigned int boostResult = blkSpin ^ redSpin;
 
  buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
  
  redSpin  = (blkSpin + numFlats++) & 0xffffffff;
  halfBuff = numBuff/2 - halfBuff;
  
  return ( blkSpin * twoToMinus_32() +            // most significant part
       (boostResult>>11) * twoToMinus_53() +  // fill in remaining bits
       nearlyTwoToMinus_54());          // non-zero
}

Referenced by flatArray(), main(), and RanshiEngine().

◆ flatArray()

void CLHEP::RanshiEngine::flatArray	(	const int	size,
		double *	vect
	)

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 136 of file RanshiEngine.cc.

                                                         {
  for (int i = 0; i < size; ++i) {
    vect[i] = flat();
  }
}

◆ get() [1/2]

bool CLHEP::RanshiEngine::get ( const std::vector< unsigned long > & v )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 371 of file RanshiEngine.cc.

                                                          {
  if ((v[0] & 0xffffffffUL) != engineIDulong<RanshiEngine>()) {
    std::cerr << 
        "\nRanshiEngine get:state vector has wrong ID word - state unchanged\n";
    return false;
  }
  return getState(v);
}

◆ get() [2/2]

std::istream & CLHEP::RanshiEngine::get ( std::istream & is )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 312 of file RanshiEngine.cc.

                                             {
  char beginMarker [MarkerLen];
  is >> std::ws;
  is.width(MarkerLen);  // causes the next read to the char* to be <=
            // that many bytes, INCLUDING A TERMINATION \0 
            // (Stroustrup, section 21.3.2)
  is >> beginMarker;
  if (strcmp(beginMarker,"RanshiEngine-begin")) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "\nInput mispositioned or"
          << "\nRanshiEngine state description missing or"
          << "\nwrong engine type found." << std::endl;
    return is;
  }
  return getState(is);
}

◆ getState() [1/2]

bool CLHEP::RanshiEngine::getState ( const std::vector< unsigned long > & v )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 380 of file RanshiEngine.cc.

                                                               {
  if (v.size() != VECTOR_STATE_SIZE ) {
    std::cerr << 
        "\nRanshiEngine get:state vector has wrong length - state unchanged\n";
    return false;
  }
  for (int i = 0; i < numBuff; ++i) {
    buffer[i] = (unsigned int)v[i+1];
  }
  redSpin  = (unsigned int)v[numBuff+1];
  numFlats = (unsigned int)v[numBuff+2]; 
  halfBuff = (unsigned int)v[numBuff+3];
  return true;
}

◆ getState() [2/2]

std::istream & CLHEP::RanshiEngine::getState ( std::istream & is )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 333 of file RanshiEngine.cc.

                                                  {
  if ( possibleKeywordInput ( is, "Uvec", theSeed ) ) {
    std::vector<unsigned long> v;
    unsigned long uu;
    for (unsigned int ivec=0; ivec < VECTOR_STATE_SIZE; ++ivec) {
      is >> uu;
      if (!is) {
        is.clear(std::ios::badbit | is.rdstate());
        std::cerr << "\nRanshiEngine state (vector) description improper."
        << "\ngetState() has failed."
           << "\nInput stream is probably mispositioned now." << std::endl;
        return is;
      }
      v.push_back(uu);
    }
    getState(v);
    return (is);
  }
 
//  is >> theSeed;  Removed, encompassed by possibleKeywordInput()
 
  char endMarker   [MarkerLen];
  for (int i = 0; i < numBuff; ++i) {
    is >> buffer[i];
  }
  is >> redSpin >> numFlats >> halfBuff;
  is >> std::ws;
  is.width(MarkerLen);  
  is >> endMarker;
  if (strcmp(endMarker,"RanshiEngine-end")) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "\nRanshiEngine state description incomplete."
          << "\nInput stream is probably mispositioned now." << std::endl;
    return is;
  }
  return is;
}

Referenced by get(), getState(), and restoreStatus().

◆ name()

std::string CLHEP::RanshiEngine::name ( ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 65 of file RanshiEngine.cc.

65{return "RanshiEngine";}

◆ operator double()

CLHEP::RanshiEngine::operator double ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 247 of file RanshiEngine.cc.

                              {
  return flat();
}

◆ operator float()

CLHEP::RanshiEngine::operator float ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 251 of file RanshiEngine.cc.

                             {
  unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
  unsigned int blkSpin  = buffer[redAngle] & 0xffffffff;
  
  buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
  
  redSpin  = (blkSpin + numFlats++) & 0xffffffff;
  halfBuff = numBuff/2 - halfBuff;
  
  return float(blkSpin * twoToMinus_32());
}

◆ operator unsigned int()

CLHEP::RanshiEngine::operator unsigned int ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 263 of file RanshiEngine.cc.

                                    {
  unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
  unsigned int blkSpin  = buffer[redAngle] & 0xffffffff;
  
  buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
  
  redSpin  = (blkSpin + numFlats++) & 0xffffffff;
  halfBuff = numBuff/2 - halfBuff;
  
  return blkSpin;
}

◆ put() [1/2]

std::vector< unsigned long > CLHEP::RanshiEngine::put ( ) const

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 300 of file RanshiEngine.cc.

                                                  {
  std::vector<unsigned long> v;
  v.push_back (engineIDulong<RanshiEngine>());
  for (int i = 0; i < numBuff; ++i) {
    v.push_back(static_cast<unsigned long>(buffer[i]));
  }
  v.push_back(static_cast<unsigned long>(redSpin));
  v.push_back(static_cast<unsigned long>(numFlats));
  v.push_back(static_cast<unsigned long>(halfBuff));  
  return v;
}

Referenced by put(), and saveStatus().

◆ put() [2/2]

std::ostream & CLHEP::RanshiEngine::put ( std::ostream & os ) const

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 275 of file RanshiEngine.cc.

                                                    {
  char beginMarker[] = "RanshiEngine-begin";
  os << beginMarker << "\nUvec\n";
  std::vector<unsigned long> v = put();
  for (unsigned int i=0; i<v.size(); ++i) {
     os <<  v[i] <<  "\n";
  }
  return os;  
#ifdef REMOVED 
  char endMarker[]   = "RanshiEngine-end";
  long pr=os.precision(20);
  os << " " << beginMarker << " ";
  
  os << theSeed  << "\n";
  for (int i = 0; i < numBuff; ++i) {
    os << buffer[i]  << "\n";
  }
  os << redSpin  << " " << numFlats << "\n" << halfBuff; 
  
  os << " " << endMarker   << "\n";
  os.precision(pr);
  return os;
#endif
}

◆ restoreStatus()

void CLHEP::RanshiEngine::restoreStatus ( const char filename[] = "RanshiEngine.conf" )

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 193 of file RanshiEngine.cc.

                                                      {
  std::ifstream inFile(filename, std::ios::in);
  if (!checkFile ( inFile, filename, engineName(), "restoreStatus" )) {
    std::cerr << "  -- Engine state remains unchanged\n";
    return;
  }
  if ( possibleKeywordInput ( inFile, "Uvec", theSeed ) ) {
    std::vector<unsigned long> v;
    unsigned long xin;
    for (unsigned int ivec=0; ivec < VECTOR_STATE_SIZE; ++ivec) {
      inFile >> xin;
           #ifdef TRACE_IO
           std::cout << "ivec = " << ivec << "  xin = " << xin << "    ";
           if (ivec%3 == 0) std::cout << "\n"; 
           #endif
      if (!inFile) {
        inFile.clear(std::ios::badbit | inFile.rdstate());
        std::cerr << "\nRanshiEngine state (vector) description improper."
           << "\nrestoreStatus has failed."
           << "\nInput stream is probably mispositioned now." << std::endl;
        return;
      }
      v.push_back(xin);
    }
    getState(v);
    return;
  }
 
  if (!inFile.bad()) {
//     inFile >> theSeed;  removed -- encompased by possibleKeywordInput
    for (int i = 0; i < numBuff; ++i) {
      inFile >> buffer[i];
    }
    inFile >> redSpin >> numFlats >> halfBuff;
  }
}

◆ saveStatus()

void CLHEP::RanshiEngine::saveStatus ( const char filename[] = "RanshiEngine.conf" ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 163 of file RanshiEngine.cc.

                                                         {
  std::ofstream outFile(filename, std::ios::out);
  if (!outFile.bad()) {
    outFile << "Uvec\n";
    std::vector<unsigned long> v = put();
             #ifdef TRACE_IO
             std::cout << "Result of v = put() is:\n"; 
             #endif
    for (unsigned int i=0; i<v.size(); ++i) {
      outFile << v[i] << "\n";
             #ifdef TRACE_IO
               std::cout << v[i] << " ";
               if (i%6==0) std::cout << "\n";
             #endif
    }
             #ifdef TRACE_IO
             std::cout << "\n";
             #endif
  }
#ifdef REMOVED
  if (!outFile.bad()) {
    outFile << std::setprecision(20) << theSeed << std::endl;
    for (int i = 0; i < numBuff; ++i) {
      outFile << buffer[i] << " ";
    }
    outFile << redSpin  << " " << numFlats << " " << halfBuff << std::endl;
  }
#endif
}

◆ setSeed()

void CLHEP::RanshiEngine::setSeed	(	long	seed,
		int
	)

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 142 of file RanshiEngine.cc.

                                         {
  *this = RanshiEngine(seed); 
}

◆ setSeeds()

void CLHEP::RanshiEngine::setSeeds	(	const long *	seeds,
		int
	)

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 146 of file RanshiEngine.cc.

                                                  {
  if (*seeds) {
    int i = 0;
    while (seeds[i] && i < numBuff) {
      buffer[i] = (unsigned int)seeds[i];
      ++i;
    }
    while (i < numBuff) {
      buffer[i] = buffer[i-1];
      ++i;
    }
    theSeed = seeds[0];
    redSpin = (unsigned int)theSeed;
  }
  theSeeds = seeds;
}

◆ showStatus()

void CLHEP::RanshiEngine::showStatus ( ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 230 of file RanshiEngine.cc.

                                    {
  std::cout << std::setprecision(20) << std::endl;
  std::cout << "----------- Ranshi engine status ----------" << std::endl;
  std::cout << "Initial seed      = " << theSeed << std::endl;
  std::cout << "Current red spin  = " << redSpin << std::endl;
  std::cout << "Values produced   = " << numFlats << std::endl;
  std::cout << "Side of buffer    = " << (halfBuff ? "upper" : "lower")
        << std::endl;
  std::cout << "Current buffer    = " << std::endl;
  for (int i = 0; i < numBuff; i+=4) {
    std::cout << std::setw(10) << std::setiosflags(std::ios::right)
          << buffer[i]     << std::setw(11) << buffer[i+1] << std::setw(11)
          << buffer[i+2]   << std::setw(11) << buffer[i+3] << std::endl;
  }
  std::cout << "-------------------------------------------" << std::endl;
}

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RanshiEngine() [1/4]

◆ RanshiEngine() [2/4]

◆ RanshiEngine() [3/4]

◆ RanshiEngine() [4/4]

◆ ~RanshiEngine()

Member Function Documentation

◆ beginTag()

◆ engineName()

◆ flat()

◆ flatArray()

◆ get() [1/2]

◆ get() [2/2]

◆ getState() [1/2]

◆ getState() [2/2]

◆ name()

◆ operator double()

◆ operator float()

◆ operator unsigned int()

◆ put() [1/2]

◆ put() [2/2]

◆ restoreStatus()

◆ saveStatus()

◆ setSeed()

◆ setSeeds()

◆ showStatus()