libgfortran/config/fpu-glibc.h

   1 /* FPU-related code for systems with GNU libc.
   2    Copyright (C) 2005-2013 Free Software Foundation, Inc.
   3    Contributed by Francois-Xavier Coudert <coudert@clipper.ens.fr>
   4
   5 This file is part of the GNU Fortran runtime library (libgfortran).
   6
   7 Libgfortran is free software; you can redistribute it and/or
   8 modify it under the terms of the GNU General Public
   9 License as published by the Free Software Foundation; either
  10 version 3 of the License, or (at your option) any later version.
  11
  12 Libgfortran is distributed in the hope that it will be useful,
  13 but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15 GNU General Public License for more details.
  16
  17 Under Section 7 of GPL version 3, you are granted additional
  18 permissions described in the GCC Runtime Library Exception, version
  19 3.1, as published by the Free Software Foundation.
  20
  21 You should have received a copy of the GNU General Public License and
  22 a copy of the GCC Runtime Library Exception along with this program;
  23 see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
  24 <http://www.gnu.org/licenses/>.  */
  25
  26 /* FPU-related code for systems with the GNU libc, providing the
  27    feenableexcept function in fenv.h to set individual exceptions
  28    (there's nothing to do that in C99).  */
  29
  30 #ifdef HAVE_FENV_H
  31 #include <fenv.h>
  32 #endif
  33
  34 void set_fpu (void)
  35 {
  36   if (FE_ALL_EXCEPT != 0)
  37     fedisableexcept (FE_ALL_EXCEPT);
  38
  39   if (options.fpe & GFC_FPE_INVALID)
  40 #ifdef FE_INVALID
  41     feenableexcept (FE_INVALID);
  42 #else
  43     estr_write ("Fortran runtime warning: IEEE 'invalid operation' "
  44                 "exception not supported.\n");
  45 #endif
  46
  47 /* glibc does never have a FE_DENORMAL.  */
  48   if (options.fpe & GFC_FPE_DENORMAL)
  49 #ifdef FE_DENORMAL
  50     feenableexcept (FE_DENORMAL);
  51 #else
  52     estr_write ("Fortran runtime warning: Floating point 'denormal operand' "
  53                 "exception not supported.\n");
  54 #endif
  55
  56   if (options.fpe & GFC_FPE_ZERO)
  57 #ifdef FE_DIVBYZERO
  58     feenableexcept (FE_DIVBYZERO);
  59 #else
  60     estr_write ("Fortran runtime warning: IEEE 'division by zero' "
  61                 "exception not supported.\n");
  62 #endif
  63
  64   if (options.fpe & GFC_FPE_OVERFLOW)
  65 #ifdef FE_OVERFLOW
  66     feenableexcept (FE_OVERFLOW);
  67 #else
  68     estr_write ("Fortran runtime warning: IEEE 'overflow' "
  69                 "exception not supported.\n");
  70 #endif
  71
  72   if (options.fpe & GFC_FPE_UNDERFLOW)
  73 #ifdef FE_UNDERFLOW
  74     feenableexcept (FE_UNDERFLOW);
  75 #else
  76     estr_write ("Fortran runtime warning: IEEE 'underflow' "
  77                 "exception not supported.\n");
  78 #endif
  79
  80   if (options.fpe & GFC_FPE_INEXACT)
  81 #ifdef FE_INEXACT
  82     feenableexcept (FE_INEXACT);
  83 #else
  84     estr_write ("Fortran runtime warning: IEEE 'inexact' "
  85                 "exception not supported.\n");
  86 #endif
  87 }
  88
  89
  90 int
  91 get_fpu_except_flags (void)
  92 {
  93   int result, set_excepts;
  94
  95   result = 0;
  96   set_excepts = fetestexcept (FE_ALL_EXCEPT);
  97
  98 #ifdef FE_INVALID
  99   if (set_excepts & FE_INVALID)
 100     result |= GFC_FPE_INVALID;
 101 #endif
 102
 103 #ifdef FE_DIVBYZERO
 104   if (set_excepts & FE_DIVBYZERO)
 105     result |= GFC_FPE_ZERO;
 106 #endif
 107
 108 #ifdef FE_OVERFLOW
 109   if (set_excepts & FE_OVERFLOW)
 110     result |= GFC_FPE_OVERFLOW;
 111 #endif
 112
 113 #ifdef FE_UNDERFLOW
 114   if (set_excepts & FE_UNDERFLOW)
 115     result |= GFC_FPE_UNDERFLOW;
 116 #endif
 117
 118 #ifdef FE_DENORMAL
 119   if (set_excepts & FE_DENORMAL)
 120     result |= GFC_FPE_DENORMAL;
 121 #endif
 122
 123 #ifdef FE_INEXACT
 124   if (set_excepts & FE_INEXACT)
 125     result |= GFC_FPE_INEXACT;
 126 #endif
 127
 128   return result;
 129 }