libgfortran/config/fpu-aix.h

   1 /* AIX FPU-related code.
   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
  27 /* FPU-related code for AIX.  */
  28 #ifdef HAVE_FPTRAP_H
  29 #include <fptrap.h>
  30 #endif
  31
  32 #ifdef HAVE_FPXCP_H
  33 #include <fpxcp.h>
  34 #endif
  35
  36 void
  37 set_fpu (void)
  38 {
  39   fptrap_t mode = 0;
  40
  41   if (options.fpe & GFC_FPE_INVALID)
  42 #ifdef TRP_INVALID
  43     mode |= TRP_INVALID;
  44 #else
  45     estr_write ("Fortran runtime warning: IEEE 'invalid operation' "
  46                 "exception not supported.\n");
  47 #endif
  48
  49   if (options.fpe & GFC_FPE_DENORMAL)
  50     estr_write ("Fortran runtime warning: Floating point 'denormal operand' "
  51                 "exception not supported.\n");
  52
  53   if (options.fpe & GFC_FPE_ZERO)
  54 #ifdef TRP_DIV_BY_ZERO
  55     mode |= TRP_DIV_BY_ZERO;
  56 #else
  57     estr_write ("Fortran runtime warning: IEEE 'division by zero' "
  58                 "exception not supported.\n");
  59 #endif
  60
  61   if (options.fpe & GFC_FPE_OVERFLOW)
  62 #ifdef TRP_OVERFLOW
  63     mode |= TRP_OVERFLOW;
  64 #else
  65     estr_write ("Fortran runtime warning: IEEE 'overflow' "
  66                 "exception not supported.\n");
  67 #endif
  68
  69   if (options.fpe & GFC_FPE_UNDERFLOW)
  70 #ifdef TRP_UNDERFLOW
  71     mode |= TRP_UNDERFLOW;
  72 #else
  73     estr_write ("Fortran runtime warning: IEEE 'underflow' "
  74                 "exception not supported.\n");
  75 #endif
  76
  77   if (options.fpe & GFC_FPE_INEXACT)
  78 #ifdef TRP_INEXACT
  79     mode |= TRP_INEXACT;
  80 #else
  81     estr_write ("Fortran runtime warning: IEEE 'inexact' "
  82                 "exception not supported.\n");
  83 #endif
  84
  85   fp_trap(FP_TRAP_SYNC);
  86   fp_enable(mode);
  87 }
  88
  89
  90 int
  91 get_fpu_except_flags (void)
  92 {
  93   int result, set_excepts;
  94
  95   result = 0;
  96
  97 #ifdef HAVE_FPXCP_H
  98   if (!fp_any_xcp ())
  99     return 0;
 100
 101   if (fp_invalid_op ())
 102     result |= GFC_FPE_INVALID;
 103
 104   if (fp_divbyzero ())
 105     result |= GFC_FPE_ZERO;
 106
 107   if (fp_overflow ())
 108     result |= GFC_FPE_OVERFLOW;
 109
 110   if (fp_underflow ())
 111     result |= GFC_FPE_UNDERFLOW;
 112
 113   if (fp_inexact ())
 114     result |= GFC_FPE_INEXACT;
 115 #endif
 116
 117   return result;
 118 }