added amorphous pressure ~1/r^2
[physik/morpheus.git] / random.c
index 0b7304f..51b310b 100644 (file)
--- a/random.c
+++ b/random.c
@@ -6,13 +6,14 @@
 #include <stdio.h>
 #include <stdlib.h>
 #include <unistd.h>
+
 #include "defines.h"
 
-/* global & extern variables */
-extern u32 sum_z_segments;
+/* global and extern variables */
+extern u32 sum_z_cells;
 extern int random_fd;
 
-/* return random integer between 0 - max */
+/* return random integer between 0 - max-1 */
 u32 rand_get(u32 max) 
 {
  u32 rand_int;
@@ -21,19 +22,21 @@ u32 rand_get(u32 max)
   puts("failed reading 4 bytes of random data");
   return -23;
  }
- /* segments numbered 0...max-1 */
+ /* cells numbered 0...max-1 */
  return((u32)(rand_int*(max*1.0/URAND_MAX)));
 }
 
 /* get z value (linear growth of probability with depths) */
-u32 rand_get_lgp(int slope_nel,int start_nel) {
+u32 rand_get_lgp(int slope_nel,int start_nel)
+{
  int z;
  u32 i;
- z=rand_get(sum_z_segments)+1; /* +1 as rand_get returns values 0...max-1 */
+ z=rand_get(sum_z_cells)+1; /* +1 as rand_get returns values 0...max-1 */
  for(i=1;;i++) {
   z-=(start_nel+i*slope_nel);
   if(z<=0) break;
  }
- return(i-1); /* return values 0...z_seg-1 */
+ return(i-1); /* return values 0...z_cell-1 */
 }
 
+