ieee754d.c
上传用户:lgb322
上传日期:2013-02-24
资源大小:30529k
文件大小:4k
源码类别:

嵌入式Linux

开发平台:

Unix_Linux

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