diff options
| author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-07-23 15:29:56 +0000 |
|---|---|---|
| committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-07-23 15:29:56 +0000 |
| commit | 10254df8c62dd286e5ba2ca68e76d993ed49889f (patch) | |
| tree | fc41d7bb1a291d6619cfdbffa557f5f35d3a4f39 /src/GRHydro_Eigenproblem_Marquina.F90 | |
| parent | c1454e6bc104a14c0f669389c5e90c1d5559c9bd (diff) | |
* remove incomplete and slow NaN checks.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@141 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Eigenproblem_Marquina.F90')
| -rw-r--r-- | src/GRHydro_Eigenproblem_Marquina.F90 | 30 |
1 files changed, 0 insertions, 30 deletions
diff --git a/src/GRHydro_Eigenproblem_Marquina.F90 b/src/GRHydro_Eigenproblem_Marquina.F90 index 806f952..ddd590e 100644 --- a/src/GRHydro_Eigenproblem_Marquina.F90 +++ b/src/GRHydro_Eigenproblem_Marquina.F90 @@ -61,8 +61,6 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& CCTK_REAL flux1l,flux2l,flux3l,flux4l,flux5l CCTK_REAL flux1,flux2,flux3,flux4,flux5 - character(len=256) NaN_WarnLine - !!$ LOCAL VARS CCTK_REAL du(5),aa(5,5),qdiff1,qdiff2,qdiff3,qdiff4,qdiff5 @@ -154,13 +152,6 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& wl = one / sqrt(one - v2l) -!!$ BEGIN: Check for NaN value (1st check) - if (wl .ne. wl) then - write(NaN_WarnLine,'(a100,2g15.6)') 'NaN produced: (one, v2l)', one, v2l - call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) - endif -!!$ END: Check for NaN value (1st check) - !!$ EIGENVALUES betainvalp = beta / alp @@ -175,13 +166,6 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& lamp_nobetal = (velxl*(one-cs2l) + sqrt(lamp_tmp))*invfacl lamm_nobetal = (velxl*(one-cs2l) - sqrt(lamp_tmp))*invfacl -!!$ BEGIN: Check for NaN value (2nd check) - if (lamp_nobetal .ne. lamp_nobetal) then - write(NaN_WarnLine,'(a150,5g15.6)') 'NaN produced: (cs2l, one, v2l, u, velxl)', cs2l, one, v2l, u, velxl - call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) - endif -!!$ END: Check for NaN value (2nd check) - lampl = lamp_nobetal - betainvalp lamml = lamm_nobetal - betainvalp @@ -220,13 +204,6 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& !!$ Assume consistent primitive data wr = one / sqrt(one - v2r) -!!$ BEGIN: Check for NaN value (3rd check) - if (wr .ne. wr) then - write(NaN_WarnLine,'(a100,2g15.6)') 'NaN produced: (one, v2r)', one, v2r - call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) - endif -!!$ END: Check for NaN value (3rd check) - !!$ EIGENVALUES lam1r = velxr - betainvalp @@ -240,13 +217,6 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& lamp_nobetar = (velxr*(one-cs2r) + sqrt(lamp_tmp))*invfacr lamm_nobetar = (velxr*(one-cs2r) - sqrt(lamp_tmp))*invfacr -!!$ BEGIN: Check for NaN value (4th check) - if (lamp_nobetar .ne. lamp_nobetar) then - write(NaN_WarnLine,'(a150,5g15.6)') 'NaN produced: (cs2r, one, v2r, u, velxr)', cs2r, one, v2r, u, velxr - call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine) - endif -!!$ END: Check for NaN value (4th check) - lampr = lamp_nobetar - betainvalp lammr = lamm_nobetar - betainvalp |