开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > 嵌入式Linux > 查看源码
sp_div.c
资源名称：SBC-2410X_kernel.tar.gz　[点击查看]
上传用户：lgb322
上传日期：2013-02-24
资源大小：30529k
文件大小：4k
源码类别：

嵌入式Linux

开发平台：

Unix_Linux

sp_div.c：源码内容
  1. /* IEEE754 floating point arithmetic
  2.  * single precision
  3.  */
  4. /*
  5.  * MIPS floating point support
  6.  * Copyright (C) 1994-2000 Algorithmics Ltd.  All rights reserved.
  7.  * http://www.algor.co.uk
  8.  *
  9.  * ########################################################################
  10.  *
  11.  *  This program is free software; you can distribute it and/or modify it
  12.  *  under the terms of the GNU General Public License (Version 2) as
  13.  *  published by the Free Software Foundation.
  14.  *
  15.  *  This program is distributed in the hope it will be useful, but WITHOUT
  16.  *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  17.  *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  18.  *  for more details.
  19.  *
  20.  *  You should have received a copy of the GNU General Public License along
  21.  *  with this program; if not, write to the Free Software Foundation, Inc.,
  22.  *  59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
  23.  *
  24.  * ########################################################################
  25.  */
  28. #include "ieee754sp.h"
  30. ieee754sp ieee754sp_div(ieee754sp x, ieee754sp y)
  31. {
  32. COMPXSP;
  33. COMPYSP;
  35. CLEARCX;
  37. EXPLODEXSP;
  38. EXPLODEYSP;
  40. switch (CLPAIR(xc, yc)) {
  41. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_SNAN):
  42. return ieee754sp_nanxcpt(ieee754sp_bestnan(x, y), "div", x,
  43.  y);
  45. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_SNAN):
  46. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_SNAN):
  47. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_SNAN):
  48. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_SNAN):
  49. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_SNAN):
  50. return ieee754sp_nanxcpt(y, "div", x, y);
  52. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_QNAN):
  53. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_ZERO):
  54. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_NORM):
  55. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_DNORM):
  56. case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_INF):
  57. return ieee754sp_nanxcpt(x, "div", x, y);
  59. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_QNAN):
  60. return ieee754sp_bestnan(x, y);
  62. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_QNAN):
  63. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_QNAN):
  64. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_QNAN):
  65. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_QNAN):
  66. return y;
  68. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_ZERO):
  69. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_NORM):
  70. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_DNORM):
  71. case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_INF):
  72. return x;
  75. /* Infinity handeling 
  76.  */
  78. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_INF):
  79. SETCX(IEEE754_INVALID_OPERATION);
  80. return ieee754sp_xcpt(ieee754sp_indef(), "div", x, y);
  82. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_INF):
  83. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_INF):
  84. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_INF):
  85. return ieee754sp_zero(xs ^ ys);
  87. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_ZERO):
  88. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_NORM):
  89. case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_DNORM):
  90. return ieee754sp_inf(xs ^ ys);
  92. /* Zero handeling 
  93.  */
  95. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_ZERO):
  96. SETCX(IEEE754_INVALID_OPERATION);
  97. return ieee754sp_xcpt(ieee754sp_indef(), "div", x, y);
  99. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
  100. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
  101. SETCX(IEEE754_ZERO_DIVIDE);
  102. return ieee754sp_xcpt(ieee754sp_inf(xs ^ ys), "div", x, y);
  104. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_NORM):
  105. case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_DNORM):
  106. return ieee754sp_zero(xs == ys ? 0 : 1);
  108. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_DNORM):
  109. SPDNORMX;
  111. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_DNORM):
  112. SPDNORMY;
  113. break;
  115. case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_NORM):
  116. SPDNORMX;
  117. break;
  119. case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_NORM):
  120. break;
  121. }
  122. assert(xm & SP_HIDDEN_BIT);
  123. assert(ym & SP_HIDDEN_BIT);
  125. /* provide rounding space */
  126. xm <<= 3;
  127. ym <<= 3;
  129. {
  130. /* now the dirty work */
  132. unsigned rm = 0;
  133. int re = xe - ye;
  134. unsigned bm;
  136. for (bm = SP_MBIT(SP_MBITS + 2); bm; bm >>= 1) {
  137. if (xm >= ym) {
  138. xm -= ym;
  139. rm |= bm;
  140. if (xm == 0)
  141. break;
  142. }
  143. xm <<= 1;
  144. }
  145. rm <<= 1;
  146. if (xm)
  147. rm |= 1; /* have remainder, set sticky */
  149. assert(rm);
  151. /* normalise rm to rounding precision ?
  152.  */
  153. while ((rm >> (SP_MBITS + 3)) == 0) {
  154. rm <<= 1;
  155. re--;
  156. }
  158. SPNORMRET2(xs == ys ? 0 : 1, re, rm, "div", x, y);
  159. }
  160. }