#include <MixMaxRng.h>

Inheritance diagram for CLHEP::MixMaxRng:

Data Structures
struct	rng_state_st

Public Member Functions
double	flat ()

void	flatArray (const int size, double *vect)

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

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

long	getSeed () const

const long *	getSeeds () const

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

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

	MixMaxRng ()

	MixMaxRng (const MixMaxRng &rng)

	MixMaxRng (long seed)

	MixMaxRng (std::istream &is)

std::string	name () const

	operator double ()

	operator float ()

	operator unsigned int ()

bool	operator!= (const HepRandomEngine &engine)

MixMaxRng &	operator= (const MixMaxRng &rng)

bool	operator== (const HepRandomEngine &engine)

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

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

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

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

void	setSeed (long seed, int dum=0)

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

void	showStatus () const

	~MixMaxRng ()

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

static std::string	engineName ()

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

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

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

static double	exponent_bit_32 ()

static double	mantissa_bit_12 ()

static double	mantissa_bit_24 ()

static double	mantissa_bit_32 ()

static double	nearlyTwoToMinus_54 ()

static double	twoToMinus_32 ()

static double	twoToMinus_48 ()

static double	twoToMinus_49 ()

static double	twoToMinus_53 ()

Protected Attributes
long	theSeed

const long *	theSeeds

Private Types
typedef struct rng_state_st	rng_state_t

Private Member Functions
myuint_t	apply_bigskip (myuint_t Vout, myuint_t Vin, myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)

MixMaxRng	Branch ()

void	BranchInplace (int id)

double	convert1double (myuint_t u)

myuint_t	fmodmulM61 (myuint_t cum, myuint_t s, myuint_t a)

double	generate (int i)

myuint_t	get_next ()

double	get_next_float ()

double	get_next_float_packbits ()

double	iterate ()

myuint_t	iterate_raw_vec (myuint_t *Y, myuint_t sumtotOld)

	MixMaxRng (myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)

myuint_t	MOD_MULSPEC (myuint_t k)

myuint_t	modadd (myuint_t foo, myuint_t bar)

myuint_t	MULWU (myuint_t k)

myuint_t	precalc ()

void	print_state () const

void	seed64 (myuint_t seedval)

void	seed_spbox (myuint_t)

void	seed_uniquestream (myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)

void	seed_vielbein (unsigned int i)

Static Private Member Functions
static constexpr int	rng_get_N ()

static constexpr long long int	rng_get_SPECIAL ()

static constexpr int	rng_get_SPECIALMUL ()

Private Attributes
rng_state_t	S

Static Private Attributes
static constexpr int	BITS =61

static constexpr double	INV_M61 =0.43368086899420177360298E-18

static constexpr myuint_t	M61 =2305843009213693951ULL

static const int	N = 17

static constexpr long long int	SPECIAL = ((N==17)? 0 : ((N==240)? 487013230256099140ULL:0) )

static constexpr long long int	SPECIALMUL = ((N==17)? 36: ((N==240)? 51 :53) )

static constexpr unsigned int	VECTOR_STATE_SIZE = 2*N+4

Detailed Description

Definition at line 49 of file MixMaxRng.h.

Member Typedef Documentation

◆ rng_state_t

typedef struct rng_state_st CLHEP::MixMaxRng::rng_state_t

private

Definition at line 182 of file MixMaxRng.h.

Constructor & Destructor Documentation

◆ MixMaxRng() [1/5]

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

Definition at line 69 of file MixMaxRng.cc.

: HepRandomEngine()
{
   get(is);
}

References get().

◆ MixMaxRng() [2/5]

CLHEP::MixMaxRng::MixMaxRng ( )

Definition at line 55 of file MixMaxRng.cc.

: HepRandomEngine()
{
   int numEngines = ++numberOfEngines;
   setSeed(static_cast<long>(numEngines));
}

References CLHEP::anonymous_namespace{MixMaxRng.cc}::numberOfEngines(), and setSeed().

◆ MixMaxRng() [3/5]

CLHEP::MixMaxRng::MixMaxRng ( long seed )

Definition at line 62 of file MixMaxRng.cc.

: HepRandomEngine()
{
   theSeed=seed;
   setSeed(seed);
}

References setSeed(), and CLHEP::HepRandomEngine::theSeed.

◆ ~MixMaxRng()

CLHEP::MixMaxRng::~MixMaxRng ( )

Definition at line 75 of file MixMaxRng.cc.

76{

77}

◆ MixMaxRng() [4/5]

CLHEP::MixMaxRng::MixMaxRng ( const MixMaxRng & rng )

Definition at line 79 of file MixMaxRng.cc.

  : HepRandomEngine(rng)
{
   S.V = rng.S.V;
   S.sumtot= rng.S.sumtot;
   S.counter= rng.S.counter;
}

References CLHEP::MixMaxRng::rng_state_st::counter, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ MixMaxRng() [5/5]

CLHEP::MixMaxRng::MixMaxRng	(	myID_t	clusterID,
		myID_t	machineID,
		myID_t	runID,
		myID_t	streamID
	)

private

Member Function Documentation

◆ apply_bigskip()

myuint_t CLHEP::MixMaxRng::apply_bigskip	(	myuint_t *	Vout,
		myuint_t *	Vin,
		myID_t	clusterID,
		myID_t	machineID,
		myID_t	runID,
		myID_t	streamID
	)

private

Definition at line 599 of file MixMaxRng.cc.

{
  /*
    makes a derived state vector, Vout, from the mother state vector Vin
    by skipping a large number of steps, determined by the given seeding ID's
             
    it is mathematically guaranteed that the substreams derived in this way from the SAME (!!!) Vin will not collide provided
    1) at least one bit of ID is different
    2) less than 10^100 numbers are drawn from the stream
    (this is good enough : a single CPU will not exceed this in the lifetime of the universe, 10^19 sec,
    even if it had a clock cycle of Planch time, 10^44 Hz )
             
    Caution: never apply this to a derived vector, just choose some mother vector Vin, for example the unit vector by seed_vielbein(X,0),
    and use it in all your runs, just change runID to get completely nonoverlapping streams of random numbers on a different day.
             
    clusterID and machineID are provided for the benefit of large organizations who wish to ensure that a simulation
    which is running in parallel on a large number of  clusters and machines will have non-colliding source of random numbers.
             
    did i repeat it enough times? the non-collision guarantee is absolute, not probabilistic
 
   */
 
   const myuint_t skipMat17[128][17] =
   #include "CLHEP/Random/mixmax_skip_N17.icc"
   ;
            
   const myuint_t* skipMat[128];
   for (int i=0; i<128; i++) { skipMat[i] = skipMat17[i];}
            
   myID_t IDvec[4] = {streamID, runID, machineID, clusterID};
   int r,i,j,  IDindex;
   myID_t id;
   myuint_t Y[N], cum[N];
   myuint_t coeff;
   myuint_t* rowPtr;
   myuint_t sumtot=0;
            
   for (i=0; i<N; i++) { Y[i] = Vin[i]; sumtot = modadd( sumtot, Vin[i]); } ;
   for (IDindex=0; IDindex<4; IDindex++)
   { // go from lower order to higher order ID
     id=IDvec[IDindex];
     //printf("now doing ID at level %d, with ID = %d\n", IDindex, id);
     r = 0;
     while (id)
     {
       if (id & 1)
       {
         rowPtr = (myuint_t*)skipMat[r + IDindex*8*sizeof(myID_t)];
         for (i=0; i<N; i++){ cum[i] = 0; }
         for (j=0; j<N; j++)
         { // j is lag, enumerates terms of the poly
           // for zero lag Y is already given
           coeff = rowPtr[j]; // same coeff for all i
           for (i =0; i<N; i++){
             cum[i] =  fmodmulM61( cum[i], coeff ,  Y[i] ) ;
           }
           sumtot = iterate_raw_vec(Y, sumtot);
         }
         sumtot=0;
         for (i=0; i<N; i++){ Y[i] = cum[i]; sumtot = modadd( sumtot, cum[i]); } ;
       }
       id = (id >> 1); r++; // bring up the r-th bit in the ID
     }
   }
   sumtot=0;
   for (i=0; i<N; i++){ Vout[i] = Y[i]; sumtot = modadd( sumtot, Y[i]); }
   // returns sumtot, and copy the vector over to Vout
   return (sumtot) ;
}

References fmodmulM61(), iterate_raw_vec(), modadd(), N, and Y().

Referenced by seed_uniquestream().

◆ beginTag()

std::string CLHEP::MixMaxRng::beginTag ( )

static

Definition at line 402 of file MixMaxRng.cc.

{ 
   return "MixMaxRng-begin"; 
}

◆ Branch()

MixMaxRng CLHEP::MixMaxRng::Branch ( )

private

Definition at line 731 of file MixMaxRng.cc.

{
   S.sumtot = iterate_raw_vec(S.V.data(), S.sumtot); S.counter = 1;
   MixMaxRng tmp=*this;
   tmp.BranchInplace(0); // daughter id
   return tmp;
}

References BranchInplace(), CLHEP::MixMaxRng::rng_state_st::counter, iterate_raw_vec(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ BranchInplace()

void CLHEP::MixMaxRng::BranchInplace ( int id )

private

Definition at line 739 of file MixMaxRng.cc.

{
   // Dont forget to iterate the mother, when branching the daughter, or else will have collisions!
   // a 64-bit LCG from Knuth line 26, is used to mangle a vector component
   constexpr myuint_t MULT64=6364136223846793005ULL;
   myuint_t tmp=S.V[id];
   S.V[1] *= MULT64; S.V[id] &= M61;
   S.sumtot = MIXMAX_MOD_MERSENNE( S.sumtot + S.V[id] - tmp + M61);
   S.sumtot = iterate_raw_vec(S.V.data(), S.sumtot);// printf("iterating!\n");
   S.counter = 1;
}

References CLHEP::MixMaxRng::rng_state_st::counter, iterate_raw_vec(), M61, MIXMAX_MOD_MERSENNE, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by Branch().

◆ checkFile()

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

staticprotectedinherited

Definition at line 47 of file RandomEngine.cc.

                                           {
  if (!file) {
    std::cerr << "Failure to find or open file " << filename <<
    " in " << classname << "::" << methodname << "()\n";
    return false;
  }  
  return true;
}                

References geant4_check_module_cycles::file.

Referenced by CLHEP::DualRand::restoreStatus(), CLHEP::HepJamesRandom::restoreStatus(), CLHEP::MTwistEngine::restoreStatus(), CLHEP::RanecuEngine::restoreStatus(), CLHEP::RanluxEngine::restoreStatus(), CLHEP::Ranlux64Engine::restoreStatus(), and CLHEP::RanshiEngine::restoreStatus().

◆ convert1double()

double CLHEP::MixMaxRng::convert1double ( myuint_t u )

inlineprivate

Definition at line 149 of file MixMaxRng.h.

  {
    const double one = 1;
    const myuint_t onemask = *(myuint_t*)&one;
    myuint_t tmp = (u>>9) | onemask; // bits between 52 and 62 dont affect the result!
    double d = *(double*)&tmp;
    return d-1.0;
  }

Referenced by generate(), and get_next_float_packbits().

◆ engineName()

std::string CLHEP::MixMaxRng::engineName ( )

static

Definition at line 252 of file MixMaxRng.cc.

{
   return "MixMaxRng";
}

◆ exponent_bit_32()

static double CLHEP::HepRandomEngine::exponent_bit_32 ( )

inlinestaticprotectedinherited

◆ flat()

double CLHEP::MixMaxRng::flat ( )

inlinevirtual

Implements CLHEP::HepRandomEngine.

Definition at line 65 of file MixMaxRng.h.

65{ return (S.counter<=(N-1)) ? generate(S.counter):iterate(); }

CLHEP::MixMaxRng::generate

double generate(int i)

Definition: MixMaxRng.cc:272

CLHEP::MixMaxRng::iterate

double iterate()

Definition: MixMaxRng.cc:296

References CLHEP::MixMaxRng::rng_state_st::counter, generate(), iterate(), N, and S.

Referenced by flatArray().

◆ flatArray()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 326 of file MixMaxRng.cc.

{
   // fill_array( S, size, arrayDbl );
   for (int i=0; i<size; ++i) { vect[i] = flat(); }
}

References flat().

◆ fmodmulM61()

myuint_t CLHEP::MixMaxRng::fmodmulM61	(	myuint_t	cum,
		myuint_t	s,
		myuint_t	a
	)

private

Definition at line 683 of file MixMaxRng.cc.

  {
    const myuint_t MASK32=0xFFFFFFFFULL;
    myuint_t o,ph,pl,ah,al;
    o=(s)*a;
    ph = ((s)>>32);
    pl = (s) & MASK32;
    ah = a>>32;
    al = a & MASK32;
    o = (o & M61) + ((ph*ah)<<3) + ((ah*pl+al*ph + ((al*pl)>>32))>>29) ;
    o += cum;
    o = (o & M61) + ((o>>61));
    return o;
  }

References G4INCL::DeuteronDensity::anonymous_namespace{G4INCLDeuteronDensity.cc}::al, M61, MASK32, and CLHEP::s.

Referenced by apply_bigskip(), and MOD_MULSPEC().

◆ generate()

double CLHEP::MixMaxRng::generate ( int i )

private

Definition at line 272 of file MixMaxRng.cc.

{
   S.counter++;
#if defined(__clang__) || defined(__llvm__)
   return INV_M61*static_cast<double>(S.V[i]);
#elif defined(__GNUC__) && (__GNUC__ < 7) && (!defined(__ICC)) && defined(__x86_64__) && defined(__SSE2_MATH__)
   int64_t Z=S.V[i];
   double F=0.0;
   //#warning Using the inline assembler
    /* using SSE inline assemly to zero the xmm register, just before int64 -> double conversion,
       not necessary in GCC-5 or better, but huge penalty on earlier compilers
    */
   __asm__ __volatile__(  "pxor %0, %0;"
                          "cvtsi2sdq %1, %0;"
                          :"=x"(F)
                          :"r"(Z)
                       );
   return F*INV_M61;
#else
  //#warning other method
   return convert1double(S.V[i]); //get_next_float_packbits();
#endif
}

References convert1double(), CLHEP::MixMaxRng::rng_state_st::counter, INV_M61, S, CLHEP::MixMaxRng::rng_state_st::V, and Z.

Referenced by flat().

◆ get() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 440 of file MixMaxRng.cc.

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

References getState().

◆ get() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

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

References getState(), and CLHEP::MarkerLen.

Referenced by MixMaxRng().

◆ get_next()

myuint_t CLHEP::MixMaxRng::get_next ( )

private

Definition at line 517 of file MixMaxRng.cc.

{
   int i;
   i=S.counter;
            
   if ((i<=(N-1)) )
   {
     S.counter++;
     return S.V[i];
   }
   else
   {
     S.sumtot = iterate_raw_vec(S.V.data(), S.sumtot);
     S.counter=2;
     return S.V[1];
   }
}

References CLHEP::MixMaxRng::rng_state_st::counter, iterate_raw_vec(), N, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by get_next_float_packbits().

◆ get_next_float()

double CLHEP::MixMaxRng::get_next_float ( )

inlineprivate

Definition at line 131 of file MixMaxRng.h.

131{ return get_next_float_packbits(); }

CLHEP::MixMaxRng::get_next_float_packbits

double get_next_float_packbits()

Definition: MixMaxRng.cc:547

References get_next_float_packbits().

◆ get_next_float_packbits()

double CLHEP::MixMaxRng::get_next_float_packbits ( )

private

Definition at line 547 of file MixMaxRng.cc.

{
   myuint_t Z=get_next();
   return convert1double(Z);
}

References convert1double(), get_next(), and Z.

Referenced by get_next_float().

◆ getSeed()

long CLHEP::HepRandomEngine::getSeed ( ) const

inlineinherited

Definition at line 112 of file RandomEngine.h.

112{ return theSeed; }

References CLHEP::HepRandomEngine::theSeed.

◆ getSeeds()

const long * CLHEP::HepRandomEngine::getSeeds ( ) const

inlineinherited

Definition at line 115 of file RandomEngine.h.

115{ return theSeeds; }

CLHEP::HepRandomEngine::theSeeds

const long * theSeeds

Definition: RandomEngine.h:143

References CLHEP::HepRandomEngine::theSeeds.

◆ getState() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 450 of file MixMaxRng.cc.

{
   if (v.size() != VECTOR_STATE_SIZE ) {
     std::cerr <<
        "\nMixMaxRng::getState(): vector has wrong length - state unchanged\n";
     return false;
   }
   for (int i=1; i<2*rng_get_N() ; i=i+2) {
     S.V[i/2]= ( (v[i] & MASK32) | ( (myuint_t)(v[i+1]) << 32 ) );
     // unpack from a data structure which is 32-bit on some platforms
   }
   S.counter = v[2*rng_get_N()+1];
   precalc();
   if ( ( (v[2*rng_get_N()+2] & MASK32)
        | ( (myuint_t)(v[2*rng_get_N()+3]) << 32 ) ) != S.sumtot) {
     std::cerr << "\nMixMaxRng::getState(): vector has wrong checksum!"
               << "\nInput vector is probably mispositioned now.\n";
     return false;
   }
   return true;
}

References CLHEP::MixMaxRng::rng_state_st::counter, MASK32, precalc(), rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, CLHEP::MixMaxRng::rng_state_st::V, and VECTOR_STATE_SIZE.

◆ getState() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 407 of file MixMaxRng.cc.

{
   char endMarker[MarkerLen];
   is >> theSeed;
   for (int i=0; i<rng_get_N(); ++i)  is >> S.V[i];
   is >> S.counter;
   myuint_t checksum;
   is >> checksum;
   is >> std::ws;
   is.width(MarkerLen);
   is >> endMarker;
   if (strcmp(endMarker,"MixMaxRng-end")) {
       is.clear(std::ios::badbit | is.rdstate());
       std::cerr << "\nMixMaxRng state description incomplete."
                 << "\nInput stream is probably mispositioned now.\n";
       return is;
   }
   if ( S.counter < 0 || S.counter > rng_get_N() ) {
       std::cerr << "\nMixMaxRng::getState(): "
                 << "vector read wrong value of counter from file!"
                 << "\nInput stream is probably mispositioned now.\n";
       return is;
   }
   precalc();
   if ( checksum != S.sumtot) {
       std::cerr << "\nMixMaxRng::getState(): "
                 << "checksum disagrees with value stored in file!"
                 << "\nInput stream is probably mispositioned now.\n";
       return is;
   }
   return is;
}

References CLHEP::MixMaxRng::rng_state_st::counter, CLHEP::MarkerLen, precalc(), rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, CLHEP::HepRandomEngine::theSeed, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by get().

◆ iterate()

double CLHEP::MixMaxRng::iterate ( )

private

Definition at line 296 of file MixMaxRng.cc.

{
   myuint_t* Y=S.V.data();
   myuint_t  tempP, tempV;
   Y[0] = ( tempV = S.sumtot);
   myuint_t sumtot = Y[0], ovflow = 0; // will keep a running sum of all new elements
   tempP = 0;              // will keep a partial sum of all old elements
   myuint_t tempPO;
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[1] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[1]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[2] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[2]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[3] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[3]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[4] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[4]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[5] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[5]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[6] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[6]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[7] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[7]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[8] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[8]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[9] ); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[9]  = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[10]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[10] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[11]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[11] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[12]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[12] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[13]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[13] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[14]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[14] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[15]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[15] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   tempPO = MULWU(tempP); tempP = modadd(tempP, Y[16]); tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); Y[16] = tempV; sumtot += tempV; if (sumtot < tempV) {ovflow++;};
   S.sumtot = MIXMAX_MOD_MERSENNE(MIXMAX_MOD_MERSENNE(sumtot) + (ovflow <<3 ));
 
   S.counter=2;
   return double(S.V[1])*INV_M61;
}

References CLHEP::MixMaxRng::rng_state_st::counter, INV_M61, MIXMAX_MOD_MERSENNE, modadd(), MULWU(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, CLHEP::MixMaxRng::rng_state_st::V, and Y().

Referenced by flat().

◆ iterate_raw_vec()

myuint_t CLHEP::MixMaxRng::iterate_raw_vec	(	myuint_t *	Y,
		myuint_t	sumtotOld
	)

private

Definition at line 497 of file MixMaxRng.cc.

{
   // operates with a raw vector, uses known sum of elements of Y
   int i;
            
   myuint_t  tempP, tempV;
   Y[0] = ( tempV = sumtotOld);
   myuint_t sumtot = Y[0], ovflow = 0; // will keep a running sum of all new elements
   tempP = 0;              // will keep a partial sum of all old elements
   for (i=1; (i<N); i++)
   {
     myuint_t tempPO = MULWU(tempP);
     tempP = modadd(tempP, Y[i]);
     tempV = MIXMAX_MOD_MERSENNE(tempV+tempP+tempPO); // new Y[i] = old Y[i] + old partial * m
     Y[i] = tempV;
     sumtot += tempV; if (sumtot < tempV) {ovflow++;}
   }
   return MIXMAX_MOD_MERSENNE(MIXMAX_MOD_MERSENNE(sumtot) + (ovflow <<3 ));
}

References MIXMAX_MOD_MERSENNE, modadd(), MULWU(), N, and Y().

Referenced by apply_bigskip(), Branch(), BranchInplace(), and get_next().

◆ mantissa_bit_12()

static double CLHEP::HepRandomEngine::mantissa_bit_12 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::RanluxEngine::flat(), and CLHEP::RanluxEngine::flatArray().

◆ mantissa_bit_24()

static double CLHEP::HepRandomEngine::mantissa_bit_24 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::RanluxEngine::flat(), CLHEP::RanluxEngine::flatArray(), CLHEP::RanluxEngine::getState(), CLHEP::RanluxEngine::put(), CLHEP::RanluxEngine::setSeed(), and CLHEP::RanluxEngine::setSeeds().

◆ mantissa_bit_32()

static double CLHEP::HepRandomEngine::mantissa_bit_32 ( )

inlinestaticprotectedinherited

◆ MOD_MULSPEC()

myuint_t CLHEP::MixMaxRng::MOD_MULSPEC ( myuint_t k )

private

Definition at line 472 of file MixMaxRng.cc.

{
   switch (N)
   {
      case 17:
          return 0;
          break;
      case 8:
          return 0;
          break;
      case 240:
          return fmodmulM61( 0, SPECIAL , (k) );
          break;
      default:
          std::cerr << "MIXMAX ERROR: " << "Disallowed value of parameter N\n";
          std::terminate();
          break;
   }
}

References fmodmulM61(), N, and SPECIAL.

◆ modadd()

myuint_t CLHEP::MixMaxRng::modadd	(	myuint_t	foo,
		myuint_t	bar
	)

private

Definition at line 699 of file MixMaxRng.cc.

{
#if defined(__x86_64__) && defined(__GNUC__) && (!defined(__ICC))
   //#warning Using assembler routine in modadd
   myuint_t out;
   /* Assembler trick suggested by Andrzej Görlich     */
   __asm__ ("addq %2, %0; "
            "btrq $61, %0; "
            "adcq $0, %0; "
            :"=r"(out)
            :"0"(foo), "r"(bar)
           );
   return out;
#else
   return MIXMAX_MOD_MERSENNE(foo+bar);
#endif
}

References CLHEP::bar, and MIXMAX_MOD_MERSENNE.

Referenced by apply_bigskip(), iterate(), and iterate_raw_vec().

◆ MULWU()

myuint_t CLHEP::MixMaxRng::MULWU ( myuint_t k )

private

Definition at line 492 of file MixMaxRng.cc.

{
   return (( (k)<<(SPECIALMUL) & M61) ^ ( (k) >> (BITS-SPECIALMUL))  );
}

References BITS, M61, and SPECIALMUL.

Referenced by iterate(), and iterate_raw_vec().

◆ name()

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

inlinevirtual

Implements CLHEP::HepRandomEngine.

Definition at line 104 of file MixMaxRng.h.

104{ return "MixMaxRng"; }

Referenced by source.g4viscp.G4Scene::create_scene(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

◆ nearlyTwoToMinus_54()

static double CLHEP::HepRandomEngine::nearlyTwoToMinus_54 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::DualRand::flat(), CLHEP::MTwistEngine::flat(), and CLHEP::RanshiEngine::flat().

◆ newEngine() [1/2]

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

staticinherited

Definition at line 96 of file RandomEngine.cc.

                                                             {
  return EngineFactory::newEngine(v);
}

References CLHEP::EngineFactory::newEngine().

◆ newEngine() [2/2]

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

staticinherited

Definition at line 91 of file RandomEngine.cc.

                                                          {
  return EngineFactory::newEngine(is);
}

References CLHEP::EngineFactory::newEngine().

Referenced by CLHEP::StaticRandomStates::restore().

◆ operator double()

CLHEP::MixMaxRng::operator double ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 332 of file MixMaxRng.cc.

{
  return flat();
}

References G4AblaRandom::flat().

◆ operator float()

CLHEP::MixMaxRng::operator float ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 337 of file MixMaxRng.cc.

{
  return float( flat() );
}

References G4AblaRandom::flat().

◆ operator unsigned int()

CLHEP::MixMaxRng::operator unsigned int ( )

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 342 of file MixMaxRng.cc.

{
   return static_cast<unsigned int>(get_next());
   // clhep_get_next returns a 64-bit integer, of which the lower 61 bits
   // are random and upper 3 bits are zero
}

◆ operator!=()

bool CLHEP::HepRandomEngine::operator!= ( const HepRandomEngine & engine )

inlineinherited

◆ operator=()

MixMaxRng & CLHEP::MixMaxRng::operator= ( const MixMaxRng & rng )

Definition at line 87 of file MixMaxRng.cc.

{
   // Check assignment to self
   //
   if (this == &rng)  { return *this; }
 
   // Copy base class data
   //
   HepRandomEngine::operator=(rng);
 
   S.V = rng.S.V;
   S.sumtot= rng.S.sumtot;
   S.counter= rng.S.counter;
 
   return *this;
}

References CLHEP::MixMaxRng::rng_state_st::counter, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ operator==()

bool CLHEP::HepRandomEngine::operator== ( const HepRandomEngine & engine )

inlineinherited

◆ precalc()

myuint_t CLHEP::MixMaxRng::precalc ( )

private

Definition at line 535 of file MixMaxRng.cc.

{
   int i;
   myuint_t temp;
   temp = 0;
   for (i=0; i < N; i++){
     temp = MIXMAX_MOD_MERSENNE(temp + S.V[i]);
   }
   S.sumtot = temp;
   return temp;
}

References MIXMAX_MOD_MERSENNE, N, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by getState(), and restoreStatus().

◆ print_state()

void CLHEP::MixMaxRng::print_state ( ) const

private

Definition at line 717 of file MixMaxRng.cc.

{
   int j;
   std::cout << "mixmax state, file version 1.0\n";
   std::cout << "N=" << rng_get_N() << "; V[N]={";
   for (j=0; (j< (rng_get_N()-1) ); j++) {
     std::cout << S.V[j] << ", ";
   }
   std::cout << S.V[rng_get_N()-1];
   std::cout << "}; ";
   std::cout << "counter= " << S.counter;
   std::cout << "sumtot= " << S.sumtot << "\n";
}

References CLHEP::MixMaxRng::rng_state_st::counter, rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by restoreStatus(), and showStatus().

◆ put() [1/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 367 of file MixMaxRng.cc.

{
   std::vector<unsigned long> v;
   v.push_back (engineIDulong<MixMaxRng>());
   for (int i=0; i<rng_get_N(); ++i)
   {
     v.push_back(static_cast<unsigned long>(S.V[i] & MASK32));
       // little-ended order on all platforms
     v.push_back(static_cast<unsigned long>(S.V[i] >> 32  ));
       // pack uint64 into a data structure which is 32-bit on some platforms
   }
   v.push_back(static_cast<unsigned long>(S.counter));
   v.push_back(static_cast<unsigned long>(S.sumtot & MASK32));
   v.push_back(static_cast<unsigned long>(S.sumtot >> 32));
   return v;
}

References CLHEP::MixMaxRng::rng_state_st::counter, MASK32, rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ put() [2/2]

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

virtual

Reimplemented from CLHEP::HepRandomEngine.

Definition at line 349 of file MixMaxRng.cc.

{
   char beginMarker[] = "MixMaxRng-begin";
   char endMarker[]   = "MixMaxRng-end";
 
   int pr = os.precision(24);
   os << beginMarker << " ";
   os << theSeed << "\n";
   for (int i=0; i<rng_get_N(); ++i) {
      os <<  S.V[i] << "\n";
   }
   os << S.counter << "\n";
   os << S.sumtot << "\n";
   os << endMarker << "\n";
   os.precision(pr);
   return os;  
}

References CLHEP::MixMaxRng::rng_state_st::counter, rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, CLHEP::HepRandomEngine::theSeed, and CLHEP::MixMaxRng::rng_state_st::V.

◆ restoreStatus()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 124 of file MixMaxRng.cc.

{
   // a function for reading the state from a file
   FILE* fin;
   if(  ( fin = fopen(filename, "r") ) )
   {
      char l=0;
      while ( l != '{' ) { // 0x7B = "{"
        l=fgetc(fin); // proceed until hitting opening bracket
      }
      ungetc(' ', fin);
   }
   else
   {
      fprintf(stderr, "mixmax -> read_state: error reading file %s\n", filename);
      throw std::runtime_error("Error in reading state file");
   }
    
   myuint_t vecVal;
   //printf("mixmax -> read_state: starting to read state from file\n");
   if (!fscanf(fin, "%llu", &S.V[0]) )
   {
     fprintf(stderr, "mixmax -> read_state: error reading file %s\n", filename);
     throw std::runtime_error("Error in reading state file");
   }
 
   int i;
   for( i = 1; i < rng_get_N(); i++)
   {
     if (!fscanf(fin, ", %llu", &vecVal) )
     {
       fprintf(stderr, "mixmax -> read_state: error reading vector component i=%d from file %s\n", i, filename);
       throw std::runtime_error("Error in reading state file");
     }
     if(  vecVal <= MixMaxRng::M61 )
     {
       S.V[i] = vecVal;
     }
     else
     {
       fprintf(stderr, "mixmax -> read_state: Invalid state vector value= %llu"
               " ( must be less than %llu ) "
               " obtained from reading file %s\n"
               , vecVal, MixMaxRng::M61, filename);
     }
   }
    
   int counter;
   if (!fscanf( fin, "}; counter=%i; ", &counter))
   {
     fprintf(stderr, "mixmax -> read_state: error reading counter from file %s\n", filename);
     throw std::runtime_error("Error in reading state file");
   }
   if( counter <= rng_get_N() )
   {
     S.counter= counter;
   }
   else
   {
     fprintf(stderr, "mixmax -> read_state: Invalid counter = %d"
             "  Must be 0 <= counter < %u\n" , counter, rng_get_N());
     print_state();
     throw std::runtime_error("Error in reading state counter");
   }
   precalc();
   myuint_t sumtot;
   if (!fscanf( fin, "sumtot=%llu\n", &sumtot))
   {
     fprintf(stderr, "mixmax -> read_state: error reading checksum from file %s\n", filename);
     throw std::runtime_error("Error in reading state file");
   }
 
   if (S.sumtot != sumtot)
   {
     fprintf(stderr, "mixmax -> checksum error while reading state from file %s - corrupted?\n", filename);
     throw std::runtime_error("Error in reading state checksum");
   }
   fclose(fin);
}

References CLHEP::MixMaxRng::rng_state_st::counter, M61, precalc(), print_state(), rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ rng_get_N()

constexpr int CLHEP::MixMaxRng::rng_get_N ( )

staticconstexprprivate

Definition at line 257 of file MixMaxRng.cc.

{
   return N;
}

References N.

Referenced by getState(), print_state(), put(), restoreStatus(), and saveStatus().

◆ rng_get_SPECIAL()

constexpr long long int CLHEP::MixMaxRng::rng_get_SPECIAL ( )

staticconstexprprivate

Definition at line 262 of file MixMaxRng.cc.

{
   return SPECIAL;
}

References SPECIAL.

◆ rng_get_SPECIALMUL()

constexpr int CLHEP::MixMaxRng::rng_get_SPECIALMUL ( )

staticconstexprprivate

Definition at line 267 of file MixMaxRng.cc.

{
   return SPECIALMUL;
}

References SPECIALMUL.

◆ saveStatus()

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

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 104 of file MixMaxRng.cc.

{
   // Create a C file-handle or an object that can accept the same output
   FILE *fh= fopen(filename, "w");
   if( fh )
   {
     int j;
     fprintf(fh, "mixmax state, file version 1.0\n" );
     fprintf(fh, "N=%u; V[N]={", rng_get_N() );
     for (j=0; (j< (rng_get_N()-1) ); j++) {
         fprintf(fh, "%llu, ", S.V[j] );
     }
     fprintf(fh, "%llu", S.V[rng_get_N()-1] );
     fprintf(fh, "}; " );
     fprintf(fh, "counter=%u; ", S.counter );
     fprintf(fh, "sumtot=%llu;\n", S.sumtot );
     fclose(fh);
   }
}

References CLHEP::MixMaxRng::rng_state_st::counter, rng_get_N(), S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

◆ seed64()

void CLHEP::MixMaxRng::seed64 ( myuint_t seedval )

inlineprivate

Definition at line 138 of file MixMaxRng.h.

138{ seed_uniquestream( 0, 0, (myID_t)(seedval>>32), (myID_t)seedval ); } // seed with one 64-bit seed

CLHEP::MixMaxRng::seed_uniquestream

void seed_uniquestream(myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)

Definition: MixMaxRng.cc:592

References seed_uniquestream().

◆ seed_spbox()

void CLHEP::MixMaxRng::seed_spbox ( myuint_t seed )

private

Definition at line 571 of file MixMaxRng.cc.

{
   // a 64-bit LCG from Knuth line 26, in combination with a bit swap is used to seed
 
   const myuint_t MULT64=6364136223846793005ULL;
   int i;
            
   myuint_t sumtot=0,ovflow=0;
   if (seed == 0) throw std::runtime_error("try seeding with nonzero seed next time");
            
   myuint_t l = seed;
            
   for (i=0; i < N; i++){
     l*=MULT64; l = (l << 32) ^ (l>>32);
     S.V[i] = l & M61;
     sumtot += S.V[(i)]; if (sumtot < S.V[(i)]) {ovflow++;}
   }
   S.counter = N;  // set the counter to N if iteration should happen right away
   S.sumtot = MIXMAX_MOD_MERSENNE(MIXMAX_MOD_MERSENNE(sumtot) + (ovflow <<3 ));
}

References CLHEP::MixMaxRng::rng_state_st::counter, M61, MIXMAX_MOD_MERSENNE, N, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by setSeed().

◆ seed_uniquestream()

void CLHEP::MixMaxRng::seed_uniquestream	(	myID_t	clusterID,
		myID_t	machineID,
		myID_t	runID,
		myID_t	streamID
	)

private

Definition at line 592 of file MixMaxRng.cc.

{
   seed_vielbein(0);
   S.sumtot = apply_bigskip(S.V.data(), S.V.data(),  clusterID,  machineID,  runID,   streamID );
   S.counter = 1;
}

References apply_bigskip(), CLHEP::MixMaxRng::rng_state_st::counter, S, seed_vielbein(), CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by seed64(), and setSeeds().

◆ seed_vielbein()

void CLHEP::MixMaxRng::seed_vielbein ( unsigned int i )

private

Definition at line 553 of file MixMaxRng.cc.

{
   int i;
   if (index<N)
   {
     for (i=0; i < N; i++){
       S.V[i] = 0;
     }
     S.V[index] = 1;
   }
   else
   {
     std::terminate();
   }
   S.counter = N;  // set the counter to N if iteration should happen right away
   S.sumtot = 1;
}

References CLHEP::MixMaxRng::rng_state_st::counter, N, S, CLHEP::MixMaxRng::rng_state_st::sumtot, and CLHEP::MixMaxRng::rng_state_st::V.

Referenced by seed_uniquestream().

◆ setSeed()

void CLHEP::MixMaxRng::setSeed	(	long	seed,
		int	dum = `0`
	)

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 214 of file MixMaxRng.cc.

{
   //seed_uniquestream(0,0,0,longSeed);
   theSeed = longSeed;
   seed_spbox(longSeed);
}

References seed_spbox(), and CLHEP::HepRandomEngine::theSeed.

Referenced by MixMaxRng().

◆ setSeeds()

void CLHEP::MixMaxRng::setSeeds	(	const long *	seeds,
		int	seedNum = `0`
	)

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 225 of file MixMaxRng.cc.

{
   myID_t seed0, seed1= 0, seed2= 0, seed3= 0;
 
   if( seedNum < 1 ) {  // Assuming at least 2 seeds in vector...
       seed0= static_cast<myID_t>(Seeds[0]) & MASK32;
       seed1= static_cast<myID_t>(Seeds[1]) & MASK32;
   } 
   else
   {
     if( seedNum < 4 ) {
       seed0= static_cast<myID_t>(Seeds[0]) & MASK32;
       if( seedNum > 1){ seed1= static_cast<myID_t>(Seeds[1]) & MASK32; }
       if( seedNum > 2){ seed2= static_cast<myID_t>(Seeds[2]) & MASK32; }
     }
     if( seedNum >= 4 ) {
       seed0= static_cast<myID_t>(Seeds[0]) & MASK32;
       seed1= static_cast<myID_t>(Seeds[1]) & MASK32;
       seed2= static_cast<myID_t>(Seeds[2]) & MASK32;
       seed3= static_cast<myID_t>(Seeds[3]) & MASK32;
     }
   }
   theSeed = Seeds[0];
   theSeeds = Seeds;
   seed_uniquestream(seed3, seed2, seed1, seed0);
}

References MASK32, seed_uniquestream(), CLHEP::HepRandomEngine::theSeed, and CLHEP::HepRandomEngine::theSeeds.

◆ showStatus()

void CLHEP::MixMaxRng::showStatus ( ) const

virtual

Implements CLHEP::HepRandomEngine.

Definition at line 204 of file MixMaxRng.cc.

{
   std::cout << std::endl;
   std::cout << "------- MixMaxRng engine status -------" << std::endl;
 
   std::cout << " Current state vector is:" << std::endl;
   print_state();
   std::cout << "---------------------------------------" << std::endl;
}

References print_state().

◆ twoToMinus_32()

static double CLHEP::HepRandomEngine::twoToMinus_32 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::DualRand::flat(), CLHEP::MTwistEngine::flat(), CLHEP::RanshiEngine::flat(), CLHEP::Ranlux64Engine::setSeed(), and CLHEP::Ranlux64Engine::setSeeds().

◆ twoToMinus_48()

static double CLHEP::HepRandomEngine::twoToMinus_48 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::Ranlux64Engine::advance(), CLHEP::Ranlux64Engine::setSeed(), CLHEP::Ranlux64Engine::setSeeds(), and CLHEP::Ranlux64Engine::update().

◆ twoToMinus_49()

static double CLHEP::HepRandomEngine::twoToMinus_49 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::Ranlux64Engine::flat().

◆ twoToMinus_53()

static double CLHEP::HepRandomEngine::twoToMinus_53 ( )

inlinestaticprotectedinherited

Referenced by CLHEP::DualRand::flat(), CLHEP::MTwistEngine::flat(), and CLHEP::RanshiEngine::flat().

Field Documentation

◆ BITS

constexpr int CLHEP::MixMaxRng::BITS =61

staticconstexprprivate

Definition at line 116 of file MixMaxRng.h.

Referenced by MULWU().

◆ INV_M61

constexpr double CLHEP::MixMaxRng::INV_M61 =0.43368086899420177360298E-18

staticconstexprprivate

Definition at line 118 of file MixMaxRng.h.

Referenced by generate(), and iterate().

◆ M61

constexpr myuint_t CLHEP::MixMaxRng::M61 =2305843009213693951ULL

staticconstexprprivate

Definition at line 117 of file MixMaxRng.h.

Referenced by BranchInplace(), fmodmulM61(), MULWU(), restoreStatus(), and seed_spbox().

◆ N

const int CLHEP::MixMaxRng::N = 17

staticprivate

Definition at line 51 of file MixMaxRng.h.

Referenced by apply_bigskip(), flat(), get_next(), iterate_raw_vec(), MOD_MULSPEC(), precalc(), rng_get_N(), seed_spbox(), and seed_vielbein().

◆ S

rng_state_t CLHEP::MixMaxRng::S

private

Definition at line 183 of file MixMaxRng.h.

Referenced by Branch(), BranchInplace(), flat(), generate(), get_next(), getState(), iterate(), MixMaxRng(), operator=(), precalc(), print_state(), put(), restoreStatus(), saveStatus(), seed_spbox(), seed_uniquestream(), and seed_vielbein().

◆ SPECIAL

constexpr long long int CLHEP::MixMaxRng::SPECIAL = ((N==17)? 0 : ((N==240)? 487013230256099140ULL:0) )

staticconstexprprivate

Definition at line 113 of file MixMaxRng.h.

Referenced by MOD_MULSPEC(), and rng_get_SPECIAL().

◆ SPECIALMUL

constexpr long long int CLHEP::MixMaxRng::SPECIALMUL = ((N==17)? 36: ((N==240)? 51 :53) )

staticconstexprprivate

Definition at line 114 of file MixMaxRng.h.

Referenced by MULWU(), and rng_get_SPECIALMUL().

◆ theSeed

long CLHEP::HepRandomEngine::theSeed

protectedinherited

Definition at line 142 of file RandomEngine.h.

Referenced by CLHEP::DualRand::DualRand(), CLHEP::HepRandomEngine::getSeed(), CLHEP::DualRand::getState(), CLHEP::HepJamesRandom::getState(), getState(), CLHEP::MTwistEngine::getState(), CLHEP::RanecuEngine::getState(), CLHEP::Ranlux64Engine::getState(), CLHEP::RanluxEngine::getState(), CLHEP::RanshiEngine::getState(), CLHEP::HepJamesRandom::HepJamesRandom(), MixMaxRng(), CLHEP::RanecuEngine::put(), put(), CLHEP::MTwistEngine::put(), CLHEP::RanecuEngine::RanecuEngine(), CLHEP::RanluxEngine::RanluxEngine(), CLHEP::RanluxppEngine::RanluxppEngine(), CLHEP::RanshiEngine::RanshiEngine(), CLHEP::DualRand::restoreStatus(), CLHEP::HepJamesRandom::restoreStatus(), CLHEP::MTwistEngine::restoreStatus(), CLHEP::RanecuEngine::restoreStatus(), CLHEP::RanluxEngine::restoreStatus(), CLHEP::Ranlux64Engine::restoreStatus(), CLHEP::RanshiEngine::restoreStatus(), CLHEP::MTwistEngine::saveStatus(), CLHEP::RanecuEngine::setIndex(), CLHEP::RanecuEngine::setSeed(), CLHEP::HepJamesRandom::setSeed(), setSeed(), CLHEP::RanluxppEngine::setSeed(), CLHEP::Ranlux64Engine::setSeed(), CLHEP::RanluxEngine::setSeed(), CLHEP::DualRand::setSeed(), CLHEP::MTwistEngine::setSeed(), CLHEP::RanecuEngine::setSeeds(), CLHEP::Ranlux64Engine::setSeeds(), CLHEP::RanluxEngine::setSeeds(), setSeeds(), CLHEP::RanshiEngine::setSeeds(), CLHEP::DualRand::showStatus(), CLHEP::HepJamesRandom::showStatus(), CLHEP::MTwistEngine::showStatus(), CLHEP::RanecuEngine::showStatus(), CLHEP::Ranlux64Engine::showStatus(), CLHEP::RanluxEngine::showStatus(), and CLHEP::RanshiEngine::showStatus().

◆ theSeeds

const long* CLHEP::HepRandomEngine::theSeeds

protectedinherited

Definition at line 143 of file RandomEngine.h.

Referenced by CLHEP::HepRandomEngine::getSeeds(), CLHEP::RanecuEngine::RanecuEngine(), CLHEP::RanecuEngine::setIndex(), CLHEP::RanecuEngine::setSeed(), CLHEP::HepJamesRandom::setSeeds(), CLHEP::RanluxppEngine::setSeeds(), CLHEP::RanecuEngine::setSeeds(), CLHEP::Ranlux64Engine::setSeeds(), CLHEP::RanluxEngine::setSeeds(), setSeeds(), CLHEP::DualRand::setSeeds(), CLHEP::MTwistEngine::setSeeds(), and CLHEP::RanshiEngine::setSeeds().

◆ VECTOR_STATE_SIZE

constexpr unsigned int CLHEP::MixMaxRng::VECTOR_STATE_SIZE = 2*N+4

staticconstexprprivate

Definition at line 119 of file MixMaxRng.h.

Referenced by getState().

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

source/externals/clhep/include/CLHEP/Random/MixMaxRng.h
source/externals/clhep/src/MixMaxRng.cc

Data Structures

Public Member Functions

Static Public Member Functions

Static Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ rng_state_t

Constructor & Destructor Documentation

◆ MixMaxRng() [1/5]

◆ MixMaxRng() [2/5]

◆ MixMaxRng() [3/5]

◆ ~MixMaxRng()

◆ MixMaxRng() [4/5]

◆ MixMaxRng() [5/5]

Member Function Documentation

◆ apply_bigskip()

◆ beginTag()

◆ Branch()

◆ BranchInplace()

◆ checkFile()

◆ convert1double()

◆ engineName()

◆ exponent_bit_32()

◆ flat()

◆ flatArray()

◆ fmodmulM61()

◆ generate()

◆ get() [1/2]

◆ get() [2/2]

◆ get_next()

◆ get_next_float()

◆ get_next_float_packbits()

◆ getSeed()

◆ getSeeds()

◆ getState() [1/2]

◆ getState() [2/2]

◆ iterate()

◆ iterate_raw_vec()

◆ mantissa_bit_12()

◆ mantissa_bit_24()

◆ mantissa_bit_32()

◆ MOD_MULSPEC()

◆ modadd()

◆ MULWU()

◆ name()

◆ nearlyTwoToMinus_54()

◆ newEngine() [1/2]

◆ newEngine() [2/2]

◆ operator double()

◆ operator float()

◆ operator unsigned int()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ precalc()

◆ print_state()

◆ put() [1/2]

◆ put() [2/2]

◆ restoreStatus()

◆ rng_get_N()

◆ rng_get_SPECIAL()

◆ rng_get_SPECIALMUL()

◆ saveStatus()

◆ seed64()

◆ seed_spbox()

◆ seed_uniquestream()

◆ seed_vielbein()

◆ setSeed()

◆ setSeeds()

◆ showStatus()

◆ twoToMinus_32()

◆ twoToMinus_48()

◆ twoToMinus_49()

◆ twoToMinus_53()