From e14b402798458759073298a0c31a7bd5f18e7774 Mon Sep 17 00:00:00 2001 From: bmundim Date: Fri, 31 Dec 2010 19:07:40 +0000 Subject: 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 --- interface.ccl | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) (limited to 'interface.ccl') 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 -- cgit v1.2.3