config/libbid/bid128_sqrt.c

36      UINT128 CX, CX1, CX2, A10, S2, T128, TP128, CS, CSM, res;
160   CSM.w[1] = (CS.w[1] << 1) | (CS.w[0] >> 63);
161   CSM.w[0] = (CS.w[0] + CS.w[0]) | 1;
162   // CSM^2
163   //__mul_128x128_to_256(M256, CSM, CSM);
164   __sqr128_to_256 (M256, CSM);
180     // M256 - 8*CSM
186     // if CSM' > C256, round up
299      UINT128 CX, CX1, CX2, A10, S2, T128, TP128, CS, CSM, res;
430   CSM.w[1] = (CS.w[1] << 1) | (CS.w[0] >> 63);
431   CSM.w[0] = (CS.w[0] + CS.w[0]) | 1;
432   // CSM^2
433   //__mul_128x128_to_256(M256, CSM, CSM);
434   __sqr128_to_256 (M256, CSM);
450     // M256 - 8*CSM
456     // if CSM' > C256, round up