removed all nan's

sf_calc_action.c

@@ -38,42 +38,6 @@ void dirichlet_boundary_conditions(int t) {
   }
 }
 
-#define _su3_nan(u) \
-   (u).c00.re=1./0.0; \
-   (u).c00.im=1./0.0; \
-   (u).c01.re=1./0.0; \
-   (u).c01.im=1./0.0; \
-   (u).c02.re=1./0.0; \
-   (u).c02.im=1./0.0; \
-   (u).c10.re=1./0.0; \
-   (u).c10.im=1./0.0; \
-   (u).c11.re=1./0.0; \
-   (u).c11.im=1./0.0; \
-   (u).c12.re=1./0.0; \
-   (u).c12.im=1./0.0; \
-   (u).c20.re=1./0.0; \
-   (u).c20.im=1./0.0; \
-   (u).c21.re=1./0.0; \
-   (u).c21.im=1./0.0; \
-   (u).c22.re=1./0.0; \
-   (u).c22.im=1./0.0;
-
-/* the next function sets
-all the gauge links in the time direction (from t on) to nan.
-Note that the rest of the links at the boundaries (spatial links) are not yet touched here */
-void nan_dirichlet_boundary_conditions(int t) {
-
-  int ix;
-
-  for (ix=0;ix<VOLUME;ix++){
-
-    if (g_t[ix] == t) {
-
-      _su3_nan(g_gauge_field[ix][0]);
-
-    }   
-  }
-}
 
 /* this function sets all the spatial links located at the time boundaries to one
    it does nothing to the time-like links */

sf_calc_action.h

@@ -14,7 +14,6 @@ void set_all_links_to_one();
 void set_all_links_to_one_with_dirichlet(int t);
 void print_su3_matrix (su3 u);
 void sf_boundary_conditions_spatially_constant_abelian_field(int t, double eta);
-void nan_dirichlet_boundary_conditions(int t);
 
 /*** MEASUREMENTS ***/