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