passagemath-flint 10.6.1rc10__cp313-cp313-manylinux_2_27_x86_64.manylinux_2_28_x86_64.whl

sage/libs/flint/acb_poly.pxd

@@ -0,0 +1,251 @@
+# sage_setup: distribution = sagemath-flint
+# distutils: libraries = flint
+# distutils: depends = flint/acb_poly.h
+
+################################################################################
+# This file is auto-generated by the script
+#   SAGE_ROOT/src/sage_setup/autogen/flint_autogen.py.
+# From the commit 3e2c3a3e091106a25ca9c6fba28e02f2cbcd654a
+# Do not modify by hand! Fix and rerun the script instead.
+################################################################################
+
+from libc.stdio cimport FILE
+from sage.libs.gmp.types cimport *
+from sage.libs.mpfr.types cimport *
+from sage.libs.flint.types cimport *
+
+cdef extern from "flint_wrap.h":
+    void acb_poly_init(acb_poly_t poly) noexcept
+    void acb_poly_clear(acb_poly_t poly) noexcept
+    void acb_poly_fit_length(acb_poly_t poly, slong len) noexcept
+    void _acb_poly_set_length(acb_poly_t poly, slong len) noexcept
+    void _acb_poly_normalise(acb_poly_t poly) noexcept
+    void acb_poly_swap(acb_poly_t poly1, acb_poly_t poly2) noexcept
+    slong acb_poly_allocated_bytes(const acb_poly_t x) noexcept
+    slong acb_poly_length(const acb_poly_t poly) noexcept
+    slong acb_poly_degree(const acb_poly_t poly) noexcept
+    bint acb_poly_is_zero(const acb_poly_t poly) noexcept
+    bint acb_poly_is_one(const acb_poly_t poly) noexcept
+    bint acb_poly_is_x(const acb_poly_t poly) noexcept
+    void acb_poly_zero(acb_poly_t poly) noexcept
+    void acb_poly_one(acb_poly_t poly) noexcept
+    void acb_poly_set(acb_poly_t dest, const acb_poly_t src) noexcept
+    void acb_poly_set_round(acb_poly_t dest, const acb_poly_t src, slong prec) noexcept
+    void acb_poly_set_trunc(acb_poly_t dest, const acb_poly_t src, slong n) noexcept
+    void acb_poly_set_trunc_round(acb_poly_t dest, const acb_poly_t src, slong n, slong prec) noexcept
+    void acb_poly_set_coeff_si(acb_poly_t poly, slong n, slong c) noexcept
+    void acb_poly_set_coeff_acb(acb_poly_t poly, slong n, const acb_t c) noexcept
+    void acb_poly_get_coeff_acb(acb_t v, const acb_poly_t poly, slong n) noexcept
+    void _acb_poly_shift_right(acb_ptr res, acb_srcptr poly, slong len, slong n) noexcept
+    void acb_poly_shift_right(acb_poly_t res, const acb_poly_t poly, slong n) noexcept
+    void _acb_poly_shift_left(acb_ptr res, acb_srcptr poly, slong len, slong n) noexcept
+    void acb_poly_shift_left(acb_poly_t res, const acb_poly_t poly, slong n) noexcept
+    void acb_poly_truncate(acb_poly_t poly, slong n) noexcept
+    slong acb_poly_valuation(const acb_poly_t poly) noexcept
+    void acb_poly_printd(const acb_poly_t poly, slong digits) noexcept
+    void acb_poly_fprintd(FILE * file, const acb_poly_t poly, slong digits) noexcept
+    void acb_poly_randtest(acb_poly_t poly, flint_rand_t state, slong len, slong prec, slong mag_bits) noexcept
+    bint acb_poly_equal(const acb_poly_t A, const acb_poly_t B) noexcept
+    bint acb_poly_contains(const acb_poly_t poly1, const acb_poly_t poly2) noexcept
+    bint acb_poly_contains_fmpz_poly(const acb_poly_t poly1, const fmpz_poly_t poly2) noexcept
+    bint acb_poly_contains_fmpq_poly(const acb_poly_t poly1, const fmpq_poly_t poly2) noexcept
+    bint _acb_poly_overlaps(acb_srcptr poly1, slong len1, acb_srcptr poly2, slong len2) noexcept
+    bint acb_poly_overlaps(const acb_poly_t poly1, const acb_poly_t poly2) noexcept
+    int acb_poly_get_unique_fmpz_poly(fmpz_poly_t z, const acb_poly_t x) noexcept
+    bint acb_poly_is_real(const acb_poly_t poly) noexcept
+    void acb_poly_set_fmpz_poly(acb_poly_t poly, const fmpz_poly_t re, slong prec) noexcept
+    void acb_poly_set2_fmpz_poly(acb_poly_t poly, const fmpz_poly_t re, const fmpz_poly_t im, slong prec) noexcept
+    void acb_poly_set_arb_poly(acb_poly_t poly, const arb_poly_t re) noexcept
+    void acb_poly_set2_arb_poly(acb_poly_t poly, const arb_poly_t re, const arb_poly_t im) noexcept
+    void acb_poly_set_fmpq_poly(acb_poly_t poly, const fmpq_poly_t re, slong prec) noexcept
+    void acb_poly_set2_fmpq_poly(acb_poly_t poly, const fmpq_poly_t re, const fmpq_poly_t im, slong prec) noexcept
+    void acb_poly_set_acb(acb_poly_t poly, const acb_t src) noexcept
+    void acb_poly_set_si(acb_poly_t poly, slong src) noexcept
+    void _acb_poly_majorant(arb_ptr res, acb_srcptr poly, slong len, slong prec) noexcept
+    void acb_poly_majorant(arb_poly_t res, const acb_poly_t poly, slong prec) noexcept
+    void _acb_poly_add(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    void acb_poly_add(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong prec) noexcept
+    void acb_poly_add_si(acb_poly_t C, const acb_poly_t A, slong B, slong prec) noexcept
+    void _acb_poly_sub(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    void acb_poly_sub(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong prec) noexcept
+    void acb_poly_add_series(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong len, slong prec) noexcept
+    void acb_poly_sub_series(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong len, slong prec) noexcept
+    void acb_poly_neg(acb_poly_t C, const acb_poly_t A) noexcept
+    void acb_poly_scalar_mul_2exp_si(acb_poly_t C, const acb_poly_t A, slong c) noexcept
+    void acb_poly_scalar_mul(acb_poly_t C, const acb_poly_t A, const acb_t c, slong prec) noexcept
+    void acb_poly_scalar_div(acb_poly_t C, const acb_poly_t A, const acb_t c, slong prec) noexcept
+    void _acb_poly_mullow_classical(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong n, slong prec) noexcept
+    void _acb_poly_mullow_transpose(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong n, slong prec) noexcept
+    void _acb_poly_mullow_transpose_gauss(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong n, slong prec) noexcept
+    void _acb_poly_mullow(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong n, slong prec) noexcept
+    void acb_poly_mullow_classical(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong n, slong prec) noexcept
+    void acb_poly_mullow_transpose(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong n, slong prec) noexcept
+    void acb_poly_mullow_transpose_gauss(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong n, slong prec) noexcept
+    void acb_poly_mullow(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong n, slong prec) noexcept
+    void _acb_poly_mul(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    void acb_poly_mul(acb_poly_t C, const acb_poly_t A1, const acb_poly_t B2, slong prec) noexcept
+    void _acb_poly_inv_series(acb_ptr Qinv, acb_srcptr Q, slong Qlen, slong len, slong prec) noexcept
+    void acb_poly_inv_series(acb_poly_t Qinv, const acb_poly_t Q, slong n, slong prec) noexcept
+    void _acb_poly_div_series(acb_ptr Q, acb_srcptr A, slong Alen, acb_srcptr B, slong Blen, slong n, slong prec) noexcept
+    void acb_poly_div_series(acb_poly_t Q, const acb_poly_t A, const acb_poly_t B, slong n, slong prec) noexcept
+    void _acb_poly_div(acb_ptr Q, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    void _acb_poly_rem(acb_ptr R, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    void _acb_poly_divrem(acb_ptr Q, acb_ptr R, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec) noexcept
+    int acb_poly_divrem(acb_poly_t Q, acb_poly_t R, const acb_poly_t A, const acb_poly_t B, slong prec) noexcept
+    void _acb_poly_div_root(acb_ptr Q, acb_t R, acb_srcptr A, slong len, const acb_t c, slong prec) noexcept
+    void _acb_poly_taylor_shift(acb_ptr g, const acb_t c, slong n, slong prec) noexcept
+    void acb_poly_taylor_shift(acb_poly_t g, const acb_poly_t f, const acb_t c, slong prec) noexcept
+    void _acb_poly_compose(acb_ptr res, acb_srcptr poly1, slong len1, acb_srcptr poly2, slong len2, slong prec) noexcept
+    void acb_poly_compose(acb_poly_t res, const acb_poly_t poly1, const acb_poly_t poly2, slong prec) noexcept
+    void _acb_poly_compose_series(acb_ptr res, acb_srcptr poly1, slong len1, acb_srcptr poly2, slong len2, slong n, slong prec) noexcept
+    void acb_poly_compose_series(acb_poly_t res, const acb_poly_t poly1, const acb_poly_t poly2, slong n, slong prec) noexcept
+    void _acb_poly_revert_series(acb_ptr h, acb_srcptr f, slong flen, slong n, slong prec) noexcept
+    void acb_poly_revert_series(acb_poly_t h, const acb_poly_t f, slong n, slong prec) noexcept
+    void _acb_poly_evaluate_horner(acb_t y, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate_horner(acb_t y, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_evaluate_rectangular(acb_t y, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate_rectangular(acb_t y, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_evaluate(acb_t y, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate(acb_t y, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_evaluate2_horner(acb_t y, acb_t z, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate2_horner(acb_t y, acb_t z, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_evaluate2_rectangular(acb_t y, acb_t z, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate2_rectangular(acb_t y, acb_t z, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_evaluate2(acb_t y, acb_t z, acb_srcptr f, slong len, const acb_t x, slong prec) noexcept
+    void acb_poly_evaluate2(acb_t y, acb_t z, const acb_poly_t f, const acb_t x, slong prec) noexcept
+    void _acb_poly_product_roots(acb_ptr poly, acb_srcptr xs, slong n, slong prec) noexcept
+    void acb_poly_product_roots(acb_poly_t poly, acb_srcptr xs, slong n, slong prec) noexcept
+    acb_ptr * _acb_poly_tree_alloc(slong len) noexcept
+    void _acb_poly_tree_free(acb_ptr * tree, slong len) noexcept
+    void _acb_poly_tree_build(acb_ptr * tree, acb_srcptr roots, slong len, slong prec) noexcept
+    void _acb_poly_evaluate_vec_iter(acb_ptr ys, acb_srcptr poly, slong plen, acb_srcptr xs, slong n, slong prec) noexcept
+    void acb_poly_evaluate_vec_iter(acb_ptr ys, const acb_poly_t poly, acb_srcptr xs, slong n, slong prec) noexcept
+    void _acb_poly_evaluate_vec_fast_precomp(acb_ptr vs, acb_srcptr poly, slong plen, acb_ptr * tree, slong len, slong prec) noexcept
+    void _acb_poly_evaluate_vec_fast(acb_ptr ys, acb_srcptr poly, slong plen, acb_srcptr xs, slong n, slong prec) noexcept
+    void acb_poly_evaluate_vec_fast(acb_ptr ys, const acb_poly_t poly, acb_srcptr xs, slong n, slong prec) noexcept
+    void _acb_poly_interpolate_newton(acb_ptr poly, acb_srcptr xs, acb_srcptr ys, slong n, slong prec) noexcept
+    void acb_poly_interpolate_newton(acb_poly_t poly, acb_srcptr xs, acb_srcptr ys, slong n, slong prec) noexcept
+    void _acb_poly_interpolate_barycentric(acb_ptr poly, acb_srcptr xs, acb_srcptr ys, slong n, slong prec) noexcept
+    void acb_poly_interpolate_barycentric(acb_poly_t poly, acb_srcptr xs, acb_srcptr ys, slong n, slong prec) noexcept
+    void _acb_poly_interpolation_weights(acb_ptr w, acb_ptr * tree, slong len, slong prec) noexcept
+    void _acb_poly_interpolate_fast_precomp(acb_ptr poly, acb_srcptr ys, acb_ptr * tree, acb_srcptr weights, slong len, slong prec) noexcept
+    void _acb_poly_interpolate_fast(acb_ptr poly, acb_srcptr xs, acb_srcptr ys, slong len, slong prec) noexcept
+    void acb_poly_interpolate_fast(acb_poly_t poly, acb_srcptr xs, acb_srcptr ys, slong n, slong prec) noexcept
+    void _acb_poly_derivative(acb_ptr res, acb_srcptr poly, slong len, slong prec) noexcept
+    void acb_poly_derivative(acb_poly_t res, const acb_poly_t poly, slong prec) noexcept
+    void _acb_poly_nth_derivative(acb_ptr res, acb_srcptr poly, ulong n, slong len, slong prec) noexcept
+    void acb_poly_nth_derivative(acb_poly_t res, const acb_poly_t poly, ulong n, slong prec) noexcept
+    void _acb_poly_integral(acb_ptr res, acb_srcptr poly, slong len, slong prec) noexcept
+    void acb_poly_integral(acb_poly_t res, const acb_poly_t poly, slong prec) noexcept
+    void _acb_poly_borel_transform(acb_ptr res, acb_srcptr poly, slong len, slong prec) noexcept
+    void acb_poly_borel_transform(acb_poly_t res, const acb_poly_t poly, slong prec) noexcept
+    void _acb_poly_inv_borel_transform(acb_ptr res, acb_srcptr poly, slong len, slong prec) noexcept
+    void acb_poly_inv_borel_transform(acb_poly_t res, const acb_poly_t poly, slong prec) noexcept
+    void _acb_poly_binomial_transform_basecase(acb_ptr b, acb_srcptr a, slong alen, slong len, slong prec) noexcept
+    void acb_poly_binomial_transform_basecase(acb_poly_t b, const acb_poly_t a, slong len, slong prec) noexcept
+    void _acb_poly_binomial_transform_convolution(acb_ptr b, acb_srcptr a, slong alen, slong len, slong prec) noexcept
+    void acb_poly_binomial_transform_convolution(acb_poly_t b, const acb_poly_t a, slong len, slong prec) noexcept
+    void _acb_poly_binomial_transform(acb_ptr b, acb_srcptr a, slong alen, slong len, slong prec) noexcept
+    void acb_poly_binomial_transform(acb_poly_t b, const acb_poly_t a, slong len, slong prec) noexcept
+    void _acb_poly_graeffe_transform(acb_ptr b, acb_srcptr a, slong len, slong prec) noexcept
+    void acb_poly_graeffe_transform(acb_poly_t b, const acb_poly_t a, slong prec) noexcept
+    void _acb_poly_pow_ui_trunc_binexp(acb_ptr res, acb_srcptr f, slong flen, ulong exp, slong len, slong prec) noexcept
+    void acb_poly_pow_ui_trunc_binexp(acb_poly_t res, const acb_poly_t poly, ulong exp, slong len, slong prec) noexcept
+    void _acb_poly_pow_ui(acb_ptr res, acb_srcptr f, slong flen, ulong exp, slong prec) noexcept
+    void acb_poly_pow_ui(acb_poly_t res, const acb_poly_t poly, ulong exp, slong prec) noexcept
+    void _acb_poly_pow_series(acb_ptr h, acb_srcptr f, slong flen, acb_srcptr g, slong glen, slong len, slong prec) noexcept
+    void acb_poly_pow_series(acb_poly_t h, const acb_poly_t f, const acb_poly_t g, slong len, slong prec) noexcept
+    void _acb_poly_pow_acb_series(acb_ptr h, acb_srcptr f, slong flen, const acb_t g, slong len, slong prec) noexcept
+    void acb_poly_pow_acb_series(acb_poly_t h, const acb_poly_t f, const acb_t g, slong len, slong prec) noexcept
+    void _acb_poly_sqrt_series(acb_ptr g, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sqrt_series(acb_poly_t g, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_rsqrt_series(acb_ptr g, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_rsqrt_series(acb_poly_t g, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_log_series(acb_ptr res, acb_srcptr f, slong flen, slong n, slong prec) noexcept
+    void acb_poly_log_series(acb_poly_t res, const acb_poly_t f, slong n, slong prec) noexcept
+    void _acb_poly_log1p_series(acb_ptr res, acb_srcptr f, slong flen, slong n, slong prec) noexcept
+    void acb_poly_log1p_series(acb_poly_t res, const acb_poly_t f, slong n, slong prec) noexcept
+    void _acb_poly_atan_series(acb_ptr res, acb_srcptr f, slong flen, slong n, slong prec) noexcept
+    void acb_poly_atan_series(acb_poly_t res, const acb_poly_t f, slong n, slong prec) noexcept
+    void _acb_poly_exp_series_basecase(acb_ptr f, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_exp_series_basecase(acb_poly_t f, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_exp_series(acb_ptr f, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_exp_series(acb_poly_t f, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_exp_pi_i_series(acb_ptr f, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_exp_pi_i_series(acb_poly_t f, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sin_cos_series(acb_ptr s, acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sin_cos_series(acb_poly_t s, acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sin_series(acb_ptr s, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sin_series(acb_poly_t s, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_cos_series(acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_cos_series(acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_tan_series(acb_ptr g, acb_srcptr h, slong hlen, slong len, slong prec) noexcept
+    void acb_poly_tan_series(acb_poly_t g, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sin_cos_pi_series(acb_ptr s, acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sin_cos_pi_series(acb_poly_t s, acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sin_pi_series(acb_ptr s, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sin_pi_series(acb_poly_t s, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_cos_pi_series(acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_cos_pi_series(acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_cot_pi_series(acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_cot_pi_series(acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sinh_cosh_series_basecase(acb_ptr s, acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sinh_cosh_series_basecase(acb_poly_t s, acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sinh_cosh_series_exponential(acb_ptr s, acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sinh_cosh_series_exponential(acb_poly_t s, acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sinh_cosh_series(acb_ptr s, acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sinh_cosh_series(acb_poly_t s, acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sinh_series(acb_ptr s, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sinh_series(acb_poly_t s, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_cosh_series(acb_ptr c, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_cosh_series(acb_poly_t c, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_sinc_series(acb_ptr s, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_sinc_series(acb_poly_t s, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_lambertw_series(acb_ptr res, acb_srcptr z, slong zlen, const fmpz_t k, int flags, slong len, slong prec) noexcept
+    void acb_poly_lambertw_series(acb_poly_t res, const acb_poly_t z, const fmpz_t k, int flags, slong len, slong prec) noexcept
+    void _acb_poly_gamma_series(acb_ptr res, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_gamma_series(acb_poly_t res, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_rgamma_series(acb_ptr res, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_rgamma_series(acb_poly_t res, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_lgamma_series(acb_ptr res, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_lgamma_series(acb_poly_t res, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_digamma_series(acb_ptr res, acb_srcptr h, slong hlen, slong n, slong prec) noexcept
+    void acb_poly_digamma_series(acb_poly_t res, const acb_poly_t h, slong n, slong prec) noexcept
+    void _acb_poly_rising_ui_series(acb_ptr res, acb_srcptr f, slong flen, ulong r, slong trunc, slong prec) noexcept
+    void acb_poly_rising_ui_series(acb_poly_t res, const acb_poly_t f, ulong r, slong trunc, slong prec) noexcept
+    void _acb_poly_powsum_series_naive(acb_ptr z, const acb_t s, const acb_t a, const acb_t q, slong n, slong len, slong prec) noexcept
+    void _acb_poly_powsum_series_naive_threaded(acb_ptr z, const acb_t s, const acb_t a, const acb_t q, slong n, slong len, slong prec) noexcept
+    void _acb_poly_powsum_one_series_sieved(acb_ptr z, const acb_t s, slong n, slong len, slong prec) noexcept
+    void _acb_poly_zeta_em_choose_param(mag_t bound, ulong * N, ulong * M, const acb_t s, const acb_t a, slong d, slong target, slong prec) noexcept
+    void _acb_poly_zeta_em_bound1(mag_t bound, const acb_t s, const acb_t a, slong N, slong M, slong d, slong wp) noexcept
+    void _acb_poly_zeta_em_bound(arb_ptr vec, const acb_t s, const acb_t a, ulong N, ulong M, slong d, slong wp) noexcept
+    void _acb_poly_zeta_em_tail_naive(acb_ptr z, const acb_t s, const acb_t Na, acb_srcptr Nasx, slong M, slong len, slong prec) noexcept
+    void _acb_poly_zeta_em_tail_bsplit(acb_ptr z, const acb_t s, const acb_t Na, acb_srcptr Nasx, slong M, slong len, slong prec) noexcept
+    void _acb_poly_zeta_em_sum(acb_ptr z, const acb_t s, const acb_t a, int deflate, ulong N, ulong M, slong d, slong prec) noexcept
+    void _acb_poly_zeta_cpx_series(acb_ptr z, const acb_t s, const acb_t a, int deflate, slong d, slong prec) noexcept
+    void _acb_poly_zeta_series(acb_ptr res, acb_srcptr h, slong hlen, const acb_t a, int deflate, slong len, slong prec) noexcept
+    void acb_poly_zeta_series(acb_poly_t res, const acb_poly_t f, const acb_t a, int deflate, slong n, slong prec) noexcept
+    void _acb_poly_polylog_cpx_small(acb_ptr w, const acb_t s, const acb_t z, slong len, slong prec) noexcept
+    void _acb_poly_polylog_cpx_zeta(acb_ptr w, const acb_t s, const acb_t z, slong len, slong prec) noexcept
+    void _acb_poly_polylog_cpx(acb_ptr w, const acb_t s, const acb_t z, slong len, slong prec) noexcept
+    void _acb_poly_polylog_series(acb_ptr w, acb_srcptr s, slong slen, const acb_t z, slong len, slong prec) noexcept
+    void acb_poly_polylog_series(acb_poly_t w, const acb_poly_t s, const acb_t z, slong len, slong prec) noexcept
+    void _acb_poly_erf_series(acb_ptr res, acb_srcptr z, slong zlen, slong n, slong prec) noexcept
+    void acb_poly_erf_series(acb_poly_t res, const acb_poly_t z, slong n, slong prec) noexcept
+    void _acb_poly_agm1_series(acb_ptr res, acb_srcptr z, slong zlen, slong len, slong prec) noexcept
+    void acb_poly_agm1_series(acb_poly_t res, const acb_poly_t z, slong n, slong prec) noexcept
+    void _acb_poly_elliptic_k_series(acb_ptr res, acb_srcptr z, slong zlen, slong len, slong prec) noexcept
+    void acb_poly_elliptic_k_series(acb_poly_t res, const acb_poly_t z, slong n, slong prec) noexcept
+    void _acb_poly_elliptic_p_series(acb_ptr res, acb_srcptr z, slong zlen, const acb_t tau, slong len, slong prec) noexcept
+    void acb_poly_elliptic_p_series(acb_poly_t res, const acb_poly_t z, const acb_t tau, slong n, slong prec) noexcept
+    void _acb_poly_root_bound_fujiwara(mag_t bound, acb_srcptr poly, slong len) noexcept
+    void acb_poly_root_bound_fujiwara(mag_t bound, acb_poly_t poly) noexcept
+    void _acb_poly_root_inclusion(acb_t r, const acb_t m, acb_srcptr poly, acb_srcptr polyder, slong len, slong prec) noexcept
+    slong _acb_poly_validate_roots(acb_ptr roots, acb_srcptr poly, slong len, slong prec) noexcept
+    void _acb_poly_refine_roots_durand_kerner(acb_ptr roots, acb_srcptr poly, slong len, slong prec) noexcept
+    slong _acb_poly_find_roots(acb_ptr roots, acb_srcptr poly, acb_srcptr initial, slong len, slong maxiter, slong prec) noexcept
+    slong acb_poly_find_roots(acb_ptr roots, const acb_poly_t poly, acb_srcptr initial, slong maxiter, slong prec) noexcept
+    int _acb_poly_validate_real_roots(acb_srcptr roots, acb_srcptr poly, slong len, slong prec) noexcept
+    int acb_poly_validate_real_roots(acb_srcptr roots, const acb_poly_t poly, slong prec) noexcept
+
+from .acb_poly_macros cimport *

sage/libs/flint/acb_poly_macros.pxd

@@ -0,0 +1,8 @@
+# sage_setup: distribution = sagemath-flint
+# Macros from acb_poly.h
+# See https://github.com/flintlib/flint/issues/1529
+
+from .types cimport *
+
+cdef extern from "flint_wrap.h":
+    acb_ptr acb_poly_get_coeff_ptr(acb_poly_t p, ulong n)

sage/libs/flint/acb_theta.pxd

@@ -0,0 +1,124 @@
+# sage_setup: distribution = sagemath-flint
+# distutils: libraries = flint
+# distutils: depends = flint/acb_theta.h
+
+################################################################################
+# This file is auto-generated by the script
+#   SAGE_ROOT/src/sage_setup/autogen/flint_autogen.py.
+# Do not modify by hand! Fix and rerun the script instead.
+################################################################################
+
+from libc.stdio cimport FILE
+from sage.libs.gmp.types cimport *
+from sage.libs.mpfr.types cimport *
+from sage.libs.flint.types cimport *
+
+cdef extern from "flint_wrap.h":
+    void acb_theta_all(acb_ptr th, acb_srcptr z, const acb_mat_t tau, int sqr, slong prec) noexcept
+    void acb_theta_naive_fixed_ab(acb_ptr th, ulong ab, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_naive_all(acb_ptr th, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_jet_all(acb_ptr dth, acb_srcptr z, const acb_mat_t tau, slong ord, slong prec) noexcept
+    void acb_theta_jet_naive_fixed_ab(acb_ptr dth, ulong ab, acb_srcptr z, const acb_mat_t tau, slong ord, slong prec) noexcept
+    void acb_theta_jet_naive_all(acb_ptr dth, acb_srcptr z, const acb_mat_t tau, slong ord, slong prec) noexcept
+    slong sp2gz_dim(const fmpz_mat_t mat) noexcept
+    void sp2gz_set_blocks(fmpz_mat_t mat, const fmpz_mat_t alpha, const fmpz_mat_t beta, const fmpz_mat_t gamma, const fmpz_mat_t delta) noexcept
+    void sp2gz_j(fmpz_mat_t mat) noexcept
+    void sp2gz_block_diag(fmpz_mat_t mat, const fmpz_mat_t U) noexcept
+    void sp2gz_trig(fmpz_mat_t mat, const fmpz_mat_t S) noexcept
+    void sp2gz_embed(fmpz_mat_t res, const fmpz_mat_t mat) noexcept
+    void sp2gz_restrict(fmpz_mat_t res, const fmpz_mat_t mat) noexcept
+    slong sp2gz_nb_fundamental(slong g) noexcept
+    void sp2gz_fundamental(fmpz_mat_t mat, slong j) noexcept
+    bint sp2gz_is_correct(const fmpz_mat_t mat) noexcept
+    bint sp2gz_is_j(const fmpz_mat_t mat) noexcept
+    bint sp2gz_is_block_diag(const fmpz_mat_t mat) noexcept
+    bint sp2gz_is_trig(const fmpz_mat_t mat) noexcept
+    bint sp2gz_is_embedded(fmpz_mat_t res, const fmpz_mat_t mat) noexcept
+    void sp2gz_inv(fmpz_mat_t inv, const fmpz_mat_t mat) noexcept
+    fmpz_mat_struct * sp2gz_decompose(slong * nb, const fmpz_mat_t mat) noexcept
+    void sp2gz_randtest(fmpz_mat_t mat, flint_rand_t state, slong bits) noexcept
+    void acb_siegel_cocycle(acb_mat_t c, const fmpz_mat_t mat, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_transform_cocycle_inv(acb_mat_t w, acb_mat_t c, acb_mat_t cinv, const fmpz_mat_t mat, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_transform(acb_mat_t w, const fmpz_mat_t mat, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_transform_z(acb_ptr r, acb_mat_t w, const fmpz_mat_t mat, acb_srcptr z, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_cho(arb_mat_t C, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_yinv(arb_mat_t Yinv, const acb_mat_t tau, slong prec) noexcept
+    void acb_siegel_reduce(fmpz_mat_t mat, const acb_mat_t tau, slong prec) noexcept
+    bint acb_siegel_is_reduced(const acb_mat_t tau, slong tol_exp, slong prec) noexcept
+    void acb_siegel_randtest(acb_mat_t tau, flint_rand_t state, slong prec, slong mag_bits) noexcept
+    void acb_siegel_randtest_reduced(acb_mat_t tau, flint_rand_t state, slong prec, slong mag_bits) noexcept
+    void acb_siegel_randtest_vec(acb_ptr z, flint_rand_t state, slong g, slong prec) noexcept
+    void acb_theta_char_get_slong(slong * n, ulong a, slong g) noexcept
+    ulong acb_theta_char_get_a(const slong * n, slong g) noexcept
+    void acb_theta_char_get_arb(arb_ptr v, ulong a, slong g) noexcept
+    void acb_theta_char_get_acb(acb_ptr v, ulong a, slong g) noexcept
+    slong acb_theta_char_dot(ulong a, ulong b, slong g) noexcept
+    slong acb_theta_char_dot_slong(ulong a, const slong * n, slong g) noexcept
+    void acb_theta_char_dot_acb(acb_t x, ulong a, acb_srcptr z, slong g, slong prec) noexcept
+    bint acb_theta_char_is_even(ulong ab, slong g) noexcept
+    bint acb_theta_char_is_goepel(ulong ch1, ulong ch2, ulong ch3, ulong ch4, slong g) noexcept
+    bint acb_theta_char_is_syzygous(ulong ch1, ulong ch2, ulong ch3, slong g) noexcept
+    void acb_theta_eld_init(acb_theta_eld_t E, slong d, slong g) noexcept
+    void acb_theta_eld_clear(acb_theta_eld_t E) noexcept
+    int acb_theta_eld_set(acb_theta_eld_t E, const arb_mat_t C, const arf_t R2, arb_srcptr v) noexcept
+    void acb_theta_eld_points(slong * pts, const acb_theta_eld_t E) noexcept
+    void acb_theta_eld_border(slong * pts, const acb_theta_eld_t E) noexcept
+    bint acb_theta_eld_contains(const acb_theta_eld_t E, slong * pt) noexcept
+    void acb_theta_eld_print(const acb_theta_eld_t E) noexcept
+    void acb_theta_naive_radius(arf_t R2, arf_t eps, const arb_mat_t C, slong ord, slong prec) noexcept
+    void acb_theta_naive_reduce(arb_ptr v, acb_ptr new_zs, arb_ptr as, acb_ptr cs, arb_ptr us, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_naive_term(acb_t res, acb_srcptr z, const acb_mat_t tau, slong * tup, slong * n, slong prec) noexcept
+    void acb_theta_naive_worker(acb_ptr th, slong len, acb_srcptr zs, slong nb, const acb_mat_t tau, const acb_theta_eld_t E, slong ord, slong prec, acb_theta_naive_worker_t worker) noexcept
+    void acb_theta_naive_00(acb_ptr th, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_naive_0b(acb_ptr th, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_naive_fixed_a(acb_ptr th, ulong a, acb_srcptr zs, slong nb, const acb_mat_t tau, slong prec) noexcept
+    slong acb_theta_jet_nb(slong ord, slong g) noexcept
+    slong acb_theta_jet_total_order(const slong * tup, slong g) noexcept
+    void acb_theta_jet_tuples(slong * tups, slong ord, slong g) noexcept
+    slong acb_theta_jet_index(const slong * tup, slong g) noexcept
+    void acb_theta_jet_mul(acb_ptr res, acb_srcptr v1, acb_srcptr v2, slong ord, slong g, slong prec) noexcept
+    void acb_theta_jet_compose(acb_ptr res, acb_srcptr v, const acb_mat_t N, slong ord, slong prec) noexcept
+    void acb_theta_jet_exp_pi_i(acb_ptr res, arb_srcptr a, slong ord, slong g, slong prec) noexcept
+    void acb_theta_jet_naive_radius(arf_t R2, arf_t eps, arb_srcptr v, const arb_mat_t C, slong ord, slong prec) noexcept
+    void acb_theta_jet_naive_00(acb_ptr dth, acb_srcptr z, const acb_mat_t tau, slong ord, slong prec) noexcept
+    void acb_theta_jet_error_bounds(arb_ptr err, acb_srcptr z, const acb_mat_t tau, acb_srcptr dth, slong ord, slong prec) noexcept
+    void acb_theta_dist_pt(arb_t d, arb_srcptr v, const arb_mat_t C, slong * n, slong prec) noexcept
+    void acb_theta_dist_lat(arb_t d, arb_srcptr v, const arb_mat_t C, slong prec) noexcept
+    void acb_theta_dist_a0(arb_ptr d, acb_srcptr z, const acb_mat_t tau, slong prec) noexcept
+    slong acb_theta_dist_addprec(const arb_t d) noexcept
+    void acb_theta_agm_hadamard(acb_ptr res, acb_srcptr a, slong g, slong prec) noexcept
+    void acb_theta_agm_sqrt(acb_ptr res, acb_srcptr a, acb_srcptr rts, slong nb, slong prec) noexcept
+    void acb_theta_agm_mul(acb_ptr res, acb_srcptr a1, acb_srcptr a2, slong g, slong prec) noexcept
+    void acb_theta_agm_mul_tight(acb_ptr res, acb_srcptr a0, acb_srcptr a, arb_srcptr d0, arb_srcptr d, slong g, slong prec) noexcept
+    int acb_theta_ql_a0_naive(acb_ptr th, acb_srcptr t, acb_srcptr z, arb_srcptr d0, arb_srcptr d, const acb_mat_t tau, slong guard, slong prec) noexcept
+    int acb_theta_ql_a0_split(acb_ptr th, acb_srcptr t, acb_srcptr z, arb_srcptr d, const acb_mat_t tau, slong s, slong guard, slong prec, acb_theta_ql_worker_t worker) noexcept
+    int acb_theta_ql_a0_steps(acb_ptr th, acb_srcptr t, acb_srcptr z, arb_srcptr d0, arb_srcptr d, const acb_mat_t tau, slong nb_steps, slong s, slong guard, slong prec, acb_theta_ql_worker_t worker) noexcept
+    slong acb_theta_ql_a0_nb_steps(const arb_mat_t C, slong s, slong prec) noexcept
+    int acb_theta_ql_a0(acb_ptr th, acb_srcptr t, acb_srcptr z, arb_srcptr d0, arb_srcptr d, const acb_mat_t tau, slong guard, slong prec) noexcept
+    slong acb_theta_ql_reduce(acb_ptr new_z, acb_t c, arb_t u, slong * n1, acb_srcptr z, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_ql_all(acb_ptr th, acb_srcptr z, const acb_mat_t tau, int sqr, slong prec) noexcept
+    void acb_theta_jet_ql_bounds(arb_t c, arb_t rho, acb_srcptr z, const acb_mat_t tau, slong ord) noexcept
+    void acb_theta_jet_ql_radius(arf_t eps, arf_t err, const arb_t c, const arb_t rho, slong ord, slong g, slong prec) noexcept
+    void acb_theta_jet_ql_finite_diff(acb_ptr dth, const arf_t eps, const arf_t err, acb_srcptr val, slong ord, slong g, slong prec) noexcept
+    void acb_theta_jet_ql_all(acb_ptr dth, acb_srcptr z, const acb_mat_t tau, slong ord, slong prec) noexcept
+    ulong acb_theta_transform_char(slong * e, const fmpz_mat_t mat, ulong ab) noexcept
+    void acb_theta_transform_sqrtdet(acb_t res, const acb_mat_t tau, slong prec) noexcept
+    slong acb_theta_transform_kappa(acb_t sqrtdet, const fmpz_mat_t mat, const acb_mat_t tau, slong prec) noexcept
+    slong acb_theta_transform_kappa2(const fmpz_mat_t mat) noexcept
+    void acb_theta_transform_proj(acb_ptr res, const fmpz_mat_t mat, acb_srcptr th, int sqr, slong prec) noexcept
+    void acb_theta_g2_jet_naive_1(acb_ptr dth, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_g2_detk_symj(acb_poly_t res, const acb_mat_t m, const acb_poly_t f, slong k, slong j, slong prec) noexcept
+    void acb_theta_g2_transvectant(acb_poly_t res, const acb_poly_t g, const acb_poly_t h, slong m, slong n, slong k, slong prec) noexcept
+    void acb_theta_g2_transvectant_lead(acb_t res, const acb_poly_t g, const acb_poly_t h, slong m, slong n, slong k, slong prec) noexcept
+    slong acb_theta_g2_character(const fmpz_mat_t mat) noexcept
+    void acb_theta_g2_psi4(acb_t res, acb_srcptr th2, slong prec) noexcept
+    void acb_theta_g2_psi6(acb_t res, acb_srcptr th2, slong prec) noexcept
+    void acb_theta_g2_chi10(acb_t res, acb_srcptr th2, slong prec) noexcept
+    void acb_theta_g2_chi12(acb_t res, acb_srcptr th2, slong prec) noexcept
+    void acb_theta_g2_chi5(acb_t res, acb_srcptr th, slong prec) noexcept
+    void acb_theta_g2_chi35(acb_t res, acb_srcptr th, slong prec) noexcept
+    void acb_theta_g2_chi3_6(acb_poly_t res, acb_srcptr dth, slong prec) noexcept
+    void acb_theta_g2_sextic(acb_poly_t res, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_g2_sextic_chi5(acb_poly_t res, acb_t chi5, const acb_mat_t tau, slong prec) noexcept
+    void acb_theta_g2_covariants(acb_poly_struct * res, const acb_poly_t f, slong prec) noexcept
+    void acb_theta_g2_covariants_lead(acb_ptr res, const acb_poly_t f, slong prec) noexcept

sage/libs/flint/acf.pxd

@@ -0,0 +1,32 @@
+# sage_setup: distribution = sagemath-flint
+# distutils: libraries = flint
+# distutils: depends = flint/acf.h
+
+################################################################################
+# This file is auto-generated by the script
+#   SAGE_ROOT/src/sage_setup/autogen/flint_autogen.py.
+# From the commit 3e2c3a3e091106a25ca9c6fba28e02f2cbcd654a
+# Do not modify by hand! Fix and rerun the script instead.
+################################################################################
+
+from libc.stdio cimport FILE
+from sage.libs.gmp.types cimport *
+from sage.libs.mpfr.types cimport *
+from sage.libs.flint.types cimport *
+
+cdef extern from "flint_wrap.h":
+    void acf_init(acf_t x) noexcept
+    void acf_clear(acf_t x) noexcept
+    void acf_swap(acf_t z, acf_t x) noexcept
+    slong acf_allocated_bytes(const acf_t x) noexcept
+    arf_ptr acf_real_ptr(acf_t z) noexcept
+    arf_ptr acf_imag_ptr(acf_t z) noexcept
+    void acf_set(acf_t z, const acf_t x) noexcept
+    bint acf_equal(const acf_t x, const acf_t y) noexcept
+    int acf_add(acf_t res, const acf_t x, const acf_t y, slong prec, arf_rnd_t rnd) noexcept
+    int acf_sub(acf_t res, const acf_t x, const acf_t y, slong prec, arf_rnd_t rnd) noexcept
+    int acf_mul(acf_t res, const acf_t x, const acf_t y, slong prec, arf_rnd_t rnd) noexcept
+    void acf_approx_inv(acf_t res, const acf_t x, slong prec, arf_rnd_t rnd) noexcept
+    void acf_approx_div(acf_t res, const acf_t x, const acf_t y, slong prec, arf_rnd_t rnd) noexcept
+    void acf_approx_sqrt(acf_t res, const acf_t x, slong prec, arf_rnd_t rnd) noexcept
+    void acf_approx_dot(acf_t res, const acf_t initial, int subtract, acf_srcptr x, slong xstep, acf_srcptr y, slong ystep, slong len, slong prec, arf_rnd_t rnd) noexcept

sage/libs/flint/aprcl.pxd

@@ -0,0 +1,84 @@
+# sage_setup: distribution = sagemath-flint
+# distutils: libraries = flint
+# distutils: depends = flint/aprcl.h
+
+################################################################################
+# This file is auto-generated by the script
+#   SAGE_ROOT/src/sage_setup/autogen/flint_autogen.py.
+# From the commit 3e2c3a3e091106a25ca9c6fba28e02f2cbcd654a
+# Do not modify by hand! Fix and rerun the script instead.
+################################################################################
+
+from libc.stdio cimport FILE
+from sage.libs.gmp.types cimport *
+from sage.libs.mpfr.types cimport *
+from sage.libs.flint.types cimport *
+
+cdef extern from "flint_wrap.h":
+    bint aprcl_is_prime(const fmpz_t n) noexcept
+    bint aprcl_is_prime_jacobi(const fmpz_t n) noexcept
+    bint aprcl_is_prime_gauss(const fmpz_t n) noexcept
+    primality_test_status _aprcl_is_prime_jacobi(const fmpz_t n, const aprcl_config config) noexcept
+    primality_test_status _aprcl_is_prime_gauss(const fmpz_t n, const aprcl_config config) noexcept
+    bint aprcl_is_prime_gauss_min_R(const fmpz_t n, ulong R) noexcept
+    bint aprcl_is_prime_final_division(const fmpz_t n, const fmpz_t s, ulong r) noexcept
+    void aprcl_config_gauss_init(aprcl_config conf, const fmpz_t n) noexcept
+    void aprcl_config_gauss_init_min_R(aprcl_config conf, const fmpz_t n, ulong R) noexcept
+    void aprcl_config_gauss_clear(aprcl_config conf) noexcept
+    ulong aprcl_R_value(const fmpz_t n) noexcept
+    void aprcl_config_jacobi_init(aprcl_config conf, const fmpz_t n) noexcept
+    void aprcl_config_jacobi_clear(aprcl_config conf) noexcept
+    void unity_zp_init(unity_zp f, ulong p, ulong exp, const fmpz_t n) noexcept
+    void unity_zp_clear(unity_zp f) noexcept
+    void unity_zp_copy(unity_zp f, const unity_zp g) noexcept
+    void unity_zp_swap(unity_zp f, unity_zp q) noexcept
+    void unity_zp_set_zero(unity_zp f) noexcept
+    slong unity_zp_is_unity(unity_zp f) noexcept
+    bint unity_zp_equal(unity_zp f, unity_zp g) noexcept
+    void unity_zp_print(const unity_zp f) noexcept
+    void unity_zp_coeff_set_fmpz(unity_zp f, ulong ind, const fmpz_t x) noexcept
+    void unity_zp_coeff_set_ui(unity_zp f, ulong ind, ulong x) noexcept
+    void unity_zp_coeff_add_fmpz(unity_zp f, ulong ind, const fmpz_t x) noexcept
+    void unity_zp_coeff_add_ui(unity_zp f, ulong ind, ulong x) noexcept
+    void unity_zp_coeff_inc(unity_zp f, ulong ind) noexcept
+    void unity_zp_coeff_dec(unity_zp f, ulong ind) noexcept
+    void unity_zp_mul_scalar_fmpz(unity_zp f, const unity_zp g, const fmpz_t s) noexcept
+    void unity_zp_mul_scalar_ui(unity_zp f, const unity_zp g, ulong s) noexcept
+    void unity_zp_add(unity_zp f, const unity_zp g, const unity_zp h) noexcept
+    void unity_zp_mul(unity_zp f, const unity_zp g, const unity_zp h) noexcept
+    void unity_zp_sqr(unity_zp f, const unity_zp g) noexcept
+    void unity_zp_mul_inplace(unity_zp f, const unity_zp g, const unity_zp h, fmpz_t * t) noexcept
+    void unity_zp_sqr_inplace(unity_zp f, const unity_zp g, fmpz_t * t) noexcept
+    void unity_zp_pow_fmpz(unity_zp f, const unity_zp g, const fmpz_t pow) noexcept
+    void unity_zp_pow_ui(unity_zp f, const unity_zp g, ulong pow) noexcept
+    ulong _unity_zp_pow_select_k(const fmpz_t n) noexcept
+    void unity_zp_pow_2k_fmpz(unity_zp f, const unity_zp g, const fmpz_t pow) noexcept
+    void unity_zp_pow_2k_ui(unity_zp f, const unity_zp g, ulong pow) noexcept
+    void unity_zp_pow_sliding_fmpz(unity_zp f, unity_zp g, const fmpz_t pow) noexcept
+    void _unity_zp_reduce_cyclotomic_divmod(unity_zp f) noexcept
+    void _unity_zp_reduce_cyclotomic(unity_zp f) noexcept
+    void unity_zp_reduce_cyclotomic(unity_zp f, const unity_zp g) noexcept
+    void unity_zp_aut(unity_zp f, const unity_zp g, ulong x) noexcept
+    void unity_zp_aut_inv(unity_zp f, const unity_zp g, ulong x) noexcept
+    void unity_zp_jacobi_sum_pq(unity_zp f, ulong q, ulong p) noexcept
+    void unity_zp_jacobi_sum_2q_one(unity_zp f, ulong q) noexcept
+    void unity_zp_jacobi_sum_2q_two(unity_zp f, ulong q) noexcept
+    void unity_zpq_init(unity_zpq f, ulong q, ulong p, const fmpz_t n) noexcept
+    void unity_zpq_clear(unity_zpq f) noexcept
+    void unity_zpq_copy(unity_zpq f, const unity_zpq g) noexcept
+    void unity_zpq_swap(unity_zpq f, unity_zpq q) noexcept
+    bint unity_zpq_equal(const unity_zpq f, const unity_zpq g) noexcept
+    slong unity_zpq_p_unity(const unity_zpq f) noexcept
+    bint unity_zpq_is_p_unity(const unity_zpq f) noexcept
+    bint unity_zpq_is_p_unity_generator(const unity_zpq f) noexcept
+    void unity_zpq_coeff_set_fmpz(unity_zpq f, slong i, slong j, const fmpz_t x) noexcept
+    void unity_zpq_coeff_set_ui(unity_zpq f, slong i, slong j, ulong x) noexcept
+    void unity_zpq_coeff_add(unity_zpq f, slong i, slong j, const fmpz_t x) noexcept
+    void unity_zpq_add(unity_zpq f, const unity_zpq g, const unity_zpq h) noexcept
+    void unity_zpq_mul(unity_zpq f, const unity_zpq g, const unity_zpq h) noexcept
+    void _unity_zpq_mul_unity_p(unity_zpq f) noexcept
+    void unity_zpq_mul_unity_p_pow(unity_zpq f, const unity_zpq g, slong k) noexcept
+    void unity_zpq_pow(unity_zpq f, const unity_zpq g, const fmpz_t p) noexcept
+    void unity_zpq_pow_ui(unity_zpq f, const unity_zpq g, ulong p) noexcept
+    void unity_zpq_gauss_sum(unity_zpq f, ulong q, ulong p) noexcept
+    void unity_zpq_gauss_sum_sigma_pow(unity_zpq f, ulong q, ulong p) noexcept