ranlxd.c File Reference

#include <float.h>
#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>

Go to the source code of this file.

Classes
struct	vec_t
struct	dble_vec_t

Macros
#define	BASE   0x1000000
#define	MASK   0xffffff
#define	STEP(pi, pj)

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

Macro Definition Documentation

BASE

#define BASE   0x1000000

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

MASK

#define MASK   0xffffff

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

STEP

#define STEP	(		pi,
			pj )

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

#  define STEP( pi, pj )                                                                      \
    d = ( *pj ).c1.c1 - ( *pi ).c1.c1 - carry.c1;                                             \
    ( *pi ).c2.c1 += ( d < 0 );                                                               \
    d += BASE;                                                                                \
    ( *pi ).c1.c1 = d & MASK;                                                                 \
    d             = ( *pj ).c1.c2 - ( *pi ).c1.c2 - carry.c2;                                 \
    ( *pi ).c2.c2 += ( d < 0 );                                                               \
    d += BASE;                                                                                \
    ( *pi ).c1.c2 = d & MASK;                                                                 \
    d             = ( *pj ).c1.c3 - ( *pi ).c1.c3 - carry.c3;                                 \
    ( *pi ).c2.c3 += ( d < 0 );                                                               \
    d += BASE;                                                                                \
    ( *pi ).c1.c3 = d & MASK;                                                                 \
    d             = ( *pj ).c1.c4 - ( *pi ).c1.c4 - carry.c4;                                 \
    ( *pi ).c2.c4 += ( d < 0 );                                                               \
    d += BASE;                                                                                \
    ( *pi ).c1.c4 = d & MASK;                                                                 \
    d             = ( *pj ).c2.c1 - ( *pi ).c2.c1;                                            \
    carry.c1      = ( d < 0 );                                                                \
    d += BASE;                                                                                \
    ( *pi ).c2.c1 = d & MASK;                                                                 \
    d             = ( *pj ).c2.c2 - ( *pi ).c2.c2;                                            \
    carry.c2      = ( d < 0 );                                                                \
    d += BASE;                                                                                \
    ( *pi ).c2.c2 = d & MASK;                                                                 \
    d             = ( *pj ).c2.c3 - ( *pi ).c2.c3;                                            \
    carry.c3      = ( d < 0 );                                                                \
    d += BASE;                                                                                \
    ( *pi ).c2.c3 = d & MASK;                                                                 \
    d             = ( *pj ).c2.c4 - ( *pi ).c2.c4;                                            \
    carry.c4      = ( d < 0 );                                                                \
    d += BASE;                                                                                \
    ( *pi ).c2.c4 = d & MASK

Function Documentation

ranlxd()

void ranlxd	(	double	r[],
		int	n )

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] ) );
  }
}

Referenced by ranlxdf_().

rlxd_get()

void rlxd_get

(

int

state[]

)

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;
}

Referenced by rlxdgetf_().

rlxd_init()

void rlxd_init	(	int	level,
		int	seed )

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;
}

Referenced by ranlxd(), and rlxdinit_().

rlxd_reset()

void rlxd_reset

(

int

state[]

)

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 );
}

Referenced by rlxdresetf_().

rlxd_size()

int rlxd_size

(

void

)

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

{ return ( 105 ); }

Referenced by rlxd_get(), rlxd_reset(), and rlxd_sizef_().

Variable Documentation

num

int num[96]

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

Referenced by EvtPDL::alias(), Dalitz::Babar_sakurai(), barrel_conv2root(), ExtSteppingAction::CalculateEmcEndCopyNb(), ExtSteppingAction::CalculateEmcEndPhiNb(), BesTofDigitizerEcV4::StripStruct::calNextN(), MucCalibMgr::CheckEvent(), BesEmcConstruction::ComputeEndCopyNb(), BesEmcSD::ComputeEndCopyNb(), EmcROOTGeo::ComputeEndCopyNb(), ExtBesEmcConstruction::ComputeEndCopyNb(), EvtDToKmPipPipPi0::decay(), EFEventLoopMgr::EF_Fragment_Address(), endcap_conv2root(), EvtSpinAmp::EvtSpinAmp(), EvtSpinAmp::EvtSpinAmp(), EvtSpinAmp::EvtSpinAmp(), EFtoTDS::execute(), ResetEtsAlg::execute(), QCMCFilter::findD0Decay(), BesEmcEndGeometry::GetCryNumInOneLayer(), ExtBesEmcEndGeometry::GetCryNumInOneLayer(), RecMucTrack::GetHitInSeg(), TofShower::getNeighbors(), EmcROOTGeo::GetPhysicalCrystal(), RootTofElecDataCnv::i_createObj(), TreeTofElecDataCnv::i_createObj(), BesMucEfficiency::Initialize(), BesGlobalTrigSvc::initialize(), QCMCFilter::initialize(), MucMark::NumInCol(), KalmanVertexFit::numTrack(), EmcSelBhaEvent::OutputMV(), EvtPdfDiv< T >::pdf(), EmcRecShowerPosLoglin::Position(), BesMucTrig::readIndex(), MdcTSF::ReadLut(), MdcTF::ReadTFLut(), Dalitz::sakurai(), CalibData::EndcapElec::setNum(), CalibData::EndcapElec::setNum(), EndcapElec::setNum(), EndcapElec::setNum(), BarrelElec::setNumEast(), BarrelElec::setNumEast(), CalibData::BarrelElec::setNumEast(), CalibData::BarrelElec::setNumEast(), BarrelElec::setNumWest(), BarrelElec::setNumWest(), CalibData::BarrelElec::setNumWest(), CalibData::BarrelElec::setNumWest(), KalmanVertexFit::setTrackIteration(), and KalmanVertexFit::setVertexIteration().

vec

dble_vec_t vec[12]

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

Referenced by EFResult::EFResult(), ElementLinkVector< DOBJ, StoragePolicy, IndexingPolicy >::ElementLinkVector(), Gam4pikp::execute(), ExtSteppingAction::GetGapID(), BesTofGeoParameter::GetVec(), Intersection(), DSemilepAlg::isGoodTrack(), DTagTool::isGoodTrack(), DTagTool::isGoodTrack_loose(), BesGeoTrack::Paint(), ranlxdf_(), BesPolygon2D::RotatePhi(), HltInf::setConData(), THltInf::setConData(), HltInf::setEmcData(), THltInf::setEmcData(), HltInf::setMdcData(), THltInf::setMdcData(), HltInf::setMucData(), THltInf::setMucData(), HltInf::setTofData(), THltInf::setTofData(), HltProcessor::Signature::Signature(), ElementLinkVector< DOBJ, StoragePolicy, IndexingPolicy >::swap(), and MdcSWire::yAxis().