Merge pull request #260 from bonachea/ieee-inexact

Fix issue #259: prif_init signals IEEE_INEXACT

Author: bonachea

src/caffeine/program_startup_s.F90

@@ -7,6 +7,7 @@
 contains
 
   module procedure prif_init
+    use ieee_arithmetic, only: ieee_inexact, ieee_set_flag
     logical, save :: prif_init_called_previously = .false.
 
     if (prif_init_called_previously) then
@@ -27,6 +28,10 @@
 
        call sync_init()
 
+       ! issue #259: Ensure we clear any IEEE FP exceptions potentially
+       ! signalled from within the C-based initialization code above
+       call ieee_set_flag(ieee_inexact, .false.)
+
        prif_init_called_previously = .true.
        stat = 0
     end if