From 0a39384311781e47154eb333d5e74dacab8297da Mon Sep 17 00:00:00 2001 From: rhaas Date: Tue, 15 Apr 2014 19:50:06 +0000 Subject: GRhydro: detect EOS errors in MHD Con2Prim pressure call this is a cut down version of 3e18d28da06966fddcaf0612afa1006ceb027d18 which also contains vector potential changes. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@639 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Con2PrimM_pt_EOSOmni.c | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/src/GRHydro_Con2PrimM_pt_EOSOmni.c b/src/GRHydro_Con2PrimM_pt_EOSOmni.c index 1fc555b..bcaa1dd 100644 --- a/src/GRHydro_Con2PrimM_pt_EOSOmni.c +++ b/src/GRHydro_Con2PrimM_pt_EOSOmni.c @@ -397,6 +397,12 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_pt) ( // eosvars.eosanyerr[0]=5; p = pressure_rho0_eps_eosomni(rho0,*epsilon,&dum1,&dum2,&eosvars) ; // EOSOMNI // p = pressure_rho0_u(rho0,u,gammaeos) ; // EOS + // error handling + if (eosvars.eosanyerr[0] != 0) { + if(eosvars.eoskeyerr[0] == 668) { + *retval = 3.0; + } + } w = rho0 + u + p; -- cgit v1.2.3