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

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