1 /* 2 Copyright (c) 1991, 1993 3 The Regents of the University of California. All rights reserved. 4 c) UNIX System Laboratories, Inc. 5 All or some portions of this file are derived from material licensed 6 to the University of California by American Telephone and Telegraph 7 Co. or Unix System Laboratories, Inc. and are reproduced herein with 8 the permission of UNIX System Laboratories, Inc. 9 10 Redistribution and use in source and binary forms, with or without 11 modification, are permitted provided that the following conditions 12 are met: 13 1. Redistributions of source code must retain the above copyright 14 notice, this list of conditions and the following disclaimer. 15 2. Redistributions in binary form must reproduce the above copyright 16 notice, this list of conditions and the following disclaimer in the 17 documentation and/or other materials provided with the distribution. 18 3. Neither the name of the University nor the names of its contributors 19 may be used to endorse or promote products derived from this software 20 without specific prior written permission. 21 22 THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND 23 ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 24 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 25 ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE 26 FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 27 DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 28 OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 29 HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 30 LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 31 OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 32 SUCH DAMAGE. 33 */ 34 #ifndef _IEEE_FP_H_ 35 #define _IEEE_FP_H_ 36 37 #include "_ansi.h" 38 39 #include <machine/ieeefp.h> 40 #include <float.h> 41 #include <stdint.h> 42 43 _BEGIN_STD_C 44 45 #ifndef _LDBL_EQ_DBL 46 47 #ifndef LDBL_MANT_DIG 48 #error "LDBL_MANT_DIG not defined - should be found in float.h" 49 50 #elif LDBL_MANT_DIG == DBL_MANT_DIG 51 #error "double and long double are the same size but LDBL_EQ_DBL is not defined" 52 53 #elif LDBL_MANT_DIG == 53 54 /* This happens when doubles are 32-bits and long doubles are 64-bits. */ 55 #define EXT_EXPBITS 11 56 #define EXT_FRACHBITS 20 57 #define EXT_FRACLBITS 32 58 #define __ieee_ext_field_type unsigned long 59 60 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 65) && LDBL_MAX_EXP == 16384 61 /* 80-bit floats, as on x86 */ 62 #define EXT_EXPBITS 15 63 #define EXT_FRACHBITS 32 64 #define EXT_FRACLBITS 32 65 #define __ieee_ext_field_type unsigned int 66 67 #elif (LDBL_MANT_DIG == 113 || LDBL_MANT_DIG == 112) && LDBL_MAX_EXP == 16384 68 /* 128-bit IEEE floats, as on risc-v */ 69 #define EXT_EXPBITS 15 70 #define EXT_FRACHBITS 48 71 #define EXT_FRACLBITS 64 72 #define __ieee_ext_field_type unsigned long long 73 74 #elif LDBL_MANT_DIG == 106 && DBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 75 /* 128-bit double-double, as on powerpc */ 76 #define EXT_EXPBITS 11 77 #define EXT_FRACHBITS DBL_MANT_DIG 78 #define EXT_FRACLBITS DBL_MANT_DIG 79 #define EXT_FRACGAP EXT_EXPBITS 80 #define __ieee_ext_field_type unsigned long long 81 82 #else 83 #error Unsupported value for LDBL_MANT_DIG 84 #endif 85 86 #define EXT_EXP_INFNAN ((1 << EXT_EXPBITS) - 1) /* 32767 */ 87 #define EXT_EXP_BIAS ((1 << (EXT_EXPBITS - 1)) - 1) /* 16383 */ 88 #define EXT_FRACBITS (EXT_FRACLBITS + EXT_FRACHBITS) 89 90 typedef struct ieee_ext 91 { 92 __ieee_ext_field_type ext_fracl : EXT_FRACLBITS; 93 #ifdef EXT_FRACGAP 94 __ieee_ext_field_type ext_gap : EXT_FRACGAP; 95 #endif 96 __ieee_ext_field_type ext_frach : EXT_FRACHBITS; 97 __ieee_ext_field_type ext_exp : EXT_EXPBITS; 98 __ieee_ext_field_type ext_sign : 1; 99 } ieee_ext; 100 101 typedef union ieee_ext_u 102 { 103 long double extu_ld; 104 struct ieee_ext extu_ext; 105 } ieee_ext_u; 106 107 #endif /* ! _LDBL_EQ_DBL */ 108 109 110 /* FLOATING ROUNDING */ 111 112 typedef int fp_rnd; 113 #define FP_RN 0 /* Round to nearest */ 114 #define FP_RM 1 /* Round down */ 115 #define FP_RP 2 /* Round up */ 116 #define FP_RZ 3 /* Round to zero (trunate) */ 117 118 fp_rnd fpgetround (void); 119 fp_rnd fpsetround (fp_rnd); 120 121 /* EXCEPTIONS */ 122 123 typedef int fp_except; 124 #define FP_X_INV 0x10 /* Invalid operation */ 125 #define FP_X_DX 0x80 /* Divide by zero */ 126 #define FP_X_OFL 0x04 /* Overflow exception */ 127 #define FP_X_UFL 0x02 /* Underflow exception */ 128 #define FP_X_IMP 0x01 /* imprecise exception */ 129 130 fp_except fpgetmask (void); 131 fp_except fpsetmask (fp_except); 132 fp_except fpgetsticky (void); 133 fp_except fpsetsticky (fp_except); 134 135 /* INTEGER ROUNDING */ 136 137 typedef int fp_rdi; 138 #define FP_RDI_TOZ 0 /* Round to Zero */ 139 #define FP_RDI_RD 1 /* Follow float mode */ 140 141 fp_rdi fpgetroundtoi (void); 142 fp_rdi fpsetroundtoi (fp_rdi); 143 144 #define __IEEE_DBL_EXPBIAS 1023 145 #define __IEEE_FLT_EXPBIAS 127 146 147 #define __IEEE_DBL_EXPLEN 11 148 #define __IEEE_FLT_EXPLEN 8 149 150 151 #define __IEEE_DBL_FRACLEN (64 - (__IEEE_DBL_EXPLEN + 1)) 152 #define __IEEE_FLT_FRACLEN (32 - (__IEEE_FLT_EXPLEN + 1)) 153 154 #define __IEEE_DBL_MAXPOWTWO ((double)(1L << 32 - 2) * (1L << (32-11) - 32 + 1)) 155 #define __IEEE_FLT_MAXPOWTWO ((float)(1L << (32-8) - 1)) 156 157 #define __IEEE_DBL_NAN_EXP 0x7ff 158 #define __IEEE_FLT_NAN_EXP 0xff 159 160 #ifdef __ieeefp_isnanf 161 #define isnanf(x) __ieeefp_isnanf(x) 162 #endif 163 164 #ifdef __ieeefp_isinff 165 #define isinff(x) __ieeefp_isinff(x) 166 #endif 167 168 #ifdef __ieeefp_finitef 169 #define finitef(x) __ieeefp_finitef(x) 170 #endif 171 172 #ifdef _DOUBLE_IS_32BITS 173 #undef __IEEE_DBL_EXPBIAS 174 #define __IEEE_DBL_EXPBIAS __IEEE_FLT_EXPBIAS 175 176 #undef __IEEE_DBL_EXPLEN 177 #define __IEEE_DBL_EXPLEN __IEEE_FLT_EXPLEN 178 179 #undef __IEEE_DBL_FRACLEN 180 #define __IEEE_DBL_FRACLEN __IEEE_FLT_FRACLEN 181 182 #undef __IEEE_DBL_MAXPOWTWO 183 #define __IEEE_DBL_MAXPOWTWO __IEEE_FLT_MAXPOWTWO 184 185 #undef __IEEE_DBL_NAN_EXP 186 #define __IEEE_DBL_NAN_EXP __IEEE_FLT_NAN_EXP 187 188 #undef __ieee_double_shape_type 189 #define __ieee_double_shape_type __ieee_float_shape_type 190 191 #endif /* _DOUBLE_IS_32BITS */ 192 193 _END_STD_C 194 195 #endif /* _IEEE_FP_H_ */ 196
