[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 289)

[bcmsma at astro.rit.edu]

User: bmundim
Date: 2011/11/01 05:22 PM

Modified:
 /trunk/src/
  GRHydro_Con2PrimM_pt_EOSOmni.c

Log:
 RIT dev: Minor change in sign convention of an equation.

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Con2PrimM_pt_EOSOmni.c
Delta lines: +7 -7
===================================================================
--- trunk/src/GRHydro_Con2PrimM_pt_EOSOmni.c	2011-11-01 22:21:25 UTC (rev 288)
+++ trunk/src/GRHydro_Con2PrimM_pt_EOSOmni.c	2011-11-01 22:22:31 UTC (rev 289)
@@ -903,12 +903,12 @@
   //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 @@
   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];