From efa9ed11cdfc52ff82fb56a4c0a841ccca4a4836 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Ole=20Sch=C3=BCtt?= <ole.schuett@cp2k.org>
Date: Tue, 29 May 2018 21:55:18 +0000
Subject: [PATCH] pexsi: ignore floating point exceptions

svn-origin-rev: 18453
---
 src/pexsi_interface.F | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)

diff --git a/src/pexsi_interface.F b/src/pexsi_interface.F
index 80de68098b..640e1dba68 100644
--- a/src/pexsi_interface.F
+++ b/src/pexsi_interface.F
@@ -20,6 +20,11 @@ MODULE pexsi_interface
       f_ppexsi_plan_initialize, &
       f_ppexsi_retrieve_real_dft_matrix, &
       f_ppexsi_set_default_options
+#endif
+#if defined (__HAS_IEEE_EXCEPTIONS)
+   USE ieee_exceptions, ONLY: ieee_get_halting_mode, &
+      ieee_set_halting_mode, &
+      ieee_all
 #endif
    USE kinds, ONLY: int_8, &
       real_8
@@ -380,11 +385,27 @@ CONTAINS
 
 #ifdef  __LIBPEXSI
       INTEGER                                  :: handle, info
+#if defined (__HAS_IEEE_EXCEPTIONS)
+      LOGICAL, DIMENSION(5)                    :: halt
+#endif
 
       CALL timeset(routineN, handle)
+
+      ! Unfortuntatelly, some PEXSI kernels raise IEEE754 exceptions.
+      ! Therefore, we disable floating point traps temporarily.
+#if defined (__HAS_IEEE_EXCEPTIONS)
+      CALL ieee_get_halting_mode(IEEE_ALL, halt)
+      CALL ieee_set_halting_mode(IEEE_ALL, .FALSE.)
+#endif
+
       CALL f_ppexsi_dft_driver(plan, pexsi_options%options, numElectronExact, muPEXSI, &
                                numElectronPEXSI, muMinInertia, muMaxInertia, &
                                numTotalInertiaIter, numTotalPEXSIIter, info)
+
+#if defined (__HAS_IEEE_EXCEPTIONS)
+      CALL ieee_set_halting_mode(IEEE_ALL, halt)
+#endif
+
       IF (info .NE. 0) &
          CPABORT("Pexsi returned an error. Consider logPEXSI0 for details.")
       CALL timestop(handle)