diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-10-25 16:37:02 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-10-25 16:37:02 +0000 |
| commit | 72a0532057c4270ab82157eb257aa9286ceb2f11 (patch) | |
| tree | 64a2b6ead0ae4e6d8c83f993cee1a97479cd00b8 /src/GRHydro_Con2PrimM_pt_EOSOmni.c | |
| parent | fc7d0ec3fb8c96476bff5e025255818eaf44ed1a (diff) | |
GRHydro: Unify function interface for GRHydro_Con2PrimM_pt.c and GRHydro_Con2PrimM_pt_EOSOmni.c
* Interestingly Cactus doesn't complain about the function interface in
interface.ccl. It has the old set of arguments for that function and declare it
as written in Fortran instead of C.
* Cylindrical explosion works again with GRHydro_Con2PrimM_pt.c routine. You
can just swap both routines in make.code.defn
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@427 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Con2PrimM_pt_EOSOmni.c')
| -rw-r--r-- | src/GRHydro_Con2PrimM_pt_EOSOmni.c | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/GRHydro_Con2PrimM_pt_EOSOmni.c b/src/GRHydro_Con2PrimM_pt_EOSOmni.c index 4a52414..8a72a6b 100644 --- a/src/GRHydro_Con2PrimM_pt_EOSOmni.c +++ b/src/GRHydro_Con2PrimM_pt_EOSOmni.c @@ -369,7 +369,12 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( CCTK_REAL dum1,dum2; p = pressure_rho0_eps_eosomni(rho0,*epsilon,&dum1,&dum2,&eosvars) ; // EOSOMNI - // p = pressure_rho0_u(rho0,u,gammaeos) ; // EOS + CCTK_REAL p2 = pressure_rho0_u(rho0,u,gammaeos) ; // EOS + + if (fabs(p-p2) > 1.0e-8){ + printf("p = %20.16g, p2 = %20.16g, delta pressure = %20.16g \n",p, p2, p-p2); + } + w = rho0 + u + p ; W_last = w*gammasq ; @@ -406,6 +411,9 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( if( ((*retval) != 0.) || (W == FAIL_VAL) ) { *retval = *retval*100.+1.; fprintf(stdout," *retval = %26.16e \n", *retval ); + fprintf(stdout," x_3d[%d] = %26.16e \n", 0, x_3d[0] ); + fprintf(stdout," x_3d[%d] = %26.16e \n", 1, x_3d[1] ); + fprintf(stdout," x_3d[%d] = %26.16e \n", 2, x_3d[2] ); return; } else{ |