File tree 1 file changed +4
-3
lines changed 1 file changed +4
-3
lines changed Original file line number Diff line number Diff line change @@ -335,14 +335,15 @@ def _integer_sqrt_of_frac_rto(n: int, m: int) -> int:
335
335
return a | (a * a * m != n )
336
336
337
337
338
- # For 53 bit precision floats, the _float_sqrt_of_frac() shift is 109.
339
- _sqrt_shift : int = 2 * sys .float_info .mant_dig + 3
338
+ # For 53 bit precision floats, the bit width used in
339
+ # _float_sqrt_of_frac() is 109.
340
+ _sqrt_bit_width : int = 2 * sys .float_info .mant_dig + 3
340
341
341
342
342
343
def _float_sqrt_of_frac (n : int , m : int ) -> float :
343
344
"""Square root of n/m as a float, correctly rounded."""
344
345
# See principle and proof sketch at: https://bugs.python.org/msg407078
345
- q = (n .bit_length () - m .bit_length () - _sqrt_shift ) // 2
346
+ q = (n .bit_length () - m .bit_length () - _sqrt_bit_width ) // 2
346
347
if q >= 0 :
347
348
numerator = _integer_sqrt_of_frac_rto (n , m << 2 * q ) << q
348
349
denominator = 1
You can’t perform that action at this time.
0 commit comments