fsultra.c File Reference

#include <math.h>
#include "others_generators.h"
#include "machine.h"
#include "sciprint.h"

Include dependency graph for fsultra.c:

Go to the source code of this file.

Defines
#define	N   37
#define	N2   24
#define	SWB(c, x, y, z)   c = (y<0) ? (((z=x-y-c) < 0) || (x>=0)) : (((z=x-y-c) < 0) && (x>=0));
Functions
void sciprint	__PARAMS ((char *fmt,...))
void	advance_state_swb ()
int	set_state_fsultra_simple (double s1, double s2)
int	set_state_fsultra (double *s)
void	get_state_fsultra (double s[])
unsigned long	fsultra ()
Variables
static int	is_init = 0
static long	swb_state [N]
static int	swb_index = N
static int	swb_flag
static unsigned long	cong_state
static double	DEFAULT_SEED1 = 1234567.0
static double	DEFAULT_SEED2 = 7654321.0

---

Define Documentation

#define N   37

Definition at line 66 of file fsultra.c.

Referenced by advance_state_swb(), C2F(), ChildrenCounter(), contourif(), creadbmat(), createvarfrom(), DrawMerge3d(), fsultra(), get_state_fsultra(), get_state_mt(), intdgeesx(), intdgesvd(), intdgetrf(), intdlassq(), intdsyev(), intreadxls(), intzgeesx(), Merge3dBuildTable(), Merge3dDimension(), mxCalloc(), mxgetscalar(), mxSetJc(), randmt(), set_state_fsultra(), set_state_fsultra_simple(), set_state_mt(), set_state_mt_simple(), set_surface_color_property(), set_tics_style_property(), set_xtics_coord_property(), and set_ytics_coord_property().

#define N2   24

Definition at line 67 of file fsultra.c.

Referenced by advance_state_swb().

#define SWB	(	c,
		x,
		y,
		z		)	c = (y<0) ? (((z=x-y-c) < 0) || (x>=0)) : (((z=x-y-c) < 0) && (x>=0));

Definition at line 103 of file fsultra.c.

Referenced by advance_state_swb().

---

Function Documentation

void sciprint __PARAMS

(

(char *fmt,...)

)

void advance_state_swb

(

)

Definition at line 107 of file fsultra.c.

References i, N, N2, SWB, swb_flag, swb_index, and swb_state.

Referenced by fsultra(), and set_state_fsultra_simple().

{ 
  int i;
  /*
   *  The following are the heart of the system and should be
   *  written is assembler to be as fast as possible. It may even make sense
   *  to unravel the loop and simply wirte 37 consecutive SWB operations!
   */
  for (i=0;  i<N2; i++) 
    SWB(swb_flag,swb_state[i+N-N2],swb_state[i],swb_state[i]);
  for (i=N2; i<N;  i++) 
    SWB(swb_flag,swb_state[i  -N2],swb_state[i],swb_state[i]);
  swb_index = 0;
}

Here is the caller graph for this function:

unsigned long fsultra

(

)

Definition at line 225 of file fsultra.c.

References advance_state_swb(), cong_state, DEFAULT_SEED1, DEFAULT_SEED2, is_init, N, set_state_fsultra_simple(), swb_index, and swb_state.

{
  if (swb_index >= N)  /* generate N words at one time */
    { 
      if ( ! is_init )
        set_state_fsultra_simple(DEFAULT_SEED1, DEFAULT_SEED2);
      else
        advance_state_swb();
    }
  return (swb_state[swb_index++] ^ (cong_state = cong_state * 69069));
}

Here is the call graph for this function:

void get_state_fsultra

(

double

s[]

)

Definition at line 211 of file fsultra.c.

References cong_state, DEFAULT_SEED1, DEFAULT_SEED2, i, is_init, N, set_state_fsultra_simple(), swb_flag, swb_index, and swb_state.

Referenced by RandI().

{
  int i;

  if ( ! is_init )
    set_state_fsultra_simple(DEFAULT_SEED1, DEFAULT_SEED2);

  s[0] = (double)  swb_index;
  s[1] = (double)  swb_flag;
  s[2] = (double)  cong_state;
  for (i = 0 ; i < N ; i++) 
    s[i+3] = (double) (unsigned long) swb_state[i];
}

Here is the call graph for this function:

Here is the caller graph for this function:

int set_state_fsultra

(

double *

s

)

Definition at line 172 of file fsultra.c.

References cong_state, i, int, is_init, long, N, sciprint(), swb_flag, swb_index, and swb_state.

{ 
  double try;
  int i;

  try = s[0];
  if ( floor(try) != try || try < 0.0  ||  try > (double) N)
    {
      sciprint("\n\r the first component of the fsultra state, must be an integer in [0, %d] \n\r",N);
      return 0;
    }
  swb_index = (int) try;

  try = s[1];
  if ( try != 0.0  &&  try != 1.0)
    {
      sciprint("\n\r the second component of the fsultra state, must be 0 or 1 \n\r");
      return 0;
    }
  swb_flag = (int) try;

  try = s[2];
  if ( floor(try) != try  ||  try <= 0 ||  try > 4294967295.0 )
    {
      sciprint("\n\r the third component of the fsultra state, must be an integer in [1, 2^32-1] \n\r");
      return 0;
    }
  cong_state = (unsigned long) try;
 
  /* no verif here ... */
  for (i = 0 ; i < N ; i++) 
    swb_state[i] = (long) (((unsigned long) s[i+3]) & 0xffffffff);

  is_init = 1;
  return 1;
}

Here is the call graph for this function:

int set_state_fsultra_simple	(	double	s1,
		double	s2
	)

Definition at line 138 of file fsultra.c.

References advance_state_swb(), cong_state, i, is_init, j, long, N, sciprint(), swb_flag, swb_index, and swb_state.

Referenced by fsultra(), get_state_fsultra(), and RandI().

{ 
  unsigned long shrgx, tidbits=0;
  int i, j;

  if (    (s1 == floor(s1) && 0.0 <= s1 && s1 <= 4294967295.0)
       && (s2 == floor(s2) && 0.0 <= s2 && s2 <= 4294967295.0) )
    {
      cong_state = ((unsigned long) s1)*2 + 1;
      shrgx = (unsigned long) s2;
      for ( i=0 ; i<N ; i++)
        {
          for ( j=32 ; j>0 ; j--)
            { 
              cong_state = cong_state * 69069;
              shrgx = shrgx ^ (shrgx >> 15);
              shrgx = shrgx ^ (shrgx << 17);
              tidbits = (tidbits>>1) | (0x80000000 & (cong_state^shrgx));
            }
          swb_state[i] = tidbits;
        }
      swb_index = 0;
      swb_flag = 0;
      advance_state_swb();  /* pour retrouver la même séquence que ds scilab V3.0 */
      is_init = 1;
      return 1;
    }
  else
    {
      sciprint("\n\r bad seed for fsultra, must be integers in [0, 2^32-1] \n\r");
      return 0;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

Variable Documentation

unsigned long cong_state [static]

Definition at line 73 of file fsultra.c.

Referenced by fsultra(), get_state_fsultra(), set_state_fsultra(), and set_state_fsultra_simple().

double DEFAULT_SEED1 = 1234567.0 [static]

Definition at line 81 of file fsultra.c.

Referenced by fsultra(), and get_state_fsultra().

double DEFAULT_SEED2 = 7654321.0 [static]

Definition at line 81 of file fsultra.c.

Referenced by fsultra(), and get_state_fsultra().

int is_init = 0 [static]

Definition at line 69 of file fsultra.c.

int swb_flag [static]

Definition at line 72 of file fsultra.c.

Referenced by advance_state_swb(), get_state_fsultra(), set_state_fsultra(), and set_state_fsultra_simple().

int swb_index = N [static]

Definition at line 71 of file fsultra.c.

Referenced by advance_state_swb(), fsultra(), get_state_fsultra(), set_state_fsultra(), and set_state_fsultra_simple().

long swb_state[N] [static]

Definition at line 70 of file fsultra.c.

Referenced by advance_state_swb(), fsultra(), get_state_fsultra(), set_state_fsultra(), and set_state_fsultra_simple().

---