29     uint32_t a = s ? -x : x;  in gglRecipQNormalized()  local
32     if (ggl_unlikely(a <= 2LU)) {  in gglRecipQNormalized()
40     const int32_t lz = gglClz(a);  in gglRecipQNormalized()
41     a <<= lz;  // 0.32  in gglRecipQNormalized()
42     uint32_t r = a;  in gglRecipQNormalized()
46     if (a != 0x80000000) {  in gglRecipQNormalized()
50         r = (((2LU<<30) - uint32_t((uint64_t(a)*r) >> 32)) * uint64_t(r)) >> 30;  in gglRecipQNormalized()
51         r = (((2LU<<30) - uint32_t((uint64_t(a)*r) >> 32)) * uint64_t(r)) >> 30;  in gglRecipQNormalized()
82     const GGLfixed a = x;  in gglSqrtRecipx()  local
84     x = ggl_sqrt_reciproc_approx_tab[(a>>(28-lz))&0x7];  in gglSqrtRecipx()
92     x = gglMulx((x>>1),(0x30000 - gglMulx(gglMulx(a,x),x)));  in gglSqrtRecipx()
93     x = gglMulx((x>>1),(0x30000 - gglMulx(gglMulx(a,x),x)));  in gglSqrtRecipx()
97 GGLfixed gglSqrtx(GGLfixed a)  in gglSqrtx()  argument
107         if (a >= temp) {  in gglSqrtx()
109             a -= temp;  in gglSqrtx()
237 static inline void swap(T& a, T& b) {  in swap()  argument
238     T t(a);  in swap()
239     a = b;  in swap()
244 int32_t slow_muldiv(uint32_t a, uint32_t b, uint32_t c)  in slow_muldiv()  argument
248     uint64_t ab = uint64_t(a)*b;  in slow_muldiv()
280 int32_t quick_muldiv(int32_t a, int32_t b, int32_t c)  in quick_muldiv()  argument
283     int leading = gglClz(a);  in quick_muldiv()
285     a <<= leading;  in quick_muldiv()
288         if (a < 0)  in quick_muldiv()
290         a <<= 1;  in quick_muldiv()
309 int32_t gglMulDivi(int32_t a, int32_t b, int32_t c)  in gglMulDivi()  argument
312 	int32_t sign = a^b^c;  in gglMulDivi()
314 	if (a < 0) a = -a;  in gglMulDivi()
318     if (a < b) {  in gglMulDivi()
319         swap(a, b);  in gglMulDivi()
322 	if (b <= c) result = quick_muldiv(a, b, c);  in gglMulDivi()
323 	else        result = slow_muldiv((uint32_t)a, (uint32_t)b, (uint32_t)c);  in gglMulDivi()