diff --git a/Processes/Sibyll/Decay.h b/Processes/Sibyll/Decay.h
index 93516c141ad13888e47b57d72d247c897c2ebdbf..4bf06c3ef022bb3744674469731c95544579742b 100644
--- a/Processes/Sibyll/Decay.h
+++ b/Processes/Sibyll/Decay.h
@@ -10,6 +10,8 @@
 
 #include <corsika/particles/ParticleProperties.h>
 
+#include <fenv.h>
+
 namespace corsika::process {
 
   namespace sibyll {
@@ -145,6 +147,9 @@ namespace corsika::process {
       void DoDecay(Particle& p, Stack& s) const {
         using corsika::geometry::Point;
         using namespace corsika::units::si;
+
+	feenableexcept(FE_INVALID);
+
         fCount++;
         SibStack ss;
         ss.Clear();
@@ -188,6 +193,8 @@ namespace corsika::process {
         }
         // empty sibyll stack
         ss.Clear();
+
+	fedisableexcept(FE_INVALID);
       }
     };
   } // namespace sibyll