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 }
 119
 120
 121 int
 122 get_fpu_rounding_mode (void)
 123 {
 124   int rnd_mode;
 125
 126   rnd_mode = fegetround ();
 127
 128   switch (rnd_mode)
 129     {
 130 #ifdef FE_TONEAREST
 131       case FE_TONEAREST:
 132         return GFC_FPE_TONEAREST;
 133 #endif
 134
 135 #ifdef FE_UPWARD
 136       case FE_UPWARD:
 137         return GFC_FPE_UPWARD;
 138 #endif
 139
 140 #ifdef FE_DOWNWARD
 141       case FE_DOWNWARD:
 142         return GFC_FPE_DOWNWARD;
 143 #endif
 144
 145 #ifdef FE_TOWARDZERO
 146       case FE_TOWARDZERO:
 147         return GFC_FPE_TOWARDZERO;
 148 #endif
 149       default:
 150         return GFC_FPE_INVALID;
 151     }
 152 }
 153
 154
 155 void
 156 set_fpu_rounding_mode (int mode)
 157 {
 158   int rnd_mode;
 159
 160   switch (mode)
 161     {
 162 #ifdef FE_TONEAREST
 163       case GFC_FPE_TONEAREST:
 164         rnd_mode = FE_TONEAREST;
 165         break;
 166 #endif
 167
 168 #ifdef FE_UPWARD
 169       case GFC_FPE_UPWARD:
 170         rnd_mode = FE_UPWARD;
 171         break;
 172 #endif
 173
 174 #ifdef FE_DOWNWARD
 175       case GFC_FPE_DOWNWARD:
 176         rnd_mode = FE_DOWNWARD;
 177         break;
 178 #endif
 179
 180 #ifdef FE_TOWARDZERO
 181       case GFC_FPE_TOWARDZERO:
 182         rnd_mode = FE_TOWARDZERO;
 183         break;
 184 #endif
 185       default:
 186         return;
 187     }
 188
 189   fesetround (rnd_mode);
 190 }