random.c

   1 /*
   2  * random.c - functions to get random values
   3  *
   4  */
   5
   6 #include <stdio.h>
   7 #include <stdlib.h>
   8 #include <unistd.h>
   9
  10 #include "defines.h"
  11
  12 /* global and extern variables */
  13 extern u32 gr;
  14 extern int random_fd;
  15 extern u32 *rand_buf,*rand_current;
  16
  17 /* return random integer between 0 - max-1 */
  18 u32 rand_get(u32 max)
  19 {
  20  if(rand_current>=rand_buf+(RAND_BUF_SIZE/sizeof(u32)))
  21  {
  22   printfd("debug: reading new random bytes\n");
  23   if(read(random_fd,rand_buf,RAND_BUF_SIZE)!=RAND_BUF_SIZE)
  24   {
  25    puts("failed reading 1 mega bytes of random data");
  26    return -23;
  27   }
  28   rand_current=rand_buf;
  29  }
  30  /* cells numbered 0...max-1 */
  31  return((u32)(*(rand_current++)*(max*1.0/URAND_MAX)));
  32 }
  33
  34 /* get z value (linear growth of probability with depths) */
  35 u32 rand_get_lgp(int slope,int start,u32 z_max)
  36 {
  37  int z;
  38  u32 i,weighted_sum_z;
  39  weighted_sum_z=z_max*start+slope*gr;
  40  z=rand_get(weighted_sum_z)+1; /* +1 as rand_get returns values 0...max-1 */
  41  for(i=1;;i++)
  42  {
  43   z-=(start+i*slope);
  44   if(z<=0) break;
  45  }
  46  return(i-1); /* return values 0...z_cell-1 */
  47 }
  48