ranlxd.cxx 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	RANLXD_C
#define	BASE   0x1000000
#define	MASK   0xffffff
#define	STEP(pi, pj)

Functions
static void	error (int no)
static void	update (void)
static void	define_constants (void)
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[])

Variables
static int	init =0
static int	pr
static int	prm
static int	ir
static int	jr
static int	is
static int	is_old
static int	next [96]
static double	one_bit
static vec_t	carry
union {
dble_vec_t   vec [12]
int   num [96]
}	x

Macro Definition Documentation

#define BASE   0x1000000

Definition at line 361 of file ranlxd.cxx.

#define MASK   0xffffff

Definition at line 362 of file ranlxd.cxx.

#define RANLXD_C

Definition at line 41 of file ranlxd.cxx.

#define STEP	(		pi,
			pj
	)

Definition at line 384 of file ranlxd.cxx.

Referenced by update().

Function Documentation

static void define_constants ( void  )

static

Definition at line 486 of file ranlxd.cxx.

Referenced by rlxd_init(), and rlxd_reset().

{
   int k;

   one_bit=ldexp(1.0,-24);

   for (k=0;k<96;k++)
   {
      next[k]=(k+1)%96;
      if ((k%4)==3)
         next[k]=(k+5)%96;
   }   
}

static void error ( int  no )

static

Definition at line 419 of file ranlxd.cxx.

References exit(), and printf().

Referenced by RhoCandidate::Correlation(), PndTrkTracking2::EliminateSpuriousSZ_bis(), PndTrkTracking2::EliminateSpuriousSZ_ter(), KFParticleBaseSIMD::GetMomentum(), KFParticleBaseSIMD::GetPt(), PndLmdSensorAligner::readPairsFromBinary(), RKTrackRep::RKutta(), rlxd_get(), rlxd_init(), and rlxd_reset().

{
   switch(no)
   {
      case 0:
         printf("Error in rlxd_init\n");
         printf("Arithmetic on this machine is not suitable for ranlxd\n");
         break;
      case 1:
         printf("Error in subroutine rlxd_init\n");
         printf("Bad choice of luxury level (should be 1 or 2)\n");
         break;
      case 2:
         printf("Error in subroutine rlxd_init\n");
         printf("Bad choice of seed (should be between 1 and 2^31-1)\n");
         break;
      case 3:
         printf("Error in rlxd_get\n");
         printf("Undefined state (ranlxd is not initialized)\n");
         break;
      case 4:
         printf("Error in rlxd_reset\n");
         printf("Arithmetic on this machine is not suitable for ranlxd\n");
         break;
      case 5:
         printf("Error in rlxd_reset\n");
         printf("Unexpected input data\n");
         break;
   }         
   printf("Program aborted\n");
   exit(0);
}

void ranlxd	(	double	r[],
		int	n
	)

Definition at line 571 of file ranlxd.cxx.

References is, n, rlxd_init(), update(), and x.

Referenced by mzrnd().

{
   int k;

   if (init==0)
      rlxd_init(1,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]));      
   }
}

void rlxd_get

(

int

state[]

)

Definition at line 594 of file ranlxd.cxx.

References vec_t::c1, vec_t::c2, vec_t::c3, vec_t::c4, error(), pr, rlxd_size(), and x.

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

void rlxd_init	(	int	level,
		int	seed
	)

Definition at line 501 of file ranlxd.cxx.

References vec_t::c1, vec_t::c2, vec_t::c3, vec_t::c4, define_constants(), error(), i, pr, seed, and x.

Referenced by mz_pp_to_leplep_vandewi_init(), mz_pp_to_pipi_vandewi_init(), and ranlxd().

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

void rlxd_reset

(

int

state[]

)

Definition at line 618 of file ranlxd.cxx.

References vec_t::c1, vec_t::c2, vec_t::c3, vec_t::c4, define_constants(), error(), ir, pr, rlxd_size(), and x.

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

int rlxd_size

(

void

)

Definition at line 588 of file ranlxd.cxx.

Referenced by rlxd_get(), and rlxd_reset().

{
   return(105);
}

static void update ( void  )

static

Definition at line 453 of file ranlxd.cxx.

References d, ir, is, jr, pi, pr, prm, STEP, and x.

Referenced by Photospp::PH_HEPEVT_Interface::get(), GFKalman::processHit(), and ranlxd().

{
   int k,kmax,d;
   dble_vec_t *pmin,*pmax,*pi,*pj;

   kmax=pr;
   pmin=&x.vec[0];
   pmax=pmin+12;
   pi=&x.vec[ir];
   pj=&x.vec[jr];
      
   for (k=0;k<kmax;k++) 
   {
      STEP(pi,pj);
      pi+=1;
      pj+=1;
      if (pi==pmax)
         pi=pmin;      
      if (pj==pmax)
         pj=pmin; 
   }

   ir+=prm;
   jr+=prm;
   if (ir>=12)
      ir-=12;
   if (jr>=12)
      jr-=12;
   is=8*ir;
   is_old=is;
}

Variable Documentation

vec_t carry

static

Definition at line 376 of file ranlxd.cxx.

int init =0

static

Definition at line 374 of file ranlxd.cxx.

Referenced by autocutx(), crosstag(), cutfinderx(), evaltrig(), PndLLbarAnaTask::Exec(), PndHypMicroWriter::Exec(), PndHypMicroIdealWriter::Exec(), PndFTSCAGBTracker::FitTracks(), KFParticleBaseSIMD::InvertSym3(), RhoError::RhoError(), simubg(), softtrigger_kin5(), softtrigger_toy12(), TMVAApply(), TMVATester(), and toy_core().

int ir

static

Definition at line 374 of file ranlxd.cxx.

Referenced by PndTrkCombiLegendreTransform::ComputeThetaR(), PndFastSim::cutAndSmear(), PndMdtParamDigi::Init(), rlxd_reset(), and update().

int is

static

Definition at line 374 of file ranlxd.cxx.

Referenced by createRootGeometry_RICH(), PndMdtDigitization::Digitize(), PndMdtParamDigi::Digitize(), PndMdtPointsToWaveform::Digitize(), PndLmdBPtestTask::Exec(), PndTrkCTGeometryCalculations::IntersectionsWithClosedbiHexagonLeft(), PndSttTrackFinderReal::IntersectionsWithClosedbiHexagonLeft(), PndTrkCTGeometryCalculations::IntersectionsWithClosedbiHexagonRight(), PndSttTrackFinderReal::IntersectionsWithClosedbiHexagonRight(), PndTrkCTGeometryCalculations::IntersectionsWithClosedPolygon(), PndSttTrackFinderReal::IntersectionsWithClosedPolygon(), PndTrkCTGeometryCalculations::IntersectionsWithOpenPolygon(), PndSttTrackFinderReal::IntersectionsWithOpenPolygon(), PndGemTrackFinderQA::PrepareMCTracks(), PndGemTrackFinderOnHits::PrintMCTrackSegments(), PndGemTrackFinderOnHitsTB::PrintMCTrackSegments(), ranlxd(), PndLmdAlignManager::readConfigFile(), PndLmdSensorAligner::readPairsFromBinary(), and update().

int is_old

static

Definition at line 374 of file ranlxd.cxx.

int jr

static

Definition at line 374 of file ranlxd.cxx.

Referenced by update().

int next[96]

static

Definition at line 374 of file ranlxd.cxx.

Referenced by PndDrcReco::CanvasSave(), check_complete(), check_pid(), comp_complete(), comp_multiFiles(), comp_pid(), comp_recoqa(), PndEmcStructure::crystal_name_analysis(), evaltrig(), eventDisplay(), PndGemDigiAna::Finish(), PndGemFindHitsAna::Finish(), PndGemFindHitsQA::Finish(), PndGemMagneticFieldVsTrackParameters::Finish(), PndGemTrackFinderQA::Finish(), get_tuples(), PndDrcAccuDigiPixelDraw::Init(), PndTrackCandDraw::Init(), LoadManySimFiles(), lutaddbar(), main(), chigen::PythiaChiGen::next(), plot_all(), plotmyhistos(), PndEmcStructure::PndEmcStructure(), PndMasterTask::PrintTaskList(), PndEmc::ProcessHits(), PndMasterTask::SetVerbose(), sizetree(), tut_makegifs(), DrawHits::WriteToFile(), and PndDrcRecoLookupMap::WriteToFile().

int num[96]

Definition at line 381 of file ranlxd.cxx.

Referenced by PndLmdStripClusterTask::AddMSErr(), ana_d0d0b(), ana_d0d0b_micro(), ana_dsdsj(), ana_dsdsj2(), ana_dsdsj2_EvtLoop(), ana_dsdsj_full(), ana_etac_3pi0(), ana_Lambda(), ana_Lambda_fit(), ana_Lambda_fit2(), ana_pid(), PndRiemannTrack::calcSZChi2(), create_lumi(), create_lumi_geanetest(), PndDrc::ProcessHits(), PndBranchBurstBuilder_eventBased::SetNEvents(), GFBookkeeping::setNumber(), PndLmdTrackQ::SetNumDoubleMChits(), PndRecoDafFit::SetNumIterations(), PndRecoKalmanFit::SetNumIterations(), PndRecoDafFit2::SetNumIterations(), PndRecoKalmanFit2::SetNumIterations(), PndRecoMultiKalmanTask2::SetNumIterations(), PndRecoKalmanTask2::SetNumIterations(), PndRecoMultiKalmanTask::SetNumIterations(), PndRecoKalmanTask::SetNumIterations(), PndLmdTrackQ::SetNumMChits(), PndRecoDafFit::SetTrackRep(), PndRecoKalmanFit::SetTrackRep(), PndRecoMultiKalmanTask::SetTrackRep(), PndRecoKalmanTask::SetTrackRep(), PndTrkBoundaryParStraws2::SttTubeList(), and PndRiemannTrack::szFit().

double one_bit

static

Definition at line 375 of file ranlxd.cxx.

int pr

static

Definition at line 374 of file ranlxd.cxx.

Referenced by ana_dsdsj_full(), rlxd_get(), rlxd_init(), rlxd_reset(), and update().

int prm

static

Definition at line 374 of file ranlxd.cxx.

Referenced by fitsb(), and update().

dble_vec_t vec[12]

Definition at line 380 of file ranlxd.cxx.

Referenced by PndTrkTools::ComputePocaToPointOnCircle3(), PndLmdAlignQA::createHist(), PndDrcReco::DetermineCherenkov(), PndSttHelixHitProducer::Exec(), PndSttTrackFitterQATask::Exec(), PndLmdTrackFinderCATask::Exec(), OnlineDBuilderTask::Exec(), PndMvdAllDataEventAna::FillHitProjHistos(), PndLmdTrackFinderTask::FindHitsI(), PndStraightLineTrackFinderTask::FindHitsI(), PndLmdTrackFinderTask::FindHitsII(), PndStraightLineTrackFinderTask::FindHitsII(), PndLmdTrackFinderTask::FindHitsIII(), PndStraightLineTrackFinderTask::FindHitsIII(), GFDetPlane::getGraphics(), PndMQTopixHitProducer::GetTransformation(), PndMvdTopixHitProducer::GetTransformation(), PndTrackCandDraw::GetVector(), PndSttHelixTrackFitter::IntersectionFinder(), PndTrkLegendreNew::IntersectionFinder(), PndTrkCombiLegendreTask::IntersectionFinder(), PndTrkTrackFinder::IntersectionFinder(), PndPidIdealFwdOnlyAssociatorTask::IsForward(), PndEmcStructure::locateIndex(), and PndDrcLutFill::ProcessPhotonHit().

union { ... } x