A change in the test to check if we are inside the mask. I was being too strict

asking for the interpolated mask to be exactly 1.0.  Round off error can in some
cases change this a tiny bit and the horizon finder would give up.  Now I just
check if the interpolated value is less than 0.99.  I also added a longer comment
on how this test works.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@120 89daf98e-ef62-4674-b946-b8ff9de2216c

1 files changed, 11 insertions, 2 deletions

diff --git a/src/AHFinder_int.F b/src/AHFinder_int.F
index ae5a202..5fd058d 100644
--- a/src/AHFinder_int.F
+++ b/src/AHFinder_int.F
@@ -596,7 +596,15 @@
      .   CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL)
 
 !     Check if we are either inside mask, or to close for comfort.
-      
+!     The way I do this is very simple, I just check the interpolated
+!     values of the mask.  If any of these values is not 1.0, it means
+!     that a point inside the mask is contaminating the interpolation,
+!     and this means in turn that we are too close to the mask.
+
+!     Notice that in practice I don't ask for the interpolated values
+!     to be exactly 1.0, since round off error can modify the value
+!     slightly without really meaning that we are close to the mask.
+
       if (.not.CCTK_EQUALS(ahf_mask,'off')) then
 
          interror3 = 0
@@ -604,8 +612,9 @@
          if (myproc.lt.npt*npp) then
             do j=0,l_nphi
                do i=0,l_ntheta
-                  if (intmask(i,j).ne.1.0D0) then
+                  if (intmask(i,j).lt.0.99D0) then
                      interror3 = 1
+                     print *,i,j,intmask(i,j)
                   end if
                end do
             end do