[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 201)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Wed Mar 27 20:48:59 CDT 2013
User: rhaas
Date: 2013/03/27 08:48 PM
Modified:
/trunk/src/
GRHydro_BondiM.c
Log:
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 at tapir.caltech.edu>
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_BondiM.c
Delta lines: +4 -12
===================================================================
--- trunk/src/GRHydro_BondiM.c 2013-03-28 01:48:57 UTC (rev 200)
+++ trunk/src/GRHydro_BondiM.c 2013-03-28 01:48:59 UTC (rev 201)
@@ -1225,21 +1225,13 @@
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]);
More information about the Commits
mailing list