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 | union ieee754dp ieee754dp_neg(union ieee754dp x) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 30 | { |
| 31 | COMPXDP; |
| 32 | |
| 33 | EXPLODEXDP; |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame] | 34 | ieee754_clearcx(); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 35 | FLUSHXDP; |
| 36 | |
Atsushi Nemoto | c5033d7 | 2006-01-02 21:59:49 +0900 | [diff] [blame] | 37 | /* |
| 38 | * Invert the sign ALWAYS to prevent an endless recursion on |
| 39 | * pow() in libc. |
| 40 | */ |
| 41 | /* quick fix up */ |
| 42 | DPSIGN(x) ^= 1; |
| 43 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 44 | if (xc == IEEE754_CLASS_SNAN) { |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 45 | union ieee754dp y = ieee754dp_indef(); |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame] | 46 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
Atsushi Nemoto | c5033d7 | 2006-01-02 21:59:49 +0900 | [diff] [blame] | 47 | DPSIGN(y) = DPSIGN(x); |
Ralf Baechle | 90efba3 | 2014-04-25 03:19:57 +0200 | [diff] [blame^] | 48 | return ieee754dp_nanxcpt(y); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 49 | } |
| 50 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 51 | return x; |
| 52 | } |
| 53 | |
| 54 | |
Ralf Baechle | 2209bcb | 2014-04-16 01:31:11 +0200 | [diff] [blame] | 55 | union ieee754dp ieee754dp_abs(union ieee754dp x) |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 56 | { |
| 57 | COMPXDP; |
| 58 | |
| 59 | EXPLODEXDP; |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame] | 60 | ieee754_clearcx(); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 61 | FLUSHXDP; |
| 62 | |
Nigel Stephens | cea2be4 | 2009-10-12 14:57:18 -0700 | [diff] [blame] | 63 | /* Clear sign ALWAYS, irrespective of NaN */ |
| 64 | DPSIGN(x) = 0; |
| 65 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 66 | if (xc == IEEE754_CLASS_SNAN) { |
Ralf Baechle | 9e8bad1 | 2014-04-19 00:36:32 +0200 | [diff] [blame] | 67 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
Ralf Baechle | 90efba3 | 2014-04-25 03:19:57 +0200 | [diff] [blame^] | 68 | return ieee754dp_nanxcpt(ieee754dp_indef()); |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 69 | } |
| 70 | |
Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 71 | return x; |
| 72 | } |