BornAgain  1.19.79
Open-source research software to simulate and fit neutron and x-ray reflectometry and grazing-incidence small-angle scattering
Classes | Macros | Typedefs | Functions
mixmax.h File Reference

Go to the source code of this file.

Classes

struct  rng_state_st
 

Macros

#define ARRAY_INDEX_OUT_OF_BOUNDS   0xFF01
 
#define BITS   61
 
#define ERROR_READING_STATE_CHECKSUM   0xFF05
 
#define ERROR_READING_STATE_COUNTER   0xFF04
 
#define ERROR_READING_STATE_FILE   0xFF03
 
#define get_next(X)   GET_BY_MACRO(X)
 
#define get_next_float(X)   get_next_float_BY_MACRO(X)
 
#define INV_MERSBASE   (0.4336808689942017736029811203479766845703E-18)
 
#define M61   2305843009213693951ULL
 
#define MERSBASE   M61
 
#define MOD_MERSENNE(k)   MOD_PAYNE(k)
 
#define MOD_MULSPEC(k)   fmodmulM61( 0, SPECIAL , (k) )
 
#define MOD_PAYNE(k)   ((((k)) & MERSBASE) + (((k)) >> BITS) )
 
#define MOD_REM(k)   ((k) % MERSBASE )
 
#define N   240
 
#define SEED_WAS_ZERO   0xFF02
 
#define SPECIAL   487013230256099140ULL
 
#define SPECIALMUL   51
 

Typedefs

typedef uint32_t myID_t
 
typedef uint64_t myuint
 
typedef struct rng_state_st rng_state_t
 

Functions

myuint apply_bigskip (myuint *Vout, myuint *Vin, myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)
 
void branch_inplace (rng_state_t *Xin, myID_t *ID)
 
void fill_array (rng_state_t *X, unsigned int n, double *array)
 
myuint fmodmulM61 (myuint cum, myuint s, myuint a)
 
myuint GET_BY_MACRO (rng_state_t *X)
 
double get_next_float_BY_MACRO (rng_state_t *X)
 
int iterate (rng_state_t *X)
 
void iterate_and_fill_array (rng_state_t *X, double *array)
 
myuint iterate_raw_vec (myuint *Y, myuint sumtotOld)
 
myuint modadd (myuint foo, myuint bar)
 
myuint precalc (rng_state_t *X)
 
void print_state (rng_state_t *X)
 
void read_state (rng_state_t *X, const char filename[])
 
rng_state_t * rng_alloc ()
 
rng_state_t * rng_copy (myuint *Y)
 
int rng_free (rng_state_t *X)
 
int rng_get_N (void)
 
void seed_spbox (rng_state_t *X, myuint seed)
 
void seed_uniquestream (rng_state_t *X, myID_t clusterID, myID_t machineID, myID_t runID, myID_t streamID)
 
void seed_vielbein (rng_state_t *X, unsigned int i)
 

Macro Definition Documentation

◆ ARRAY_INDEX_OUT_OF_BOUNDS

#define ARRAY_INDEX_OUT_OF_BOUNDS   0xFF01

Definition at line 228 of file mixmax.h.

◆ BITS

#define BITS   61

Definition at line 110 of file mixmax.h.

◆ ERROR_READING_STATE_CHECKSUM

#define ERROR_READING_STATE_CHECKSUM   0xFF05

Definition at line 232 of file mixmax.h.

◆ ERROR_READING_STATE_COUNTER

#define ERROR_READING_STATE_COUNTER   0xFF04

Definition at line 231 of file mixmax.h.

◆ ERROR_READING_STATE_FILE

#define ERROR_READING_STATE_FILE   0xFF03

Definition at line 230 of file mixmax.h.

◆ get_next

#define get_next (   X)    GET_BY_MACRO(X)

Definition at line 189 of file mixmax.h.

◆ get_next_float

#define get_next_float (   X)    get_next_float_BY_MACRO(X)

Definition at line 190 of file mixmax.h.

◆ INV_MERSBASE

#define INV_MERSBASE   (0.4336808689942017736029811203479766845703E-18)

Definition at line 126 of file mixmax.h.

◆ M61

#define M61   2305843009213693951ULL

Definition at line 114 of file mixmax.h.

◆ MERSBASE

#define MERSBASE   M61

Definition at line 120 of file mixmax.h.

◆ MOD_MERSENNE

#define MOD_MERSENNE (   k)    MOD_PAYNE(k)

Definition at line 123 of file mixmax.h.

◆ MOD_MULSPEC

#define MOD_MULSPEC (   k)    fmodmulM61( 0, SPECIAL , (k) )

Definition at line 171 of file mixmax.h.

◆ MOD_PAYNE

#define MOD_PAYNE (   k)    ((((k)) & MERSBASE) + (((k)) >> BITS) )

Definition at line 121 of file mixmax.h.

◆ MOD_REM

#define MOD_REM (   k)    ((k) % MERSBASE )

Definition at line 122 of file mixmax.h.

◆ N

#define N   240

Definition at line 31 of file mixmax.h.

◆ SEED_WAS_ZERO

#define SEED_WAS_ZERO   0xFF02

Definition at line 229 of file mixmax.h.

◆ SPECIAL

#define SPECIAL   487013230256099140ULL

Definition at line 170 of file mixmax.h.

◆ SPECIALMUL

#define SPECIALMUL   51

Definition at line 169 of file mixmax.h.

Typedef Documentation

◆ myID_t

typedef uint32_t myID_t

Definition at line 76 of file mixmax.h.

◆ myuint

typedef uint64_t myuint

Definition at line 45 of file mixmax.h.

◆ rng_state_t

typedef struct rng_state_st rng_state_t

Definition at line 45 of file mixmax.h.

Function Documentation

◆ apply_bigskip()

myuint apply_bigskip ( myuint *  Vout,
myuint *  Vin,
myID_t  clusterID,
myID_t  machineID,
myID_t  runID,
myID_t  streamID 
)

◆ branch_inplace()

void branch_inplace ( rng_state_t *  Xin,
myID_t *  ID 
)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::Branch().

◆ fill_array()

void fill_array ( rng_state_t *  X,
unsigned int  n,
double *  array 
)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::RndmArray().

◆ fmodmulM61()

myuint fmodmulM61 ( myuint  cum,
myuint  s,
myuint  a 
)

◆ GET_BY_MACRO()

myuint GET_BY_MACRO ( rng_state_t *  X)
inline

Definition at line 195 of file mixmax.h.

195  {
196  int i;
197  i=X->counter;
198 
199  if (i<=(N-1) ){
200  X->counter++;
201  return X->V[i];
202  }else{
203  X->sumtot = iterate_raw_vec(X->V, X->sumtot);
204  X->counter=2;
205  return X->V[1];
206  }
207 }
N
#define N
Definition: mixmax.h:31
iterate_raw_vec
myuint iterate_raw_vec(myuint *Y, myuint sumtotOld)
rng_state_st::V
myuint V[N]
Definition: mixmax.h:54
rng_state_st::counter
int counter
Definition: mixmax.h:56
rng_state_st::sumtot
myuint sumtot
Definition: mixmax.h:55

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

Referenced by get_next_float_BY_MACRO().

Here is the call graph for this function:

◆ get_next_float_BY_MACRO()

double get_next_float_BY_MACRO ( rng_state_t *  X)
inline

Definition at line 209 of file mixmax.h.

209  {
210  /* cast to signed int trick suggested by Andrzej Gòˆrlich */
211  int64_t Z=(int64_t)GET_BY_MACRO(X);
212  double F;
213 #if defined(__GNUC__) && (__GNUC__ < 5) && (!defined(__ICC)) && defined(__x86_64__) && defined(__SSE2_MATH__) && defined(USE_INLINE_ASM)
214 //#warning Using the inline assembler
215 /* using SSE inline assemly to zero the xmm register, just before int64 -> double conversion,
216  not really necessary in GCC-5 or better, but huge penalty on earlier compilers
217  */
218  __asm__ __volatile__("pxor %0, %0; "
219  :"=x"(F)
220  );
221 #endif
222  F=Z;
223  return F*INV_MERSBASE;
224 }
INV_MERSBASE
#define INV_MERSBASE
Definition: mixmax.h:126
GET_BY_MACRO
myuint GET_BY_MACRO(rng_state_t *X)
Definition: mixmax.h:195

References GET_BY_MACRO(), and INV_MERSBASE.

Here is the call graph for this function:

◆ iterate()

int iterate ( rng_state_t *  X)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::Iterate().

◆ iterate_and_fill_array()

void iterate_and_fill_array ( rng_state_t *  X,
double *  array 
)

◆ iterate_raw_vec()

myuint iterate_raw_vec ( myuint *  Y,
myuint  sumtotOld 
)

Referenced by GET_BY_MACRO().

◆ modadd()

myuint modadd ( myuint  foo,
myuint  bar 
)

◆ precalc()

myuint precalc ( rng_state_t *  X)

◆ print_state()

void print_state ( rng_state_t *  X)

◆ read_state()

void read_state ( rng_state_t *  X,
const char  filename[] 
)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::ReadState().

◆ rng_alloc()

rng_state_t* rng_alloc ( )

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::MixMaxEngineImpl().

◆ rng_copy()

rng_state_t* rng_copy ( myuint *  Y)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::SetState().

◆ rng_free()

int rng_free ( rng_state_t *  X)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::~MixMaxEngineImpl(), and ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::SetState().

◆ rng_get_N()

int rng_get_N ( void  )

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::GetState(), and ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::Size().

◆ seed_spbox()

void seed_spbox ( rng_state_t *  X,
myuint  seed 
)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::SetSeedFast().

◆ seed_uniquestream()

void seed_uniquestream ( rng_state_t *  X,
myID_t  clusterID,
myID_t  machineID,
myID_t  runID,
myID_t  streamID 
)

Referenced by ROOT::Math::MixMaxEngineImpl< ROOT_MM_N >::SetSeed().

◆ seed_vielbein()

void seed_vielbein ( rng_state_t *  X,
unsigned int  i 
)