def mpi_sqrt(s, prec): sa, sb = s # sqrt is monotonous a = mpf_sqrt(sa, prec, round_floor) b = mpf_sqrt(sb, prec, round_ceiling) return a, b
def sqrt_fixed(x, prec): return to_fixed(mpf_sqrt(from_man_exp(x, -prec, prec), prec), prec)