diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-11 06:29:27 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-11 06:29:27 +0000 |
commit | d73ce95f867e1a041cc51b82a954b6578e24610c (patch) | |
tree | 7ddc171f63a7bed8ba60fed9cfddee787d454029 | |
parent | c0c3fc9a92ff79ec8f5fa30d362a05924a31f585 (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.c | 2 |
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; |