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 #include "defines.h"
  10
  11 /* global & extern variables */
  12 extern u32 sum_z_segments;
  13 extern int random_fd;
  14
  15 /* return random integer between 0 - max */
  16 u32 rand_get(u32 max)
  17 {
  18  u32 rand_int;
  19  if(read(random_fd,&rand_int,4)!=4)
  20  {
  21   puts("failed reading 4 bytes of random data");
  22   return -23;
  23  }
  24  /* segments numbered 0...max-1 */
  25  return((u32)(rand_int*(max*1.0/URAND_MAX)));
  26 }
  27
  28 /* get z value (linear growth of probability with depths) */
  29 u32 rand_get_lgp(int slope_nel,int start_nel) {
  30  int z;
  31  u32 i;
  32  z=rand_get(sum_z_segments)+1; /* +1 as rand_get returns values 0...max-1 */
  33  for(i=1;;i++) {
  34   z-=(start_nel+i*slope_nel);
  35   if(z<=0) break;
  36  }
  37  return(i-1); /* return values 0...z_seg-1 */
  38 }
  39