FLINT 3.7 (in development)

FLINT 3.7.0 (in development)#

Replace GMP's mpz_gcd, mpz_gcdext and mpz_invert in fmpz_gcd, fmpz_xgcd, fmpz_invmod and related functions by ports of GMP's half-GCD code using FLINT's multiplication and FFTs. GCDs are typically 1.1x faster than with GMP 6.3.0 up to 1000 limbs and 2x faster at $10^6$ limbs, extended GCDs up to 3x faster. Adds flint_mpz_gcd, flint_mpz_gcdext, flint_mpz_invert, flint_mpn_gcdext, _flint_mpn_gcd_inplace and _flint_mpn_gcdext_inplace [FJ].
flint_mpn_gcd_full is renamed flint_mpn_gcd and handles small operands faster. flint_mpn_gcd_full and flint_mpn_gcd_full2 remain as deprecated macros [FJ].
Faster contents of integer vectors (_fmpz_vec_content and the contents and canonical forms of fmpz_poly, fmpq_poly, fmpz_mat, fmpz_mpoly, fmpq_mpoly and fmpz_mpoly_q), using divisibility tests instead of chained GCDs: 10x to 75x faster for long vectors of single-limb entries with a content other than 1. Adds flint_mpn_gcd_vec_1 and flint_mpn_gcd_vec [FJ].
Faster LCMs of integer vectors and clearing of denominators (_fmpz_vec_lcm, _fmpq_vec_get_fmpz_vec_fmpz, the fmpq_mat_get_fmpz_mat functions and their uses in fmpq_poly and gr): for 1000 denominators that are 1, equal or divisors of a common denominator, the LCM is 5x to 18x faster and the clearing of denominators 2x to 6x faster. Adds _fmpz_vec_lcm_strided and _fmpq_vec_get_fmpz_vec_scaled [FJ].

Fix the modular inverse in the curve selection of fmpz_factor_ecm, which was wrong when its cofactor had fewer limbs than the modulus (wasting the curve), and which could call GMP's mpn_gcdext with its operands in the wrong order [FJ].
Fix crashes in finite-field multivariate Wang and Zippel factorization when an evaluation gives incompatible leading coefficient divisors [#2905].
Fix arb_pow losing containment when the output aliases an inexact exponent and the base interval has a zero endpoint [#2879].
Fix arb_pow returning nan for a small exact integer exponent when the base has a zero midpoint and a nonzero radius [#2908].