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

knarf at cct.lsu.edu knarf at cct.lsu.edu
Tue Apr 13 17:40:33 CDT 2010


User: knarf
Date: 2010/04/13 05:40 PM

Modified:
 /trunk/src/
  Whisky_Con2Prim.F90, Whisky_ENOReconstruct.F90, Whisky_PPM.F90, Whisky_TVDReconstruct.F90, Whisky_UpdateMask.F90

Log:
 correct logical/integer expressions

File Changes:

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

File [modified]: Whisky_Con2Prim.F90
Delta lines: +8 -8
===================================================================
--- trunk/src/Whisky_Con2Prim.F90	2010-04-11 03:32:45 UTC (rev 105)
+++ trunk/src/Whisky_Con2Prim.F90	2010-04-13 22:40:33 UTC (rev 106)
@@ -604,7 +604,7 @@
 
         !do not compute if in atmosphere or in an excised region
         if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-            whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+            whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
          
         gxxl = 0.5d0 * (gxx(i,j,k) + gxx(i-xoffset,j-yoffset,k-zoffset))
         gxyl = 0.5d0 * (gxy(i,j,k) + gxy(i-xoffset,j-yoffset,k-zoffset))
@@ -818,7 +818,7 @@
 
         !do not compute if in atmosphere or in an excised region
         if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-            whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+            whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
 
         if (conformal_state .eq. 0) then
           call SpatialDeterminant(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),&
@@ -1146,7 +1146,7 @@
         
         !do not compute if in atmosphere or in an excised region
         if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-            whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+            whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
         
         gxxl = 0.5d0 * (gxx(i,j,k) + gxx(i-xoffset,j-yoffset,k-zoffset))
         gxyl = 0.5d0 * (gxy(i,j,k) + gxy(i-xoffset,j-yoffset,k-zoffset))
@@ -1629,7 +1629,7 @@
 
           !do not compute if in atmosphere or in an excised region
           if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-             whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) &
+             whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) &
                then
              press_new(i,j,k) = press(i,j,k)
              press_old(i,j,k) = press(i,j,k)
@@ -1733,7 +1733,7 @@
 
           !do not compute if in atmosphere or in an excised region
           if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-              whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+              whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
  
           f = press_new(i,j,k) - press(i,j,k)
 
@@ -1787,7 +1787,7 @@
         
         !do not compute if in atmosphere or in an excised region
         if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-            whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+            whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
 
 !!$         Note that in the following we enforce eps > 0
 !!$         This makes the later check redundant
@@ -2017,7 +2017,7 @@
 
           !do not compute if in atmosphere or in an excised region
           if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-              whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+              whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
 
           enthalpy = 1.d0 + eps(i,j,k) + press(i,j,k) / press_new(i,j,k)
 
@@ -2070,7 +2070,7 @@
 
         !do not compute if in atmosphere or in an excised region
         if (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, type_bits, atmosphere) .or. &
-            whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
+            whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0)) cycle
  
         udens = dens(i,j,k) / sqrt(Whisky_det(i,j,k))
 

File [modified]: Whisky_ENOReconstruct.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/Whisky_ENOReconstruct.F90	2010-04-11 03:32:45 UTC (rev 105)
+++ trunk/src/Whisky_ENOReconstruct.F90	2010-04-13 22:40:33 UTC (rev 106)
@@ -181,7 +181,7 @@
 
 !!$    Initialize excision
   do i = 1, nx
-    if (whisky_enable_internal_excision .and. (hydro_excision_mask(i) .ne. 0)) then
+    if (whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i) .ne. 0)) then
       vdiff(1, i) = large * (-1.d0*(1.d0+mod(i,10)))**i
       trivial_rp(i) = .true.
       excise(i) = .true.

File [modified]: Whisky_PPM.F90
Delta lines: +5 -7
===================================================================
--- trunk/src/Whisky_PPM.F90	2010-04-11 03:32:45 UTC (rev 105)
+++ trunk/src/Whisky_PPM.F90	2010-04-13 22:40:33 UTC (rev 106)
@@ -430,9 +430,8 @@
   end if
 
   !!$ excision
- if (hydro_excision) then
   do i = 1, nx
-    if (whisky_enable_internal_excision .and. &
+    if (whisky_enable_internal_excision /= 0 .and. &
         (hydro_excision_mask(i) .ne. 0)) then
       if (i .gt. 1) then
         trivial_rp(i-1)=.true.
@@ -443,7 +442,7 @@
       !!$  have to follow the order of sub-expressions given here and might
       !!$  access outside the array range
       cond = i .gt. 1
-      if (cond .and. whisky_enable_internal_excision) then
+      if (cond .and. whisky_enable_internal_excision /= 0) then
         cond = hydro_excision_mask(i-1) .ne. 0
       end if
       if (cond) then
@@ -471,7 +470,7 @@
         end if
       else
         cond = (i.gt.2) .and. (i.lt.nx)
-        if (cond .and. whisky_enable_internal_excision) then
+        if (cond .and. whisky_enable_internal_excision /= 0) then
             cond = (ppm_mppm .eq. 0) .and. (hydro_excision_mask(i-2) .ne. 0)
         end if
         if (cond) then
@@ -485,7 +484,7 @@
         end if
       end if
       cond = i.lt.nx
-      if (cond .and. whisky_enable_internal_excision) then
+      if (cond .and. whisky_enable_internal_excision /= 0) then
           cond = hydro_excision_mask(i+1) .ne. 0
       end if
       if (cond) then
@@ -513,7 +512,7 @@
         endif
       else
         cond = (i.lt.nx-1) .and. (i.gt.1)
-        if (cond .and. whisky_enable_internal_excision) then
+        if (cond .and. whisky_enable_internal_excision /= 0) then
             cond = (ppm_mppm .eq. 0) .and. (hydro_excision_mask(i+2) .ne. 0)
         end if
         if (cond) then
@@ -528,7 +527,6 @@
       end if
     end if
   end do
- end if
 return
 
 end subroutine SimplePPM_1d

File [modified]: Whisky_TVDReconstruct.F90
Delta lines: +3 -3
===================================================================
--- trunk/src/Whisky_TVDReconstruct.F90	2010-04-11 03:32:45 UTC (rev 105)
+++ trunk/src/Whisky_TVDReconstruct.F90	2010-04-13 22:40:33 UTC (rev 106)
@@ -56,18 +56,18 @@
     do k = whisky_stencil, nz-whisky_stencil+1
       do j = whisky_stencil, ny-whisky_stencil+1
         do i = whisky_stencil, nx-whisky_stencil+1
-          if (whisky_enable_internal_excision .and. &
+          if (whisky_enable_internal_excision /= 0 .and. &
               (hydro_excision_mask(i,j,k) .ne. 0)) then
             trivial_rp(i-xoffset, j-yoffset, k-zoffset) = .true.
             trivial_rp(i, j, k) = .true.
             bextm(i, j, k) = orig(i, j, k)
             bextp(i, j, k) = orig(i, j, k)
-            if (whisky_enable_internal_excision .and. &
+            if (whisky_enable_internal_excision /= 0 .and. &
                 (hydro_excision_mask(i+xoffset,j+yoffset,k+zoffset) .eq. 0)) then
               bextm(i, j, k) = orig(i+xoffset, j+yoffset, k+zoffset)
               bextp(i, j, k) = orig(i+xoffset, j+yoffset, k+zoffset)
             end if
-          else if (whisky_enable_internal_excision .and. &
+          else if (whisky_enable_internal_excision /= 0 .and. &
                    ((hydro_excision_mask(i-xoffset,j-yoffset,k-zoffset) .ne. 0) .or. &
                     (hydro_excision_mask(i+xoffset,j+yoffset,k+zoffset) .ne. 0))) then
             bextm(i, j, k) = orig(i, j, k)

File [modified]: Whisky_UpdateMask.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/Whisky_UpdateMask.F90	2010-04-11 03:32:45 UTC (rev 105)
+++ trunk/src/Whisky_UpdateMask.F90	2010-04-13 22:40:33 UTC (rev 106)
@@ -38,7 +38,7 @@
   do k = 1, cctk_lsh(3)
     do j = 1, cctk_lsh(2)
       do i = 1, cctk_lsh(1)
-        if ( whisky_enable_internal_excision .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. &
+        if ( whisky_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. &
              (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, \
                                           type_bits, atmosphere)) .or. &
              (tau(i,j,k) + frac * taurhs(i,j,k) .le. 0.d0) .or. &



More information about the Commits mailing list