From d73ce95f867e1a041cc51b82a954b6578e24610c Mon Sep 17 00:00:00 2001 From: rhaas Date: Sat, 11 May 2013 06:29:27 +0000 Subject: GRhydro: Pressure in polytypeM_pt was off by factor inv_sqrt_detg From: srichers git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@510 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Con2PrimM_polytype_pt.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/GRHydro_Con2PrimM_polytype_pt.c b/src/GRHydro_Con2PrimM_polytype_pt.c index 4931415..108d700 100644 --- a/src/GRHydro_Con2PrimM_polytype_pt.c +++ b/src/GRHydro_Con2PrimM_polytype_pt.c @@ -406,7 +406,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt) ( // printf("doublecheck - S, rho, gamma: %e %e %e\n",*sc_in, rho_gm1,gamma); - p = (*sc_in) * rho_gm1 / gamma; + p = (*sc_in) * inv_sqrt_detg * rho_gm1 / gamma; u = p / gm1; -- cgit v1.2.3