28 #define _FP_FRAC_DECL_2(X)	_FP_W_TYPE X##_f0 = 0, X##_f1 = 0  argument
30 #define _FP_FRAC_SET_2(X,I)	__FP_FRAC_SET_2(X, I)  argument
31 #define _FP_FRAC_HIGH_2(X)	(X##_f1)  argument
32 #define _FP_FRAC_LOW_2(X)	(X##_f0)  argument
33 #define _FP_FRAC_WORD_2(X,w)	(X##_f##w)  argument
34 #define _FP_FRAC_SLL_2(X, N) (						       \  argument
38 			X##_f1 = X##_f1 + X##_f1 +			       \
39 				(((_FP_WS_TYPE) (X##_f0)) < 0);		       \
40 			X##_f0 += X##_f0;				       \
42 			X##_f1 = X##_f1 << (N) | X##_f0 >>		       \
44 			X##_f0 <<= (N);					       \
49 	      X##_f1 = X##_f0 << ((N) - _FP_W_TYPE_SIZE);		       \
50 	      X##_f0 = 0;						       \
54 #define _FP_FRAC_SRL_2(X, N) (						       \  argument
57 	      X##_f0 = X##_f0 >> (N) | X##_f1 << (_FP_W_TYPE_SIZE - (N));      \
58 	      X##_f1 >>= (N);						       \
61 	      X##_f0 = X##_f1 >> ((N) - _FP_W_TYPE_SIZE);		       \
62 	      X##_f1 = 0;						       \
67 #define _FP_FRAC_SRS_2(X, N, sz) (					       \  argument
70 	      X##_f0 = (X##_f1 << (_FP_W_TYPE_SIZE - (N)) | X##_f0 >> (N)      \
72 			   ? X##_f0 & 1					       \
73 			   : (X##_f0 << (_FP_W_TYPE_SIZE - (N))) != 0));       \
74 		X##_f1 >>= (N);						       \
77 	      X##_f0 = (X##_f1 >> ((N) - _FP_W_TYPE_SIZE)		       \
80 			     : (X##_f1 << (2*_FP_W_TYPE_SIZE - (N))))          \
81 			    | X##_f0) != 0));				       \
82 	      X##_f1 = 0;						       \
85 #define _FP_FRAC_ADDI_2(X,I)	\  argument
86   __FP_FRAC_ADDI_2(X##_f1, X##_f0, I)
88 #define _FP_FRAC_ADD_2(R,X,Y)	\  argument
89   __FP_FRAC_ADD_2(R##_f1, R##_f0, X##_f1, X##_f0, Y##_f1, Y##_f0)
91 #define _FP_FRAC_SUB_2(R,X,Y)	\  argument
92   __FP_FRAC_SUB_2(R##_f1, R##_f0, X##_f1, X##_f0, Y##_f1, Y##_f0)
94 #define _FP_FRAC_DEC_2(X,Y)	\  argument
95   __FP_FRAC_DEC_2(X##_f1, X##_f0, Y##_f1, Y##_f0)
97 #define _FP_FRAC_CLZ_2(R,X)	\  argument
99     if (X##_f1)			\
100       __FP_CLZ(R,X##_f1);	\
103       __FP_CLZ(R,X##_f0);	\
109 #define _FP_FRAC_NEGP_2(X)	((_FP_WS_TYPE)X##_f1 < 0)  argument
110 #define _FP_FRAC_ZEROP_2(X)	((X##_f1 | X##_f0) == 0)  argument
111 #define _FP_FRAC_OVERP_2(fs,X)	(_FP_FRAC_HIGH_##fs(X) & _FP_OVERFLOW_##fs)  argument
112 #define _FP_FRAC_CLEAR_OVERP_2(fs,X)	(_FP_FRAC_HIGH_##fs(X) &= ~_FP_OVERFLOW_##fs)  argument
113 #define _FP_FRAC_EQ_2(X, Y)	(X##_f1 == Y##_f1 && X##_f0 == Y##_f0)  argument
114 #define _FP_FRAC_GT_2(X, Y)	\  argument
115   (X##_f1 > Y##_f1 || (X##_f1 == Y##_f1 && X##_f0 > Y##_f0))
116 #define _FP_FRAC_GE_2(X, Y)	\  argument
117   (X##_f1 > Y##_f1 || (X##_f1 == Y##_f1 && X##_f0 >= Y##_f0))
127 #define __FP_FRAC_SET_2(X,I1,I0)	(X##_f0 = I0, X##_f1 = I1)  argument
180 #define _FP_UNPACK_RAW_2(fs, X, val)			\  argument
184     X##_f0 = _flo.bits.frac0;				\
185     X##_f1 = _flo.bits.frac1;				\
186     X##_e  = _flo.bits.exp;				\
187     X##_s  = _flo.bits.sign;				\
190 #define _FP_UNPACK_RAW_2_P(fs, X, val)			\  argument
195     X##_f0 = _flo->bits.frac0;				\
196     X##_f1 = _flo->bits.frac1;				\
197     X##_e  = _flo->bits.exp;				\
198     X##_s  = _flo->bits.sign;				\
206 #define _FP_PACK_RAW_2(fs, val, X)			\  argument
210     _flo.bits.frac0 = X##_f0;				\
211     _flo.bits.frac1 = X##_f1;				\
212     _flo.bits.exp   = X##_e;				\
213     _flo.bits.sign  = X##_s;				\
218 #define _FP_PACK_RAW_2_P(fs, val, X)			\  argument
223     _flo->bits.frac0 = X##_f0;				\
224     _flo->bits.frac1 = X##_f1;				\
225     _flo->bits.exp   = X##_e;				\
226     _flo->bits.sign  = X##_s;				\
236 #define _FP_MUL_MEAT_2_wide(wfracbits, R, X, Y, doit)			\  argument
240     doit(_FP_FRAC_WORD_4(_z,1), _FP_FRAC_WORD_4(_z,0), X##_f0, Y##_f0);	\
241     doit(_b_f1, _b_f0, X##_f0, Y##_f1);					\
242     doit(_c_f1, _c_f0, X##_f1, Y##_f0);					\
243     doit(_FP_FRAC_WORD_4(_z,3), _FP_FRAC_WORD_4(_z,2), X##_f1, Y##_f1);	\
266 #define _FP_MUL_MEAT_2_wide_3mul(wfracbits, R, X, Y, doit)		\  argument
272     _b_f0 = X##_f0 + X##_f1;						\
273     _c1 = _b_f0 < X##_f0;						\
276     doit(_d, _FP_FRAC_WORD_4(_z,0), X##_f0, Y##_f0);			\
278     doit(_c_f1, _c_f0, X##_f1, Y##_f1);					\
306 #define _FP_MUL_MEAT_2_gmp(wfracbits, R, X, Y)				\  argument
310     _x[0] = X##_f0; _x[1] = X##_f1;					\
332 #define _FP_MUL_MEAT_2_120_240_double(wfracbits, R, X, Y, setfetz, resetfe)	\  argument
357     _e240 = (double)(long)(X##_f0 & 0xffffff);					\
359     _d240 = (double)(long)((X##_f0 >> 24) & 0xffffff);				\
361     _c240 = (double)(long)(((X##_f1 << 16) & 0xffffff) | (X##_f0 >> 48));	\
363     _b240 = (double)(long)((X##_f1 >> 8) & 0xffffff);				\
365     _a240 = (double)(long)(X##_f1 >> 32);					\
423 #define _FP_DIV_MEAT_2_udiv(fs, R, X, Y)				\  argument
426     if (_FP_FRAC_GT_2(X, Y))						\
428 	_n_f2 = X##_f1 >> 1;						\
429 	_n_f1 = X##_f1 << (_FP_W_TYPE_SIZE - 1) | X##_f0 >> 1;		\
430 	_n_f0 = X##_f0 << (_FP_W_TYPE_SIZE - 1);			\
435 	_n_f2 = X##_f1;							\
436 	_n_f1 = X##_f0;							\
491 #define _FP_DIV_MEAT_2_gmp(fs, R, X, Y)					\  argument
496     if (_FP_FRAC_GT_2(X, Y))						\
499 	_x[1] = (X##_f0 << (_FP_WFRACBITS_##fs-1 - _FP_W_TYPE_SIZE) |	\
500 		 X##_f1 >> (_FP_W_TYPE_SIZE -				\
502 	_x[2] = X##_f1 << (_FP_WFRACBITS_##fs-1 - _FP_W_TYPE_SIZE);	\
506 	_x[1] = (X##_f0 << (_FP_WFRACBITS_##fs - _FP_W_TYPE_SIZE) |	\
507 		 X##_f1 >> (_FP_W_TYPE_SIZE -				\
509 	_x[2] = X##_f1 << (_FP_WFRACBITS_##fs - _FP_W_TYPE_SIZE);	\
524 #define _FP_SQRT_MEAT_2(R, S, T, X, q)			\  argument
529 	if (T##_f1 <= X##_f1)				\
532 	    X##_f1 -= T##_f1;				\
535 	_FP_FRAC_SLL_2(X, 1);				\
543 	if (T##_f1 < X##_f1 || 				\
544 	    (T##_f1 == X##_f1 && T##_f0 <= X##_f0))	\
548 	    _FP_FRAC_DEC_2(X, T);			\
551 	_FP_FRAC_SLL_2(X, 1);				\
554     if (X##_f0 | X##_f1)				\
556 	if (S##_f1 < X##_f1 || 				\
557 	    (S##_f1 == X##_f1 && S##_f0 < X##_f0))	\
569 #define _FP_FRAC_ASSEMBLE_2(r, X, rsize)	\  argument
571 		? ({ (r) = X##_f0; })		\
573 		     (r) = X##_f1;		\
575 		     (r) += X##_f0;		\
578 #define _FP_FRAC_DISASSEMBLE_2(X, r, rsize)				\  argument
580     X##_f0 = r;								\
581     X##_f1 = (rsize <= _FP_W_TYPE_SIZE ? 0 : r >> _FP_W_TYPE_SIZE);	\