diff options
author | bmundim <bmundim@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2010-12-31 19:07:40 +0000 |
---|---|---|
committer | bmundim <bmundim@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2010-12-31 19:07:40 +0000 |
commit | e14b402798458759073298a0c31a7bd5f18e7774 (patch) | |
tree | 165f8e4fdda1ff7a1b3dca94962aa86f19d09dfa /interface.ccl | |
parent | 85cb2e6f9c6191fc9889a1055017236c9a071192 (diff) |
RIT MHD development:
Add #include "cctk_Functions.h" and DECLARE_CCTK_FUNCTIONS
to GRHydro_C2P2CM.F90.
Solve several issues of rank mismatch in the EOS Omni
interface.
Initialize the velocities at GRHydro_C2P2CM.F90 since
GRHydro_Con2PrimM_pt.c uses past values as initial guess
for its Newton-Raphson procedure. These unitialized values
were causing C2P2CM to fail from time to time, since it
would access memory positions with very large numerical
values, i.e. v>>1. C2P2CM test always passes now.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@118 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'interface.ccl')
-rw-r--r-- | interface.ccl | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/interface.ccl b/interface.ccl index b8f7af7..4f3f7b5 100644 --- a/interface.ccl +++ b/interface.ccl @@ -79,12 +79,29 @@ SUBROUTINE Con2PrimPoly(CCTK_INT IN handle, CCTK_REAL OUT dens, \ CCTK_INT IN GRHydro_reflevel, CCTK_REAL IN GRHydro_C2P_failed) +void FUNCTION Con2PrimGenM(CCTK_INT INOUT handle, CCTK_REAL INOUT gamma_eos, CCTK_REAL INOUT dens, \ + CCTK_REAL INOUT sx, CCTK_REAL INOUT sy, CCTK_REAL INOUT sz, \ + CCTK_REAL INOUT tau, CCTK_REAL INOUT rho, \ + CCTK_REAL INOUT velx, CCTK_REAL INOUT vely, CCTK_REAL INOUT velz, \ + CCTK_REAL INOUT epsilon, CCTK_REAL INOUT pressure, \ + CCTK_REAL INOUT w_lorentz, \ + CCTK_REAL INOUT gxx, CCTK_REAL INOUT gxy, CCTK_REAL INOUT gxz, \ + CCTK_REAL INOUT gyy, CCTK_REAL INOUT gyz, CCTK_REAL INOUT gzz, \ + CCTK_REAL INOUT uxx, CCTK_REAL INOUT uxy, CCTK_REAL INOUT uxz, \ + CCTK_REAL INOUT uyy, CCTK_REAL INOUT uyz, CCTK_REAL INOUT uzz, \ + CCTK_REAL INOUT det, \ + CCTK_REAL INOUT Bvecx, CCTK_REAL INOUT Bvecy, CCTK_REAL INOUT Bvecz, \ + CCTK_REAL INOUT Bvecsq, \ + CCTK_INT OUT epsnegative, \ + CCTK_REAL OUT retval) + USES FUNCTION SpatialDet USES FUNCTION Prim2ConPoly USES FUNCTION Prim2ConGen USES FUNCTION Prim2ConPolyM USES FUNCTION Prim2ConGenM USES FUNCTION Con2PrimPoly +USES FUNCTION Con2PrimGenM CCTK_INT FUNCTION EOS_Omni_GetHandle(CCTK_STRING IN name) USES FUNCTION EOS_Omni_GetHandle |