X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fnlsop.git;a=blobdiff_plain;f=nlsop.c;h=deb65c9f208c6cb334a19ee2dd414e44020d7d06;hp=ae629256ef3fb40ee25b56e95d16345d4e0b38c7;hb=6ca951e0c48b5b9372a34d4ea5c9e9e9704336f0;hpb=2cdf781ce8f65b1b008efd254325ae597b37d0d6

diff --git a/nlsop.c b/nlsop.c
index ae62925..deb65c9 100644
--- 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);