| /* ieee754 floating point arithmetic |
| * single and double precision |
| * |
| * BUGS |
| * not much dp done |
| * doesn't generate IEEE754_INEXACT |
| * |
| */ |
| /* |
| * MIPS floating point support |
| * Copyright (C) 1994-2000 Algorithmics Ltd. |
| * |
| * ######################################################################## |
| * |
| * This program is free software; you can distribute it and/or modify it |
| * under the terms of the GNU General Public License (Version 2) as |
| * published by the Free Software Foundation. |
| * |
| * This program is distributed in the hope it will be useful, but WITHOUT |
| * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| * for more details. |
| * |
| * You should have received a copy of the GNU General Public License along |
| * with this program; if not, write to the Free Software Foundation, Inc., |
| * 59 Temple Place - Suite 330, Boston MA 02111-1307, USA. |
| * |
| * ######################################################################## |
| */ |
| |
| #include <linux/compiler.h> |
| |
| #include "ieee754.h" |
| #include "ieee754sp.h" |
| #include "ieee754dp.h" |
| |
| /* special constants |
| */ |
| |
| #define DPCNST(s, b, m) \ |
| { \ |
| .sign = (s), \ |
| .bexp = (b) + DP_EBIAS, \ |
| .mant = (m) \ |
| } |
| |
| const union ieee754dp __ieee754dp_spcvals[] = { |
| DPCNST(0, DP_EMIN - 1, 0x0000000000000ULL), /* + zero */ |
| DPCNST(1, DP_EMIN - 1, 0x0000000000000ULL), /* - zero */ |
| DPCNST(0, 0, 0x0000000000000ULL), /* + 1.0 */ |
| DPCNST(1, 0, 0x0000000000000ULL), /* - 1.0 */ |
| DPCNST(0, 3, 0x4000000000000ULL), /* + 10.0 */ |
| DPCNST(1, 3, 0x4000000000000ULL), /* - 10.0 */ |
| DPCNST(0, DP_EMAX + 1, 0x0000000000000ULL), /* + infinity */ |
| DPCNST(1, DP_EMAX + 1, 0x0000000000000ULL), /* - infinity */ |
| DPCNST(0, DP_EMAX + 1, 0x7FFFFFFFFFFFFULL), /* + indef quiet Nan */ |
| DPCNST(0, DP_EMAX, 0xFFFFFFFFFFFFFULL), /* + max */ |
| DPCNST(1, DP_EMAX, 0xFFFFFFFFFFFFFULL), /* - max */ |
| DPCNST(0, DP_EMIN, 0x0000000000000ULL), /* + min normal */ |
| DPCNST(1, DP_EMIN, 0x0000000000000ULL), /* - min normal */ |
| DPCNST(0, DP_EMIN - 1, 0x0000000000001ULL), /* + min denormal */ |
| DPCNST(1, DP_EMIN - 1, 0x0000000000001ULL), /* - min denormal */ |
| DPCNST(0, 31, 0x0000000000000ULL), /* + 1.0e31 */ |
| DPCNST(0, 63, 0x0000000000000ULL), /* + 1.0e63 */ |
| }; |
| |
| #define SPCNST(s, b, m) \ |
| { \ |
| .sign = (s), \ |
| .bexp = (b) + SP_EBIAS, \ |
| .mant = (m) \ |
| } |
| |
| const union ieee754sp __ieee754sp_spcvals[] = { |
| SPCNST(0, SP_EMIN - 1, 0x000000), /* + zero */ |
| SPCNST(1, SP_EMIN - 1, 0x000000), /* - zero */ |
| SPCNST(0, 0, 0x000000), /* + 1.0 */ |
| SPCNST(1, 0, 0x000000), /* - 1.0 */ |
| SPCNST(0, 3, 0x200000), /* + 10.0 */ |
| SPCNST(1, 3, 0x200000), /* - 10.0 */ |
| SPCNST(0, SP_EMAX + 1, 0x000000), /* + infinity */ |
| SPCNST(1, SP_EMAX + 1, 0x000000), /* - infinity */ |
| SPCNST(0, SP_EMAX + 1, 0x3FFFFF), /* + indef quiet Nan */ |
| SPCNST(0, SP_EMAX, 0x7FFFFF), /* + max normal */ |
| SPCNST(1, SP_EMAX, 0x7FFFFF), /* - max normal */ |
| SPCNST(0, SP_EMIN, 0x000000), /* + min normal */ |
| SPCNST(1, SP_EMIN, 0x000000), /* - min normal */ |
| SPCNST(0, SP_EMIN - 1, 0x000001), /* + min denormal */ |
| SPCNST(1, SP_EMIN - 1, 0x000001), /* - min denormal */ |
| SPCNST(0, 31, 0x000000), /* + 1.0e31 */ |
| SPCNST(0, 63, 0x000000), /* + 1.0e63 */ |
| }; |