fixed get_reject method
[physik/nlsop.git] / nlsop.c
diff --git a/nlsop.c b/nlsop.c
index 26ba461..0cdecc4 100644 (file)
--- a/nlsop.c
+++ b/nlsop.c
@@ -817,7 +817,7 @@ u32 get_reject_graph(info *my_info,d3_lattice *d3_l,char *file,u32 *graph) {
    if(a>d3_l->max_z*CELL_LENGTH) k=0;
    else 
    {
-    graph[(int)(a/CELL_LENGTH)]=(int)(URAND_MAX*b);
+    graph[(int)(a/CELL_LENGTH)]=(int)(URAND_MAX/100*b);
     flag[(int)(a/CELL_LENGTH)]=1;
    }
   }
@@ -838,7 +838,7 @@ u32 get_reject_graph(info *my_info,d3_lattice *d3_l,char *file,u32 *graph) {
  free(flag);
 
 #ifdef DEBUG_INTERPOL_PROFILE
- printf("debug: (interpolated profile)\n");
+ printf("debug: (interpolated profile)\n");
  for(i=0;i<d3_l->max_z;i++) printf("%d %d\n",i,graph[i]);
 #endif