ranlxs.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	RANLXS_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	rlxs_init (int level, int seed)
void	ranlxs (float r[], int n)
int	rlxs_size (void)
void	rlxs_get (int state[])
void	rlxs_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 float	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 359 of file ranlxs.cxx.

#define MASK   0xffffff

Definition at line 360 of file ranlxs.cxx.

#define RANLXS_C

Definition at line 41 of file ranlxs.cxx.

#define STEP	(		pi,
			pj
	)

Definition at line 382 of file ranlxs.cxx.

Referenced by update().

Function Documentation

static void define_constants ( void  )

static

Definition at line 484 of file ranlxs.cxx.

Referenced by rlxs_init(), and rlxs_reset().

{
   int k;

   one_bit=(float)(ldexp(1.0,-24));

   for (k=0;k<96;k++)
      next[k]=(k+1)%96;
}

static void error ( int  no )

static

Definition at line 417 of file ranlxs.cxx.

References exit(), and printf().

Referenced by rlxs_get(), rlxs_init(), and rlxs_reset().

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

void ranlxs	(	float	r[],
		int	n
	)

Author

Manuel Zambrana zambr.nosp@m.ana@.nosp@m.kph.u.nosp@m.ni-m.nosp@m.ainz..nosp@m.de and Dmitry Khaneftd khane.nosp@m.ftd@.nosp@m.kph.u.nosp@m.ni-m.nosp@m.ainz..nosp@m.de

Mainz, May 2011

Definition at line 566 of file ranlxs.cxx.

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

{
   int k;

   if (init==0)
      rlxs_init(0,1);

   for (k=0;k<n;k++) 
   {
      is=next[is];
      if (is==is_old)
         update();
      r[k]=one_bit*(float)(x.num[is]);      
   }
}

void rlxs_get

(

int

state[]

)

Definition at line 589 of file ranlxs.cxx.

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

{
   int k;

   if (init==0)
      error(3);

   state[0]=rlxs_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 rlxs_init	(	int	level,
		int	seed
	)

Definition at line 495 of file ranlxs.cxx.

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

Referenced by ranlxs().

{
   int i,k,l;
   int ibit,jbit,xbit[31];
   int ix,iy;

   if ((INT_MAX<2147483647)||(FLT_RADIX!=2)||(FLT_MANT_DIG<24))
      error(0);

   define_constants();
   
   if (level==0)
      pr=109;
   else 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=95;
   is_old=0;
   prm=pr%12;
   init=1;
}

void rlxs_reset

(

int

state[]

)

Definition at line 613 of file ranlxs.cxx.

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

{
   int k;

   if ((INT_MAX<2147483647)||(FLT_RADIX!=2)||(FLT_MANT_DIG<24))
      error(4);

   define_constants();

   if (state[0]!=rlxs_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!=109)&&(pr!=202)&&(pr!=397))||
       (ir<0)||(ir>11)||(jr<0)||(jr>11)||(jr!=((ir+7)%12))||
       (is<0)||(is>95))
      error(5);
}

int rlxs_size

(

void

)

Definition at line 583 of file ranlxs.cxx.

Referenced by rlxs_get(), and rlxs_reset().

{
   return(105);
}

static void update ( void  )

static

Definition at line 451 of file ranlxs.cxx.

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

Referenced by ranlxs().

{
   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 374 of file ranlxs.cxx.

int init =0

static

Definition at line 372 of file ranlxs.cxx.

int ir

static

Definition at line 372 of file ranlxs.cxx.

Referenced by rlxs_reset(), and update().

int is

static

Definition at line 372 of file ranlxs.cxx.

Referenced by ranlxs(), and update().

int is_old

static

Definition at line 372 of file ranlxs.cxx.

int jr

static

Definition at line 372 of file ranlxs.cxx.

Referenced by update().

int next[96]

static

Definition at line 372 of file ranlxs.cxx.

int num[96]

Definition at line 379 of file ranlxs.cxx.

float one_bit

static

Definition at line 373 of file ranlxs.cxx.

int pr

static

Definition at line 372 of file ranlxs.cxx.

Referenced by rlxs_get(), rlxs_init(), rlxs_reset(), and update().

int prm

static

Definition at line 372 of file ranlxs.cxx.

Referenced by update().

dble_vec_t vec[12]

Definition at line 378 of file ranlxs.cxx.

union { ... } x