CLHEP::DualRand Class Reference

#include <DualRand.h>

Inheritance diagram for CLHEP::DualRand:

Public Member Functions
	DualRand ()
	DualRand (long seed)
	DualRand (std::istream &is)
	DualRand (int rowIndex, int colIndex)
virtual	~DualRand ()
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[]="DualRand.conf") const
void	restoreStatus (const char filename[]="DualRand.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)

Static Public Attributes
static const unsigned int	VECTOR_STATE_SIZE = 9

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 51 of file DualRand.h.

Constructor & Destructor Documentation

DualRand() [1/4]

CLHEP::DualRand::DualRand

(

)

Definition at line 81 of file DualRand.cc.

: HepRandomEngine(),
  numEngines(numberOfEngines++),
  tausworthe (1234567 + numEngines + 175321),
  integerCong(69607 * tausworthe + 54329, numEngines)
{
  theSeed = 1234567;
}

DualRand() [2/4]

CLHEP::DualRand::DualRand

(

long

seed

)

Definition at line 90 of file DualRand.cc.

: HepRandomEngine(),
  numEngines(0),
  tausworthe ((unsigned int)seed + 175321),
  integerCong(69607 * tausworthe + 54329, 8043) // MF - not numEngines
{
  theSeed = seed;
}

DualRand() [3/4]

CLHEP::DualRand::DualRand

(

std::istream &

is

)

Definition at line 99 of file DualRand.cc.

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

DualRand() [4/4]

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

Definition at line 106 of file DualRand.cc.

: HepRandomEngine(),
  numEngines(0),
  tausworthe (rowIndex + 1000 * colIndex + 85329),
  integerCong(69607 * tausworthe + 54329, 1123) // MF - not numengines
{
  theSeed = rowIndex;
}

~DualRand()

CLHEP::DualRand::~DualRand ( )

virtual

Definition at line 115 of file DualRand.cc.

{ }

Member Function Documentation

beginTag()

std::string CLHEP::DualRand::beginTag ( )

static

Definition at line 282 of file DualRand.cc.

                                { 
  return "DualRand-begin"; 
}

engineName()

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

inlinestatic

Definition at line 98 of file DualRand.h.

{return "DualRand";}

Referenced by restoreStatus().

flat()

double CLHEP::DualRand::flat ( )

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 117 of file DualRand.cc.

                      {
  unsigned int ic ( integerCong );
  unsigned int t  ( tausworthe  );
  return ( (t  ^ ic) * twoToMinus_32() +              // most significant part
           (t >> 11) * twoToMinus_53() +              // fill in remaining bits
           nearlyTwoToMinus_54()                // make sure non-zero
         );
}

Referenced by flatArray(), and main().

flatArray()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 126 of file DualRand.cc.

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

get() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 322 of file DualRand.cc.

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

get() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 265 of file DualRand.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,"DualRand-begin")) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "\nInput mispositioned or"
          << "\nDualRand state description missing or"
          << "\nwrong engine type found." << std::endl;
    return is;
  }
  return getState(is);
}

getState() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 336 of file DualRand.cc.

                                                           {
  std::vector<unsigned long>::const_iterator iv = v.begin()+1;
  if (!tausworthe.get(iv)) return false;
  if (!integerCong.get(iv)) return false;
  if (iv != v.end()) {
    std::cerr << 
        "\nDualRand get:state vector has wrong size: " << v.size() 
    << "\n         Apparently " << iv-v.begin() << " words were consumed\n";
    return false;
  }
  return true;
}

getState() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 286 of file DualRand.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 << "\nDualRand 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];
  tausworthe.get(is);
  integerCong.get(is);
  is >> std::ws;
  is.width(MarkerLen);  
   is >> endMarker;
  if (strcmp(endMarker,"DualRand-end")) {
    is.clear(std::ios::badbit | is.rdstate());
    std::cerr << "DualRand 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::DualRand::name ( ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 73 of file DualRand.cc.

{return "DualRand";}

operator double()

CLHEP::DualRand::operator double ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 222 of file DualRand.cc.

                          {
   return flat();
}

operator float()

CLHEP::DualRand::operator float ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 226 of file DualRand.cc.

                         {
  return (float) ( (integerCong ^ tausworthe) * twoToMinus_32() 
                        + nearlyTwoToMinus_54()    ); 
                    // add this so that zero never happens
}

operator unsigned int()

CLHEP::DualRand::operator unsigned int ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 232 of file DualRand.cc.

                                {
  return (integerCong ^ tausworthe) & 0xffffffff;
}

put() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 257 of file DualRand.cc.

                                              {
  std::vector<unsigned long> v;
  v.push_back (engineIDulong<DualRand>());
  tausworthe.put(v);
  integerCong.put(v);
  return v;
}

Referenced by put(), and saveStatus().

put() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 236 of file DualRand.cc.

                                                {
  char beginMarker[] = "DualRand-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[]   = "DualRand-end";
  long pr=os.precision(20);
  os << " " << beginMarker << " ";
  os << theSeed << " ";
  tausworthe.put(os);
  integerCong.put(os);
  os << " " <<  endMarker  << "\n";
  os.precision(pr);
  return os;
#endif
}

restoreStatus()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 172 of file DualRand.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 << "\nDualRand 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
    tausworthe.get(inFile);
    integerCong.get(inFile);
  }
}

saveStatus()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 143 of file DualRand.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
    long pr=outFile.precision(20);
    outFile << theSeed << std::endl;
    tausworthe.put(outFile);
    integerCong.put(outFile);
    outFile << std::endl; // This is superfluous but harmless
    outFile.precision(pr);
#endif
}

setSeed()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 132 of file DualRand.cc.

                                     {
  theSeed = seed;
  tausworthe  = Tausworthe((unsigned int)seed + 175321);
  integerCong = IntegerCong(69607 * tausworthe + 54329, 8043);
}

Referenced by setSeeds().

setSeeds()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 138 of file DualRand.cc.

                                               {
  setSeed(seeds ? *seeds : 1234567, 0);
  theSeeds = seeds;
}

showStatus()

void CLHEP::DualRand::showStatus ( ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 207 of file DualRand.cc.

                                {
  long pr=std::cout.precision(20);
  std::cout << std::endl;
  std::cout <<         "-------- DualRand engine status ---------"
        << std::endl;
  std::cout << "Initial seed          = " << theSeed << std::endl;
  std::cout << "Tausworthe generator  = " << std::endl;
  tausworthe.put(std::cout);
  std::cout << "\nIntegerCong generator = " << std::endl;
  integerCong.put(std::cout);
  std::cout << std::endl << "-----------------------------------------"
        << std::endl;
  std::cout.precision(pr);
}

Member Data Documentation

VECTOR_STATE_SIZE

const unsigned int CLHEP::DualRand::VECTOR_STATE_SIZE = 9

static

Definition at line 104 of file DualRand.h.

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

---

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

• DualRand.h
• DualRand.cc