VisionFive2 Linux kernel

StarFive Tech Linux Kernel for VisionFive (JH7110) boards (mirror)

More than 9999 Commits   32 Branches   54 Tags
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   1) /* Software floating-point emulation.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   2)    Basic two-word fraction declaration and manipulation.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   3)    Copyright (C) 1997,1998,1999 Free Software Foundation, Inc.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   4)    This file is part of the GNU C Library.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   5)    Contributed by Richard Henderson (rth@cygnus.com),
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   6) 		  Jakub Jelinek (jj@ultra.linux.cz),
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   7) 		  David S. Miller (davem@redhat.com) and
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   8) 		  Peter Maydell (pmaydell@chiark.greenend.org.uk).
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700   9) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  10)    The GNU C Library is free software; you can redistribute it and/or
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  11)    modify it under the terms of the GNU Library General Public License as
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  12)    published by the Free Software Foundation; either version 2 of the
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  13)    License, or (at your option) any later version.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  14) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  15)    The GNU C Library is distributed in the hope that it will be useful,
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  16)    but WITHOUT ANY WARRANTY; without even the implied warranty of
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  17)    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  18)    Library General Public License for more details.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  19) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  20)    You should have received a copy of the GNU Library General Public
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  21)    License along with the GNU C Library; see the file COPYING.LIB.  If
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  22)    not, write to the Free Software Foundation, Inc.,
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  23)    59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  24) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  25) #ifndef __MATH_EMU_OP_2_H__
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  26) #define __MATH_EMU_OP_2_H__
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  27) 
40d3057ac036f (Kumar Gala     2008-06-27 09:33:59 -0500  28) #define _FP_FRAC_DECL_2(X)	_FP_W_TYPE X##_f0 = 0, X##_f1 = 0
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  29) #define _FP_FRAC_COPY_2(D,S)	(D##_f0 = S##_f0, D##_f1 = S##_f1)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  30) #define _FP_FRAC_SET_2(X,I)	__FP_FRAC_SET_2(X, I)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  31) #define _FP_FRAC_HIGH_2(X)	(X##_f1)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  32) #define _FP_FRAC_LOW_2(X)	(X##_f0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  33) #define _FP_FRAC_WORD_2(X,w)	(X##_f##w)
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  34) #define _FP_FRAC_SLL_2(X, N) (						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  35) 	(void) (((N) < _FP_W_TYPE_SIZE)					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  36) 	  ? ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  37) 		if (__builtin_constant_p(N) && (N) == 1) {		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  38) 			X##_f1 = X##_f1 + X##_f1 +			       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  39) 				(((_FP_WS_TYPE) (X##_f0)) < 0);		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  40) 			X##_f0 += X##_f0;				       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  41) 		} else {						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  42) 			X##_f1 = X##_f1 << (N) | X##_f0 >>		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  43) 						(_FP_W_TYPE_SIZE - (N));       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  44) 			X##_f0 <<= (N);					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  45) 		}							       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  46) 		0;							       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  47) 	    })								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  48) 	  : ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  49) 	      X##_f1 = X##_f0 << ((N) - _FP_W_TYPE_SIZE);		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  50) 	      X##_f0 = 0;						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  51) 	  })))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  52) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  53) 
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  54) #define _FP_FRAC_SRL_2(X, N) (						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  55) 	(void) (((N) < _FP_W_TYPE_SIZE)					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  56) 	  ? ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  57) 	      X##_f0 = X##_f0 >> (N) | X##_f1 << (_FP_W_TYPE_SIZE - (N));      \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  58) 	      X##_f1 >>= (N);						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  59) 	    })								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  60) 	  : ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  61) 	      X##_f0 = X##_f1 >> ((N) - _FP_W_TYPE_SIZE);		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  62) 	      X##_f1 = 0;						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  63) 	    })))
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  64) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  65) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  66) /* Right shift with sticky-lsb.  */
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  67) #define _FP_FRAC_SRS_2(X, N, sz) (					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  68) 	(void) (((N) < _FP_W_TYPE_SIZE)					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  69) 	  ? ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  70) 	      X##_f0 = (X##_f1 << (_FP_W_TYPE_SIZE - (N)) | X##_f0 >> (N)      \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  71) 			| (__builtin_constant_p(N) && (N) == 1		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  72) 			   ? X##_f0 & 1					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  73) 			   : (X##_f0 << (_FP_W_TYPE_SIZE - (N))) != 0));       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  74) 		X##_f1 >>= (N);						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  75) 	    })								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  76) 	  : ({								       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  77) 	      X##_f0 = (X##_f1 >> ((N) - _FP_W_TYPE_SIZE)		       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  78) 			| ((((N) == _FP_W_TYPE_SIZE			       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  79) 			     ? 0					       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  80) 			     : (X##_f1 << (2*_FP_W_TYPE_SIZE - (N))))          \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  81) 			    | X##_f0) != 0));				       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  82) 	      X##_f1 = 0;						       \
7adb3e998f5be (Vincent Chen   2018-11-22 11:14:37 +0800  83) 	    })))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  84) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  85) #define _FP_FRAC_ADDI_2(X,I)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  86)   __FP_FRAC_ADDI_2(X##_f1, X##_f0, I)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  87) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  88) #define _FP_FRAC_ADD_2(R,X,Y)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  89)   __FP_FRAC_ADD_2(R##_f1, R##_f0, X##_f1, X##_f0, Y##_f1, Y##_f0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  90) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  91) #define _FP_FRAC_SUB_2(R,X,Y)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  92)   __FP_FRAC_SUB_2(R##_f1, R##_f0, X##_f1, X##_f0, Y##_f1, Y##_f0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  93) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  94) #define _FP_FRAC_DEC_2(X,Y)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  95)   __FP_FRAC_DEC_2(X##_f1, X##_f0, Y##_f1, Y##_f0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  96) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  97) #define _FP_FRAC_CLZ_2(R,X)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  98)   do {				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700  99)     if (X##_f1)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 100)       __FP_CLZ(R,X##_f1);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 101)     else 			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 102)     {				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 103)       __FP_CLZ(R,X##_f0);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 104)       R += _FP_W_TYPE_SIZE;	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 105)     }				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 106)   } while(0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 107) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 108) /* Predicates */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 109) #define _FP_FRAC_NEGP_2(X)	((_FP_WS_TYPE)X##_f1 < 0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 110) #define _FP_FRAC_ZEROP_2(X)	((X##_f1 | X##_f0) == 0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 111) #define _FP_FRAC_OVERP_2(fs,X)	(_FP_FRAC_HIGH_##fs(X) & _FP_OVERFLOW_##fs)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 112) #define _FP_FRAC_CLEAR_OVERP_2(fs,X)	(_FP_FRAC_HIGH_##fs(X) &= ~_FP_OVERFLOW_##fs)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 113) #define _FP_FRAC_EQ_2(X, Y)	(X##_f1 == Y##_f1 && X##_f0 == Y##_f0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 114) #define _FP_FRAC_GT_2(X, Y)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 115)   (X##_f1 > Y##_f1 || (X##_f1 == Y##_f1 && X##_f0 > Y##_f0))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 116) #define _FP_FRAC_GE_2(X, Y)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 117)   (X##_f1 > Y##_f1 || (X##_f1 == Y##_f1 && X##_f0 >= Y##_f0))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 118) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 119) #define _FP_ZEROFRAC_2		0, 0
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 120) #define _FP_MINFRAC_2		0, 1
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 121) #define _FP_MAXFRAC_2		(~(_FP_WS_TYPE)0), (~(_FP_WS_TYPE)0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 122) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 123) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 124)  * Internals 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 125)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 126) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 127) #define __FP_FRAC_SET_2(X,I1,I0)	(X##_f0 = I0, X##_f1 = I1)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 128) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 129) #define __FP_CLZ_2(R, xh, xl)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 130)   do {				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 131)     if (xh)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 132)       __FP_CLZ(R,xh);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 133)     else 			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 134)     {				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 135)       __FP_CLZ(R,xl);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 136)       R += _FP_W_TYPE_SIZE;	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 137)     }				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 138)   } while(0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 139) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 140) #if 0
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 141) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 142) #ifndef __FP_FRAC_ADDI_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 143) #define __FP_FRAC_ADDI_2(xh, xl, i)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 144)   (xh += ((xl += i) < i))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 145) #endif
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 146) #ifndef __FP_FRAC_ADD_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 147) #define __FP_FRAC_ADD_2(rh, rl, xh, xl, yh, yl)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 148)   (rh = xh + yh + ((rl = xl + yl) < xl))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 149) #endif
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 150) #ifndef __FP_FRAC_SUB_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 151) #define __FP_FRAC_SUB_2(rh, rl, xh, xl, yh, yl)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 152)   (rh = xh - yh - ((rl = xl - yl) > xl))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 153) #endif
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 154) #ifndef __FP_FRAC_DEC_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 155) #define __FP_FRAC_DEC_2(xh, xl, yh, yl)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 156)   do {					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 157)     UWtype _t = xl;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 158)     xh -= yh + ((xl -= yl) > _t);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 159)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 160) #endif
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 161) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 162) #else
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 163) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 164) #undef __FP_FRAC_ADDI_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 165) #define __FP_FRAC_ADDI_2(xh, xl, i)	add_ssaaaa(xh, xl, xh, xl, 0, i)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 166) #undef __FP_FRAC_ADD_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 167) #define __FP_FRAC_ADD_2			add_ssaaaa
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 168) #undef __FP_FRAC_SUB_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 169) #define __FP_FRAC_SUB_2			sub_ddmmss
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 170) #undef __FP_FRAC_DEC_2
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 171) #define __FP_FRAC_DEC_2(xh, xl, yh, yl)	sub_ddmmss(xh, xl, xh, xl, yh, yl)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 172) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 173) #endif
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 174) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 175) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 176)  * Unpack the raw bits of a native fp value.  Do not classify or
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 177)  * normalize the data.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 178)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 179) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 180) #define _FP_UNPACK_RAW_2(fs, X, val)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 181)   do {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 182)     union _FP_UNION_##fs _flo; _flo.flt = (val);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 183) 							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 184)     X##_f0 = _flo.bits.frac0;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 185)     X##_f1 = _flo.bits.frac1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 186)     X##_e  = _flo.bits.exp;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 187)     X##_s  = _flo.bits.sign;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 188)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 189) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 190) #define _FP_UNPACK_RAW_2_P(fs, X, val)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 191)   do {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 192)     union _FP_UNION_##fs *_flo =			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 193)       (union _FP_UNION_##fs *)(val);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 194) 							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 195)     X##_f0 = _flo->bits.frac0;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 196)     X##_f1 = _flo->bits.frac1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 197)     X##_e  = _flo->bits.exp;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 198)     X##_s  = _flo->bits.sign;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 199)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 200) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 201) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 202) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 203)  * Repack the raw bits of a native fp value.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 204)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 205) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 206) #define _FP_PACK_RAW_2(fs, val, X)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 207)   do {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 208)     union _FP_UNION_##fs _flo;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 209) 							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 210)     _flo.bits.frac0 = X##_f0;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 211)     _flo.bits.frac1 = X##_f1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 212)     _flo.bits.exp   = X##_e;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 213)     _flo.bits.sign  = X##_s;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 214) 							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 215)     (val) = _flo.flt;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 216)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 217) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 218) #define _FP_PACK_RAW_2_P(fs, val, X)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 219)   do {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 220)     union _FP_UNION_##fs *_flo =			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 221)       (union _FP_UNION_##fs *)(val);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 222) 							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 223)     _flo->bits.frac0 = X##_f0;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 224)     _flo->bits.frac1 = X##_f1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 225)     _flo->bits.exp   = X##_e;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 226)     _flo->bits.sign  = X##_s;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 227)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 228) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 229) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 230) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 231)  * Multiplication algorithms:
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 232)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 233) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 234) /* Given a 1W * 1W => 2W primitive, do the extended multiplication.  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 235) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 236) #define _FP_MUL_MEAT_2_wide(wfracbits, R, X, Y, doit)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 237)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 238)     _FP_FRAC_DECL_4(_z); _FP_FRAC_DECL_2(_b); _FP_FRAC_DECL_2(_c);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 239) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 240)     doit(_FP_FRAC_WORD_4(_z,1), _FP_FRAC_WORD_4(_z,0), X##_f0, Y##_f0);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 241)     doit(_b_f1, _b_f0, X##_f0, Y##_f1);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 242)     doit(_c_f1, _c_f0, X##_f1, Y##_f0);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 243)     doit(_FP_FRAC_WORD_4(_z,3), _FP_FRAC_WORD_4(_z,2), X##_f1, Y##_f1);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 244) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 245)     __FP_FRAC_ADD_3(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 246) 		    _FP_FRAC_WORD_4(_z,1), 0, _b_f1, _b_f0,		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 247) 		    _FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 248) 		    _FP_FRAC_WORD_4(_z,1));				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 249)     __FP_FRAC_ADD_3(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 250) 		    _FP_FRAC_WORD_4(_z,1), 0, _c_f1, _c_f0,		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 251) 		    _FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 252) 		    _FP_FRAC_WORD_4(_z,1));				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 253) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 254)     /* Normalize since we know where the msb of the multiplicands	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 255)        were (bit B), we know that the msb of the of the product is	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 256)        at either 2B or 2B-1.  */					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 257)     _FP_FRAC_SRS_4(_z, wfracbits-1, 2*wfracbits);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 258)     R##_f0 = _FP_FRAC_WORD_4(_z,0);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 259)     R##_f1 = _FP_FRAC_WORD_4(_z,1);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 260)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 261) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 262) /* Given a 1W * 1W => 2W primitive, do the extended multiplication.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 263)    Do only 3 multiplications instead of four. This one is for machines
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 264)    where multiplication is much more expensive than subtraction.  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 265) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 266) #define _FP_MUL_MEAT_2_wide_3mul(wfracbits, R, X, Y, doit)		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 267)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 268)     _FP_FRAC_DECL_4(_z); _FP_FRAC_DECL_2(_b); _FP_FRAC_DECL_2(_c);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 269)     _FP_W_TYPE _d;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 270)     int _c1, _c2;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 271) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 272)     _b_f0 = X##_f0 + X##_f1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 273)     _c1 = _b_f0 < X##_f0;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 274)     _b_f1 = Y##_f0 + Y##_f1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 275)     _c2 = _b_f1 < Y##_f0;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 276)     doit(_d, _FP_FRAC_WORD_4(_z,0), X##_f0, Y##_f0);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 277)     doit(_FP_FRAC_WORD_4(_z,2), _FP_FRAC_WORD_4(_z,1), _b_f0, _b_f1);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 278)     doit(_c_f1, _c_f0, X##_f1, Y##_f1);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 279) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 280)     _b_f0 &= -_c2;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 281)     _b_f1 &= -_c1;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 282)     __FP_FRAC_ADD_3(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 283) 		    _FP_FRAC_WORD_4(_z,1), (_c1 & _c2), 0, _d,		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 284) 		    0, _FP_FRAC_WORD_4(_z,2), _FP_FRAC_WORD_4(_z,1));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 285)     __FP_FRAC_ADDI_2(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 286) 		     _b_f0);						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 287)     __FP_FRAC_ADDI_2(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 288) 		     _b_f1);						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 289)     __FP_FRAC_DEC_3(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 290) 		    _FP_FRAC_WORD_4(_z,1),				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 291) 		    0, _d, _FP_FRAC_WORD_4(_z,0));			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 292)     __FP_FRAC_DEC_3(_FP_FRAC_WORD_4(_z,3),_FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 293) 		    _FP_FRAC_WORD_4(_z,1), 0, _c_f1, _c_f0);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 294)     __FP_FRAC_ADD_2(_FP_FRAC_WORD_4(_z,3), _FP_FRAC_WORD_4(_z,2),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 295) 		    _c_f1, _c_f0,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 296) 		    _FP_FRAC_WORD_4(_z,3), _FP_FRAC_WORD_4(_z,2));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 297) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 298)     /* Normalize since we know where the msb of the multiplicands	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 299)        were (bit B), we know that the msb of the of the product is	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 300)        at either 2B or 2B-1.  */					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 301)     _FP_FRAC_SRS_4(_z, wfracbits-1, 2*wfracbits);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 302)     R##_f0 = _FP_FRAC_WORD_4(_z,0);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 303)     R##_f1 = _FP_FRAC_WORD_4(_z,1);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 304)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 305) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 306) #define _FP_MUL_MEAT_2_gmp(wfracbits, R, X, Y)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 307)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 308)     _FP_FRAC_DECL_4(_z);						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 309)     _FP_W_TYPE _x[2], _y[2];						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 310)     _x[0] = X##_f0; _x[1] = X##_f1;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 311)     _y[0] = Y##_f0; _y[1] = Y##_f1;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 312) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 313)     mpn_mul_n(_z_f, _x, _y, 2);						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 314) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 315)     /* Normalize since we know where the msb of the multiplicands	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 316)        were (bit B), we know that the msb of the of the product is	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 317)        at either 2B or 2B-1.  */					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 318)     _FP_FRAC_SRS_4(_z, wfracbits-1, 2*wfracbits);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 319)     R##_f0 = _z_f[0];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 320)     R##_f1 = _z_f[1];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 321)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 322) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 323) /* Do at most 120x120=240 bits multiplication using double floating
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 324)    point multiplication.  This is useful if floating point
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 325)    multiplication has much bigger throughput than integer multiply.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 326)    It is supposed to work for _FP_W_TYPE_SIZE 64 and wfracbits
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 327)    between 106 and 120 only.  
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 328)    Caller guarantees that X and Y has (1LLL << (wfracbits - 1)) set.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 329)    SETFETZ is a macro which will disable all FPU exceptions and set rounding
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 330)    towards zero,  RESETFE should optionally reset it back.  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 331) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 332) #define _FP_MUL_MEAT_2_120_240_double(wfracbits, R, X, Y, setfetz, resetfe)	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 333)   do {										\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 334)     static const double _const[] = {						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 335)       /* 2^-24 */ 5.9604644775390625e-08,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 336)       /* 2^-48 */ 3.5527136788005009e-15,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 337)       /* 2^-72 */ 2.1175823681357508e-22,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 338)       /* 2^-96 */ 1.2621774483536189e-29,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 339)       /* 2^28 */ 2.68435456e+08,						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 340)       /* 2^4 */ 1.600000e+01,							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 341)       /* 2^-20 */ 9.5367431640625e-07,						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 342)       /* 2^-44 */ 5.6843418860808015e-14,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 343)       /* 2^-68 */ 3.3881317890172014e-21,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 344)       /* 2^-92 */ 2.0194839173657902e-28,					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 345)       /* 2^-116 */ 1.2037062152420224e-35};					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 346)     double _a240, _b240, _c240, _d240, _e240, _f240, 				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 347) 	   _g240, _h240, _i240, _j240, _k240;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 348)     union { double d; UDItype i; } _l240, _m240, _n240, _o240,			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 349) 				   _p240, _q240, _r240, _s240;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 350)     UDItype _t240, _u240, _v240, _w240, _x240, _y240 = 0;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 351) 										\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 352)     if (wfracbits < 106 || wfracbits > 120)					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 353)       abort();									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 354) 										\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 355)     setfetz;									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 356) 										\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 357)     _e240 = (double)(long)(X##_f0 & 0xffffff);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 358)     _j240 = (double)(long)(Y##_f0 & 0xffffff);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 359)     _d240 = (double)(long)((X##_f0 >> 24) & 0xffffff);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 360)     _i240 = (double)(long)((Y##_f0 >> 24) & 0xffffff);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 361)     _c240 = (double)(long)(((X##_f1 << 16) & 0xffffff) | (X##_f0 >> 48));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 362)     _h240 = (double)(long)(((Y##_f1 << 16) & 0xffffff) | (Y##_f0 >> 48));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 363)     _b240 = (double)(long)((X##_f1 >> 8) & 0xffffff);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 364)     _g240 = (double)(long)((Y##_f1 >> 8) & 0xffffff);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 365)     _a240 = (double)(long)(X##_f1 >> 32);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 366)     _f240 = (double)(long)(Y##_f1 >> 32);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 367)     _e240 *= _const[3];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 368)     _j240 *= _const[3];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 369)     _d240 *= _const[2];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 370)     _i240 *= _const[2];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 371)     _c240 *= _const[1];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 372)     _h240 *= _const[1];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 373)     _b240 *= _const[0];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 374)     _g240 *= _const[0];								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 375)     _s240.d =							      _e240*_j240;\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 376)     _r240.d =						_d240*_j240 + _e240*_i240;\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 377)     _q240.d =				  _c240*_j240 + _d240*_i240 + _e240*_h240;\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 378)     _p240.d =		    _b240*_j240 + _c240*_i240 + _d240*_h240 + _e240*_g240;\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 379)     _o240.d = _a240*_j240 + _b240*_i240 + _c240*_h240 + _d240*_g240 + _e240*_f240;\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 380)     _n240.d = _a240*_i240 + _b240*_h240 + _c240*_g240 + _d240*_f240;		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 381)     _m240.d = _a240*_h240 + _b240*_g240 + _c240*_f240;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 382)     _l240.d = _a240*_g240 + _b240*_f240;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 383)     _k240 =   _a240*_f240;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 384)     _r240.d += _s240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 385)     _q240.d += _r240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 386)     _p240.d += _q240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 387)     _o240.d += _p240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 388)     _n240.d += _o240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 389)     _m240.d += _n240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 390)     _l240.d += _m240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 391)     _k240 += _l240.d;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 392)     _s240.d -= ((_const[10]+_s240.d)-_const[10]);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 393)     _r240.d -= ((_const[9]+_r240.d)-_const[9]);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 394)     _q240.d -= ((_const[8]+_q240.d)-_const[8]);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 395)     _p240.d -= ((_const[7]+_p240.d)-_const[7]);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 396)     _o240.d += _const[7];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 397)     _n240.d += _const[6];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 398)     _m240.d += _const[5];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 399)     _l240.d += _const[4];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 400)     if (_s240.d != 0.0) _y240 = 1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 401)     if (_r240.d != 0.0) _y240 = 1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 402)     if (_q240.d != 0.0) _y240 = 1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 403)     if (_p240.d != 0.0) _y240 = 1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 404)     _t240 = (DItype)_k240;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 405)     _u240 = _l240.i;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 406)     _v240 = _m240.i;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 407)     _w240 = _n240.i;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 408)     _x240 = _o240.i;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 409)     R##_f1 = (_t240 << (128 - (wfracbits - 1)))					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 410) 	     | ((_u240 & 0xffffff) >> ((wfracbits - 1) - 104));			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 411)     R##_f0 = ((_u240 & 0xffffff) << (168 - (wfracbits - 1)))			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 412)     	     | ((_v240 & 0xffffff) << (144 - (wfracbits - 1)))			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 413)     	     | ((_w240 & 0xffffff) << (120 - (wfracbits - 1)))			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 414)     	     | ((_x240 & 0xffffff) >> ((wfracbits - 1) - 96))			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 415)     	     | _y240;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 416)     resetfe;									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 417)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 418) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 419) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 420)  * Division algorithms:
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 421)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 422) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 423) #define _FP_DIV_MEAT_2_udiv(fs, R, X, Y)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 424)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 425)     _FP_W_TYPE _n_f2, _n_f1, _n_f0, _r_f1, _r_f0, _m_f1, _m_f0;		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 426)     if (_FP_FRAC_GT_2(X, Y))						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 427)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 428) 	_n_f2 = X##_f1 >> 1;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 429) 	_n_f1 = X##_f1 << (_FP_W_TYPE_SIZE - 1) | X##_f0 >> 1;		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 430) 	_n_f0 = X##_f0 << (_FP_W_TYPE_SIZE - 1);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 431)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 432)     else								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 433)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 434) 	R##_e--;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 435) 	_n_f2 = X##_f1;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 436) 	_n_f1 = X##_f0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 437) 	_n_f0 = 0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 438)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 439) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 440)     /* Normalize, i.e. make the most significant bit of the 		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 441)        denominator set. */						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 442)     _FP_FRAC_SLL_2(Y, _FP_WFRACXBITS_##fs);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 443) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 444)     udiv_qrnnd(R##_f1, _r_f1, _n_f2, _n_f1, Y##_f1);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 445)     umul_ppmm(_m_f1, _m_f0, R##_f1, Y##_f0);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 446)     _r_f0 = _n_f0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 447)     if (_FP_FRAC_GT_2(_m, _r))						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 448)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 449) 	R##_f1--;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 450) 	_FP_FRAC_ADD_2(_r, Y, _r);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 451) 	if (_FP_FRAC_GE_2(_r, Y) && _FP_FRAC_GT_2(_m, _r))		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 452) 	  {								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 453) 	    R##_f1--;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 454) 	    _FP_FRAC_ADD_2(_r, Y, _r);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 455) 	  }								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 456)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 457)     _FP_FRAC_DEC_2(_r, _m);						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 458) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 459)     if (_r_f1 == Y##_f1)						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 460)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 461) 	/* This is a special case, not an optimization			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 462) 	   (_r/Y##_f1 would not fit into UWtype).			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 463) 	   As _r is guaranteed to be < Y,  R##_f0 can be either		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 464) 	   (UWtype)-1 or (UWtype)-2.  But as we know what kind		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 465) 	   of bits it is (sticky, guard, round),  we don't care.	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 466) 	   We also don't care what the reminder is,  because the	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 467) 	   guard bit will be set anyway.  -jj */			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 468) 	R##_f0 = -1;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 469)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 470)     else								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 471)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 472) 	udiv_qrnnd(R##_f0, _r_f1, _r_f1, _r_f0, Y##_f1);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 473) 	umul_ppmm(_m_f1, _m_f0, R##_f0, Y##_f0);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 474) 	_r_f0 = 0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 475) 	if (_FP_FRAC_GT_2(_m, _r))					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 476) 	  {								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 477) 	    R##_f0--;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 478) 	    _FP_FRAC_ADD_2(_r, Y, _r);					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 479) 	    if (_FP_FRAC_GE_2(_r, Y) && _FP_FRAC_GT_2(_m, _r))		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 480) 	      {								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 481) 		R##_f0--;						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 482) 		_FP_FRAC_ADD_2(_r, Y, _r);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 483) 	      }								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 484) 	  }								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 485) 	if (!_FP_FRAC_EQ_2(_r, _m))					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 486) 	  R##_f0 |= _FP_WORK_STICKY;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 487)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 488)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 489) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 490) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 491) #define _FP_DIV_MEAT_2_gmp(fs, R, X, Y)					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 492)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 493)     _FP_W_TYPE _x[4], _y[2], _z[4];					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 494)     _y[0] = Y##_f0; _y[1] = Y##_f1;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 495)     _x[0] = _x[3] = 0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 496)     if (_FP_FRAC_GT_2(X, Y))						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 497)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 498) 	R##_e++;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 499) 	_x[1] = (X##_f0 << (_FP_WFRACBITS_##fs-1 - _FP_W_TYPE_SIZE) |	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 500) 		 X##_f1 >> (_FP_W_TYPE_SIZE -				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 501) 			    (_FP_WFRACBITS_##fs-1 - _FP_W_TYPE_SIZE)));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 502) 	_x[2] = X##_f1 << (_FP_WFRACBITS_##fs-1 - _FP_W_TYPE_SIZE);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 503)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 504)     else								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 505)       {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 506) 	_x[1] = (X##_f0 << (_FP_WFRACBITS_##fs - _FP_W_TYPE_SIZE) |	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 507) 		 X##_f1 >> (_FP_W_TYPE_SIZE -				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 508) 			    (_FP_WFRACBITS_##fs - _FP_W_TYPE_SIZE)));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 509) 	_x[2] = X##_f1 << (_FP_WFRACBITS_##fs - _FP_W_TYPE_SIZE);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 510)       }									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 511) 									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 512)     (void) mpn_divrem (_z, 0, _x, 4, _y, 2);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 513)     R##_f1 = _z[1];							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 514)     R##_f0 = _z[0] | ((_x[0] | _x[1]) != 0);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 515)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 516) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 517) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 518) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 519)  * Square root algorithms:
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 520)  * We have just one right now, maybe Newton approximation
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 521)  * should be added for those machines where division is fast.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 522)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 523)  
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 524) #define _FP_SQRT_MEAT_2(R, S, T, X, q)			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 525)   do {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 526)     while (q)						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 527)       {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 528) 	T##_f1 = S##_f1 + q;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 529) 	if (T##_f1 <= X##_f1)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 530) 	  {						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 531) 	    S##_f1 = T##_f1 + q;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 532) 	    X##_f1 -= T##_f1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 533) 	    R##_f1 += q;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 534) 	  }						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 535) 	_FP_FRAC_SLL_2(X, 1);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 536) 	q >>= 1;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 537)       }							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 538)     q = (_FP_W_TYPE)1 << (_FP_W_TYPE_SIZE - 1);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 539)     while (q != _FP_WORK_ROUND)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 540)       {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 541) 	T##_f0 = S##_f0 + q;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 542) 	T##_f1 = S##_f1;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 543) 	if (T##_f1 < X##_f1 || 				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 544) 	    (T##_f1 == X##_f1 && T##_f0 <= X##_f0))	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 545) 	  {						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 546) 	    S##_f0 = T##_f0 + q;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 547) 	    S##_f1 += (T##_f0 > S##_f0);		\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 548) 	    _FP_FRAC_DEC_2(X, T);			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 549) 	    R##_f0 += q;				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 550) 	  }						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 551) 	_FP_FRAC_SLL_2(X, 1);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 552) 	q >>= 1;					\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 553)       }							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 554)     if (X##_f0 | X##_f1)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 555)       {							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 556) 	if (S##_f1 < X##_f1 || 				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 557) 	    (S##_f1 == X##_f1 && S##_f0 < X##_f0))	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 558) 	  R##_f0 |= _FP_WORK_ROUND;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 559) 	R##_f0 |= _FP_WORK_STICKY;			\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 560)       }							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 561)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 562) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 563) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 564) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 565)  * Assembly/disassembly for converting to/from integral types.  
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 566)  * No shifting or overflow handled here.
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 567)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 568) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 569) #define _FP_FRAC_ASSEMBLE_2(r, X, rsize)	\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 570) 	(void) (((rsize) <= _FP_W_TYPE_SIZE)	\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 571) 		? ({ (r) = X##_f0; })		\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 572) 		: ({				\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 573) 		     (r) = X##_f1;		\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 574) 		     (r) <<= _FP_W_TYPE_SIZE;	\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 575) 		     (r) += X##_f0;		\
8183db10db637 (Vincent Chen   2019-05-27 14:17:21 +0800 576) 		    }))
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 577) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 578) #define _FP_FRAC_DISASSEMBLE_2(X, r, rsize)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 579)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 580)     X##_f0 = r;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 581)     X##_f1 = (rsize <= _FP_W_TYPE_SIZE ? 0 : r >> _FP_W_TYPE_SIZE);	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 582)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 583) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 584) /*
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 585)  * Convert FP values between word sizes
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 586)  */
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 587) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 588) #define _FP_FRAC_CONV_1_2(dfs, sfs, D, S)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 589)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 590)     if (S##_c != FP_CLS_NAN)						\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 591)       _FP_FRAC_SRS_2(S, (_FP_WFRACBITS_##sfs - _FP_WFRACBITS_##dfs),	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 592) 		     _FP_WFRACBITS_##sfs);				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 593)     else								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 594)       _FP_FRAC_SRL_2(S, (_FP_WFRACBITS_##sfs - _FP_WFRACBITS_##dfs));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 595)     D##_f = S##_f0;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 596)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 597) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 598) #define _FP_FRAC_CONV_2_1(dfs, sfs, D, S)				\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 599)   do {									\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 600)     D##_f0 = S##_f;							\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 601)     D##_f1 = 0;								\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 602)     _FP_FRAC_SLL_2(D, (_FP_WFRACBITS_##dfs - _FP_WFRACBITS_##sfs));	\
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 603)   } while (0)
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 604) 
^1da177e4c3f4 (Linus Torvalds 2005-04-16 15:20:36 -0700 605) #endif