Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 1 | /* |
| 2 | * MIPS floating point support |
| 3 | * Copyright (C) 1994-2000 Algorithmics Ltd. |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 4 | * |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 5 | * This program is free software; you can distribute it and/or modify it |
| 6 | * under the terms of the GNU General Public License (Version 2) as |
| 7 | * published by the Free Software Foundation. |
| 8 | * |
| 9 | * This program is distributed in the hope it will be useful, but WITHOUT |
| 10 | * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| 11 | * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| 12 | * for more details. |
| 13 | * |
| 14 | * You should have received a copy of the GNU General Public License along |
| 15 | * with this program; if not, write to the Free Software Foundation, Inc., |
| 16 | * 59 Temple Place - Suite 330, Boston MA 02111-1307, USA. |
| 17 | * |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 18 | * Nov 7, 2000 |
| 19 | * Modification to allow integration with Linux kernel |
| 20 | * |
| 21 | * Kevin D. Kissell, kevink@mips.com and Carsten Langgard, carstenl@mips.com |
| 22 | * Copyright (C) 2000 MIPS Technologies, Inc. All rights reserved. |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 23 | */ |
Ralf Baechle | a663bf9 | 2005-10-23 13:44:31 +0100 | [diff] [blame] | 24 | #ifndef __ARCH_MIPS_MATH_EMU_IEEE754_H |
| 25 | #define __ARCH_MIPS_MATH_EMU_IEEE754_H |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 26 | |
Ralf Baechle | cae5506 | 2014-04-16 00:47:59 +0200 | [diff] [blame] | 27 | #include <linux/compiler.h> |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 28 | #include <asm/byteorder.h> |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 29 | #include <linux/types.h> |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 30 | #include <linux/sched.h> |
Ralf Baechle | f80cc08 | 2014-04-16 00:40:02 +0200 | [diff] [blame] | 31 | #include <asm/bitfield.h> |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 32 | |
| 33 | /* |
| 34 | * Not very pretty, but the Linux kernel's normal va_list definition |
| 35 | * does not allow it to be used as a structure element, as it is here. |
| 36 | */ |
| 37 | #ifndef _STDARG_H |
| 38 | #include <stdarg.h> |
| 39 | #endif |
| 40 | |
Ralf Baechle | 2370881 | 2014-04-16 01:16:02 +0200 | [diff] [blame] | 41 | struct ieee754dp_const { |
Ralf Baechle | f80cc08 | 2014-04-16 00:40:02 +0200 | [diff] [blame] | 42 | __BITFIELD_FIELD(unsigned sign:1, |
| 43 | __BITFIELD_FIELD(unsigned bexp:11, |
| 44 | __BITFIELD_FIELD(unsigned manthi:20, |
| 45 | __BITFIELD_FIELD(unsigned mantlo:32, |
| 46 | ;)))) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 47 | }; |
| 48 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 49 | union ieee754dp { |
Ralf Baechle | 2370881 | 2014-04-16 01:16:02 +0200 | [diff] [blame] | 50 | struct ieee754dp_const oparts; |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 51 | struct { |
Ralf Baechle | f80cc08 | 2014-04-16 00:40:02 +0200 | [diff] [blame] | 52 | __BITFIELD_FIELD(unsigned int sign:1, |
| 53 | __BITFIELD_FIELD(unsigned int bexp:11, |
| 54 | __BITFIELD_FIELD(u64 mant:52, |
| 55 | ;))) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 56 | } parts; |
| 57 | double d; |
| 58 | u64 bits; |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 59 | }; |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 60 | |
Ralf Baechle | 2370881 | 2014-04-16 01:16:02 +0200 | [diff] [blame] | 61 | struct ieee754sp_const { |
Ralf Baechle | f80cc08 | 2014-04-16 00:40:02 +0200 | [diff] [blame] | 62 | __BITFIELD_FIELD(unsigned sign:1, |
| 63 | __BITFIELD_FIELD(unsigned bexp:8, |
| 64 | __BITFIELD_FIELD(unsigned mant:23, |
| 65 | ;))) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 66 | }; |
| 67 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 68 | union ieee754sp { |
Ralf Baechle | 2370881 | 2014-04-16 01:16:02 +0200 | [diff] [blame] | 69 | struct ieee754sp_const parts; |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 70 | float f; |
| 71 | u32 bits; |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 72 | }; |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 73 | |
| 74 | /* |
| 75 | * single precision (often aka float) |
| 76 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 77 | int ieee754sp_finite(union ieee754sp x); |
| 78 | int ieee754sp_class(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 79 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 80 | union ieee754sp ieee754sp_abs(union ieee754sp x); |
| 81 | union ieee754sp ieee754sp_neg(union ieee754sp x); |
| 82 | union ieee754sp ieee754sp_scalb(union ieee754sp x, int); |
| 83 | union ieee754sp ieee754sp_logb(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 84 | |
| 85 | /* x with sign of y */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 86 | union ieee754sp ieee754sp_copysign(union ieee754sp x, union ieee754sp y); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 87 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 88 | union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y); |
| 89 | union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y); |
| 90 | union ieee754sp ieee754sp_mul(union ieee754sp x, union ieee754sp y); |
| 91 | union ieee754sp ieee754sp_div(union ieee754sp x, union ieee754sp y); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 92 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 93 | union ieee754sp ieee754sp_fint(int x); |
| 94 | union ieee754sp ieee754sp_funs(unsigned x); |
| 95 | union ieee754sp ieee754sp_flong(s64 x); |
| 96 | union ieee754sp ieee754sp_fulong(u64 x); |
| 97 | union ieee754sp ieee754sp_fdp(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 98 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 99 | int ieee754sp_tint(union ieee754sp x); |
| 100 | unsigned int ieee754sp_tuns(union ieee754sp x); |
| 101 | s64 ieee754sp_tlong(union ieee754sp x); |
| 102 | u64 ieee754sp_tulong(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 103 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 104 | int ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cop, int sig); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 105 | /* |
| 106 | * basic sp math |
| 107 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 108 | union ieee754sp ieee754sp_modf(union ieee754sp x, union ieee754sp * ip); |
| 109 | union ieee754sp ieee754sp_frexp(union ieee754sp x, int *exp); |
| 110 | union ieee754sp ieee754sp_ldexp(union ieee754sp x, int exp); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 111 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 112 | union ieee754sp ieee754sp_ceil(union ieee754sp x); |
| 113 | union ieee754sp ieee754sp_floor(union ieee754sp x); |
| 114 | union ieee754sp ieee754sp_trunc(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 115 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 116 | union ieee754sp ieee754sp_sqrt(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 117 | |
| 118 | /* |
| 119 | * double precision (often aka double) |
| 120 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 121 | int ieee754dp_finite(union ieee754dp x); |
| 122 | int ieee754dp_class(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 123 | |
| 124 | /* x with sign of y */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 125 | union ieee754dp ieee754dp_copysign(union ieee754dp x, union ieee754dp y); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 126 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 127 | union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y); |
| 128 | union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y); |
| 129 | union ieee754dp ieee754dp_mul(union ieee754dp x, union ieee754dp y); |
| 130 | union ieee754dp ieee754dp_div(union ieee754dp x, union ieee754dp y); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 131 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 132 | union ieee754dp ieee754dp_abs(union ieee754dp x); |
| 133 | union ieee754dp ieee754dp_neg(union ieee754dp x); |
| 134 | union ieee754dp ieee754dp_scalb(union ieee754dp x, int); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 135 | |
| 136 | /* return exponent as integer in floating point format |
| 137 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 138 | union ieee754dp ieee754dp_logb(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 139 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 140 | union ieee754dp ieee754dp_fint(int x); |
| 141 | union ieee754dp ieee754dp_funs(unsigned x); |
| 142 | union ieee754dp ieee754dp_flong(s64 x); |
| 143 | union ieee754dp ieee754dp_fulong(u64 x); |
| 144 | union ieee754dp ieee754dp_fsp(union ieee754sp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 145 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 146 | union ieee754dp ieee754dp_ceil(union ieee754dp x); |
| 147 | union ieee754dp ieee754dp_floor(union ieee754dp x); |
| 148 | union ieee754dp ieee754dp_trunc(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 149 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 150 | int ieee754dp_tint(union ieee754dp x); |
| 151 | unsigned int ieee754dp_tuns(union ieee754dp x); |
| 152 | s64 ieee754dp_tlong(union ieee754dp x); |
| 153 | u64 ieee754dp_tulong(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 154 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 155 | int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cop, int sig); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 156 | /* |
| 157 | * basic sp math |
| 158 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 159 | union ieee754dp ieee754dp_modf(union ieee754dp x, union ieee754dp * ip); |
| 160 | union ieee754dp ieee754dp_frexp(union ieee754dp x, int *exp); |
| 161 | union ieee754dp ieee754dp_ldexp(union ieee754dp x, int exp); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 162 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 163 | union ieee754dp ieee754dp_ceil(union ieee754dp x); |
| 164 | union ieee754dp ieee754dp_floor(union ieee754dp x); |
| 165 | union ieee754dp ieee754dp_trunc(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 166 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 167 | union ieee754dp ieee754dp_sqrt(union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 168 | |
| 169 | |
| 170 | |
| 171 | /* 5 types of floating point number |
| 172 | */ |
| 173 | #define IEEE754_CLASS_NORM 0x00 |
| 174 | #define IEEE754_CLASS_ZERO 0x01 |
| 175 | #define IEEE754_CLASS_DNORM 0x02 |
| 176 | #define IEEE754_CLASS_INF 0x03 |
| 177 | #define IEEE754_CLASS_SNAN 0x04 |
| 178 | #define IEEE754_CLASS_QNAN 0x05 |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 179 | |
| 180 | /* exception numbers */ |
| 181 | #define IEEE754_INEXACT 0x01 |
| 182 | #define IEEE754_UNDERFLOW 0x02 |
| 183 | #define IEEE754_OVERFLOW 0x04 |
| 184 | #define IEEE754_ZERO_DIVIDE 0x08 |
| 185 | #define IEEE754_INVALID_OPERATION 0x10 |
| 186 | |
| 187 | /* cmp operators |
| 188 | */ |
| 189 | #define IEEE754_CLT 0x01 |
| 190 | #define IEEE754_CEQ 0x02 |
| 191 | #define IEEE754_CGT 0x04 |
| 192 | #define IEEE754_CUN 0x08 |
| 193 | |
| 194 | /* rounding mode |
| 195 | */ |
| 196 | #define IEEE754_RN 0 /* round to nearest */ |
| 197 | #define IEEE754_RZ 1 /* round toward zero */ |
| 198 | #define IEEE754_RD 2 /* round toward -Infinity */ |
| 199 | #define IEEE754_RU 3 /* round toward +Infinity */ |
| 200 | |
| 201 | /* other naming */ |
| 202 | #define IEEE754_RM IEEE754_RD |
| 203 | #define IEEE754_RP IEEE754_RU |
| 204 | |
| 205 | /* "normal" comparisons |
| 206 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 207 | static inline int ieee754sp_eq(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 208 | { |
| 209 | return ieee754sp_cmp(x, y, IEEE754_CEQ, 0); |
| 210 | } |
| 211 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 212 | static inline int ieee754sp_ne(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 213 | { |
| 214 | return ieee754sp_cmp(x, y, |
| 215 | IEEE754_CLT | IEEE754_CGT | IEEE754_CUN, 0); |
| 216 | } |
| 217 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 218 | static inline int ieee754sp_lt(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 219 | { |
| 220 | return ieee754sp_cmp(x, y, IEEE754_CLT, 0); |
| 221 | } |
| 222 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 223 | static inline int ieee754sp_le(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 224 | { |
| 225 | return ieee754sp_cmp(x, y, IEEE754_CLT | IEEE754_CEQ, 0); |
| 226 | } |
| 227 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 228 | static inline int ieee754sp_gt(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 229 | { |
| 230 | return ieee754sp_cmp(x, y, IEEE754_CGT, 0); |
| 231 | } |
| 232 | |
| 233 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 234 | static inline int ieee754sp_ge(union ieee754sp x, union ieee754sp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 235 | { |
| 236 | return ieee754sp_cmp(x, y, IEEE754_CGT | IEEE754_CEQ, 0); |
| 237 | } |
| 238 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 239 | static inline int ieee754dp_eq(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 240 | { |
| 241 | return ieee754dp_cmp(x, y, IEEE754_CEQ, 0); |
| 242 | } |
| 243 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 244 | static inline int ieee754dp_ne(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 245 | { |
| 246 | return ieee754dp_cmp(x, y, |
| 247 | IEEE754_CLT | IEEE754_CGT | IEEE754_CUN, 0); |
| 248 | } |
| 249 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 250 | static inline int ieee754dp_lt(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 251 | { |
| 252 | return ieee754dp_cmp(x, y, IEEE754_CLT, 0); |
| 253 | } |
| 254 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 255 | static inline int ieee754dp_le(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 256 | { |
| 257 | return ieee754dp_cmp(x, y, IEEE754_CLT | IEEE754_CEQ, 0); |
| 258 | } |
| 259 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 260 | static inline int ieee754dp_gt(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 261 | { |
| 262 | return ieee754dp_cmp(x, y, IEEE754_CGT, 0); |
| 263 | } |
| 264 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 265 | static inline int ieee754dp_ge(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 266 | { |
| 267 | return ieee754dp_cmp(x, y, IEEE754_CGT | IEEE754_CEQ, 0); |
| 268 | } |
| 269 | |
| 270 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 271 | /* |
| 272 | * Like strtod |
| 273 | */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 274 | union ieee754dp ieee754dp_fstr(const char *s, char **endp); |
| 275 | char *ieee754dp_tstr(union ieee754dp x, int prec, int fmt, int af); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 276 | |
| 277 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 278 | /* |
| 279 | * The control status register |
| 280 | */ |
| 281 | struct _ieee754_csr { |
Ralf Baechle | f80cc08 | 2014-04-16 00:40:02 +0200 | [diff] [blame] | 282 | __BITFIELD_FIELD(unsigned pad0:7, |
| 283 | __BITFIELD_FIELD(unsigned nod:1, /* set 1 for no denormalised numbers */ |
| 284 | __BITFIELD_FIELD(unsigned c:1, /* condition */ |
| 285 | __BITFIELD_FIELD(unsigned pad1:5, |
| 286 | __BITFIELD_FIELD(unsigned cx:6, /* exceptions this operation */ |
| 287 | __BITFIELD_FIELD(unsigned mx:5, /* exception enable mask */ |
| 288 | __BITFIELD_FIELD(unsigned sx:5, /* exceptions total */ |
| 289 | __BITFIELD_FIELD(unsigned rm:2, /* current rounding mode */ |
| 290 | ;)))))))) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 291 | }; |
Atsushi Nemoto | eae8907 | 2006-05-16 01:26:03 +0900 | [diff] [blame] | 292 | #define ieee754_csr (*(struct _ieee754_csr *)(¤t->thread.fpu.fcr31)) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 293 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 294 | static inline unsigned ieee754_getrm(void) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 295 | { |
| 296 | return (ieee754_csr.rm); |
| 297 | } |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 298 | static inline unsigned ieee754_setrm(unsigned rm) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 299 | { |
| 300 | return (ieee754_csr.rm = rm); |
| 301 | } |
| 302 | |
| 303 | /* |
| 304 | * get current exceptions |
| 305 | */ |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 306 | static inline unsigned ieee754_getcx(void) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 307 | { |
| 308 | return (ieee754_csr.cx); |
| 309 | } |
| 310 | |
| 311 | /* test for current exception condition |
| 312 | */ |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 313 | static inline int ieee754_cxtest(unsigned n) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 314 | { |
| 315 | return (ieee754_csr.cx & n); |
| 316 | } |
| 317 | |
| 318 | /* |
| 319 | * get sticky exceptions |
| 320 | */ |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 321 | static inline unsigned ieee754_getsx(void) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 322 | { |
| 323 | return (ieee754_csr.sx); |
| 324 | } |
| 325 | |
| 326 | /* clear sticky conditions |
| 327 | */ |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 328 | static inline unsigned ieee754_clrsx(void) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 329 | { |
| 330 | return (ieee754_csr.sx = 0); |
| 331 | } |
| 332 | |
| 333 | /* test for sticky exception condition |
| 334 | */ |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 335 | static inline int ieee754_sxtest(unsigned n) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 336 | { |
| 337 | return (ieee754_csr.sx & n); |
| 338 | } |
| 339 | |
| 340 | /* debugging */ |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 341 | union ieee754sp ieee754sp_dump(char *s, union ieee754sp x); |
| 342 | union ieee754dp ieee754dp_dump(char *s, union ieee754dp x); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 343 | |
| 344 | #define IEEE754_SPCVAL_PZERO 0 |
| 345 | #define IEEE754_SPCVAL_NZERO 1 |
| 346 | #define IEEE754_SPCVAL_PONE 2 |
| 347 | #define IEEE754_SPCVAL_NONE 3 |
| 348 | #define IEEE754_SPCVAL_PTEN 4 |
| 349 | #define IEEE754_SPCVAL_NTEN 5 |
| 350 | #define IEEE754_SPCVAL_PINFINITY 6 |
| 351 | #define IEEE754_SPCVAL_NINFINITY 7 |
| 352 | #define IEEE754_SPCVAL_INDEF 8 |
| 353 | #define IEEE754_SPCVAL_PMAX 9 /* +max norm */ |
| 354 | #define IEEE754_SPCVAL_NMAX 10 /* -max norm */ |
| 355 | #define IEEE754_SPCVAL_PMIN 11 /* +min norm */ |
| 356 | #define IEEE754_SPCVAL_NMIN 12 /* +min norm */ |
| 357 | #define IEEE754_SPCVAL_PMIND 13 /* +min denorm */ |
| 358 | #define IEEE754_SPCVAL_NMIND 14 /* +min denorm */ |
| 359 | #define IEEE754_SPCVAL_P1E31 15 /* + 1.0e31 */ |
| 360 | #define IEEE754_SPCVAL_P1E63 16 /* + 1.0e63 */ |
| 361 | |
Ralf Baechle | 2370881 | 2014-04-16 01:16:02 +0200 | [diff] [blame] | 362 | extern const struct ieee754dp_const __ieee754dp_spcvals[]; |
| 363 | extern const struct ieee754sp_const __ieee754sp_spcvals[]; |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 364 | #define ieee754dp_spcvals ((const union ieee754dp *)__ieee754dp_spcvals) |
| 365 | #define ieee754sp_spcvals ((const union ieee754sp *)__ieee754sp_spcvals) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 366 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 367 | /* |
| 368 | * Return infinity with given sign |
| 369 | */ |
| 370 | #define ieee754dp_inf(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)]) |
| 371 | #define ieee754dp_zero(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PZERO+(sn)]) |
| 372 | #define ieee754dp_one(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PONE+(sn)]) |
| 373 | #define ieee754dp_ten(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PTEN+(sn)]) |
| 374 | #define ieee754dp_indef() (ieee754dp_spcvals[IEEE754_SPCVAL_INDEF]) |
| 375 | #define ieee754dp_max(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMAX+(sn)]) |
| 376 | #define ieee754dp_min(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIN+(sn)]) |
| 377 | #define ieee754dp_mind(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIND+(sn)]) |
| 378 | #define ieee754dp_1e31() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E31]) |
| 379 | #define ieee754dp_1e63() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E63]) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 380 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 381 | #define ieee754sp_inf(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)]) |
| 382 | #define ieee754sp_zero(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PZERO+(sn)]) |
| 383 | #define ieee754sp_one(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PONE+(sn)]) |
| 384 | #define ieee754sp_ten(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PTEN+(sn)]) |
| 385 | #define ieee754sp_indef() (ieee754sp_spcvals[IEEE754_SPCVAL_INDEF]) |
| 386 | #define ieee754sp_max(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMAX+(sn)]) |
| 387 | #define ieee754sp_min(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIN+(sn)]) |
| 388 | #define ieee754sp_mind(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIND+(sn)]) |
| 389 | #define ieee754sp_1e31() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E31]) |
| 390 | #define ieee754sp_1e63() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E63]) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 391 | |
Ralf Baechle | cd21dfc | 2005-04-28 13:39:10 +0000 | [diff] [blame] | 392 | /* |
| 393 | * Indefinite integer value |
| 394 | */ |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 395 | #define ieee754si_indef() INT_MAX |
| 396 | #ifdef LONG_LONG_MAX |
| 397 | #define ieee754di_indef() LONG_LONG_MAX |
| 398 | #else |
| 399 | #define ieee754di_indef() ((s64)(~0ULL>>1)) |
| 400 | #endif |
| 401 | |
| 402 | /* IEEE exception context, passed to handler */ |
| 403 | struct ieee754xctx { |
| 404 | const char *op; /* operation name */ |
| 405 | int rt; /* result type */ |
| 406 | union { |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 407 | union ieee754sp sp; /* single precision */ |
| 408 | union ieee754dp dp; /* double precision */ |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 409 | #ifdef IEEE854_XP |
| 410 | ieee754xp xp; /* extended precision */ |
| 411 | #endif |
| 412 | int si; /* standard signed integer (32bits) */ |
| 413 | s64 di; /* extended signed integer (64bits) */ |
| 414 | } rv; /* default result format implied by op */ |
| 415 | va_list ap; |
| 416 | }; |
| 417 | |
| 418 | /* result types for xctx.rt */ |
| 419 | #define IEEE754_RT_SP 0 |
| 420 | #define IEEE754_RT_DP 1 |
| 421 | #define IEEE754_RT_XP 2 |
| 422 | #define IEEE754_RT_SI 3 |
| 423 | #define IEEE754_RT_DI 4 |
| 424 | |
Ralf Baechle | cae5506 | 2014-04-16 00:47:59 +0200 | [diff] [blame] | 425 | extern void __cold ieee754_xcpt(struct ieee754xctx *xcp); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 426 | |
| 427 | /* compat */ |
| 428 | #define ieee754dp_fix(x) ieee754dp_tint(x) |
| 429 | #define ieee754sp_fix(x) ieee754sp_tint(x) |
Ralf Baechle | a663bf9 | 2005-10-23 13:44:31 +0100 | [diff] [blame] | 430 | |
| 431 | #endif /* __ARCH_MIPS_MATH_EMU_IEEE754_H */ |