parisc/math-emu/dbl_float.h

13 #define Dbl_firstword(value) Dallp1(value)
24 #define Dbl_allp1(object) Dallp1(object)
30     Dallp1(src1dst) = (Dallp1(src2)|~((unsigned int)1<<31)) & Dallp1(src1dst)
32     Dallp1(src1dst) = (Dallp1(src2)&((unsigned int)1<<31)) | Dallp1(src1dst)
38 #define Dbl_clear_sign(srcdst) Dallp1(srcdst) &= ~((unsigned int)1<<31)
40     Dallp1(srcdst) &= Dmantissap1((unsigned int)-1)
53         Dallp2(srcdstB) = Dallp1(srcdstA) >> (varamount-32);		\
54         Dallp1(srcdstA)=0;						\
57 	Variable_shift_double(Dallp1(srcdstA), Dallp2(srcdstB), 	\
59 	Dallp1(srcdstA) >>= varamount;					\
65 	Dallp1(srcdstA) &= ((unsigned int)1<<31);  /* clear expmant field */ \
76 	Dallp1(srcdstA) = Dallp2(srcdstB) << (varamount-32);		\
81 	    Dallp1(srcdstA) = (Dallp1(srcdstA) << (varamount)) |	\
87     Shiftdouble(Dallp1(lefta), Dallp2(leftb), 31, Dallp1(resulta));	\
95     Shiftdouble(Dallp1(srcdstA),Dallp2(srcdstB),1,Dallp2(srcdstB));\
96     Dallp1(srcdstA) = (int)Dallp1(srcdstA) >> 1
105     if( (Dallp2(dbl_valueB) += 1) == 0 )  Dallp1(dbl_valueA) += 1
110     if( Dallp2(dbl_valueB) == 0 )  Dallp1(dbl_valueA) -= 1; \
120     (Dallp1(dbl_valueA) || Dallp2(dbl_valueB))
131 #define Dbl_iszero(dbl_valueA,dbl_valueB) (Dallp1(dbl_valueA)==0 && \
133 #define Dbl_iszero_allp1(dbl_value) (Dallp1(dbl_value)==0)
162     (Dallp1(dbl_op1a) < Dallp1(dbl_op2a) ||			\
163      (Dallp1(dbl_op1a) == Dallp1(dbl_op2a) &&			\
166     (Dallp1(dbl_op1a) > Dallp1(dbl_op2a) ||			\
167      (Dallp1(dbl_op1a) == Dallp1(dbl_op2a) &&			\
170     (Dallp1(dbl_op1a) > Dallp1(dbl_op2a) ||			\
171      (Dallp1(dbl_op1a) == Dallp1(dbl_op2a) &&			\
174     (Dallp1(dbl_op1a) < Dallp1(dbl_op2a) ||			\
175      (Dallp1(dbl_op1a) == Dallp1(dbl_op2a) &&			\
178      ((Dallp1(dbl_op1a) == Dallp1(dbl_op2a)) &&			\
182     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),24,Dallp1(dbl_valueA)); \
185     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),25,Dallp1(dbl_valueA)); \
188     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),28,Dallp1(dbl_valueA)); \
191     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),29,Dallp1(dbl_valueA)); \
194     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),30,Dallp1(dbl_valueA)); \
197     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),31,Dallp1(dbl_valueA)); \
201     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),8,Dallp2(dbl_valueB)); \
202     Dallp1(dbl_valueA) >>= 8
204     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),4,Dallp2(dbl_valueB)); \
205     Dallp1(dbl_valueA) >>= 4
207     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),2,Dallp2(dbl_valueB)); \
208     Dallp1(dbl_valueA) >>= 2
210     Shiftdouble(Dallp1(dbl_valueA),Dallp2(dbl_valueB),1,Dallp2(dbl_valueB)); \
211     Dallp1(dbl_valueA) >>= 1
250     Dallp1(desta) = src->wd0;		\
253     dest->wd0 = Dallp1(srca);		\
264     Dallp1(dbl_valueA) 						\
268     Dallp1(dbl_valueA) = ((unsigned int)1<<31) |		\
272     Dallp1(dbl_valueA) = ((unsigned int)sign << 31) | 		\
281 #define Dbl_setzero_sign(dbl_value) Dallp1(dbl_value) &= 0x7fffffff
283     Dallp1(dbl_value) &= 0x800fffff
285     Dallp1(dbl_valueA) &= 0xfff00000; 			\
287 #define Dbl_setzero_mantissap1(dbl_value) Dallp1(dbl_value) &= 0xfff00000
290     Dallp1(dbl_valueA) &= 0x80000000;		\
293     Dallp1(dbl_valueA) &= 0x80000000
295     Dallp1(dbl_valueA) = 0; Dallp2(dbl_valueB) = 0
296 #define Dbl_setzerop1(dbl_value) Dallp1(dbl_value) = 0
299     Dallp1(dbl_value) = (unsigned int)1 << 31; Dallp2(dbl_value) = 0
300 #define Dbl_setnegativezerop1(dbl_value) Dallp1(dbl_value) = (unsigned int)1<<31
309     Dallp1(dbl_valueA) = ((DBL_EMAX+DBL_BIAS) << (32-(1+DBL_EXP_LENGTH))) \
313     Dallp1(dbl_valueA) = ((DBL_EMAX+DBL_BIAS) << (32-(1+DBL_EXP_LENGTH))) \
324     Dallp1(dbl_valueA) = ((1<<DBL_EXP_LENGTH) | DBL_INFINITY_EXPONENT) 	\
328     Dallp1(dbl_valueA) = ((unsigned int)sign << 31) |			\
344 	        Variable_shift_double(Dallp1(srcdstA),Dallp2(srcdstB),	\
349 	    Dallp2(srcdstB) = Dallp1(srcdstA) >> (shift - 32);		\
353 	    Extall(extent) = Dallp1(srcdstA);				\
357 	Dallp1(srcdstA) = 0;						\
365 	    Variable_shift_double(Dallp1(srcdstA),Dallp2(srcdstB),shift, \
367 	    Dallp1(srcdstA) >>= shift;					\
378 	    Dallp2(srcdstB) = (Dallp1(srcdstA) << 32 - (shift)) |	\
380 	    Dallp1(srcdstA) = Dallp1(srcdstA) >> shift
388     if( Dallp2(rightb) > Dallp2(leftb) ) Dallp1(lefta)--;	\
390     Dallp1(resulta) = Dallp1(lefta) - Dallp1(righta)
398         if((Dallp2(resultb)--) == 0) Dallp1(resulta)--;			\
404     Dallp1(resulta) = Dallp1(lefta) + Dallp1(righta);			\
406 	Dallp1(resulta)++
409     result = Dallp1(left) XOR Dallp1(right)
412     Dallp1(result) = left XOR Dallp1(right)
421     Dallp1(desta) = ((DBL_EMAX+DBL_BIAS)+1)<< (32-(1+DBL_EXP_LENGTH))	\
425     Dallp1(desta) = ((DBL_EMAX+DBL_BIAS)+1)<< (32-(1+DBL_EXP_LENGTH))	\
449 	if ((src1dstB) + (src2B) < (src1dstB)) Dallp1(src1dstA)++; \
450 	Dallp1(src1dstA) += (src2A);				\
459 	if ((src1dstB) < (src2B)) Dallp1(src1dstA)--;		\
460 	Dallp1(src1dstA) -= (src2A);				\
497 		Dallp1(opndp1) >>= 1-exponent;				\
500 		Dallp2(opndp2) = Dallp1(opndp1);			\
505 	    guard = (Dallp1(opndp1) >> -32-exponent) & 1;		\
507 	    else sticky |= (Dallp2(opndp2) | Dallp1(opndp1) << 64+exponent); \
508 	    Dallp2(opndp2) = Dallp1(opndp1) >> -31-exponent;		\
515 	sticky |= (Dallp1(opndp1) | Dallp2(opndp2));			\
719 	Dextallp1(dest1) = Dallp1(src1); Dextallp2(dest2) = Dallp2(src2); \