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