User: rhaas Date: 2013/05/13 12:27 AM

Modified: /branches/ET_2013_05/ /branches/ET_2013_05/src/ GRHydro_Con2PrimM_polytype_pt.c

Log: GRhydro: Pressure in polytypeM_pt was off by factor inv_sqrt_detg

Directory Changes:

Directory: /branches/ET_2013_05/ ================================

Property changes on: branches/ET_2013_05 ___________________________________________________________________ Directory: /svn:mergeinfo/ ==========================

- /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 + /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510

File Changes:

Directory: /branches/ET_2013_05/src/ ====================================

File [modified]: GRHydro_Con2PrimM_polytype_pt.c Delta lines: +1 -1 =================================================================== --- branches/ET_2013_05/src/GRHydro_Con2PrimM_polytype_pt.c 2013-05-11 06:29:37 UTC (rev 512) +++ branches/ET_2013_05/src/GRHydro_Con2PrimM_polytype_pt.c 2013-05-13 05:27:57 UTC (rev 513) @@ -406,7 +406,7 @@

// 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;