MIPS: math-emu: Switch to using the MIPS rounding modes.
Previously math-emu was using the IEEE-754 constants internally. These were differing by having the constants for rounding to +/- infinity switched, so a conversion was necessary. This would be entirely avoidable if the MIPS constants were used throughout, so get rid of the bloat. Signed-off-by: Ralf Baechle <ralf@linux-mips.org>
This commit is contained in:
		
					parent
					
						
							
								aef3fb76aa
							
						
					
				
			
			
				commit
				
					
						56a6473339
					
				
			
		
					 15 changed files with 63 additions and 91 deletions
				
			
		| 
						 | 
					@ -67,21 +67,6 @@ static int fpux_emu(struct pt_regs *,
 | 
				
			||||||
/* Determine rounding mode from the RM bits of the FCSR */
 | 
					/* Determine rounding mode from the RM bits of the FCSR */
 | 
				
			||||||
#define modeindex(v) ((v) & FPU_CSR_RM)
 | 
					#define modeindex(v) ((v) & FPU_CSR_RM)
 | 
				
			||||||
 | 
					
 | 
				
			||||||
/* Convert MIPS rounding mode (0..3) to IEEE library modes. */
 | 
					 | 
				
			||||||
static const unsigned char ieee_rm[4] = {
 | 
					 | 
				
			||||||
	[FPU_CSR_RN] = IEEE754_RN,
 | 
					 | 
				
			||||||
	[FPU_CSR_RZ] = IEEE754_RZ,
 | 
					 | 
				
			||||||
	[FPU_CSR_RU] = IEEE754_RU,
 | 
					 | 
				
			||||||
	[FPU_CSR_RD] = IEEE754_RD,
 | 
					 | 
				
			||||||
};
 | 
					 | 
				
			||||||
/* Convert IEEE library modes to MIPS rounding mode (0..3). */
 | 
					 | 
				
			||||||
static const unsigned char mips_rm[4] = {
 | 
					 | 
				
			||||||
	[IEEE754_RN] = FPU_CSR_RN,
 | 
					 | 
				
			||||||
	[IEEE754_RZ] = FPU_CSR_RZ,
 | 
					 | 
				
			||||||
	[IEEE754_RD] = FPU_CSR_RD,
 | 
					 | 
				
			||||||
	[IEEE754_RU] = FPU_CSR_RU,
 | 
					 | 
				
			||||||
};
 | 
					 | 
				
			||||||
 | 
					 | 
				
			||||||
/* convert condition code register number to csr bit */
 | 
					/* convert condition code register number to csr bit */
 | 
				
			||||||
static const unsigned int fpucondbit[8] = {
 | 
					static const unsigned int fpucondbit[8] = {
 | 
				
			||||||
	FPU_CSR_COND0,
 | 
						FPU_CSR_COND0,
 | 
				
			||||||
| 
						 | 
					@ -907,8 +892,7 @@ emul:
 | 
				
			||||||
			/* cop control register rd -> gpr[rt] */
 | 
								/* cop control register rd -> gpr[rt] */
 | 
				
			||||||
			if (MIPSInst_RD(ir) == FPCREG_CSR) {
 | 
								if (MIPSInst_RD(ir) == FPCREG_CSR) {
 | 
				
			||||||
				value = ctx->fcr31;
 | 
									value = ctx->fcr31;
 | 
				
			||||||
				value = (value & ~FPU_CSR_RM) |
 | 
									value = (value & ~FPU_CSR_RM) | modeindex(value);
 | 
				
			||||||
					mips_rm[modeindex(value)];
 | 
					 | 
				
			||||||
				pr_debug("%p gpr[%d]<-csr=%08x\n",
 | 
									pr_debug("%p gpr[%d]<-csr=%08x\n",
 | 
				
			||||||
					 (void *) (xcp->cp0_epc),
 | 
										 (void *) (xcp->cp0_epc),
 | 
				
			||||||
					 MIPSInst_RT(ir), value);
 | 
										 MIPSInst_RT(ir), value);
 | 
				
			||||||
| 
						 | 
					@ -939,9 +923,8 @@ emul:
 | 
				
			||||||
				 * Don't write reserved bits,
 | 
									 * Don't write reserved bits,
 | 
				
			||||||
				 * and convert to ieee library modes
 | 
									 * and convert to ieee library modes
 | 
				
			||||||
				 */
 | 
									 */
 | 
				
			||||||
				ctx->fcr31 = (value &
 | 
									ctx->fcr31 = (value & ~(FPU_CSR_RSVD | FPU_CSR_RM)) |
 | 
				
			||||||
						~(FPU_CSR_RSVD | FPU_CSR_RM)) |
 | 
										     modeindex(value);
 | 
				
			||||||
						ieee_rm[modeindex(value)];
 | 
					 | 
				
			||||||
			}
 | 
								}
 | 
				
			||||||
			if ((ctx->fcr31 >> 5) & ctx->fcr31 & FPU_CSR_ALL_E) {
 | 
								if ((ctx->fcr31 >> 5) & ctx->fcr31 & FPU_CSR_ALL_E) {
 | 
				
			||||||
				return SIGFPE;
 | 
									return SIGFPE;
 | 
				
			||||||
| 
						 | 
					@ -1515,7 +1498,7 @@ copcsr:
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			oldrm = ieee754_csr.rm;
 | 
								oldrm = ieee754_csr.rm;
 | 
				
			||||||
			SPFROMREG(fs, MIPSInst_FS(ir));
 | 
								SPFROMREG(fs, MIPSInst_FS(ir));
 | 
				
			||||||
			ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
 | 
								ieee754_csr.rm = modeindex(MIPSInst_FUNC(ir));
 | 
				
			||||||
			rv.w = ieee754sp_tint(fs);
 | 
								rv.w = ieee754sp_tint(fs);
 | 
				
			||||||
			ieee754_csr.rm = oldrm;
 | 
								ieee754_csr.rm = oldrm;
 | 
				
			||||||
			rfmt = w_fmt;
 | 
								rfmt = w_fmt;
 | 
				
			||||||
| 
						 | 
					@ -1539,7 +1522,7 @@ copcsr:
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			oldrm = ieee754_csr.rm;
 | 
								oldrm = ieee754_csr.rm;
 | 
				
			||||||
			SPFROMREG(fs, MIPSInst_FS(ir));
 | 
								SPFROMREG(fs, MIPSInst_FS(ir));
 | 
				
			||||||
			ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
 | 
								ieee754_csr.rm = modeindex(MIPSInst_FUNC(ir));
 | 
				
			||||||
			rv.l = ieee754sp_tlong(fs);
 | 
								rv.l = ieee754sp_tlong(fs);
 | 
				
			||||||
			ieee754_csr.rm = oldrm;
 | 
								ieee754_csr.rm = oldrm;
 | 
				
			||||||
			rfmt = l_fmt;
 | 
								rfmt = l_fmt;
 | 
				
			||||||
| 
						 | 
					@ -1692,7 +1675,7 @@ dcopuop:
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			oldrm = ieee754_csr.rm;
 | 
								oldrm = ieee754_csr.rm;
 | 
				
			||||||
			DPFROMREG(fs, MIPSInst_FS(ir));
 | 
								DPFROMREG(fs, MIPSInst_FS(ir));
 | 
				
			||||||
			ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
 | 
								ieee754_csr.rm = modeindex(MIPSInst_FUNC(ir));
 | 
				
			||||||
			rv.w = ieee754dp_tint(fs);
 | 
								rv.w = ieee754dp_tint(fs);
 | 
				
			||||||
			ieee754_csr.rm = oldrm;
 | 
								ieee754_csr.rm = oldrm;
 | 
				
			||||||
			rfmt = w_fmt;
 | 
								rfmt = w_fmt;
 | 
				
			||||||
| 
						 | 
					@ -1716,7 +1699,7 @@ dcopuop:
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			oldrm = ieee754_csr.rm;
 | 
								oldrm = ieee754_csr.rm;
 | 
				
			||||||
			DPFROMREG(fs, MIPSInst_FS(ir));
 | 
								DPFROMREG(fs, MIPSInst_FS(ir));
 | 
				
			||||||
			ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
 | 
								ieee754_csr.rm = modeindex(MIPSInst_FUNC(ir));
 | 
				
			||||||
			rv.l = ieee754dp_tlong(fs);
 | 
								rv.l = ieee754dp_tlong(fs);
 | 
				
			||||||
			ieee754_csr.rm = oldrm;
 | 
								ieee754_csr.rm = oldrm;
 | 
				
			||||||
			rfmt = l_fmt;
 | 
								rfmt = l_fmt;
 | 
				
			||||||
| 
						 | 
					@ -1926,11 +1909,7 @@ int fpu_emulator_cop1Handler(struct pt_regs *xcp, struct mips_fpu_struct *ctx,
 | 
				
			||||||
			 * ieee754_csr.	 But ieee754_csr.rm is ieee
 | 
								 * ieee754_csr.	 But ieee754_csr.rm is ieee
 | 
				
			||||||
			 * library modes. (not mips rounding mode)
 | 
								 * library modes. (not mips rounding mode)
 | 
				
			||||||
			 */
 | 
								 */
 | 
				
			||||||
			/* convert to ieee library modes */
 | 
					 | 
				
			||||||
			ieee754_csr.rm = ieee_rm[ieee754_csr.rm];
 | 
					 | 
				
			||||||
			sig = cop1Emulate(xcp, ctx, dec_insn, fault_addr);
 | 
								sig = cop1Emulate(xcp, ctx, dec_insn, fault_addr);
 | 
				
			||||||
			/* revert to mips rounding mode */
 | 
					 | 
				
			||||||
			ieee754_csr.rm = mips_rm[ieee754_csr.rm];
 | 
					 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
 | 
					
 | 
				
			||||||
		if (has_fpu)
 | 
							if (has_fpu)
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -91,7 +91,7 @@ union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y)
 | 
				
			||||||
		if (xs == ys)
 | 
							if (xs == ys)
 | 
				
			||||||
			return x;
 | 
								return x;
 | 
				
			||||||
		else
 | 
							else
 | 
				
			||||||
			return ieee754dp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754dp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
| 
						 | 
					@ -168,7 +168,7 @@ union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y)
 | 
				
			||||||
			xs = ys;
 | 
								xs = ys;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		if (xm == 0)
 | 
							if (xm == 0)
 | 
				
			||||||
			return ieee754dp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754dp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
		/*
 | 
							/*
 | 
				
			||||||
		 * Normalize to rounding precision.
 | 
							 * Normalize to rounding precision.
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -80,7 +80,7 @@ union ieee754dp ieee754dp_sqrt(union ieee754dp x)
 | 
				
			||||||
	oldcsr = ieee754_csr;
 | 
						oldcsr = ieee754_csr;
 | 
				
			||||||
	ieee754_csr.mx &= ~IEEE754_INEXACT;
 | 
						ieee754_csr.mx &= ~IEEE754_INEXACT;
 | 
				
			||||||
	ieee754_csr.sx &= ~IEEE754_INEXACT;
 | 
						ieee754_csr.sx &= ~IEEE754_INEXACT;
 | 
				
			||||||
	ieee754_csr.rm = IEEE754_RN;
 | 
						ieee754_csr.rm = FPU_CSR_RN;
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	/* adjust exponent to prevent overflow */
 | 
						/* adjust exponent to prevent overflow */
 | 
				
			||||||
	scalx = 0;
 | 
						scalx = 0;
 | 
				
			||||||
| 
						 | 
					@ -122,7 +122,7 @@ union ieee754dp ieee754dp_sqrt(union ieee754dp x)
 | 
				
			||||||
	/* twiddle last bit to force y correctly rounded */
 | 
						/* twiddle last bit to force y correctly rounded */
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	/* set RZ, clear INEX flag */
 | 
						/* set RZ, clear INEX flag */
 | 
				
			||||||
	ieee754_csr.rm = IEEE754_RZ;
 | 
						ieee754_csr.rm = FPU_CSR_RZ;
 | 
				
			||||||
	ieee754_csr.sx &= ~IEEE754_INEXACT;
 | 
						ieee754_csr.sx &= ~IEEE754_INEXACT;
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	/* t=x/y; ...chopped quotient, possibly inexact */
 | 
						/* t=x/y; ...chopped quotient, possibly inexact */
 | 
				
			||||||
| 
						 | 
					@ -139,10 +139,10 @@ union ieee754dp ieee754dp_sqrt(union ieee754dp x)
 | 
				
			||||||
		oldcsr.sx |= IEEE754_INEXACT;
 | 
							oldcsr.sx |= IEEE754_INEXACT;
 | 
				
			||||||
 | 
					
 | 
				
			||||||
		switch (oldcsr.rm) {
 | 
							switch (oldcsr.rm) {
 | 
				
			||||||
		case IEEE754_RU:
 | 
							case FPU_CSR_RU:
 | 
				
			||||||
			y.bits += 1;
 | 
								y.bits += 1;
 | 
				
			||||||
			/* drop through */
 | 
								/* drop through */
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			t.bits += 1;
 | 
								t.bits += 1;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -91,7 +91,7 @@ union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y)
 | 
				
			||||||
		if (xs != ys)
 | 
							if (xs != ys)
 | 
				
			||||||
			return x;
 | 
								return x;
 | 
				
			||||||
		else
 | 
							else
 | 
				
			||||||
			return ieee754dp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754dp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
| 
						 | 
					@ -171,7 +171,7 @@ union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y)
 | 
				
			||||||
			xs = ys;
 | 
								xs = ys;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		if (xm == 0) {
 | 
							if (xm == 0) {
 | 
				
			||||||
			if (ieee754_csr.rm == IEEE754_RD)
 | 
								if (ieee754_csr.rm == FPU_CSR_RD)
 | 
				
			||||||
				return ieee754dp_zero(1);	/* round negative inf. => sign = -1 */
 | 
									return ieee754dp_zero(1);	/* round negative inf. => sign = -1 */
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
				return ieee754dp_zero(0);	/* other round modes   => sign = 1 */
 | 
									return ieee754dp_zero(0);	/* other round modes   => sign = 1 */
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -74,17 +74,17 @@ int ieee754dp_tint(union ieee754dp x)
 | 
				
			||||||
		   to be zero */
 | 
							   to be zero */
 | 
				
			||||||
		odd = (xm & 0x1) != 0x0;
 | 
							odd = (xm & 0x1) != 0x0;
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			if (round && (sticky || odd))
 | 
								if (round && (sticky || odd))
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if ((round || sticky) && !xs)
 | 
								if ((round || sticky) && !xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if ((round || sticky) && xs)
 | 
								if ((round || sticky) && xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -79,17 +79,17 @@ s64 ieee754dp_tlong(union ieee754dp x)
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		odd = (xm & 0x1) != 0x0;
 | 
							odd = (xm & 0x1) != 0x0;
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			if (round && (sticky || odd))
 | 
								if (round && (sticky || odd))
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if ((round || sticky) && !xs)
 | 
								if ((round || sticky) && !xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if ((round || sticky) && xs)
 | 
								if ((round || sticky) && xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -126,13 +126,6 @@ enum {
 | 
				
			||||||
#define IEEE754_CGT	0x04
 | 
					#define IEEE754_CGT	0x04
 | 
				
			||||||
#define IEEE754_CUN	0x08
 | 
					#define IEEE754_CUN	0x08
 | 
				
			||||||
 | 
					
 | 
				
			||||||
/* rounding mode
 | 
					 | 
				
			||||||
*/
 | 
					 | 
				
			||||||
#define IEEE754_RN	0	/* round to nearest */
 | 
					 | 
				
			||||||
#define IEEE754_RZ	1	/* round toward zero  */
 | 
					 | 
				
			||||||
#define IEEE754_RD	2	/* round toward -Infinity */
 | 
					 | 
				
			||||||
#define IEEE754_RU	3	/* round toward +Infinity */
 | 
					 | 
				
			||||||
 | 
					 | 
				
			||||||
/* "normal" comparisons
 | 
					/* "normal" comparisons
 | 
				
			||||||
*/
 | 
					*/
 | 
				
			||||||
static inline int ieee754sp_eq(union ieee754sp x, union ieee754sp y)
 | 
					static inline int ieee754sp_eq(union ieee754sp x, union ieee754sp y)
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -67,17 +67,17 @@ static u64 ieee754dp_get_rounding(int sn, u64 xm)
 | 
				
			||||||
	 */
 | 
						 */
 | 
				
			||||||
	if (xm & (DP_MBIT(3) - 1)) {
 | 
						if (xm & (DP_MBIT(3) - 1)) {
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			xm += 0x3 + ((xm >> 3) & 1);
 | 
								xm += 0x3 + ((xm >> 3) & 1);
 | 
				
			||||||
			/* xm += (xm&0x8)?0x4:0x3 */
 | 
								/* xm += (xm&0x8)?0x4:0x3 */
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if (!sn)	/* ?? */
 | 
								if (!sn)	/* ?? */
 | 
				
			||||||
				xm += 0x8;
 | 
									xm += 0x8;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if (sn) /* ?? */
 | 
								if (sn) /* ?? */
 | 
				
			||||||
				xm += 0x8;
 | 
									xm += 0x8;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					@ -108,15 +108,15 @@ union ieee754dp ieee754dp_format(int sn, int xe, u64 xm)
 | 
				
			||||||
			ieee754_setcx(IEEE754_INEXACT);
 | 
								ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			switch(ieee754_csr.rm) {
 | 
								switch(ieee754_csr.rm) {
 | 
				
			||||||
			case IEEE754_RN:
 | 
								case FPU_CSR_RN:
 | 
				
			||||||
			case IEEE754_RZ:
 | 
								case FPU_CSR_RZ:
 | 
				
			||||||
				return ieee754dp_zero(sn);
 | 
									return ieee754dp_zero(sn);
 | 
				
			||||||
			case IEEE754_RU:    /* toward +Infinity */
 | 
								case FPU_CSR_RU:    /* toward +Infinity */
 | 
				
			||||||
				if (sn == 0)
 | 
									if (sn == 0)
 | 
				
			||||||
					return ieee754dp_min(0);
 | 
										return ieee754dp_min(0);
 | 
				
			||||||
				else
 | 
									else
 | 
				
			||||||
					return ieee754dp_zero(1);
 | 
										return ieee754dp_zero(1);
 | 
				
			||||||
			case IEEE754_RD:    /* toward -Infinity */
 | 
								case FPU_CSR_RD:    /* toward -Infinity */
 | 
				
			||||||
				if (sn == 0)
 | 
									if (sn == 0)
 | 
				
			||||||
					return ieee754dp_zero(0);
 | 
										return ieee754dp_zero(0);
 | 
				
			||||||
				else
 | 
									else
 | 
				
			||||||
| 
						 | 
					@ -172,16 +172,16 @@ union ieee754dp ieee754dp_format(int sn, int xe, u64 xm)
 | 
				
			||||||
		ieee754_setcx(IEEE754_INEXACT);
 | 
							ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
		/* -O can be table indexed by (rm,sn) */
 | 
							/* -O can be table indexed by (rm,sn) */
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			return ieee754dp_inf(sn);
 | 
								return ieee754dp_inf(sn);
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			return ieee754dp_max(sn);
 | 
								return ieee754dp_max(sn);
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if (sn == 0)
 | 
								if (sn == 0)
 | 
				
			||||||
				return ieee754dp_inf(0);
 | 
									return ieee754dp_inf(0);
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
				return ieee754dp_max(1);
 | 
									return ieee754dp_max(1);
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if (sn == 0)
 | 
								if (sn == 0)
 | 
				
			||||||
				return ieee754dp_max(0);
 | 
									return ieee754dp_max(0);
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -67,17 +67,17 @@ static unsigned ieee754sp_get_rounding(int sn, unsigned xm)
 | 
				
			||||||
	 */
 | 
						 */
 | 
				
			||||||
	if (xm & (SP_MBIT(3) - 1)) {
 | 
						if (xm & (SP_MBIT(3) - 1)) {
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			xm += 0x3 + ((xm >> 3) & 1);
 | 
								xm += 0x3 + ((xm >> 3) & 1);
 | 
				
			||||||
			/* xm += (xm&0x8)?0x4:0x3 */
 | 
								/* xm += (xm&0x8)?0x4:0x3 */
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if (!sn)	/* ?? */
 | 
								if (!sn)	/* ?? */
 | 
				
			||||||
				xm += 0x8;
 | 
									xm += 0x8;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if (sn) /* ?? */
 | 
								if (sn) /* ?? */
 | 
				
			||||||
				xm += 0x8;
 | 
									xm += 0x8;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					@ -108,15 +108,15 @@ union ieee754sp ieee754sp_format(int sn, int xe, unsigned xm)
 | 
				
			||||||
			ieee754_setcx(IEEE754_INEXACT);
 | 
								ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
			switch(ieee754_csr.rm) {
 | 
								switch(ieee754_csr.rm) {
 | 
				
			||||||
			case IEEE754_RN:
 | 
								case FPU_CSR_RN:
 | 
				
			||||||
			case IEEE754_RZ:
 | 
								case FPU_CSR_RZ:
 | 
				
			||||||
				return ieee754sp_zero(sn);
 | 
									return ieee754sp_zero(sn);
 | 
				
			||||||
			case IEEE754_RU:      /* toward +Infinity */
 | 
								case FPU_CSR_RU:      /* toward +Infinity */
 | 
				
			||||||
				if (sn == 0)
 | 
									if (sn == 0)
 | 
				
			||||||
					return ieee754sp_min(0);
 | 
										return ieee754sp_min(0);
 | 
				
			||||||
				else
 | 
									else
 | 
				
			||||||
					return ieee754sp_zero(1);
 | 
										return ieee754sp_zero(1);
 | 
				
			||||||
			case IEEE754_RD:      /* toward -Infinity */
 | 
								case FPU_CSR_RD:      /* toward -Infinity */
 | 
				
			||||||
				if (sn == 0)
 | 
									if (sn == 0)
 | 
				
			||||||
					return ieee754sp_zero(0);
 | 
										return ieee754sp_zero(0);
 | 
				
			||||||
				else
 | 
									else
 | 
				
			||||||
| 
						 | 
					@ -170,16 +170,16 @@ union ieee754sp ieee754sp_format(int sn, int xe, unsigned xm)
 | 
				
			||||||
		ieee754_setcx(IEEE754_INEXACT);
 | 
							ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
		/* -O can be table indexed by (rm,sn) */
 | 
							/* -O can be table indexed by (rm,sn) */
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			return ieee754sp_inf(sn);
 | 
								return ieee754sp_inf(sn);
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			return ieee754sp_max(sn);
 | 
								return ieee754sp_max(sn);
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if (sn == 0)
 | 
								if (sn == 0)
 | 
				
			||||||
				return ieee754sp_inf(0);
 | 
									return ieee754sp_inf(0);
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
				return ieee754sp_max(1);
 | 
									return ieee754sp_max(1);
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if (sn == 0)
 | 
								if (sn == 0)
 | 
				
			||||||
				return ieee754sp_max(0);
 | 
									return ieee754sp_max(0);
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -91,7 +91,7 @@ union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y)
 | 
				
			||||||
		if (xs == ys)
 | 
							if (xs == ys)
 | 
				
			||||||
			return x;
 | 
								return x;
 | 
				
			||||||
		else
 | 
							else
 | 
				
			||||||
			return ieee754sp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754sp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
| 
						 | 
					@ -165,7 +165,7 @@ union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y)
 | 
				
			||||||
			xs = ys;
 | 
								xs = ys;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		if (xm == 0)
 | 
							if (xm == 0)
 | 
				
			||||||
			return ieee754sp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754sp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
		/*
 | 
							/*
 | 
				
			||||||
		 * Normalize in extended single precision
 | 
							 * Normalize in extended single precision
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -57,8 +57,8 @@ union ieee754sp ieee754sp_fdp(union ieee754dp x)
 | 
				
			||||||
		/* can't possibly be sp representable */
 | 
							/* can't possibly be sp representable */
 | 
				
			||||||
		ieee754_setcx(IEEE754_UNDERFLOW);
 | 
							ieee754_setcx(IEEE754_UNDERFLOW);
 | 
				
			||||||
		ieee754_setcx(IEEE754_INEXACT);
 | 
							ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
		if ((ieee754_csr.rm == IEEE754_RU && !xs) ||
 | 
							if ((ieee754_csr.rm == FPU_CSR_RU && !xs) ||
 | 
				
			||||||
				(ieee754_csr.rm == IEEE754_RD && xs))
 | 
									(ieee754_csr.rm == FPU_CSR_RD && xs))
 | 
				
			||||||
			return ieee754sp_mind(xs);
 | 
								return ieee754sp_mind(xs);
 | 
				
			||||||
		return ieee754sp_zero(xs);
 | 
							return ieee754sp_zero(xs);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -100,10 +100,10 @@ union ieee754sp ieee754sp_sqrt(union ieee754sp x)
 | 
				
			||||||
	if (ix != 0) {
 | 
						if (ix != 0) {
 | 
				
			||||||
		ieee754_setcx(IEEE754_INEXACT);
 | 
							ieee754_setcx(IEEE754_INEXACT);
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RU:
 | 
							case FPU_CSR_RU:
 | 
				
			||||||
			q += 2;
 | 
								q += 2;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			q += (q & 1);
 | 
								q += (q & 1);
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -91,7 +91,7 @@ union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y)
 | 
				
			||||||
		if (xs != ys)
 | 
							if (xs != ys)
 | 
				
			||||||
			return x;
 | 
								return x;
 | 
				
			||||||
		else
 | 
							else
 | 
				
			||||||
			return ieee754sp_zero(ieee754_csr.rm == IEEE754_RD);
 | 
								return ieee754sp_zero(ieee754_csr.rm == FPU_CSR_RD);
 | 
				
			||||||
 | 
					
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
	case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
						case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 | 
				
			||||||
| 
						 | 
					@ -165,7 +165,7 @@ union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y)
 | 
				
			||||||
			xs = ys;
 | 
								xs = ys;
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		if (xm == 0) {
 | 
							if (xm == 0) {
 | 
				
			||||||
			if (ieee754_csr.rm == IEEE754_RD)
 | 
								if (ieee754_csr.rm == FPU_CSR_RD)
 | 
				
			||||||
				return ieee754sp_zero(1);	/* round negative inf. => sign = -1 */
 | 
									return ieee754sp_zero(1);	/* round negative inf. => sign = -1 */
 | 
				
			||||||
			else
 | 
								else
 | 
				
			||||||
				return ieee754sp_zero(0);	/* other round modes   => sign = 1 */
 | 
									return ieee754sp_zero(0);	/* other round modes   => sign = 1 */
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -79,17 +79,17 @@ int ieee754sp_tint(union ieee754sp x)
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		odd = (xm & 0x1) != 0x0;
 | 
							odd = (xm & 0x1) != 0x0;
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			if (round && (sticky || odd))
 | 
								if (round && (sticky || odd))
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if ((round || sticky) && !xs)
 | 
								if ((round || sticky) && !xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if ((round || sticky) && xs)
 | 
								if ((round || sticky) && xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
| 
						 | 
					@ -76,17 +76,17 @@ s64 ieee754sp_tlong(union ieee754sp x)
 | 
				
			||||||
		}
 | 
							}
 | 
				
			||||||
		odd = (xm & 0x1) != 0x0;
 | 
							odd = (xm & 0x1) != 0x0;
 | 
				
			||||||
		switch (ieee754_csr.rm) {
 | 
							switch (ieee754_csr.rm) {
 | 
				
			||||||
		case IEEE754_RN:
 | 
							case FPU_CSR_RN:
 | 
				
			||||||
			if (round && (sticky || odd))
 | 
								if (round && (sticky || odd))
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RZ:
 | 
							case FPU_CSR_RZ:
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RU:	/* toward +Infinity */
 | 
							case FPU_CSR_RU:	/* toward +Infinity */
 | 
				
			||||||
			if ((round || sticky) && !xs)
 | 
								if ((round || sticky) && !xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
		case IEEE754_RD:	/* toward -Infinity */
 | 
							case FPU_CSR_RD:	/* toward -Infinity */
 | 
				
			||||||
			if ((round || sticky) && xs)
 | 
								if ((round || sticky) && xs)
 | 
				
			||||||
				xm++;
 | 
									xm++;
 | 
				
			||||||
			break;
 | 
								break;
 | 
				
			||||||
| 
						 | 
					
 | 
				
			||||||
		Loading…
	
	Add table
		Add a link
		
	
		Reference in a new issue