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

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Wed Jun 20 11:27:10 CDT 2012


User: rhaas
Date: 2012/06/20 11:27 AM

Modified:
 /trunk/src/
  GRHydro_Eigenproblem.F90

Log:
 really improve handling of cs2 < 0 in the "standard" (non-MHD, non-hot EOS)
 eigenvalue routine.
 
 the previous patch with this same commit message (r359) has the wrong commit message.
 It should instead read:
 --8<--
 Introduce real-valued atmo mask that can be interpolated. Synchronize
 real-valued atmo mask and update integer-valued atmo mask based on synchronized
 value. This fix ensures that the results are independent of the number of
 processors used. It also gets rid of some extra primitive syncs that are
 otherwise necessary when sync_conserved_only = yes.
 
 Does not cause any new test failures (Roland).
 
 Patch by Christian Reisswig.
 --8<--
 
 Patch by Christian Ott.

File Changes:

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

File [modified]: GRHydro_Eigenproblem.F90
Delta lines: +5 -2
===================================================================
--- trunk/src/GRHydro_Eigenproblem.F90	2012-06-20 16:19:10 UTC (rev 360)
+++ trunk/src/GRHydro_Eigenproblem.F90	2012-06-20 16:27:09 UTC (rev 361)
@@ -86,13 +86,16 @@
        (1.0d0 + eps + press/rho)
 
   if(cs2.lt.0.0d0) then
+     !$OMP CRITICAL
      if (abs(cs2) .gt. 1.0d-4) then
         write(warnline,'(a50,6g16.7)') 'abs(cs2), rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho
         call CCTK_WARN(1,warnline)
-        call CCTK_WARN(0,"cs2 < 0! Check speed of sound calculation!")
+        call CCTK_WARN(1,"cs2 < 0! Check speed of sound calculation!")
+        cs2 = 0.0d0
      else
-        cs2 = 0
+        cs2 = 0.0d0
      endif
+     !$OMP END CRITICAL
   endif
 
 



More information about the Commits mailing list