potential bug fixed

Meinersbur · Jul 24, 2012 · 8f44aae · 8f44aae
1 parent 2b3bd2f
commit 8f44aae
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/geometry_eo.c b/geometry_eo.c
@@ -1541,7 +1541,7 @@ void Hopping_Matrix_Indices(){
   int ix;
   int ioff = (VOLUME+RAND)/2;
   /**************** loop over all lattice sites ****************/
-  for (int icx = 0, icy = VOLUME/2; icx < VOLUME/2; icx++, icy++)
+  for (int icx = 0, icy = (VOLUME+RAND)/2; icx < VOLUME/2; icx++, icy++)
   {
     ix=g_eo2lexic[icx];
     /*********************** direction +0 ************************/
@@ -1599,8 +1599,8 @@ void Hopping_Matrix_Indices(){
     /************************ end of loop ************************/
 
   }
-  g_hi[(16*VOLUME)] = 0;
-  g_hi[(16*VOLUME)+1] = 0;
+  g_hi[(16*(VOLUME+RAND))] = 0;
+  g_hi[(16*(VOLUME+RAND))+1] = 0;
   return;
 }
 
diff --git a/init_geometry_indices.c b/init_geometry_indices.c
@@ -150,7 +150,7 @@ int init_geometry_indices(const int V) {
   for(i = 1; i < (T+4); i++){
     g_ipt[i] = g_ipt[i-1]+(LX+4);
   }
-  g_hi = (int*)calloc(16*V+2,sizeof(int));
+  g_hi = (int*)calloc(16*(VOLUME+RAND)+2,sizeof(int));
   if((void*) g_hi == NULL) return(40);
 
 #ifdef WITHLAPH