48  *	v = 2^n v1  *  v0
52  *	uv = 2^2n u1 v1  +  2^n u1 v0  +  2^n v1 u0  +  u0 v0
53  *	   = 2^2n u1 v1  +     2^n (u1 v0 + v1 u0)   +  u0 v0
55  * Now add 2^n u1 v1 to the first term and subtract it from the middle,
59  *	uv = (2^2n + 2^n) (u1 v1)  +
60  *	         (2^n)    (u1 v0 - u1 v1 + u0 v1 - u0 v0)  +
65  *	uv = (2^2n + 2^n) (u1 v1)  +			[u1v1 = high]
66  *		 (2^n)    (u1 - u0) (v0 - v1)  +	[(u1-u0)... = mid]
69  * The terms (u1 v1), (u1 - u0) (v0 - v1), and (u0 v0) can all be done
71  * of (u1 - u0) or (v0 - v1) may be negative.)
105 #define	v1	v.ul[H]  in __muldi3()  macro
110 	 * u1, u0, v1, and v0 will be directly accessible through the  in __muldi3()
122 	if (u1 == 0 && v1 == 0) {  in __muldi3()
124 		 * An (I hope) important optimization occurs when u1 and v1  in __muldi3()
142 		if (v0 >= v1)  in __muldi3()
143 			vdiff = v0 - v1;  in __muldi3()
145 			vdiff = v1 - v0, negmid ^= 1;  in __muldi3()
148 		high = u1 * v1;  in __muldi3()
160 #undef v1  in __muldi3()
184 	u_long u1, u0, v1, v0, udiff, vdiff, high, mid, low;  in __lmulq()  local
191 	v1 = HHALF(v);  in __lmulq()
197 	if (u1 == 0 && v1 == 0)  in __lmulq()
204 	if (v0 >= v1)  in __lmulq()
205 		vdiff = v0 - v1;  in __lmulq()
207 		vdiff = v1 - v0, neg ^= 1;  in __lmulq()
210 	high = u1 * v1;  in __lmulq()