Merge pull request #402 from oscarbenjamin/pr_flint35_bindings

Update Cython binding declarations for FLINT 3.5

Author: oscarbenjamin

src/flint/flintlib/functions/acb_poly.pxd

@@ -78,6 +78,12 @@ cdef extern from "flint/acb_poly.h":
     void acb_poly_mullow(acb_poly_t C, const acb_poly_t A, const acb_poly_t B, slong n, slong prec)
     void _acb_poly_mul(acb_ptr C, acb_srcptr A, slong lenA, acb_srcptr B, slong lenB, slong prec)
     void acb_poly_mul(acb_poly_t C, const acb_poly_t A1, const acb_poly_t B2, slong prec)
+    void _acb_poly_mulmid_transpose(acb_ptr z, acb_srcptr x, slong xlen, acb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void acb_poly_mulmid_transpose(acb_poly_t res, const acb_poly_t poly1, const acb_poly_t poly2, slong nlo, slong nhi, slong prec)
+    void _acb_poly_mulmid_classical(acb_ptr z, acb_srcptr x, slong xlen, acb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void acb_poly_mulmid_classical(acb_poly_t res, const acb_poly_t poly1, const acb_poly_t poly2, slong nlo, slong nhi, slong prec)
+    void _acb_poly_mulmid(acb_ptr z, acb_srcptr x, slong xlen, acb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void acb_poly_mulmid(acb_poly_t res, const acb_poly_t poly1, const acb_poly_t poly2, slong nlo, slong nhi, slong prec)
     void _acb_poly_inv_series(acb_ptr Qinv, acb_srcptr Q, slong Qlen, slong len, slong prec)
     void acb_poly_inv_series(acb_poly_t Qinv, const acb_poly_t Q, slong n, slong prec)
     void _acb_poly_div_series(acb_ptr Q, acb_srcptr A, slong Alen, acb_srcptr B, slong Blen, slong n, slong prec)
@@ -240,5 +246,9 @@ cdef extern from "flint/acb_poly.h":
     slong _acb_poly_validate_roots(acb_ptr roots, acb_srcptr poly, slong len, slong prec)
     slong _acb_poly_find_roots(acb_ptr roots, acb_srcptr poly, acb_srcptr initial, slong len, slong maxiter, slong prec)
     slong acb_poly_find_roots(acb_ptr roots, const acb_poly_t poly, acb_srcptr initial, slong maxiter, slong prec)
+    double _acb_poly_find_roots_double(acb_ptr roots, acb_srcptr poly, acb_srcptr initial, slong len, slong maxiter, slong prec)
+    double cd_poly_find_roots(double * roots, const double * poly, const double * initial, slong len, slong maxiter, double reltol)
+    double cd_poly_refine_roots(double * z, const double * poly, slong len, double stepsize_bound)
+    double cd_poly_refine_roots_with_pivot(double * z, const double * poly, slong len, double stepsize_bound)
     int _acb_poly_validate_real_roots(acb_srcptr roots, acb_srcptr poly, slong len, slong prec)
     int acb_poly_validate_real_roots(acb_srcptr roots, const acb_poly_t poly, slong prec)

src/flint/flintlib/functions/arb_poly.pxd

@@ -70,6 +70,12 @@ cdef extern from "flint/arb_poly.h":
     void arb_poly_mullow(arb_poly_t C, const arb_poly_t A, const arb_poly_t B, slong n, slong prec)
     void _arb_poly_mul(arb_ptr C, arb_srcptr A, slong lenA, arb_srcptr B, slong lenB, slong prec)
     void arb_poly_mul(arb_poly_t C, const arb_poly_t A, const arb_poly_t B, slong prec)
+    void _arb_poly_mulmid_block(arb_ptr z, arb_srcptr x, slong xlen, arb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void arb_poly_mulmid_block(arb_poly_t res, const arb_poly_t poly1, const arb_poly_t poly2, slong nlo, slong nhi, slong prec)
+    void _arb_poly_mulmid_classical(arb_ptr z, arb_srcptr x, slong xlen, arb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void arb_poly_mulmid_classical(arb_poly_t res, const arb_poly_t poly1, const arb_poly_t poly2, slong nlo, slong nhi, slong prec)
+    void _arb_poly_mulmid(arb_ptr z, arb_srcptr x, slong xlen, arb_srcptr y, slong ylen, slong nlo, slong nhi, slong prec)
+    void arb_poly_mulmid(arb_poly_t res, const arb_poly_t poly1, const arb_poly_t poly2, slong nlo, slong nhi, slong prec)
     void _arb_poly_inv_series(arb_ptr Q, arb_srcptr A, slong Alen, slong len, slong prec)
     void arb_poly_inv_series(arb_poly_t Q, const arb_poly_t A, slong n, slong prec)
     void _arb_poly_div_series(arb_ptr Q, arb_srcptr A, slong Alen, arb_srcptr B, slong Blen, slong n, slong prec)

src/flint/flintlib/functions/arith.pxd

@@ -15,9 +15,6 @@ cdef extern from "flint/arith.h":
     void arith_stirling_number_1u_vec(fmpz_struct * row, ulong n, slong klen)
     void arith_stirling_number_1_vec(fmpz_struct * row, ulong n, slong klen)
     void arith_stirling_number_2_vec(fmpz_struct * row, ulong n, slong klen)
-    void arith_stirling_number_1u_vec_next(fmpz_struct * row, const fmpz_struct * prev, slong n, slong klen)
-    void arith_stirling_number_1_vec_next(fmpz_struct * row, const fmpz_struct * prev, slong n, slong klen)
-    void arith_stirling_number_2_vec_next(fmpz_struct * row, const fmpz_struct * prev, slong n, slong klen)
     void arith_stirling_matrix_1u(fmpz_mat_t mat)
     void arith_stirling_matrix_1(fmpz_mat_t mat)
     void arith_stirling_matrix_2(fmpz_mat_t mat)

src/flint/flintlib/functions/flint.pxd

@@ -34,9 +34,12 @@ cdef extern from "flint/flint.h":
     void flint_reset_num_workers(int num_workers)
     # int flint_printf(const char * format, ...)
     # int flint_fprintf(FILE * fs, const char * format, ...)
+    # int flint_snprintf(char * s, size_t n, const char * str, ...)
+    # int flint_sprintf(char * s, const char * str, ...)
     # int flint_vprintf(const char * format, va_list vlist)
     # int flint_vfprintf(FILE * fs, const char * format, va_list vlist)
-    # int flint_sprintf(char * s, const char * str, ...)
+    # int flint_vsnprintf(char * s, size_t n, const char * str, va_list vlist)
+    # int flint_vsprintf(char * s, const char * str, va_list vlist)
     # int flint_scanf(const char * str, ...)
     # int flint_fscanf(FILE * f, const char * str, ...)
     # int flint_sscanf(const char * s, const char * str, ...)

src/flint/flintlib/functions/fmpz_mat.pxd

@@ -163,6 +163,9 @@ cdef extern from "flint/fmpz_mat.h":
     void fmpz_mat_snf_kannan_bachem(fmpz_mat_t S, const fmpz_mat_t A)
     void fmpz_mat_snf_iliopoulos(fmpz_mat_t S, const fmpz_mat_t A, const fmpz_t mod)
     int fmpz_mat_is_in_snf(const fmpz_mat_t A)
+    int fmpz_mat_is_diagonal(const fmpz_mat_t A)
+    void fmpz_mat_snf_transform(fmpz_mat_t S, fmpz_mat_t U, fmpz_mat_t V, const fmpz_mat_t A)
+    slong fmpz_mat_elementary_divisors(fmpz_struct * ed, const fmpz_mat_t A)
     void fmpz_mat_gram(fmpz_mat_t B, const fmpz_mat_t A)
     int fmpz_mat_is_hadamard(const fmpz_mat_t H)
     int fmpz_mat_hadamard(fmpz_mat_t H)

src/flint/flintlib/functions/fmpz_mod_poly.pxd

@@ -78,6 +78,8 @@ cdef extern from "flint/fmpz_mod_poly.h":
     void fmpz_mod_poly_mul(fmpz_mod_poly_t res, const fmpz_mod_poly_t poly1, const fmpz_mod_poly_t poly2, const fmpz_mod_ctx_t ctx)
     void _fmpz_mod_poly_mullow(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong n, const fmpz_mod_ctx_t ctx)
     void fmpz_mod_poly_mullow(fmpz_mod_poly_t res, const fmpz_mod_poly_t poly1, const fmpz_mod_poly_t poly2, slong n, const fmpz_mod_ctx_t ctx)
+    void _fmpz_mod_poly_mulmid(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong nlo, slong nhi, const fmpz_mod_ctx_t ctx)
+    void fmpz_mod_poly_mulmid(fmpz_mod_poly_t res, const fmpz_mod_poly_t poly1, const fmpz_mod_poly_t poly2, slong nlo, slong nhi, const fmpz_mod_ctx_t ctx)
     void _fmpz_mod_poly_sqr(fmpz_struct * res, const fmpz_struct * poly, slong len, const fmpz_mod_ctx_t ctx)
     void fmpz_mod_poly_sqr(fmpz_mod_poly_t res, const fmpz_mod_poly_t poly, const fmpz_mod_ctx_t ctx)
     void fmpz_mod_poly_mulhigh(fmpz_mod_poly_t res, const fmpz_mod_poly_t poly1, const fmpz_mod_poly_t poly2, slong start, const fmpz_mod_ctx_t ctx)

src/flint/flintlib/functions/fmpz_poly.pxd

@@ -91,8 +91,8 @@ cdef extern from "flint/fmpz_poly.h":
     slong _fmpz_poly_remove_content_2exp(fmpz_struct * pol, slong len)
     void _fmpz_poly_scale_2exp(fmpz_struct * pol, slong len, slong k)
     void _fmpz_poly_bit_pack(nn_ptr arr, const fmpz_struct * poly, slong len, flint_bitcnt_t bit_size, int negate)
-    int _fmpz_poly_bit_unpack(fmpz_struct * poly, slong len, nn_srcptr arr, flint_bitcnt_t bit_size, int negate)
-    void _fmpz_poly_bit_unpack_unsigned(fmpz_struct * poly, slong len, nn_srcptr arr, flint_bitcnt_t bit_size)
+    int _fmpz_poly_bit_unpack(fmpz_struct * poly, slong nlo, slong nhi, nn_srcptr arr, flint_bitcnt_t bit_size, int negate)
+    void _fmpz_poly_bit_unpack_unsigned(fmpz_struct * poly, slong nlo, slong nhi, nn_srcptr arr, flint_bitcnt_t bit_size)
     void fmpz_poly_bit_pack(fmpz_t f, const fmpz_poly_t poly, flint_bitcnt_t bit_size)
     void fmpz_poly_bit_unpack(fmpz_poly_t poly, const fmpz_t f, flint_bitcnt_t bit_size)
     void fmpz_poly_bit_unpack_unsigned(fmpz_poly_t poly, const fmpz_t f, flint_bitcnt_t bit_size)
@@ -102,8 +102,14 @@ cdef extern from "flint/fmpz_poly.h":
     void fmpz_poly_mullow_classical(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2, slong n)
     void _fmpz_poly_mulhigh_classical(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong start)
     void fmpz_poly_mulhigh_classical(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2, slong start)
-    void _fmpz_poly_mulmid_classical(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2)
+    void _fmpz_poly_mulmid_classical(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong nlo, slong nhi)
+    void _fmpz_poly_mulmid_KS(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong nlo, slong nhi)
+    void _fmpz_poly_mulmid_SS(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong nlo, slong nhi)
+    void _fmpz_poly_mulmid(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong nlo, slong nhi)
     void fmpz_poly_mulmid_classical(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
+    void fmpz_poly_mulmid_KS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
+    void fmpz_poly_mulmid_SS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
+    void fmpz_poly_mulmid(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
     void _fmpz_poly_mul_karatsuba(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2)
     void fmpz_poly_mul_karatsuba(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
     void _fmpz_poly_mullow_karatsuba_n(fmpz_struct * res, const fmpz_struct * poly1, const fmpz_struct * poly2, slong n)
@@ -114,9 +120,9 @@ cdef extern from "flint/fmpz_poly.h":
     void fmpz_poly_mul_KS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
     void _fmpz_poly_mullow_KS(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong n)
     void fmpz_poly_mullow_KS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2, slong n)
-    void _fmpz_poly_mul_SS(fmpz_struct * output, const fmpz_struct * input1, slong length1, const fmpz_struct * input2, slong length2)
+    void _fmpz_poly_mul_SS(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2)
     void fmpz_poly_mul_SS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)
-    void _fmpz_poly_mullow_SS(fmpz_struct * output, const fmpz_struct * input1, slong length1, const fmpz_struct * input2, slong length2, slong n)
+    void _fmpz_poly_mullow_SS(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2, slong n)
     void fmpz_poly_mullow_SS(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2, slong n)
     void _fmpz_poly_mul(fmpz_struct * res, const fmpz_struct * poly1, slong len1, const fmpz_struct * poly2, slong len2)
     void fmpz_poly_mul(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2)

src/flint/flintlib/functions/gr.pxd

@@ -15,6 +15,7 @@ from flint.flintlib.types.gr cimport gr_ctx_t, gr_ptr, gr_srcptr, gr_stream_t, g
 # .. macro:: GR_TMP_CLEAR(x1, ctx)
 
 cdef extern from "flint/gr.h":
+    void gr_ctx_uninitialized(gr_ctx_t ctx)
     truth_t truth_and(truth_t x, truth_t y)
     truth_t truth_or(truth_t x, truth_t y)
     truth_t truth_not(truth_t x)
@@ -31,6 +32,7 @@ cdef extern from "flint/gr.h":
     int gr_ctx_set_gen_names(gr_ctx_t ctx, const char ** s)
     int gr_ctx_ngens(slong * ngens, gr_ctx_t ctx)
     int gr_ctx_gen_name(char ** name, slong i, gr_ctx_t ctx)
+    gr_ptr gr_ctx_base(gr_ctx_t ctx)
     void gr_init(gr_ptr res, gr_ctx_t ctx)
     void gr_clear(gr_ptr res, gr_ctx_t ctx)
     void gr_swap(gr_ptr x, gr_ptr y, gr_ctx_t ctx)

src/flint/flintlib/functions/gr_domains.pxd

@@ -4,6 +4,7 @@ from flint.flintlib.types.fmpz cimport fmpz_poly_t
 from flint.flintlib.types.gr cimport gr_ctx_t, truth_t
 from flint.flintlib.types.mpoly cimport ordering_t
 
+# unknown type ore_algebra_t
 
 
 cdef extern from "flint/gr.h":
@@ -22,6 +23,7 @@ cdef extern from "flint/gr.h":
     truth_t gr_ctx_is_finite_characteristic(gr_ctx_t ctx)
     truth_t gr_ctx_is_ordered_ring(gr_ctx_t ctx)
     truth_t gr_ctx_is_zero_ring(gr_ctx_t ctx)
+    truth_t gr_ctx_is_approx_commutative_ring(gr_ctx_t ctx)
     truth_t gr_ctx_is_exact(gr_ctx_t ctx)
     truth_t gr_ctx_is_canonical(gr_ctx_t ctx)
     truth_t gr_ctx_has_real_prec(gr_ctx_t ctx)
@@ -32,13 +34,12 @@ cdef extern from "flint/gr.h":
     void gr_ctx_init_psl2z(gr_ctx_t ctx)
     int gr_ctx_init_dirichlet_group(gr_ctx_t ctx, ulong q)
     void gr_ctx_init_random(gr_ctx_t ctx, flint_rand_t state)
+    void gr_ctx_init_random_commutative_ring(gr_ctx_t ctx, flint_rand_t state)
+    void gr_ctx_init_random_field(gr_ctx_t ctx, flint_rand_t state)
     void gr_ctx_init_fmpz(gr_ctx_t ctx)
     void gr_ctx_init_fmpq(gr_ctx_t ctx)
     void gr_ctx_init_fmpzi(gr_ctx_t ctx)
     int gr_ctx_set_is_field(gr_ctx_t ctx, truth_t is_field)
-    void gr_ctx_init_nmod(gr_ctx_t ctx, ulong n)
-    void gr_ctx_init_nmod8(gr_ctx_t ctx, unsigned char n)
-    void gr_ctx_init_nmod32(gr_ctx_t ctx, unsigned int n)
     void gr_ctx_init_fmpz_mod(gr_ctx_t ctx, const fmpz_t n)
     void gr_ctx_init_fq(gr_ctx_t ctx, const fmpz_t p, slong d, const char * var)
     void gr_ctx_init_fq_nmod(gr_ctx_t ctx, ulong p, slong d, const char * var)
@@ -58,6 +59,7 @@ cdef extern from "flint/gr.h":
     void gr_ctx_init_complex_algebraic_ca(gr_ctx_t ctx)
     void gr_ctx_ca_set_option(gr_ctx_t ctx, slong option, slong value)
     slong gr_ctx_ca_get_option(gr_ctx_t ctx, slong option)
+    void gr_ctx_init_gr_complex(gr_ctx_t ctx, gr_ctx_t real_ctx)
     void gr_ctx_init_complex_extended_ca(gr_ctx_t ctx)
     void gr_ctx_init_real_float_arf(gr_ctx_t ctx, slong prec)
     void gr_ctx_init_complex_float_acf(gr_ctx_t ctx, slong prec)
@@ -69,8 +71,13 @@ cdef extern from "flint/gr.h":
     void gr_ctx_init_fmpz_poly(gr_ctx_t ctx)
     void gr_ctx_init_fmpq_poly(gr_ctx_t ctx)
     void gr_ctx_init_gr_poly(gr_ctx_t ctx, gr_ctx_t base_ring)
+    void gr_ctx_init_random_poly(gr_ctx_t ctx, flint_rand_t state)
     void gr_ctx_init_fmpz_mpoly(gr_ctx_t ctx, slong nvars, const ordering_t ord)
+    void gr_ctx_init_fmpq_mpoly(gr_ctx_t ctx, slong nvars, const ordering_t ord)
     void gr_ctx_init_gr_mpoly(gr_ctx_t ctx, gr_ctx_t base_ring, slong nvars, const ordering_t ord)
+    void gr_ctx_init_random_mpoly(gr_ctx_t ctx, flint_rand_t state)
+    # void gr_ctx_init_gr_ore_poly(gr_ctx_t ctx, gr_ctx_t base_ring, slong base_var, const ore_algebra_t which_algebra)
+    void gr_ctx_init_random_series(gr_ctx_t ctx, flint_rand_t state)
     void gr_ctx_init_gr_fraction(gr_ctx_t ctx, gr_ctx_t domain, int flags)
     void gr_ctx_init_fmpz_mpoly_q(gr_ctx_t ctx, slong nvars, const ordering_t ord)
     void gr_ctx_init_fmpz_mod_mpoly_q(gr_ctx_t ctx, slong nvars, const ordering_t ord, const fmpz_t mod)

src/flint/flintlib/functions/gr_mpoly.pxd

@@ -45,8 +45,8 @@ cdef extern from "flint/gr_mpoly.h":
     # int gr_mpoly_set_other(gr_mpoly_t res, gr_srcptr A, gr_ctx_t A_ctx, gr_mpoly_ctx_t ctx)
     # truth_t gr_mpoly_equal(const gr_mpoly_t A, const gr_mpoly_t B, gr_mpoly_ctx_t ctx)
     # int gr_mpoly_randtest_bits(gr_mpoly_t A, flint_rand_t state, slong length, flint_bitcnt_t exp_bits, gr_mpoly_ctx_t ctx)
-    # int gr_mpoly_write_pretty(gr_stream_t out, const gr_mpoly_t A, const char ** x, gr_mpoly_ctx_t ctx)
-    # int gr_mpoly_print_pretty(const gr_mpoly_t A, const char ** x, gr_mpoly_ctx_t ctx)
+    # int gr_mpoly_write_pretty(gr_stream_t out, const gr_mpoly_t A, gr_mpoly_ctx_t ctx)
+    # int gr_mpoly_print_pretty(const gr_mpoly_t A, gr_mpoly_ctx_t ctx)
     # int gr_mpoly_get_coeff_scalar_fmpz(gr_ptr c, const gr_mpoly_t A, const fmpz_struct * exp, gr_mpoly_ctx_t ctx)
     # int gr_mpoly_get_coeff_scalar_ui(gr_ptr c, const gr_mpoly_t A, const ulong * exp, gr_mpoly_ctx_t ctx)
     # int gr_mpoly_set_coeff_scalar_fmpz(gr_mpoly_t A, gr_srcptr c, const fmpz_struct * exp, gr_mpoly_ctx_t ctx)