* remove incomplete and slow NaN checks.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@141 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 0 insertions, 23 deletions

diff --git a/src/GRHydro_Prim2Con.F90 b/src/GRHydro_Prim2Con.F90
index 8c994a5..875fa50 100644
--- a/src/GRHydro_Prim2Con.F90
+++ b/src/GRHydro_Prim2Con.F90
@@ -128,7 +128,6 @@ subroutine prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, ddens, &
   CCTK_REAL :: ddens, dsx, dsy, dsz, dtau, drho, dvelx, dvely, dvelz,&
        deps, dpress, w, vlowx, vlowy, vlowz   
   CCTK_INT :: handle
-  character(len=256) NaN_WarnLine
 
 #include "EOS_Base.inc"
   
@@ -151,15 +150,6 @@ subroutine prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, ddens, &
        *dvelz*dvelz + 2*gxy*dvelx*dvely + 2*gxz*dvelx *dvelz + 2*gyz&
        *dvely*dvelz))  
 
-!!$ BEGIN: Check for NaN value
-  if (w .ne. w) then
-    !$OMP CRITICAL
-    write(NaN_WarnLine,'(a100,3g15.6)') 'NaN produced in sqrt(): (dvelx,dvely,dvelz)', dvelx, dvely, dvelz
-    call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine)
-    !$OMP END CRITICAL
-  endif
-!!$ END: Check for NaN value
-
 #if USE_EOS_OMNI
   call EOS_Omni_press(eoskey,keytemp,n,&
        drho,deps,xtemp,xye,dpress,keyerr,anyerr)  
@@ -316,14 +306,6 @@ subroutine Prim2ConservativePolytype(CCTK_ARGUMENTS)
              rhoplus(i,j,k),velxplus(i,j,k),velyplus(i,j,k),&
              velzplus(i,j,k),epsplus(i,j,k),pressplus(i,j,k),&
              w_lorentzplus(i,j,k)) 
-        if (densminus(i,j,k) .ne. densminus(i,j,k)) then
-        !$OMP CRITICAL
-          call CCTK_WARN(1, "NaN in densminus(i,j,k) (Prim2Con)")
-          write(warnline,'  (a27,3g16.7)') 'avg_detl, velx, rhominus: ', &
-          avg_detl, velx(i,j,k), rhominus(i,j,k)
-          call CCTK_WARN(1, warnline)
-        !$OMP END CRITICAL
-        endif
       end do
     end do
   end do
@@ -643,11 +625,6 @@ subroutine primitive2conservativegeneral(CCTK_ARGUMENTS)
              rhoplus(i,j,k),velxplus(i,j,k),velyplus(i,j,k),&
              velzplus(i,j,k),epsplus(i,j,k),pressplus(i,j,k),&
              w_lorentzplus(i,j,k))
-        if (densminus(i,j,k) .ne. densminus(i,j,k)) then
-        !$OMP CRITICAL
-          call CCTK_WARN(1, "NaN in densminus(i,j,k) (Prim2Con)")
-        !$OMP END CRITICAL
-        endif
       end do
     end do
   end do