Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 1 | /* IEEE754 floating point arithmetic |
| 2 | * double precision: common utilities |
| 3 | */ |
| 4 | /* |
| 5 | * MIPS floating point support |
| 6 | * Copyright (C) 1994-2000 Algorithmics Ltd. |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 7 | * |
| 8 | * ######################################################################## |
| 9 | * |
| 10 | * This program is free software; you can distribute it and/or modify it |
| 11 | * under the terms of the GNU General Public License (Version 2) as |
| 12 | * published by the Free Software Foundation. |
| 13 | * |
| 14 | * This program is distributed in the hope it will be useful, but WITHOUT |
| 15 | * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| 16 | * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| 17 | * for more details. |
| 18 | * |
| 19 | * You should have received a copy of the GNU General Public License along |
| 20 | * with this program; if not, write to the Free Software Foundation, Inc., |
| 21 | * 59 Temple Place - Suite 330, Boston MA 02111-1307, USA. |
| 22 | * |
| 23 | * ######################################################################## |
| 24 | */ |
| 25 | |
| 26 | |
| 27 | #include "ieee754dp.h" |
| 28 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 29 | int ieee754dp_finite(union ieee754dp x) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 30 | { |
| 31 | return DPBEXP(x) != DP_EMAX + 1 + DP_EBIAS; |
| 32 | } |
| 33 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 34 | union ieee754dp ieee754dp_copysign(union ieee754dp x, union ieee754dp y) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 35 | { |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame^] | 36 | ieee754_clearcx(); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 37 | DPSIGN(x) = DPSIGN(y); |
| 38 | return x; |
| 39 | } |
| 40 | |
| 41 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 42 | union ieee754dp ieee754dp_neg(union ieee754dp x) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 43 | { |
| 44 | COMPXDP; |
| 45 | |
| 46 | EXPLODEXDP; |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame^] | 47 | ieee754_clearcx(); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 48 | FLUSHXDP; |
| 49 | |
Atsushi Nemoto | c5033d7 | 2006-01-02 21:59:49 +0900 | [diff] [blame] | 50 | /* |
| 51 | * Invert the sign ALWAYS to prevent an endless recursion on |
| 52 | * pow() in libc. |
| 53 | */ |
| 54 | /* quick fix up */ |
| 55 | DPSIGN(x) ^= 1; |
| 56 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 57 | if (xc == IEEE754_CLASS_SNAN) { |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 58 | union ieee754dp y = ieee754dp_indef(); |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame^] | 59 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
Atsushi Nemoto | c5033d7 | 2006-01-02 21:59:49 +0900 | [diff] [blame] | 60 | DPSIGN(y) = DPSIGN(x); |
| 61 | return ieee754dp_nanxcpt(y, "neg"); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 62 | } |
| 63 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 64 | return x; |
| 65 | } |
| 66 | |
| 67 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 68 | union ieee754dp ieee754dp_abs(union ieee754dp x) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 69 | { |
| 70 | COMPXDP; |
| 71 | |
| 72 | EXPLODEXDP; |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame^] | 73 | ieee754_clearcx(); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 74 | FLUSHXDP; |
| 75 | |
Nigel Stephens | cea2be4 | 2009-10-12 14:57:18 -0700 | [diff] [blame] | 76 | /* Clear sign ALWAYS, irrespective of NaN */ |
| 77 | DPSIGN(x) = 0; |
| 78 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 79 | if (xc == IEEE754_CLASS_SNAN) { |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame^] | 80 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
Nigel Stephens | cea2be4 | 2009-10-12 14:57:18 -0700 | [diff] [blame] | 81 | return ieee754dp_nanxcpt(ieee754dp_indef(), "abs"); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 82 | } |
| 83 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 84 | return x; |
| 85 | } |