1 #define _FP_DECL(wc, X)			\  argument
3   _FP_FRAC_DECL_##wc(X)
10 #define _FP_UNPACK_CANONICAL(fs, wc, X)					\  argument
15     _FP_FRAC_HIGH_##wc(X) |= _FP_IMPLBIT_##fs;				\
16     _FP_FRAC_SLL_##wc(X, _FP_WORKBITS);					\
22     if (_FP_FRAC_ZEROP_##wc(X))						\
28 	_FP_FRAC_CLZ_##wc(_shift, X);					\
30 	_FP_FRAC_SLL_##wc(X, (_shift+_FP_WORKBITS));			\
37     if (_FP_FRAC_ZEROP_##wc(X))						\
54 #define _FP_PACK_CANONICAL(fs, wc, X)				\  argument
62 	__ret |= _FP_ROUND(wc, X);				\
63 	if (_FP_FRAC_OVERP_##wc(fs, X))				\
65 	    _FP_FRAC_SRL_##wc(X, (_FP_WORKBITS+1));		\
69 	  _FP_FRAC_SRL_##wc(X, _FP_WORKBITS);			\
74 	    _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);		\
84 	    _FP_FRAC_SRS_##wc(X, X##_e, _FP_WFRACBITS_##fs);	\
85 	    __ret |= _FP_ROUND(wc, X);				\
86 	    _FP_FRAC_SLL_##wc(X, 1);				\
87 	    if (_FP_FRAC_OVERP_##wc(fs, X))			\
90 	        _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);	\
95 		_FP_FRAC_SRL_##wc(X, _FP_WORKBITS+1);		\
103 	    _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);		\
111     _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);			\
116     _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);			\
123 	_FP_FRAC_SET_##wc(X, _FP_NANFRAC_##fs);			\
127       _FP_FRAC_HIGH_##wc(X) |= _FP_QNANBIT_##fs;		\
138 #define _FP_ADD(fs, wc, R, X, Y)					     \  argument
151 	    _FP_FRAC_SRS_##wc(X, diff, _FP_WFRACBITS_##fs);		     \
152 	  else if (!_FP_FRAC_ZEROP_##wc(X))				     \
153 	    _FP_FRAC_SET_##wc(X, _FP_MINFRAC_##wc);			     \
155 	    _FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);			     \
163 	        _FP_FRAC_SRS_##wc(Y, diff, _FP_WFRACBITS_##fs);		     \
164 	      else if (!_FP_FRAC_ZEROP_##wc(Y))				     \
165 	        _FP_FRAC_SET_##wc(Y, _FP_MINFRAC_##wc);			     \
167 	        _FP_FRAC_SET_##wc(Y, _FP_ZEROFRAC_##wc);		     \
177 	  _FP_FRAC_ADD_##wc(R, X, Y);					     \
178 	  if (_FP_FRAC_OVERP_##wc(fs, R))				     \
180 	      _FP_FRAC_SRS_##wc(R, 1, _FP_WFRACBITS_##fs);		     \
187 	  _FP_FRAC_SUB_##wc(R, X, Y);					     \
188 	  if (_FP_FRAC_ZEROP_##wc(R))					     \
199 	      if (_FP_FRAC_NEGP_##wc(R))				     \
201 		  _FP_FRAC_SUB_##wc(R, Y, X);				     \
206 	      _FP_FRAC_CLZ_##wc(diff, R);				     \
211 		  _FP_FRAC_SLL_##wc(R, diff);				     \
219     _FP_CHOOSENAN(fs, wc, R, X, Y);					     \
227     _FP_FRAC_COPY_##wc(R, X);						     \
237     _FP_FRAC_COPY_##wc(R, Y);						     \
246 	_FP_FRAC_SET_##wc(R, _FP_NANFRAC_##fs);				     \
285 #define _FP_NEG(fs, wc, R, X)		\  argument
287     _FP_FRAC_COPY_##wc(R, X);		\
298 #define _FP_MUL(fs, wc, R, X, Y)			\  argument
309     if (_FP_FRAC_OVERP_##wc(fs, R))			\
310       _FP_FRAC_SRS_##wc(R, 1, _FP_WFRACBITS_##fs);	\
316     _FP_CHOOSENAN(fs, wc, R, X, Y);			\
328     _FP_FRAC_COPY_##wc(R, X);				\
339     _FP_FRAC_COPY_##wc(R, Y);				\
346     _FP_FRAC_SET_##wc(R, _FP_NANFRAC_##fs);		\
359 #define _FP_DIV(fs, wc, R, X, Y)			\  argument
372     _FP_CHOOSENAN(fs, wc, R, X, Y);			\
379     _FP_FRAC_COPY_##wc(R, X);				\
387     _FP_FRAC_COPY_##wc(R, Y);				\
406     _FP_FRAC_SET_##wc(R, _FP_NANFRAC_##fs);		\
420 #define _FP_CMP(fs, wc, ret, X, Y, un)					\  argument
423     if ((X##_e == _FP_EXPMAX_##fs && !_FP_FRAC_ZEROP_##wc(X))		\
424 	|| (Y##_e == _FP_EXPMAX_##fs && !_FP_FRAC_ZEROP_##wc(Y)))	\
430         int __x_zero = (!X##_e && _FP_FRAC_ZEROP_##wc(X)) ? 1 : 0;	\
431         int __y_zero = (!Y##_e && _FP_FRAC_ZEROP_##wc(Y)) ? 1 : 0;	\
445 	else if (_FP_FRAC_GT_##wc(X, Y))				\
447 	else if (_FP_FRAC_GT_##wc(Y, X))				\
457 #define _FP_CMP_EQ(fs, wc, ret, X, Y)					  \  argument
460     if ((X##_e == _FP_EXPMAX_##fs && !_FP_FRAC_ZEROP_##wc(X))		  \
461 	|| (Y##_e == _FP_EXPMAX_##fs && !_FP_FRAC_ZEROP_##wc(Y)))	  \
468 		&& _FP_FRAC_EQ_##wc(X, Y)				  \
469 		&& (X##_s == Y##_s || !X##_e && _FP_FRAC_ZEROP_##wc(X))); \
477 #define _FP_SQRT(fs, wc, R, X)						\  argument
479     _FP_FRAC_DECL_##wc(T); _FP_FRAC_DECL_##wc(S);			\
486     	_FP_FRAC_SET_##wc(X, _FP_ZEROFRAC_##wc);			\
513           _FP_FRAC_SLL_##wc(X, 1);					\
515         _FP_FRAC_SET_##wc(S, _FP_ZEROFRAC_##wc);			\
516         _FP_FRAC_SET_##wc(R, _FP_ZEROFRAC_##wc);			\
518         _FP_FRAC_SLL_##wc(X, 1);					\
519         _FP_SQRT_MEAT_##wc(R, S, T, X, q);				\
520         _FP_FRAC_SRL_##wc(R, 1);					\
551 #define _FP_TO_INT(fs, wc, r, X, rsize, rsigned)			\  argument
581 	    if (_FP_W_TYPE_SIZE*wc < rsize)				\
583 		_FP_FRAC_ASSEMBLE_##wc(r, X, rsize);			\
589 		  _FP_FRAC_SLL_##wc(X, (X##_e - _FP_WFRACBITS_##fs + 1));\
591 		  _FP_FRAC_SRL_##wc(X, (_FP_WFRACBITS_##fs - X##_e - 1));\
592 		_FP_FRAC_ASSEMBLE_##wc(r, X, rsize);			\
601 #define _FP_FROM_INT(fs, wc, X, r, rsize, rtype)			\  argument
624 	_FP_FRAC_DISASSEMBLE_##wc(X, ((unsigned rtype)r), rsize);	\
625 	_FP_FRAC_SLL_##wc(X, (_FP_WFRACBITS_##fs - X##_e - 1));		\