| David Brazdil | 0f672f6 | 2019-12-10 10:32:29 +0000 | [diff] [blame^] | 1 | // SPDX-License-Identifier: GPL-2.0-only |
| Andrew Scull | b4b6d4a | 2019-01-02 15:54:55 +0000 | [diff] [blame] | 2 | /* ieee754 floating point arithmetic |
| 3 | * single and double precision |
| 4 | * |
| 5 | * BUGS |
| 6 | * not much dp done |
| 7 | * doesn't generate IEEE754_INEXACT |
| Andrew Scull | b4b6d4a | 2019-01-02 15:54:55 +0000 | [diff] [blame] | 8 | */ |
| 9 | /* |
| 10 | * MIPS floating point support |
| 11 | * Copyright (C) 1994-2000 Algorithmics Ltd. |
| Andrew Scull | b4b6d4a | 2019-01-02 15:54:55 +0000 | [diff] [blame] | 12 | */ |
| 13 | |
| 14 | #include <linux/compiler.h> |
| 15 | |
| 16 | #include "ieee754.h" |
| 17 | #include "ieee754sp.h" |
| 18 | #include "ieee754dp.h" |
| 19 | |
| 20 | /* |
| 21 | * Special constants |
| 22 | */ |
| 23 | |
| 24 | /* |
| 25 | * Older GCC requires the inner braces for initialization of union ieee754dp's |
| 26 | * anonymous struct member. Without an error will result. |
| 27 | */ |
| 28 | #define xPCNST(s, b, m, ebias) \ |
| 29 | { \ |
| 30 | { \ |
| 31 | .sign = (s), \ |
| 32 | .bexp = (b) + ebias, \ |
| 33 | .mant = (m) \ |
| 34 | } \ |
| 35 | } |
| 36 | |
| 37 | #define DPCNST(s, b, m) \ |
| 38 | xPCNST(s, b, m, DP_EBIAS) |
| 39 | |
| 40 | const union ieee754dp __ieee754dp_spcvals[] = { |
| 41 | DPCNST(0, DP_EMIN - 1, 0x0000000000000ULL), /* + zero */ |
| 42 | DPCNST(1, DP_EMIN - 1, 0x0000000000000ULL), /* - zero */ |
| 43 | DPCNST(0, 0, 0x0000000000000ULL), /* + 1.0 */ |
| 44 | DPCNST(1, 0, 0x0000000000000ULL), /* - 1.0 */ |
| 45 | DPCNST(0, 3, 0x4000000000000ULL), /* + 10.0 */ |
| 46 | DPCNST(1, 3, 0x4000000000000ULL), /* - 10.0 */ |
| 47 | DPCNST(0, DP_EMAX + 1, 0x0000000000000ULL), /* + infinity */ |
| 48 | DPCNST(1, DP_EMAX + 1, 0x0000000000000ULL), /* - infinity */ |
| 49 | DPCNST(0, DP_EMAX + 1, 0x7FFFFFFFFFFFFULL), /* + ind legacy qNaN */ |
| 50 | DPCNST(0, DP_EMAX + 1, 0x8000000000000ULL), /* + indef 2008 qNaN */ |
| 51 | DPCNST(0, DP_EMAX, 0xFFFFFFFFFFFFFULL), /* + max */ |
| 52 | DPCNST(1, DP_EMAX, 0xFFFFFFFFFFFFFULL), /* - max */ |
| 53 | DPCNST(0, DP_EMIN, 0x0000000000000ULL), /* + min normal */ |
| 54 | DPCNST(1, DP_EMIN, 0x0000000000000ULL), /* - min normal */ |
| 55 | DPCNST(0, DP_EMIN - 1, 0x0000000000001ULL), /* + min denormal */ |
| 56 | DPCNST(1, DP_EMIN - 1, 0x0000000000001ULL), /* - min denormal */ |
| 57 | DPCNST(0, 31, 0x0000000000000ULL), /* + 1.0e31 */ |
| 58 | DPCNST(0, 63, 0x0000000000000ULL), /* + 1.0e63 */ |
| 59 | }; |
| 60 | |
| 61 | #define SPCNST(s, b, m) \ |
| 62 | xPCNST(s, b, m, SP_EBIAS) |
| 63 | |
| 64 | const union ieee754sp __ieee754sp_spcvals[] = { |
| 65 | SPCNST(0, SP_EMIN - 1, 0x000000), /* + zero */ |
| 66 | SPCNST(1, SP_EMIN - 1, 0x000000), /* - zero */ |
| 67 | SPCNST(0, 0, 0x000000), /* + 1.0 */ |
| 68 | SPCNST(1, 0, 0x000000), /* - 1.0 */ |
| 69 | SPCNST(0, 3, 0x200000), /* + 10.0 */ |
| 70 | SPCNST(1, 3, 0x200000), /* - 10.0 */ |
| 71 | SPCNST(0, SP_EMAX + 1, 0x000000), /* + infinity */ |
| 72 | SPCNST(1, SP_EMAX + 1, 0x000000), /* - infinity */ |
| 73 | SPCNST(0, SP_EMAX + 1, 0x3FFFFF), /* + indef legacy quiet NaN */ |
| 74 | SPCNST(0, SP_EMAX + 1, 0x400000), /* + indef 2008 quiet NaN */ |
| 75 | SPCNST(0, SP_EMAX, 0x7FFFFF), /* + max normal */ |
| 76 | SPCNST(1, SP_EMAX, 0x7FFFFF), /* - max normal */ |
| 77 | SPCNST(0, SP_EMIN, 0x000000), /* + min normal */ |
| 78 | SPCNST(1, SP_EMIN, 0x000000), /* - min normal */ |
| 79 | SPCNST(0, SP_EMIN - 1, 0x000001), /* + min denormal */ |
| 80 | SPCNST(1, SP_EMIN - 1, 0x000001), /* - min denormal */ |
| 81 | SPCNST(0, 31, 0x000000), /* + 1.0e31 */ |
| 82 | SPCNST(0, 63, 0x000000), /* + 1.0e63 */ |
| 83 | }; |