[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 403)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/07/17 12:08 PM

Modified:
 /trunk/src/
  GRHydro_UpdateMask.F90

Log:
 GRHydro: only loop over points where rhs is defined for UpdateMask
 
 From: Roland Haas <roland.haas at physics.gatech.edu>

File Changes:

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

File [modified]: GRHydro_UpdateMask.F90
Delta lines: +6 -6
===================================================================
--- trunk/src/GRHydro_UpdateMask.F90	2012-07-17 17:08:35 UTC (rev 402)
+++ trunk/src/GRHydro_UpdateMask.F90	2012-07-17 17:08:37 UTC (rev 403)
@@ -38,9 +38,9 @@
 
   if(evolve_temper.ne.1.and.evolve_Y_e.ne.1) then
      !$OMP PARALLEL DO PRIVATE(k,j,i)
-     do k = 1, cctk_lsh(3)
-        do j = 1, cctk_lsh(2)
-           do i = 1, cctk_lsh(1)
+     do k = 1+GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil
+        do j = 1+GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil
+           do i = 1+GRHydro_stencil, cctk_lsh(1)-GRHydro_stencil
               if ( GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. &
                    (atmosphere_mask(i,j,k) .ne. 0) .or. &
                       (tau(i,j,k) + frac * taurhs(i,j,k) .le. 0.d0) .or. &
@@ -59,9 +59,9 @@
   else
      ! allow negative total energy
      !$OMP PARALLEL DO PRIVATE(k,j,i)
-     do k = 1, cctk_lsh(3)
-        do j = 1, cctk_lsh(2)
-           do i = 1, cctk_lsh(1)
+     do k = 1+GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil
+        do j = 1+GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil
+           do i = 1+GRHydro_stencil, cctk_lsh(1)-GRHydro_stencil
               if ( GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. &
                    (atmosphere_mask(i,j,k) .ne. 0) .or. &
                       (dens(i,j,k) + frac * densrhs(i,j,k) .le. 0.d0) ) then