arch/mips/math-emu/ieee754d.c

   1 /*
   2  * Some debug functions
   3  *
   4  * MIPS floating point support
   5  *
   6  * Copyright (C) 1994-2000 Algorithmics Ltd.
   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 }