| Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 1 | /* | 
|  | 2 | * Some debug functions | 
|  | 3 | * | 
|  | 4 | * MIPS floating point support | 
|  | 5 | * | 
|  | 6 | * Copyright (C) 1994-2000 Algorithmics Ltd. | 
| Linus Torvalds | 1da177e | 2005-04-16 15:20:36 -0700 | [diff] [blame] | 7 | * | 
|  | 8 | *  This program is free software; you can distribute it and/or modify it | 
|  | 9 | *  under the terms of the GNU General Public License (Version 2) as | 
|  | 10 | *  published by the Free Software Foundation. | 
|  | 11 | * | 
|  | 12 | *  This program is distributed in the hope it will be useful, but WITHOUT | 
|  | 13 | *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | 
|  | 14 | *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License | 
|  | 15 | *  for more details. | 
|  | 16 | * | 
|  | 17 | *  You should have received a copy of the GNU General Public License along | 
|  | 18 | *  with this program; if not, write to the Free Software Foundation, Inc., | 
|  | 19 | *  59 Temple Place - Suite 330, Boston MA 02111-1307, USA. | 
|  | 20 | * | 
|  | 21 | *  Nov 7, 2000 | 
|  | 22 | *  Modified to build and operate in Linux kernel environment. | 
|  | 23 | * | 
|  | 24 | *  Kevin D. Kissell, kevink@mips.com and Carsten Langgaard, carstenl@mips.com | 
|  | 25 | *  Copyright (C) 2000 MIPS Technologies, Inc. All rights reserved. | 
|  | 26 | */ | 
|  | 27 |  | 
|  | 28 | #include <linux/kernel.h> | 
|  | 29 | #include "ieee754.h" | 
|  | 30 |  | 
|  | 31 | #define DP_EBIAS	1023 | 
|  | 32 | #define DP_EMIN		(-1022) | 
|  | 33 | #define DP_EMAX		1023 | 
|  | 34 | #define DP_FBITS	52 | 
|  | 35 |  | 
|  | 36 | #define SP_EBIAS	127 | 
|  | 37 | #define SP_EMIN		(-126) | 
|  | 38 | #define SP_EMAX		127 | 
|  | 39 | #define SP_FBITS	23 | 
|  | 40 |  | 
|  | 41 | #define DP_MBIT(x)	((u64)1 << (x)) | 
|  | 42 | #define DP_HIDDEN_BIT	DP_MBIT(DP_FBITS) | 
|  | 43 | #define DP_SIGN_BIT	DP_MBIT(63) | 
|  | 44 |  | 
|  | 45 |  | 
|  | 46 | #define SP_MBIT(x)	((u32)1 << (x)) | 
|  | 47 | #define SP_HIDDEN_BIT	SP_MBIT(SP_FBITS) | 
|  | 48 | #define SP_SIGN_BIT	SP_MBIT(31) | 
|  | 49 |  | 
|  | 50 |  | 
|  | 51 | #define SPSIGN(sp)	(sp.parts.sign) | 
|  | 52 | #define SPBEXP(sp)	(sp.parts.bexp) | 
|  | 53 | #define SPMANT(sp)	(sp.parts.mant) | 
|  | 54 |  | 
|  | 55 | #define DPSIGN(dp)	(dp.parts.sign) | 
|  | 56 | #define DPBEXP(dp)	(dp.parts.bexp) | 
|  | 57 | #define DPMANT(dp)	(dp.parts.mant) | 
|  | 58 |  | 
|  | 59 | ieee754dp ieee754dp_dump(char *m, ieee754dp x) | 
|  | 60 | { | 
|  | 61 | int i; | 
|  | 62 |  | 
|  | 63 | printk("%s", m); | 
|  | 64 | printk("<%08x,%08x>\n", (unsigned) (x.bits >> 32), | 
|  | 65 | (unsigned) x.bits); | 
|  | 66 | printk("\t="); | 
|  | 67 | switch (ieee754dp_class(x)) { | 
|  | 68 | case IEEE754_CLASS_QNAN: | 
|  | 69 | case IEEE754_CLASS_SNAN: | 
|  | 70 | printk("Nan %c", DPSIGN(x) ? '-' : '+'); | 
|  | 71 | for (i = DP_FBITS - 1; i >= 0; i--) | 
|  | 72 | printk("%c", DPMANT(x) & DP_MBIT(i) ? '1' : '0'); | 
|  | 73 | break; | 
|  | 74 | case IEEE754_CLASS_INF: | 
|  | 75 | printk("%cInfinity", DPSIGN(x) ? '-' : '+'); | 
|  | 76 | break; | 
|  | 77 | case IEEE754_CLASS_ZERO: | 
|  | 78 | printk("%cZero", DPSIGN(x) ? '-' : '+'); | 
|  | 79 | break; | 
|  | 80 | case IEEE754_CLASS_DNORM: | 
|  | 81 | printk("%c0.", DPSIGN(x) ? '-' : '+'); | 
|  | 82 | for (i = DP_FBITS - 1; i >= 0; i--) | 
|  | 83 | printk("%c", DPMANT(x) & DP_MBIT(i) ? '1' : '0'); | 
|  | 84 | printk("e%d", DPBEXP(x) - DP_EBIAS); | 
|  | 85 | break; | 
|  | 86 | case IEEE754_CLASS_NORM: | 
|  | 87 | printk("%c1.", DPSIGN(x) ? '-' : '+'); | 
|  | 88 | for (i = DP_FBITS - 1; i >= 0; i--) | 
|  | 89 | printk("%c", DPMANT(x) & DP_MBIT(i) ? '1' : '0'); | 
|  | 90 | printk("e%d", DPBEXP(x) - DP_EBIAS); | 
|  | 91 | break; | 
|  | 92 | default: | 
|  | 93 | printk("Illegal/Unknown IEEE754 value class"); | 
|  | 94 | } | 
|  | 95 | printk("\n"); | 
|  | 96 | return x; | 
|  | 97 | } | 
|  | 98 |  | 
|  | 99 | ieee754sp ieee754sp_dump(char *m, ieee754sp x) | 
|  | 100 | { | 
|  | 101 | int i; | 
|  | 102 |  | 
|  | 103 | printk("%s=", m); | 
|  | 104 | printk("<%08x>\n", (unsigned) x.bits); | 
|  | 105 | printk("\t="); | 
|  | 106 | switch (ieee754sp_class(x)) { | 
|  | 107 | case IEEE754_CLASS_QNAN: | 
|  | 108 | case IEEE754_CLASS_SNAN: | 
|  | 109 | printk("Nan %c", SPSIGN(x) ? '-' : '+'); | 
|  | 110 | for (i = SP_FBITS - 1; i >= 0; i--) | 
|  | 111 | printk("%c", SPMANT(x) & SP_MBIT(i) ? '1' : '0'); | 
|  | 112 | break; | 
|  | 113 | case IEEE754_CLASS_INF: | 
|  | 114 | printk("%cInfinity", SPSIGN(x) ? '-' : '+'); | 
|  | 115 | break; | 
|  | 116 | case IEEE754_CLASS_ZERO: | 
|  | 117 | printk("%cZero", SPSIGN(x) ? '-' : '+'); | 
|  | 118 | break; | 
|  | 119 | case IEEE754_CLASS_DNORM: | 
|  | 120 | printk("%c0.", SPSIGN(x) ? '-' : '+'); | 
|  | 121 | for (i = SP_FBITS - 1; i >= 0; i--) | 
|  | 122 | printk("%c", SPMANT(x) & SP_MBIT(i) ? '1' : '0'); | 
|  | 123 | printk("e%d", SPBEXP(x) - SP_EBIAS); | 
|  | 124 | break; | 
|  | 125 | case IEEE754_CLASS_NORM: | 
|  | 126 | printk("%c1.", SPSIGN(x) ? '-' : '+'); | 
|  | 127 | for (i = SP_FBITS - 1; i >= 0; i--) | 
|  | 128 | printk("%c", SPMANT(x) & SP_MBIT(i) ? '1' : '0'); | 
|  | 129 | printk("e%d", SPBEXP(x) - SP_EBIAS); | 
|  | 130 | break; | 
|  | 131 | default: | 
|  | 132 | printk("Illegal/Unknown IEEE754 value class"); | 
|  | 133 | } | 
|  | 134 | printk("\n"); | 
|  | 135 | return x; | 
|  | 136 | } |