static void arb_root_arf(arb_t z, const arf_t x, ulong k, slong prec) { int inexact = arf_root(arb_midref(z), x, k, prec, ARB_RND); if (inexact) arf_mag_set_ulp(arb_radref(z), arb_midref(z), prec); else mag_zero(arb_radref(z)); }
slong hypgeom_root_bound(const mag_t z, int r) { if (r == 0) { return 0; } else { arf_t t; slong v; arf_init(t); arf_set_mag(t, z); arf_root(t, t, r, MAG_BITS, ARF_RND_UP); arf_add_ui(t, t, 1, MAG_BITS, ARF_RND_UP); v = arf_get_si(t, ARF_RND_UP); arf_clear(t); return v; } }