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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2014/04/15 02:50 PM

Modified:
 /trunk/src/
  GRHydro_Prim2ConM.F90

Log:
 GRHydro: no longer compute speed of sound twice in F90 MHD code

File Changes:

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

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +0 -11
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2014-04-15 19:50:08 UTC (rev 640)
+++ trunk/src/GRHydro_Prim2ConM.F90	2014-04-15 19:50:11 UTC (rev 641)
@@ -85,17 +85,6 @@
     pBprim = loc(Bvec)
   end if
   
-  ! currently only the C++ code uses the speed-of-sound grid function, however
-  ! in the MHD code, both C++ and F90 code share the same prim2con routine (this
-  ! one). GRHydro_SpeedOfSound is a misnomer since it computes both cs2 and
-  ! pressure so we call it first so that the actual prim2con below overwrites
-  ! the pressur again. This is certainly not the "right" way of doing this but
-  ! it lets us have a working version for now.
-  ! TODO: extend prim2con in GRHydro_Prim2Con.c to handle MHD.
-  if (use_cxx_code .ne. 0) then
-    call GRHydro_SpeedOfSound(cctkGH)
-  end if
-
   ! constraint transport needs to be able to average fluxes in the directions
   ! other that flux_direction, which in turn need the primitives on interfaces