From c493c3fd7b27c831678b67077d4b5a830f4f5ee0 Mon Sep 17 00:00:00 2001 From: rhaas Date: Tue, 15 Apr 2014 19:50:11 +0000 Subject: GRHydro: no longer compute speed of sound twice in F90 MHD code git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@641 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Prim2ConM.F90 | 11 ----------- 1 file changed, 11 deletions(-) diff --git a/src/GRHydro_Prim2ConM.F90 b/src/GRHydro_Prim2ConM.F90 index e330ef0..978fede 100644 --- a/src/GRHydro_Prim2ConM.F90 +++ b/src/GRHydro_Prim2ConM.F90 @@ -85,17 +85,6 @@ subroutine primitive2conservativeM(CCTK_ARGUMENTS) 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 -- cgit v1.2.3