new run

[physik/nlsop.git] / nlsop.c

diff --git a/nlsop.c b/nlsop.c
index ae62925..deb65c9 100644 (file)
--- a/nlsop.c
+++ b/nlsop.c
@@ -119,6 +119,7 @@ int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info,u32 nel_z)
   j+=(*(d3_l->status+x+((y-1+d3_l->max_y)%d3_l->max_y)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
   j+=(*(d3_l->status+x+y*d3_l->max_x+((z+1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
   j+=(*(d3_l->status+x+y*d3_l->max_x+((z-1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
+
   p+=((q/6)*j);
   if(get_rand(URAND_MAX)>p) MAKE_CRYST(thiz);
  }
@@ -1168,7 +1169,7 @@ int main(int argc,char **argv)
     y_c=get_rand(d3_l.max_y);
     // z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
     z_c=get_rand(d3_l.max_z);
-    nel_z=(URAND_MAX*n_e_loss[z_c])/ne_max;
+    nel_z=URAND_MAX*(1.0*n_e_loss[z_c]/ne_max);
     process_cell(&d3_l,x_c,y_c,z_c,&my_info,nel_z);
    }
    distrib_c(&d3_l,&my_info,i,c_ratio,ip_max,c_profile);