[Commits] [svn:einsteintoolkit] NullInterp/trunk/src/ (Rev. 17)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2014/04/29 07:52 PM

Modified:
 /trunk/src/
  NullInterp_MaskInit.F90

Log:
 Fix check for guard cells
 
 before it was using unitialized memory.
 
 Fix provided by Bela Szilagyi

File Changes:

Directory: /trunk/src/
======================

File [modified]: NullInterp_MaskInit.F90
Delta lines: +3 -3
===================================================================
--- trunk/src/NullInterp_MaskInit.F90	2014-04-25 19:58:03 UTC (rev 16)
+++ trunk/src/NullInterp_MaskInit.F90	2014-04-30 00:52:32 UTC (rev 17)
@@ -99,9 +99,9 @@
     ! for guard points near the patch boundaries
 
     if (((lbnd(1).eq.0).and.any(EG_mask(1:N_ang_stencil_size,:).ne.0)) .or. &
-         ((ubnd(1).eq.gsh(1)-1).and.any(guard_mask(lsh(1)-N_ang_stencil_size+1:,:).ne.0)) .or. &
-         ((lbnd(2).eq.0).and.any(guard_mask(:,1:N_ang_stencil_size).ne.0)) .or. &
-         ((ubnd(2).eq.gsh(2)-1).and.any(guard_mask(:,lsh(2)-N_ang_stencil_size+1:).ne.0))) &
+         ((ubnd(1).eq.gsh(1)-1).and.any(EG_mask(lsh(1)-N_ang_stencil_size+1:,:).ne.0)) .or. &
+         ((lbnd(2).eq.0).and.any(EG_mask(:,1:N_ang_stencil_size).ne.0)) .or. &
+         ((ubnd(2).eq.gsh(2)-1).and.any(EG_mask(:,lsh(2)-N_ang_stencil_size+1:).ne.0))) &
          call CCTK_WARN(0, "mask setup error");
 
     ! mark guard points