contrib/gdtoa/strtordd.c

36 ULtodd(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
38 ULtodd(ULong *L, ULong *bits, Long exp, int k)
50 		L[_1] = (bits[1] >> 21 | bits[2] << 11) & (ULong)0xffffffffL;
51 		L[_0] = (bits[2] >> 21) | (bits[3] << 11 & 0xfffff)
54 		if (bits[1] &= 0x1fffff) {
55 			i = hi0bits(bits[1]) - 11;
63 				bits[1] = bits[1] << i | bits[0] >> (32-i);
64 				bits[0] = bits[0] << i & (ULong)0xffffffffL;
67 		else if (bits[0]) {
68 			i = hi0bits(bits[0]) + 21;
76 				bits[1] = bits[0] >> (32 - i);
77 				bits[0] = bits[0] << i & (ULong)0xffffffffL;
80 				bits[1] = bits[0] << (i - 32);
81 				bits[0] = 0;
88 		L[2+_1] = bits[0];
89 		L[2+_0] = (bits[1] & 0xfffff) | (exp << 20);
93 		if (bits[3])
95 		if (bits[2])
97 		if (bits[1] & 0xffe00000)
101 		L[_1] = bits[0];
102 		L[_0] = bits[1];
106 		i = hi0bits(bits[3]) - 11;	/* i >= 12 */
108 		L[_0] = ((bits[3] << i | bits[2] >> j) & 0xfffff)
110 		L[_1] = (bits[2] << i | bits[1] >> j) & 0xffffffffL;
111 		L[2+_0] = bits[1] & (((ULong)1L << j) - 1);
112 		L[2+_1] = bits[0];
116 		i = hi0bits(bits[2]) - 11;
120 			L[_0] = (bits[2] >> j & 0xfffff) | ((33 + j) << 20);
121 			L[_1] = (bits[2] << i | bits[1] >> j) & 0xffffffffL;
122 			L[2+_0] = bits[1] & (((ULong)1L << j) - 1);
123 			L[2+_1] = bits[0];
127 			L[_0] = (bits[2] & 0xfffff) | (33 << 20);
128 			L[_1] = bits[1];
130 			L[2+_1] = bits[0];
134 		L[_0] = (((bits[2] << i) | (bits[1] >> j)) & 0xfffff)
136 		L[_1] = (bits[1] << i | bits[0] >> j) & 0xffffffffL;
138 		L[2+_1] = bits[0] & ((1L << j) - 1);
142 		j = 11 - hi0bits(bits[1]);
144 		L[_0] = (bits[1] >> j & 0xfffff) | ((j + 1) << 20);
145 		L[_1] = (bits[1] << i | bits[0] >> j) & 0xffffffffL;
147 		L[2+_1] = bits[0] & (((ULong)1L << j) - 1);
161 		L[_1] = (bits[1] >> 21 | bits[2] << 11) & (ULong)0xffffffffL;
162 		L[_0] = bits[2] >> 21 | bits[3] << 11
164 		L[2+_1] = bits[0];
165 		L[2+_0] = bits[1] | (ULong)0x7ff00000L;
186 	ULong bits[4];
196 	k = strtodg(s, sp, fpi, &exp, bits);
197 	ULtodd((ULong*)dd, bits, exp, k);