libgcc/config/rs6000/float128-hw.c

   1 /* Automatic switching between software and hardware IEEE 128-bit
   2    floating-point emulation for PowerPC.
   3
   4    Copyright (C) 2016-2024 Free Software Foundation, Inc.
   5    This file is part of the GNU C Library.
   6    Contributed by Michael Meissner (meissner@linux.vnet.ibm.com)
   7    Code is based on the main soft-fp library written by:
   8         Richard Henderson (rth@cygnus.com) and
   9         Jakub Jelinek (jj@ultra.linux.cz).
  10
  11    The GNU C Library is free software; you can redistribute it and/or
  12    modify it under the terms of the GNU Lesser General Public
  13    License as published by the Free Software Foundation; either
  14    version 2.1 of the License, or (at your option) any later version.
  15
  16    In addition to the permissions in the GNU Lesser General Public
  17    License, the Free Software Foundation gives you unlimited
  18    permission to link the compiled version of this file into
  19    combinations with other programs, and to distribute those
  20    combinations without any restriction coming from the use of this
  21    file.  (The Lesser General Public License restrictions do apply in
  22    other respects; for example, they cover modification of the file,
  23    and distribution when not linked into a combine executable.)
  24
  25    The GNU C Library is distributed in the hope that it will be useful,
  26    but WITHOUT ANY WARRANTY; without even the implied warranty of
  27    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  28    Lesser General Public License for more details.
  29
  30    You should have received a copy of the GNU Lesser General Public
  31    License along with the GNU C Library; if not, see
  32    <http://www.gnu.org/licenses/>.  */
  33
  34 #include <soft-fp.h>
  35 #include <quad-float128.h>
  36
  37 #ifndef __FLOAT128_HARDWARE__
  38 #error "This module must be compiled with IEEE 128-bit hardware support"
  39 #endif
  40
  41 TFtype
  42 __addkf3_hw (TFtype a, TFtype b)
  43 {
  44   return a + b;
  45 }
  46
  47 TFtype
  48 __subkf3_hw (TFtype a, TFtype b)
  49 {
  50   return a - b;
  51 }
  52
  53 TFtype
  54 __mulkf3_hw (TFtype a, TFtype b)
  55 {
  56   return a * b;
  57 }
  58
  59 TFtype
  60 __divkf3_hw (TFtype a, TFtype b)
  61 {
  62   return a / b;
  63 }
  64
  65 TFtype
  66 __negkf2_hw (TFtype a)
  67 {
  68   return -a;
  69 }
  70
  71 TFtype
  72 __floatsikf_hw (SItype_ppc a)
  73 {
  74   return (TFtype) a;
  75 }
  76
  77 TFtype
  78 __floatunsikf_hw (USItype_ppc a)
  79 {
  80   return (TFtype) a;
  81 }
  82
  83 TFtype
  84 __floatdikf_hw (DItype_ppc a)
  85 {
  86   return (TFtype) a;
  87 }
  88
  89 TFtype
  90 __floatundikf_hw (UDItype_ppc a)
  91 {
  92   return (TFtype) a;
  93 }
  94
  95 SItype_ppc
  96 __fixkfsi_hw (TFtype a)
  97 {
  98   return (SItype_ppc) a;
  99 }
 100
 101 USItype_ppc
 102 __fixunskfsi_hw (TFtype a)
 103 {
 104   return (USItype_ppc) a;
 105 }
 106
 107 DItype_ppc
 108 __fixkfdi_hw (TFtype a)
 109 {
 110   return (DItype_ppc) a;
 111 }
 112
 113 UDItype_ppc
 114 __fixunskfdi_hw (TFtype a)
 115 {
 116   return (UDItype_ppc) a;
 117 }
 118
 119 TFtype
 120 __extendsfkf2_hw (float a)
 121 {
 122   return (TFtype) a;
 123 }
 124
 125 TFtype
 126 __extenddfkf2_hw (double a)
 127 {
 128   return (TFtype) a;
 129 }
 130
 131 float
 132 __trunckfsf2_hw (TFtype a)
 133 {
 134   return (float) a;
 135 }
 136
 137 double
 138 __trunckfdf2_hw (TFtype a)
 139 {
 140   return (double) a;
 141 }
 142
 143 /* __eqkf2 returns 0 if equal, or 1 if not equal or NaN.  */
 144 CMPtype
 145 __eqkf2_hw (TFtype a, TFtype b)
 146 {
 147   return (a != b);
 148 }
 149
 150 /* __gekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or -2 if NaN.  */
 151 CMPtype
 152 __gekf2_hw (TFtype a, TFtype b)
 153 {
 154   if (a < b)
 155     return -1;
 156
 157   else if (__builtin_isunordered (a, b))
 158     return -2;
 159
 160   else if (a == b)
 161     return 0;
 162
 163   return 1;
 164 }
 165
 166 /* __lekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or +2 if NaN.  */
 167 CMPtype
 168 __lekf2_hw (TFtype a, TFtype b)
 169 {
 170   if (a < b)
 171     return -1;
 172
 173   else if (__builtin_isunordered (a, b))
 174     return 2;
 175
 176   else if (a == b)
 177     return 0;
 178
 179   return 1;
 180 }
 181
 182 /* __unordkf2 returns 1 if NaN or 0 otherwise.  */
 183 CMPtype
 184 __unordkf2_hw (TFtype a, TFtype b)
 185 {
 186   return (__builtin_isunordered (a, b)) ? 1 : 0;
 187 }
 188
 189 /* Convert __float128 to __ibm128.  */
 190 IBM128_TYPE
 191 __extendkftf2_hw (TFtype value)
 192 {
 193   IBM128_TYPE ret;
 194
 195   CVT_FLOAT128_TO_IBM128 (ret, value);
 196   return ret;
 197 }
 198
 199 /* Convert __ibm128 to __float128.  */
 200 TFtype
 201 __trunctfkf2_hw (IBM128_TYPE value)
 202 {
 203   TFtype ret;
 204
 205   CVT_IBM128_TO_FLOAT128 (ret, value);
 206   return ret;
 207 }