/*SGPC: Simple Genetic Programming in C(c) 1993 by Walter Alden Tackett and Aviram Carmi  This code and documentation is copyrighted and is not in the public domain. All rights reserved.   - This notice may not be removed or altered.  - You may not try to make money by distributing the package or by using the   process that the code creates.  - You may not distribute modified versions without clearly documenting your   changes and notifying the principal author.  - The origin of this software must not be misrepresented, either by   explicit claim or by omission.  Since few users ever read sources,   credits must appear in the documentation.  - Altered versions must be plainly marked as such, and must not be   misrepresented as being the original software.  Since few users ever read   sources, credits must appear in the documentation.  - The authors are not responsible for the consequences of use of this    software, no matter how awful, even if they arise from flaws in it. If you make changes to the code, or have suggestions for changes,let us know!  (gpc@ipld01.hac.com)*/#ifndef lintstatic char random_c_rcsid[]="$Id: random.c,v 2.9 1993/04/22 07:39:12 gpc-avc Exp gpc-avc $";#endif/* * * $Log: random.c,v $ * Revision 2.9  1993/04/22  07:39:12  gpc-avc * Removed old log messages * * Revision 2.8  1993/04/14  04:53:06  gpc-avc * Finished modes for checkpoining * * */#include <math.h>#include "random.h"/* definitions for the random number generator (seed * GEN % MOD) */#define GEN 16807#define MOD 2147483647#define QUO 127773#define RES 2836static long	seed;#ifdef ANSI_FUNCvoid set_seed(  unsigned long s  )#elsevoid set_seed(s)  unsigned long	s;#endif{  seed = (long)s;}#ifdef ANSI_FUNCunsigned long get_seed(void)#elseunsigned long get_seed()#endif{  return seed;}#ifdef ANSI_FUNCfloat park_miller_randomizer(void)#elsefloat park_miller_randomizer()#endif{  float retval;  retval = (float) (((double) PMrand(&seed)) / ((double) MOD));  return retval;}/* Park-Miller "minimal standard" pseudo random number generator *//* Implementation by Jan Jannink (c) 1992 */#ifdef ANSI_FUNClong PMrand(	long *s	)#elselong PMrand(s)long	*s;#endif{  long   tmp;  tmp = (*s >> 16) * GEN;  *s = (*s & 65535) * GEN + ((tmp & 32767) << 16) + (tmp >> 15);  return (*s -= (*s < 0 ? MOD: 0));}#ifdef ANSI_FUNCfloat gaussian_noise(  float mean,  float sigma  )#elsefloat gaussian_noise(mean, sigma)   float mean;  float sigma;#endif{  float		gauss;  if (gaussian_noise_toggle) {    gaussian_noise_uniform1 = (random_float(1.0));    gaussian_noise_uniform2 = (random_float(1.0));    gaussian_noise_temp = sqrt(-2.0*log((double)gaussian_noise_uniform1));    gauss = gaussian_noise_temp *       (float)cos((double)(2.0*M_PI*gaussian_noise_uniform2));  }  else {    gauss = gaussian_noise_temp *       (float)sin((double)(2.0*M_PI*gaussian_noise_uniform2));  }  gaussian_noise_toggle = ! gaussian_noise_toggle;  return mean + (sigma * gauss);}#ifdef ANSI_FUNCfloat random_float( /* 0 <=  random_float() < f */  float f  )#elsefloat random_float(f) /* 0 <=  random_float() < f */  float 	f;#endif{  float	retval;  retval = (f*(1.0-park_miller_randomizer()));  return retval;}#ifdef ANSI_FUNCint random_int( /* 0 <= random_int() < i */  int i  )#elseint random_int(i) /* 0 <= random_int() < i */   int	i;#endif{  int	retval;  retval = (int) (i * random_float(1.0));  return retval;}