diff options
Diffstat (limited to 'src/GRHydro_Con2PrimM_pt_EOSOmni.c')
| -rw-r--r-- | src/GRHydro_Con2PrimM_pt_EOSOmni.c | 13 |
1 files changed, 6 insertions, 7 deletions
diff --git a/src/GRHydro_Con2PrimM_pt_EOSOmni.c b/src/GRHydro_Con2PrimM_pt_EOSOmni.c index 96e2bd8..3700b5b 100644 --- a/src/GRHydro_Con2PrimM_pt_EOSOmni.c +++ b/src/GRHydro_Con2PrimM_pt_EOSOmni.c @@ -141,7 +141,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( CCTK_REAL *gyy, CCTK_REAL *gyz, CCTK_REAL *gzz, CCTK_REAL *uxx, CCTK_REAL *uxy, CCTK_REAL *uxz, CCTK_REAL *uyy, CCTK_REAL *uyz, CCTK_REAL *uzz, - CCTK_REAL *det, + CCTK_REAL *sdet, CCTK_INT *epsnegative, CCTK_REAL *retval); @@ -225,7 +225,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( CCTK_REAL *gyy, CCTK_REAL *gyz, CCTK_REAL *gzz, CCTK_REAL *uxx, CCTK_REAL *uxy, CCTK_REAL *uxz, CCTK_REAL *uyy, CCTK_REAL *uyz, CCTK_REAL *uzz, - CCTK_REAL *det, + CCTK_REAL *sdet, CCTK_INT *epsnegative, CCTK_REAL *retval) @@ -237,8 +237,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( CCTK_REAL QdotB,VdotB; CCTK_REAL rho0,u,p,w,gammasq,gamma,gtmp,W_last,W,vsq; CCTK_REAL g_o_WBsq, QdB_o_W; - CCTK_REAL detg = (*det); - CCTK_REAL sqrt_detg = sqrt(detg); + CCTK_REAL sqrt_detg = *sdet; CCTK_REAL inv_sqrt_detg = 1./sqrt_detg; CCTK_INT i,j, i_increase ; CCTK_INT nf, nfudgemax, failwarnmode, failinfomode; @@ -313,7 +312,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( fprintf(stdout," *uyy = %26.16e \n", *uyy ); fprintf(stdout," *uyz = %26.16e \n", *uyz ); fprintf(stdout," *uzz = %26.16e \n", *uzz ); - fprintf(stdout," *det = %26.16e \n", *det ); + fprintf(stdout," *sdet = %26.16e \n", *sdet ); fprintf(stdout," *epsnegative = %10d \n", *epsnegative ); fprintf(stdout," *retval = %26.16e \n", *retval ); fflush(stdout); @@ -539,7 +538,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( fprintf(stdout," *uyy = %26.16e \n", *uyy ); fprintf(stdout," *uyz = %26.16e \n", *uyz ); fprintf(stdout," *uzz = %26.16e \n", *uzz ); - fprintf(stdout," *det = %26.16e \n", *det ); + fprintf(stdout," *sdet = %26.16e \n", *sdet ); fprintf(stdout," *epsnegative = %10d \n", *epsnegative ); fprintf(stdout," *retval = %26.16e \n", *retval ); fflush(stdout); @@ -723,7 +722,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( fprintf(stdout," *uyyo = %26.16e \n", *uyy ); fprintf(stdout," *uyzo = %26.16e \n", *uyz ); fprintf(stdout," *uzzo = %26.16e \n", *uzz ); - fprintf(stdout," *deto = %26.16e \n", *det ); + fprintf(stdout," *sdeto = %26.16e \n", *sdet ); fprintf(stdout," *epsnegativeo = %10d \n", *epsnegative ); |