From 10254df8c62dd286e5ba2ca68e76d993ed49889f Mon Sep 17 00:00:00 2001 From: cott Date: Fri, 23 Jul 2010 15:29:56 +0000 Subject: * remove incomplete and slow NaN checks. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@141 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Prim2Con.F90 | 23 ----------------------- 1 file changed, 23 deletions(-) (limited to 'src/GRHydro_Prim2Con.F90') 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 -- cgit v1.2.3