diff options
author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-01 22:22:31 +0000 |
---|---|---|
committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-01 22:22:31 +0000 |
commit | 14f73b9cd6d91a8cfae663006ffc245b28a2b725 (patch) | |
tree | b7291805b748fd10ca58dd8104cab2e79416b1ab /src | |
parent | bc12d89fb47eba72f114e8b331f5503d86fa155a (diff) |
RIT dev: Minor change in sign convention of an equation.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@289 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_Con2PrimM_pt_EOSOmni.c | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/GRHydro_Con2PrimM_pt_EOSOmni.c b/src/GRHydro_Con2PrimM_pt_EOSOmni.c index 9630e65..051db33 100644 --- a/src/GRHydro_Con2PrimM_pt_EOSOmni.c +++ b/src/GRHydro_Con2PrimM_pt_EOSOmni.c @@ -903,12 +903,12 @@ static void func_vsq_eosomni(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[], //Eq. 5: Qdotn + B^2(1+vsq)/2 - QdotBsq/2/W^2 + vsq W + rho_0(vsq) + u = 0 //rho0 = D * sqrt(1-vsq) - resid[1] = lgp->Qdotn + lgp->half_Bsq*(1.0+vsq) - 0.5*QB2Winv2 + vsq*W + rho0 + u; + resid[1] = -lgp->Qdotn - lgp->half_Bsq*(1.0+vsq) + 0.5*QB2Winv2 - vsq*W - rho0 - u; drho0_dv = -0.5*lgp->D / LorInv; - j10 = vsq+QB2Winv3; - j11 = lgp->half_Bsq + W + drho0_dv; + j10 = -vsq-QB2Winv3; + j11 = -lgp->half_Bsq - W - drho0_dv; //Eq. 6: u+ p - W(1-vsq) + rho0 = 0 => p-W = -W vsq - rho0 - u @@ -927,18 +927,18 @@ static void func_vsq_eosomni(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[], jac[0][2] = 0.0; jac[1][0] = j10; jac[1][1] = j11; - jac[1][2] = 1.0; + jac[1][2] = -1.0; jac[2][0] = vsq-1.0; jac[2][1] = dpress_dv + W + drho0_dv; jac[2][2] = dpress_du + 1.0; - c00 = j11*jac[2][2]-jac[2][1]; + c00 = j11*jac[2][2]+jac[2][1]; c01 = -1*j01*jac[2][2]; c02 = j01*jac[1][2]; - c10 = jac[2][0]-j10*jac[2][2]; + c10 = -jac[2][0]-j10*jac[2][2]; c11 = j00*jac[2][2]; - c12 = -j00; + c12 = j00; c20 = j10*jac[2][1]-j11*jac[2][0]; c21 = j01*jac[2][0]-j00*jac[2][1]; |