ieee754int.h
上传用户:lgb322
上传日期:2013-02-24
资源大小:30529k
文件大小:3k
- /*
- * IEEE754 floating point
- * common internal header file
- */
- /*
- * MIPS floating point support
- * Copyright (C) 1994-2000 Algorithmics Ltd. All rights reserved.
- * http://www.algor.co.uk
- *
- * ########################################################################
- *
- * 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 "ieee754.h"
- #define DP_EBIAS 1023
- #define DP_EMIN (-1022)
- #define DP_EMAX 1023
- #define DP_MBITS 52
- #define SP_EBIAS 127
- #define SP_EMIN (-126)
- #define SP_EMAX 127
- #define SP_MBITS 23
- #define DP_MBIT(x) ((unsigned long long)1 << (x))
- #define DP_HIDDEN_BIT DP_MBIT(DP_MBITS)
- #define DP_SIGN_BIT DP_MBIT(63)
- #define SP_MBIT(x) ((unsigned long)1 << (x))
- #define SP_HIDDEN_BIT SP_MBIT(SP_MBITS)
- #define SP_SIGN_BIT SP_MBIT(31)
- #define SPSIGN(sp) (sp.parts.sign)
- #define SPBEXP(sp) (sp.parts.bexp)
- #define SPMANT(sp) (sp.parts.mant)
- #define DPSIGN(dp) (dp.parts.sign)
- #define DPBEXP(dp) (dp.parts.bexp)
- #define DPMANT(dp) (dp.parts.mant)
- #define CLPAIR(x,y) ((x)*6+(y))
- #define CLEARCX
- (ieee754_csr.cx = 0)
- #define SETCX(x)
- (ieee754_csr.cx |= (x),ieee754_csr.sx |= (x),ieee754_csr.mx & (x))
- #define TSTX()
- (ieee754_csr.cx & ieee754_csr.mx)
- #define COMPXSP
- unsigned xm; int xe; int xs; int xc
- #define COMPYSP
- unsigned ym; int ye; int ys; int yc
- #define EXPLODESP(v,vc,vs,ve,vm)
- {
- vs = SPSIGN(v);
- ve = SPBEXP(v);
- vm = SPMANT(v);
- if(ve == SP_EMAX+1+SP_EBIAS){
- if(vm == 0)
- vc = IEEE754_CLASS_INF;
- else if(vm & SP_MBIT(SP_MBITS-1))
- vc = IEEE754_CLASS_QNAN;
- else
- vc = IEEE754_CLASS_SNAN;
- } else if(ve == SP_EMIN-1+SP_EBIAS) {
- if(vm) {
- ve = SP_EMIN;
- vc = IEEE754_CLASS_DNORM;
- } else
- vc = IEEE754_CLASS_ZERO;
- } else {
- ve -= SP_EBIAS;
- vm |= SP_HIDDEN_BIT;
- vc = IEEE754_CLASS_NORM;
- }
- }
- #define EXPLODEXSP EXPLODESP(x,xc,xs,xe,xm)
- #define EXPLODEYSP EXPLODESP(y,yc,ys,ye,ym)
- #define COMPXDP
- unsigned long long xm; int xe; int xs; int xc
- #define COMPYDP
- unsigned long long ym; int ye; int ys; int yc
- #define EXPLODEDP(v,vc,vs,ve,vm)
- {
- vm = DPMANT(v);
- vs = DPSIGN(v);
- ve = DPBEXP(v);
- if(ve == DP_EMAX+1+DP_EBIAS){
- if(vm == 0)
- vc = IEEE754_CLASS_INF;
- else if(vm & DP_MBIT(DP_MBITS-1))
- vc = IEEE754_CLASS_QNAN;
- else
- vc = IEEE754_CLASS_SNAN;
- } else if(ve == DP_EMIN-1+DP_EBIAS) {
- if(vm) {
- ve = DP_EMIN;
- vc = IEEE754_CLASS_DNORM;
- } else
- vc = IEEE754_CLASS_ZERO;
- } else {
- ve -= DP_EBIAS;
- vm |= DP_HIDDEN_BIT;
- vc = IEEE754_CLASS_NORM;
- }
- }
- #define EXPLODEXDP EXPLODEDP(x,xc,xs,xe,xm)
- #define EXPLODEYDP EXPLODEDP(y,yc,ys,ye,ym)