/dports/math/flint2/flint-2.8.4/mpoly/ |
H A D | gen_shift_left_right.c | 22 void _mpoly_gen_shift_right(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_right() argument 48 mpoly_monomial_msub(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_right() 49 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_right() 57 ulong * Aexp, in _mpoly_gen_shift_right_fmpz() argument 83 mpoly_monomial_sub_mp(Aexp + N*i, Aexp + N*i, gen, N); in _mpoly_gen_shift_right_fmpz() 91 mpoly_monomial_sub(Aexp + N*i, Aexp + N*i, gen, N); in _mpoly_gen_shift_right_fmpz() 98 void _mpoly_gen_shift_left(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_left() argument 125 mpoly_monomial_madd(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_left() 126 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_left()
|
H A D | cmp_general.c | 16 int mpoly_monomial_cmp_general(ulong * Aexp, flint_bitcnt_t Abits, in mpoly_monomial_cmp_general() argument 29 return mpoly_monomial_cmp_nomask(Aexp, Bexp, N); in mpoly_monomial_cmp_general() 40 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 42 if ((Aexp[i]^himask) > (Bexp[i]^himask)) in mpoly_monomial_cmp_general() 55 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 57 if (Aexp[i] > Bexp[i]) in mpoly_monomial_cmp_general() 68 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 70 if (Aexp[i] < Bexp[i]) in mpoly_monomial_cmp_general() 92 newAexp = Aexp; in mpoly_monomial_cmp_general() 103 mpoly_repack_monomials(newAexp, newbits, Aexp, Abits, 1, mctx); in mpoly_monomial_cmp_general()
|
H A D | monomial_cofactors.c | 32 fmpz * Aexp, * Bexp, * minAexp, * minBexp; in mpoly_monomial_cofactors() local 43 Aexp = (fmpz *) TMP_ALLOC(4*nvars*sizeof(fmpz)); in mpoly_monomial_cofactors() 44 Bexp = Aexp + 1*nvars; in mpoly_monomial_cofactors() 45 minAexp = Aexp + 2*nvars; in mpoly_monomial_cofactors() 46 minBexp = Aexp + 3*nvars; in mpoly_monomial_cofactors() 50 fmpz_init(Aexp + j); in mpoly_monomial_cofactors() 64 mpoly_get_monomial_ffmpz(Aexp, Aexps + NA*i, Abits, mctx); in mpoly_monomial_cofactors() 66 _fmpz_vec_min_inplace(minAexp, Aexp, nvars); in mpoly_monomial_cofactors() 71 fmpz_add(t2, Bbarexps + j, Aexp + j); in mpoly_monomial_cofactors() 89 fmpz_clear(Aexp + j); in mpoly_monomial_cofactors()
|
H A D | reverse.c | 15 void mpoly_reverse(ulong * Aexp, const ulong * Bexp, slong len, slong N) in mpoly_reverse() argument 19 if (Aexp == Bexp) in mpoly_reverse() 23 mpoly_monomial_swap(Aexp + N*i, Aexp + N*(len - i - 1), N); in mpoly_reverse() 30 mpoly_monomial_set(Aexp + N*i, Bexp + N*(len - i - 1), N); in mpoly_reverse()
|
H A D | main_variable_split.c | 15 void mpoly_main_variable_split_LEX(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_LEX() argument 23 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_LEX() 28 e = (e * mults[j]) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_LEX() 37 void mpoly_main_variable_split_DEG(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_DEG() argument 45 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_DEG() 50 e = (e * deg) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_DEG()
|
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/mpoly/ |
H A D | cmp_general.c | 16 int mpoly_monomial_cmp_general(ulong * Aexp, flint_bitcnt_t Abits, in mpoly_monomial_cmp_general() argument 29 return mpoly_monomial_cmp_nomask(Aexp, Bexp, N); in mpoly_monomial_cmp_general() 40 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 42 if ((Aexp[i]^himask) > (Bexp[i]^himask)) in mpoly_monomial_cmp_general() 55 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 57 if (Aexp[i] > Bexp[i]) in mpoly_monomial_cmp_general() 68 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 70 if (Aexp[i] < Bexp[i]) in mpoly_monomial_cmp_general() 92 newAexp = Aexp; in mpoly_monomial_cmp_general() 103 mpoly_repack_monomials(newAexp, newbits, Aexp, Abits, 1, mctx); in mpoly_monomial_cmp_general()
|
H A D | monomial_cofactors.c | 32 fmpz * Aexp, * Bexp, * minAexp, * minBexp; in mpoly_monomial_cofactors() local 43 Aexp = (fmpz *) TMP_ALLOC(4*nvars*sizeof(fmpz)); in mpoly_monomial_cofactors() 44 Bexp = Aexp + 1*nvars; in mpoly_monomial_cofactors() 45 minAexp = Aexp + 2*nvars; in mpoly_monomial_cofactors() 46 minBexp = Aexp + 3*nvars; in mpoly_monomial_cofactors() 50 fmpz_init(Aexp + j); in mpoly_monomial_cofactors() 64 mpoly_get_monomial_ffmpz(Aexp, Aexps + NA*i, Abits, mctx); in mpoly_monomial_cofactors() 66 _fmpz_vec_min_inplace(minAexp, Aexp, nvars); in mpoly_monomial_cofactors() 71 fmpz_add(t2, Bbarexps + j, Aexp + j); in mpoly_monomial_cofactors() 89 fmpz_clear(Aexp + j); in mpoly_monomial_cofactors()
|
H A D | gen_shift_left_right.c | 22 void _mpoly_gen_shift_right(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_right() argument 48 mpoly_monomial_msub(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_right() 49 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_right() 55 void _mpoly_gen_shift_left(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_left() argument 82 mpoly_monomial_madd(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_left() 83 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_left()
|
H A D | reverse.c | 15 void mpoly_reverse(ulong * Aexp, const ulong * Bexp, slong len, slong N) in mpoly_reverse() argument 19 if (Aexp == Bexp) in mpoly_reverse() 23 mpoly_monomial_swap(Aexp + N*i, Aexp + N*(len - i - 1), N); in mpoly_reverse() 30 mpoly_monomial_set(Aexp + N*i, Bexp + N*(len - i - 1), N); in mpoly_reverse()
|
H A D | main_variable_split.c | 15 void mpoly_main_variable_split_LEX(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_LEX() argument 23 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_LEX() 28 e = (e * mults[j]) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_LEX() 37 void mpoly_main_variable_split_DEG(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_DEG() argument 45 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_DEG() 50 e = (e * deg) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_DEG()
|
/dports/math/e-antic/flint2-ae7ec89/mpoly/ |
H A D | cmp_general.c | 16 int mpoly_monomial_cmp_general(ulong * Aexp, flint_bitcnt_t Abits, in mpoly_monomial_cmp_general() argument 29 return mpoly_monomial_cmp_nomask(Aexp, Bexp, N); in mpoly_monomial_cmp_general() 40 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 42 if ((Aexp[i]^himask) > (Bexp[i]^himask)) in mpoly_monomial_cmp_general() 55 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 57 if (Aexp[i] > Bexp[i]) in mpoly_monomial_cmp_general() 68 if (Aexp[i] != Bexp[i]) in mpoly_monomial_cmp_general() 70 if (Aexp[i] < Bexp[i]) in mpoly_monomial_cmp_general() 92 newAexp = Aexp; in mpoly_monomial_cmp_general() 103 mpoly_repack_monomials(newAexp, newbits, Aexp, Abits, 1, mctx); in mpoly_monomial_cmp_general()
|
H A D | monomial_cofactors.c | 32 fmpz * Aexp, * Bexp, * minAexp, * minBexp; in mpoly_monomial_cofactors() local 43 Aexp = (fmpz *) TMP_ALLOC(4*nvars*sizeof(fmpz)); in mpoly_monomial_cofactors() 44 Bexp = Aexp + 1*nvars; in mpoly_monomial_cofactors() 45 minAexp = Aexp + 2*nvars; in mpoly_monomial_cofactors() 46 minBexp = Aexp + 3*nvars; in mpoly_monomial_cofactors() 50 fmpz_init(Aexp + j); in mpoly_monomial_cofactors() 64 mpoly_get_monomial_ffmpz(Aexp, Aexps + NA*i, Abits, mctx); in mpoly_monomial_cofactors() 66 _fmpz_vec_min_inplace(minAexp, Aexp, nvars); in mpoly_monomial_cofactors() 71 fmpz_add(t2, Bbarexps + j, Aexp + j); in mpoly_monomial_cofactors() 89 fmpz_clear(Aexp + j); in mpoly_monomial_cofactors()
|
H A D | gen_shift_left_right.c | 22 void _mpoly_gen_shift_right(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_right() argument 48 mpoly_monomial_msub(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_right() 49 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_right() 55 void _mpoly_gen_shift_left(ulong * Aexp, flint_bitcnt_t Abits, slong Alength, in _mpoly_gen_shift_left() argument 82 mpoly_monomial_madd(Aexp + N*i, Aexp + N*i, amount, one, N); in _mpoly_gen_shift_left() 83 FLINT_ASSERT(!mpoly_monomial_overflows(Aexp + N*i, N, mask)); in _mpoly_gen_shift_left()
|
H A D | reverse.c | 15 void mpoly_reverse(ulong * Aexp, const ulong * Bexp, slong len, slong N) in mpoly_reverse() argument 19 if (Aexp == Bexp) in mpoly_reverse() 23 mpoly_monomial_swap(Aexp + N*i, Aexp + N*(len - i - 1), N); in mpoly_reverse() 30 mpoly_monomial_set(Aexp + N*i, Bexp + N*(len - i - 1), N); in mpoly_reverse()
|
H A D | main_variable_split.c | 15 void mpoly_main_variable_split_LEX(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_LEX() argument 23 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_LEX() 28 e = (e * mults[j]) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_LEX() 37 void mpoly_main_variable_split_DEG(slong * ind, ulong * pexp, const ulong * Aexp, in mpoly_main_variable_split_DEG() argument 45 slong top = Aexp[i] >> (Abits*num); in mpoly_main_variable_split_DEG() 50 e = (e * deg) + ((Aexp[i] >> (j*Abits)) & mask); in mpoly_main_variable_split_DEG()
|
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/fq_nmod_mpoly/ |
H A D | interp.c | 41 ulong * Aexp; in nmod_mpolyn_interp_reduce_lg_poly() local 50 Aexp = A->exps; in nmod_mpolyn_interp_reduce_lg_poly() 82 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_poly() local 92 Aexp = A->exps; in nmod_mpolyn_interp_lift_lg_poly() 360 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_mpolyn() local 826 ulong * Aexp; in nmod_mpolyun_interp_crt_lg_mpolyu() local 1021 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_sm_poly() local 1059 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_poly() local 1318 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_mpolyn() local 1531 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_lg_poly() local [all …]
|
/dports/math/e-antic/flint2-ae7ec89/fq_nmod_mpoly/ |
H A D | interp.c | 41 ulong * Aexp; in nmod_mpolyn_interp_reduce_lg_poly() local 50 Aexp = A->exps; in nmod_mpolyn_interp_reduce_lg_poly() 82 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_poly() local 92 Aexp = A->exps; in nmod_mpolyn_interp_lift_lg_poly() 360 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_mpolyn() local 826 ulong * Aexp; in nmod_mpolyun_interp_crt_lg_mpolyu() local 1021 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_sm_poly() local 1059 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_poly() local 1318 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_mpolyn() local 1531 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_lg_poly() local [all …]
|
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/nmod_mpoly/ |
H A D | interp.c | 88 ulong * Aexp; in nmod_mpolyn_interp_reduce_2sm_poly() local 95 Aexp = A->exps; in nmod_mpolyn_interp_reduce_2sm_poly() 143 e = Aexp; in nmod_mpolyn_interp_lift_2sm_poly() 179 if (e == Aexp) in nmod_mpolyn_interp_lift_2sm_poly() 183 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_lift_2sm_poly() 258 if (Aexp >= 0) in nmod_mpolyn_interp_crt_2sm_poly() 329 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_crt_2sm_poly() 439 ulong * Aexp; in nmod_mpolyn_interp_lift_sm_mpolyn() local 444 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpolyn() 1141 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpoly() [all …]
|
H A D | mpolyu_divides.c | 33 ulong * Aexp = *A_exp; in _nmod_mpoly_mulsub1() local 71 Aexp[Alen] = Dexp[Di]; in _nmod_mpoly_mulsub1() 79 Aexp[Alen] = exp; in _nmod_mpoly_mulsub1() 161 *A_exp = Aexp; in _nmod_mpoly_mulsub1() 189 ulong * Aexp = *A_exp; in _nmod_mpoly_mulsub() local 260 mpoly_monomial_set(Aexp + N*Alen, exp, N); in _nmod_mpoly_mulsub() 359 *A_exp = Aexp; in _nmod_mpoly_mulsub() 385 ulong * Aexp = A->exps; in nmod_mpolyuu_divides() local 400 ulong mask, exp, maxexp = Aexp[Alen - 1]; in nmod_mpolyuu_divides() 440 HEAP_ASSIGN(heap[1], Aexp[0], x); in nmod_mpolyuu_divides() [all …]
|
/dports/math/e-antic/flint2-ae7ec89/nmod_mpoly/ |
H A D | interp.c | 88 ulong * Aexp; in nmod_mpolyn_interp_reduce_2sm_poly() local 95 Aexp = A->exps; in nmod_mpolyn_interp_reduce_2sm_poly() 143 e = Aexp; in nmod_mpolyn_interp_lift_2sm_poly() 179 if (e == Aexp) in nmod_mpolyn_interp_lift_2sm_poly() 183 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_lift_2sm_poly() 258 if (Aexp >= 0) in nmod_mpolyn_interp_crt_2sm_poly() 329 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_crt_2sm_poly() 439 ulong * Aexp; in nmod_mpolyn_interp_lift_sm_mpolyn() local 444 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpolyn() 1141 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpoly() [all …]
|
/dports/math/flint2/flint-2.8.4/n_poly/ |
H A D | n_polyu1n_gcd.c | 44 slong Fi, Aexp, Bexp; in n_polyu1n_mod_interp_lift_2sm_poly() local 60 e = Aexp; in n_polyu1n_mod_interp_lift_2sm_poly() 63 if (Aexp == Bexp) in n_polyu1n_mod_interp_lift_2sm_poly() 96 if (e == Aexp) in n_polyu1n_mod_interp_lift_2sm_poly() 99 Aexp--; in n_polyu1n_mod_interp_lift_2sm_poly() 100 } while (Aexp >= 0 && Acoeffs[Aexp] == 0); in n_polyu1n_mod_interp_lift_2sm_poly() 170 if (Aexp >= 0) in n_polyu1n_mod_interp_crt_2sm_poly() 196 if (e == Aexp) in n_polyu1n_mod_interp_crt_2sm_poly() 217 if (e == Aexp) in n_polyu1n_mod_interp_crt_2sm_poly() 220 Aexp--; in n_polyu1n_mod_interp_crt_2sm_poly() [all …]
|
/dports/math/flint2/flint-2.8.4/nmod_mpoly/ |
H A D | interp.c | 88 ulong * Aexp; in nmod_mpolyn_interp_reduce_2sm_poly() local 95 Aexp = A->exps; in nmod_mpolyn_interp_reduce_2sm_poly() 143 e = Aexp; in nmod_mpolyn_interp_lift_2sm_poly() 179 if (e == Aexp) in nmod_mpolyn_interp_lift_2sm_poly() 183 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_lift_2sm_poly() 258 if (Aexp >= 0) in nmod_mpolyn_interp_crt_2sm_poly() 329 } while (Aexp >= 0 && Acoeff[Aexp] == 0); in nmod_mpolyn_interp_crt_2sm_poly() 611 ulong * Aexp; in nmod_mpolyn_interp_lift_sm_mpolyn() local 616 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpolyn() 1301 Aexp = A->exps; in nmod_mpolyn_interp_lift_sm_mpoly() [all …]
|
/dports/math/flint2/flint-2.8.4/fmpz_mod_mpoly_factor/ |
H A D | interp.c | 245 e = Aexp; in fmpz_mod_polyu1n_interp_lift_2sm_poly() 281 if (e == Aexp) in fmpz_mod_polyu1n_interp_lift_2sm_poly() 285 } while (Aexp >= 0 && fmpz_is_zero(Acoeff + Aexp)); in fmpz_mod_polyu1n_interp_lift_2sm_poly() 374 if (Aexp >= 0) in fmpz_mod_polyu1n_interp_crt_2sm_poly() 400 if (e == Aexp) in fmpz_mod_polyu1n_interp_crt_2sm_poly() 425 } while (Aexp >= 0 && fmpz_is_zero(Acoeff + Aexp)); in fmpz_mod_polyu1n_interp_crt_2sm_poly() 466 ulong * Aexp; in fmpz_mod_mpolyn_interp_reduce_sm_poly() local 475 Aexp = A->exps; in fmpz_mod_mpolyn_interp_reduce_sm_poly() 502 ulong * Aexp; in fmpz_mod_mpolyn_interp_lift_sm_poly() local 511 Aexp = A->exps; in fmpz_mod_mpolyn_interp_lift_sm_poly() [all …]
|
/dports/math/flint2/flint-2.8.4/fq_nmod_mpoly/ |
H A D | interp.c | 57 ulong * Aexp; in nmod_mpolyn_interp_reduce_lg_poly() local 66 Aexp = A->exps; in nmod_mpolyn_interp_reduce_lg_poly() 99 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_poly() local 109 Aexp = A->exps; in nmod_mpolyn_interp_lift_lg_poly() 581 ulong * Aexp; in nmod_mpolyn_interp_lift_lg_mpolyn() local 1079 ulong * Aexp; in nmod_mpolyun_interp_crt_lg_mpolyu() local 1331 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_sm_poly() local 1369 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_poly() local 1817 ulong * Aexp; in fq_nmod_mpolyn_interp_lift_sm_mpolyn() local 2091 ulong * Aexp; in fq_nmod_mpolyn_interp_reduce_lg_poly() local [all …]
|
/dports/math/flint2/flint-2.8.4/fmpz_mpoly_factor/ |
H A D | interp.c | 129 ulong * Aexp = A->exps; in fmpz_mpoly_interp_reduce_p_mpolyn() local 146 k = ((Aexp + N*Ai)[offset] >> shift) & mask; in fmpz_mpoly_interp_reduce_p_mpolyn() 153 Aexp + N*Ai, N, offset, -(k << shift))) in fmpz_mpoly_interp_reduce_p_mpolyn() 160 Eexp + N*(Ei - 1), Aexp + N*Ai, N, offset, -(k << shift))); in fmpz_mpoly_interp_reduce_p_mpolyn() 198 ulong * Aexp; in fmpz_mpoly_interp_lift_p_mpolyn() local 206 Aexp = A->exps; in fmpz_mpoly_interp_lift_p_mpolyn() 217 Aexp = A->exps; in fmpz_mpoly_interp_lift_p_mpolyn() 257 ulong * Aexp = A->exps; in fmpz_mpoly_interp_crt_p_mpolyn() local 297 Aexp + N*Ai, N, offset, vi << shift)) in fmpz_mpoly_interp_crt_p_mpolyn() 317 Fexp + N*Fi, Aexp + N*Ai, N, offset, vi << shift))) in fmpz_mpoly_interp_crt_p_mpolyn() [all …]
|