aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2013-05-11 06:29:27 +0000
committerrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2013-05-11 06:29:27 +0000
commitd73ce95f867e1a041cc51b82a954b6578e24610c (patch)
tree7ddc171f63a7bed8ba60fed9cfddee787d454029
parentc0c3fc9a92ff79ec8f5fa30d362a05924a31f585 (diff)
GRhydro: Pressure in polytypeM_pt was off by factor inv_sqrt_detg
From: srichers <srichers@caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@510 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r--src/GRHydro_Con2PrimM_polytype_pt.c2
1 files changed, 1 insertions, 1 deletions
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;