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

1 files changed, 6 insertions, 0 deletions

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;