ranlxd.c File Reference

#include <limits.h>
#include <float.h>
#include <stdlib.h>
#include <stdio.h>
#include <math.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 353 of file ranlxd.c.

MASK

#define MASK   0xffffff

Definition at line 354 of file ranlxd.c.

STEP

#define STEP	(		pi,
			pj )

Definition at line 376 of file 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 563 of file 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 586 of file 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 493 of file 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 610 of file 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 580 of file ranlxd.c.

{
   return(105);
}

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

Variable Documentation

num

int num[96]

Definition at line 373 of file ranlxd.c.

Referenced by EvtPDL::alias(), barrel_conv2root(), ExtSteppingAction::CalculateEmcEndCopyNb(), ExtSteppingAction::CalculateEmcEndPhiNb(), MucCalibMgr::CheckEvent(), BesEmcConstruction::ComputeEndCopyNb(), BesEmcSD::ComputeEndCopyNb(), EmcROOTGeo::ComputeEndCopyNb(), ExtBesEmcConstruction::ComputeEndCopyNb(), EFEventLoopMgr::EF_Fragment_Address(), endcap_conv2root(), EvtSpinAmp::EvtSpinAmp(), EvtSpinAmp::EvtSpinAmp(), EvtSpinAmp::EvtSpinAmp(), CgemSegmentRecAlg::exe_v1(), EFtoTDS::execute(), BesEmcEndGeometry::GetCryNumInOneLayer(), ExtBesEmcEndGeometry::GetCryNumInOneLayer(), RecMucTrack::GetHitInSeg(), TofShower::getNeighbors(), RootTofElecDataCnv::i_createObj(), TreeTofElecDataCnv::i_createObj(), BesMucEfficiency::Initialize(), BesGlobalTrigSvc::initialize(), StringParse::intpiece(), StringParse::longpiece(), MucMark::NumInCol(), StringParse::numpiece(), KalmanVertexFit::numTrack(), EmcSelBhaEvent::OutputMV(), EvtPdfDiv< T >::pdf(), StringParse::piece(), BesMucTrig::readIndex(), MdcTSF::ReadLut(), MdcTF::ReadTFLut(), InductionGar2::setMultiElectrons(), 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 372 of file ranlxd.c.

Referenced by EFResult::EFResult(), Gam4pikp::execute(), ExtSteppingAction::GetGapID(), KalFitGemHit::getVecPos(), Intersection(), DSemilepAlg::isGoodTrack(), DTagTool::isGoodTrack(), BesGeoTrack::Paint(), ranlxdf_(), BesPolygon2D::RotatePhi(), HltInf::setCgemData(), THltInf::setCgemData(), HltInf::setConData(), THltInf::setConData(), DotsHelixFitter::setDChitsFitFlag(), DotsHelixFitter::setDChitsPhiAmbi(), Event::CgemMcHit::SetDigiIdxVec(), TCgemMc::SetDigiIdxVec(), HltInf::setEmcData(), THltInf::setEmcData(), HltInf::setMdcData(), THltInf::setMdcData(), HltInf::setMucData(), THltInf::setMucData(), HltInf::setTofData(), THltInf::setTofData(), Event::CgemMcHit::SetVclusterIdxVec(), TCgemMc::SetVclusterIdxVec(), Event::CgemMcHit::SetXclusterIdxVec(), TCgemMc::SetXclusterIdxVec(), StraightLine::StraightLine(), and MdcSWire::yAxis().