123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207 |
- /* Automatic switching between software and hardware IEEE 128-bit
- floating-point emulation for PowerPC.
- Copyright (C) 2016-2022 Free Software Foundation, Inc.
- This file is part of the GNU C Library.
- Contributed by Michael Meissner (meissner@linux.vnet.ibm.com)
- Code is based on the main soft-fp library written by:
- Richard Henderson (rth@cygnus.com) and
- Jakub Jelinek (jj@ultra.linux.cz).
- The GNU C Library is free software; you can redistribute it and/or
- modify it under the terms of the GNU Lesser General Public
- License as published by the Free Software Foundation; either
- version 2.1 of the License, or (at your option) any later version.
- In addition to the permissions in the GNU Lesser General Public
- License, the Free Software Foundation gives you unlimited
- permission to link the compiled version of this file into
- combinations with other programs, and to distribute those
- combinations without any restriction coming from the use of this
- file. (The Lesser General Public License restrictions do apply in
- other respects; for example, they cover modification of the file,
- and distribution when not linked into a combine executable.)
- The GNU C Library is distributed in the hope that it will be useful,
- but WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- Lesser General Public License for more details.
- You should have received a copy of the GNU Lesser General Public
- License along with the GNU C Library; if not, see
- <http://www.gnu.org/licenses/>. */
- #include <soft-fp.h>
- #include <quad-float128.h>
- #ifndef __FLOAT128_HARDWARE__
- #error "This module must be compiled with IEEE 128-bit hardware support"
- #endif
- TFtype
- __addkf3_hw (TFtype a, TFtype b)
- {
- return a + b;
- }
- TFtype
- __subkf3_hw (TFtype a, TFtype b)
- {
- return a - b;
- }
- TFtype
- __mulkf3_hw (TFtype a, TFtype b)
- {
- return a * b;
- }
- TFtype
- __divkf3_hw (TFtype a, TFtype b)
- {
- return a / b;
- }
- TFtype
- __negkf2_hw (TFtype a)
- {
- return -a;
- }
- TFtype
- __floatsikf_hw (SItype_ppc a)
- {
- return (TFtype) a;
- }
- TFtype
- __floatunsikf_hw (USItype_ppc a)
- {
- return (TFtype) a;
- }
- TFtype
- __floatdikf_hw (DItype_ppc a)
- {
- return (TFtype) a;
- }
- TFtype
- __floatundikf_hw (UDItype_ppc a)
- {
- return (TFtype) a;
- }
- SItype_ppc
- __fixkfsi_hw (TFtype a)
- {
- return (SItype_ppc) a;
- }
- USItype_ppc
- __fixunskfsi_hw (TFtype a)
- {
- return (USItype_ppc) a;
- }
- DItype_ppc
- __fixkfdi_hw (TFtype a)
- {
- return (DItype_ppc) a;
- }
- UDItype_ppc
- __fixunskfdi_hw (TFtype a)
- {
- return (UDItype_ppc) a;
- }
- TFtype
- __extendsfkf2_hw (float a)
- {
- return (TFtype) a;
- }
- TFtype
- __extenddfkf2_hw (double a)
- {
- return (TFtype) a;
- }
- float
- __trunckfsf2_hw (TFtype a)
- {
- return (float) a;
- }
- double
- __trunckfdf2_hw (TFtype a)
- {
- return (double) a;
- }
- /* __eqkf2 returns 0 if equal, or 1 if not equal or NaN. */
- CMPtype
- __eqkf2_hw (TFtype a, TFtype b)
- {
- return (a != b);
- }
- /* __gekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or -2 if NaN. */
- CMPtype
- __gekf2_hw (TFtype a, TFtype b)
- {
- if (a < b)
- return -1;
- else if (__builtin_isunordered (a, b))
- return -2;
- else if (a == b)
- return 0;
- return 1;
- }
- /* __lekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or +2 if NaN. */
- CMPtype
- __lekf2_hw (TFtype a, TFtype b)
- {
- if (a < b)
- return -1;
- else if (__builtin_isunordered (a, b))
- return 2;
- else if (a == b)
- return 0;
- return 1;
- }
- /* __unordkf2 returns 1 if NaN or 0 otherwise. */
- CMPtype
- __unordkf2_hw (TFtype a, TFtype b)
- {
- return (__builtin_isunordered (a, b)) ? 1 : 0;
- }
- /* Convert __float128 to __ibm128. */
- IBM128_TYPE
- __extendkftf2_hw (TFtype value)
- {
- IBM128_TYPE ret;
- CVT_FLOAT128_TO_IBM128 (ret, value);
- return ret;
- }
- /* Convert __ibm128 to __float128. */
- TFtype
- __trunctfkf2_hw (IBM128_TYPE value)
- {
- TFtype ret;
- CVT_IBM128_TO_FLOAT128 (ret, value);
- return ret;
- }
|