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

[reisswig at tapir.caltech.edu]

User: reisswig
Date: 2012/02/16 05:18 PM

Modified:
 /trunk/src/
  GRHydro_Eigenproblem.F90, GRHydro_PPM.F90

Log:
 Cosmetic change to ePPM3.
 
 Improved handling of cs2 < 0.

File Changes:

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

File [modified]: GRHydro_Eigenproblem.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_Eigenproblem.F90	2012-02-09 01:03:16 UTC (rev 313)
+++ trunk/src/GRHydro_Eigenproblem.F90	2012-02-16 23:18:45 UTC (rev 314)
@@ -87,7 +87,7 @@
 
   if(cs2.lt.0.0d0) then
      if (abs(cs2) .gt. 1.0d-4) then
-        write(warnline,'(a50,6g16.7)') 'rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho
+        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!")
      else

File [modified]: GRHydro_PPM.F90
Delta lines: +5 -5
===================================================================
--- trunk/src/GRHydro_PPM.F90	2012-02-09 01:03:16 UTC (rev 313)
+++ trunk/src/GRHydro_PPM.F90	2012-02-16 23:18:45 UTC (rev 314)
@@ -667,16 +667,16 @@
 		  aplus(i)  = a(i) &&\
 	       endif &&\
 	    endif &&\
-	 endif      &&\
+	 endif &&\
 	 trivial_rp(i-1) = .false.                                 &&\
 	 trivial_rp(i) = .false.                                   &&\
       else                                                         &&\
 	 if (abs(daplus) .le. abs(a(i)) .and. abs(daminus) .le. abs(a(i))) then &&\
 	    if (abs(daplus) .ge. 2.0d0*abs(daminus)) then     &&\
-	       aplus(i) = a(i) - 2.0d0*(aminus(i) - a(i))               &&\
+	       aplus(i) = a(i) + 2.0d0*daminus               &&\
 	    else if (abs(daminus) .ge. 2.0d0*abs(daplus)) then &&\
-	       aminus(i) = a(i) - 2.0d0*(aplus(i) - a(i))               &&\
-	    end if  &&\
+	       aminus(i) = a(i) - 2.0d0*daplus               &&\
+	    end if                                                    &&\
 	 else &&\
 	    if (abs(daplus) .ge. 2.0d0*abs(daminus)) then     &&\
 	       aplus(i) = a(i)               &&\
@@ -686,7 +686,7 @@
 	 endif &&\
 	 trivial_rp(i-1) = .false.  &&\
 	 trivial_rp(i) = .false.    &&\
-      endif
+      endif