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

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 Phokhara/Phokhara-00-01-00/src/ranlxd.c.

◆ MASK

#define MASK 0xffffff

Definition at line 354 of file Phokhara/Phokhara-00-01-00/src/ranlxd.c.

◆ STEP

#define STEP	(		pi,
			pj )

Definition at line 376 of file Phokhara/Phokhara-00-01-00/src/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 Phokhara/Phokhara-00-01-00/src/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[] )

Definition at line 586 of file Phokhara/Phokhara-00-01-00/src/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 )

Definition at line 493 of file Phokhara/Phokhara-00-01-00/src/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().

◆ rlxd_reset()

void rlxd_reset ( int state[] )

Definition at line 610 of file Phokhara/Phokhara-00-01-00/src/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);
}