diff options
author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:48:59 +0000 |
---|---|---|
committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:48:59 +0000 |
commit | 17c962fe5711750f28560798f4213988eb470e9e (patch) | |
tree | 5d73c44a1a5ef75b2c0d395e03c8bd94f0224328 | |
parent | 7964bfddef9723adbbf7b10dd1163a806668127e (diff) |
GRHydro_Init_Data: compute w_lorentz from uT and alp instead of vel
later on prim2con will of course compute the grid function using vel
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@201 ac85fae7-cede-4708-beff-ae01c7fa1c26
-rw-r--r-- | src/GRHydro_BondiM.c | 16 |
1 files changed, 4 insertions, 12 deletions
diff --git a/src/GRHydro_BondiM.c b/src/GRHydro_BondiM.c index 0b2e6f7..2fd4184 100644 --- a/src/GRHydro_BondiM.c +++ b/src/GRHydro_BondiM.c @@ -1225,21 +1225,13 @@ static void GRHydro_BondiM_Internal(CCTK_ARGUMENTS, CCTK_REAL range_min, CCTK_RE vely(i) = (ucon[YY]/ucon[TT] + betay[i]) * inv_alp; velz(i) = (ucon[ZZ]/ucon[TT] + betaz[i]) * inv_alp; - w_lorentz[i] = 1./sqrt(1-(gxx[i]*SQR(velx(i))+gyy[i]*SQR(vely(i))+gzz[i]*SQR(velz(i))+ - 2*(gxy[i]*velx(i)*vely(i)+gxz[i]*velx(i)*velz(i)+gyz[i]*vely(i)*velz(i))) ); - if(isnan(w_lorentz[i])) { - CCTK_VWarn(CCTK_WARN_ALERT, __LINE__, __FILE__, CCTK_THORNSTRING, - "nan in w_lorentz at (%g,%g,%g): v2 = %g", - x[i],y[i],z[i],(gxx[i]*SQR(velx(i))+gyy[i]*SQR(vely(i))+gzz[i]*SQR(velz(i))+ - 2*(gxy[i]*velx(i)*vely(i)+gxz[i]*velx(i)*velz(i)+gyz[i]*vely(i)*velz(i)))); - } - SpatialDet(gxx[i],gxy[i],gxz[i],gyy[i],gyz[i],gzz[i],&det); if(bvec_method == BVEC_TRANSFORM) { - Bvecx(i) = w_lorentz[i]*bcon[XX] - alp[i]*bcon[TT]*ucon[XX]; - Bvecy(i) = w_lorentz[i]*bcon[YY] - alp[i]*bcon[TT]*ucon[YY]; - Bvecz(i) = w_lorentz[i]*bcon[ZZ] - alp[i]*bcon[TT]*ucon[ZZ]; + const CCTK_REAL w = ucon[TT] * alp[i]; + Bvecx(i) = w*bcon[XX] - alp[i]*bcon[TT]*ucon[XX]; + Bvecy(i) = w*bcon[YY] - alp[i]*bcon[TT]*ucon[YY]; + Bvecz(i) = w*bcon[ZZ] - alp[i]*bcon[TT]*ucon[ZZ]; } else { Bvecx(i) = bondi_bmag*SQR(M)*x[i]/sqrt(det)/CUBE(r[i]); Bvecy(i) = bondi_bmag*SQR(M)*y[i]/sqrt(det)/CUBE(r[i]); |