ranlxd.h File Reference

Go to the source code of this file.

Functions
void	ranlxd (double r[], int n)
void	rlxd_init (int level, int seed)
int	rlxd_size (void)
void	rlxd_get (int state[])
void	rlxd_reset (int state[])

Function Documentation

ranlxd()

void ranlxd ( double r[], int n )

extern

Definition at line 502 of file BesEvtGen/phokhara/PHOKHARA/ranlxd.c.

                                 {
  int k;
 
  if ( init == 0 ) rlxd_init( 0, 1 );
 
  for ( k = 0; k < n; k++ )
  {
    is = next[is];
    if ( is == is_old ) update();
    r[k] = one_bit * ( (double)( x.num[is + 4] ) + one_bit * (double)( x.num[is] ) );
  }
}

rlxd_get()

void rlxd_get ( int state[] )

extern

Definition at line 526 of file BesEvtGen/phokhara/PHOKHARA/ranlxd.c.

                             {
  int k;
 
  if ( init == 0 ) error( 3 );
 
  state[0] = rlxd_size();
 
  for ( k = 0; k < 96; k++ ) state[k + 1] = x.num[k];
 
  state[97]  = carry.c1;
  state[98]  = carry.c2;
  state[99]  = carry.c3;
  state[100] = carry.c4;
 
  state[101] = pr;
  state[102] = ir;
  state[103] = jr;
  state[104] = is;
}

rlxd_init()

void rlxd_init ( int level, int seed )

extern

Definition at line 439 of file BesEvtGen/phokhara/PHOKHARA/ranlxd.c.

                                      {
  int i, k, l;
  int ibit, jbit, xbit[31];
  int ix, iy;
 
  if ( ( INT_MAX < 2147483647 ) || ( FLT_RADIX != 2 ) || ( FLT_MANT_DIG < 24 ) ||
       ( DBL_MANT_DIG < 48 ) )
    error( 0 );
 
  define_constants();
 
  if ( level == 1 ) pr = 202;
  else if ( level == 2 ) pr = 397;
  else error( 1 );
 
  i = seed;
 
  for ( k = 0; k < 31; k++ )
  {
    xbit[k] = i % 2;
    i /= 2;
  }
 
  if ( ( seed <= 0 ) || ( i != 0 ) ) error( 2 );
 
  ibit = 0;
  jbit = 18;
 
  for ( i = 0; i < 4; i++ )
  {
    for ( k = 0; k < 24; k++ )
    {
      ix = 0;
 
      for ( l = 0; l < 24; l++ )
      {
        iy = xbit[ibit];
        ix = 2 * ix + iy;
 
        xbit[ibit] = ( xbit[ibit] + xbit[jbit] ) % 2;
        ibit       = ( ibit + 1 ) % 31;
        jbit       = ( jbit + 1 ) % 31;
      }
 
      if ( ( k % 4 ) != i ) ix = 16777215 - ix;
 
      x.num[4 * k + i] = ix;
    }
  }
 
  carry.c1 = 0;
  carry.c2 = 0;
  carry.c3 = 0;
  carry.c4 = 0;
 
  ir     = 0;
  jr     = 7;
  is     = 91;
  is_old = 0;
  prm    = pr % 12;
  init   = 1;
}

rlxd_reset()

void rlxd_reset ( int state[] )

extern

Definition at line 546 of file BesEvtGen/phokhara/PHOKHARA/ranlxd.c.

                               {
  int k;
 
  if ( ( INT_MAX < 2147483647 ) || ( FLT_RADIX != 2 ) || ( FLT_MANT_DIG < 24 ) ||
       ( DBL_MANT_DIG < 48 ) )
    error( 4 );
 
  define_constants();
 
  if ( state[0] != rlxd_size() ) error( 5 );
 
  for ( k = 0; k < 96; k++ )
  {
    if ( ( state[k + 1] < 0 ) || ( state[k + 1] >= 167777216 ) ) error( 5 );
 
    x.num[k] = state[k + 1];
  }
 
  if ( ( ( state[97] != 0 ) && ( state[97] != 1 ) ) ||
       ( ( state[98] != 0 ) && ( state[98] != 1 ) ) ||
       ( ( state[99] != 0 ) && ( state[99] != 1 ) ) ||
       ( ( state[100] != 0 ) && ( state[100] != 1 ) ) )
    error( 5 );
 
  carry.c1 = state[97];
  carry.c2 = state[98];
  carry.c3 = state[99];
  carry.c4 = state[100];
 
  pr     = state[101];
  ir     = state[102];
  jr     = state[103];
  is     = state[104];
  is_old = 8 * ir;
  prm    = pr % 12;
  init   = 1;
 
  if ( ( ( pr != 202 ) && ( pr != 397 ) ) || ( ir < 0 ) || ( ir > 11 ) || ( jr < 0 ) ||
       ( jr > 11 ) || ( jr != ( ( ir + 7 ) % 12 ) ) || ( is < 0 ) || ( is > 91 ) )
    error( 5 );
}

rlxd_size()

int rlxd_size ( void )

extern

Definition at line 524 of file BesEvtGen/phokhara/PHOKHARA/ranlxd.c.

{ return ( 105 ); }