config/libbid/bid128_sqrt.c

35      UINT256 M256, C256, C4, C8;
139   __mul_128x128_to_256 (C256, CX1, TP128);
142   __mul_128x128_to_256 (C256, CX, T128);
146   // 4*C256
147 C4.w[3] = (C256.w[3] << 2) | (C256.w[2] >> 62);
148 C4.w[2] = (C256.w[2] << 2) | (C256.w[1] >> 62);
149 C4.w[1] = (C256.w[1] << 2) | (C256.w[0] >> 62);
150 C4.w[0] = C256.w[0] << 2;
152 long_sqrt128 (&CS, C256);
186     // if CSM' > C256, round up
206   if (M256.w[3] > C256.w[3]
207       || (M256.w[3] == C256.w[3]
208 	  && (M256.w[2] > C256.w[2]
209 	      || (M256.w[2] == C256.w[2]
210 		  && (M256.w[1] > C256.w[1]
211 		      || (M256.w[1] == C256.w[1]
212 			  && M256.w[0] > C256.w[0])))))) {
231     if (M256.w[3] > C256.w[3]
232 	|| (M256.w[3] == C256.w[3]
233 	    && (M256.w[2] > C256.w[2]
234 		|| (M256.w[2] == C256.w[2]
235 		    && (M256.w[1] > C256.w[1]
236 			|| (M256.w[1] == C256.w[1]
237 			    && M256.w[0] > C256.w[0])))))) {
259     if (M256.w[3] < C256.w[3]
260 	|| (M256.w[3] == C256.w[3]
261 	    && (M256.w[2] < C256.w[2]
262 		|| (M256.w[2] == C256.w[2]
263 		    && (M256.w[1] < C256.w[1]
264 			|| (M256.w[1] == C256.w[1]
265 			    && M256.w[0] <= C256.w[0])))))) {
298      UINT256 M256, C256, C4, C8;
409   __mul_128x128_to_256 (C256, CX1, TP128);
412   __mul_128x128_to_256 (C256, CX, T128);
416 	   // 4*C256
417 C4.w[3] = (C256.w[3] << 2) | (C256.w[2] >> 62);
418 C4.w[2] = (C256.w[2] << 2) | (C256.w[1] >> 62);
419 C4.w[1] = (C256.w[1] << 2) | (C256.w[0] >> 62);
420 C4.w[0] = C256.w[0] << 2;
422 long_sqrt128 (&CS, C256);
456     // if CSM' > C256, round up
476   if (M256.w[3] > C256.w[3]
477       || (M256.w[3] == C256.w[3]
478 	  && (M256.w[2] > C256.w[2]
479 	      || (M256.w[2] == C256.w[2]
480 		  && (M256.w[1] > C256.w[1]
481 		      || (M256.w[1] == C256.w[1]
482 			  && M256.w[0] > C256.w[0])))))) {
501     if (M256.w[3] > C256.w[3]
502 	|| (M256.w[3] == C256.w[3]
503 	    && (M256.w[2] > C256.w[2]
504 		|| (M256.w[2] == C256.w[2]
505 		    && (M256.w[1] > C256.w[1]
506 			|| (M256.w[1] == C256.w[1]
507 			    && M256.w[0] > C256.w[0])))))) {
529     if (M256.w[3] < C256.w[3]
530 	|| (M256.w[3] == C256.w[3]
531 	    && (M256.w[2] < C256.w[2]
532 		|| (M256.w[2] == C256.w[2]
533 		    && (M256.w[1] < C256.w[1]
534 			|| (M256.w[1] == C256.w[1]
535 			    && M256.w[0] <= C256.w[0])))))) {