/dports/math/arb/arb-2.21.1/arb/ |
H A D | exp_sum_bs_simple.c | 15 bsplit(fmpz_t P, fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, const fmpz_t x, in bsplit() argument 22 *Qexp = r; in bsplit() 37 bsplit(P, T, Q, Qexp, x, r, a, m, 1); in bsplit() 45 *Qexp = *Qexp + *Q2exp; in bsplit() 57 _arb_exp_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in _arb_exp_sum_bs_simple() argument 62 bsplit(P, T, Q, Qexp, x, r, 0, N, 0); in _arb_exp_sum_bs_simple()
|
H A D | atan_sum_bs_simple.c | 15 bsplit(fmpz_t P, fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in bsplit() argument 28 *Qexp = 2 * r; in bsplit() 43 bsplit(P, T, Q, Qexp, x, r, a, m); in bsplit() 52 *Qexp = *Qexp + *Q2exp; in bsplit() 61 _arb_atan_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in _arb_atan_sum_bs_simple() argument 66 bsplit(P, T, Q, Qexp, x, r, 0, N); in _arb_atan_sum_bs_simple()
|
H A D | sin_cos_arf_bb.c | 18 bsplit(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in bsplit() argument 29 *Qexp = 2 * r + cc; in bsplit() 43 *Qexp = 2 * r + cc; in bsplit() 48 *Qexp += 2 * r + cc; in bsplit() 74 *Qexp = *Qexp + *Q2exp; in bsplit() 123 bsplit(T, Q, Qexp, xexp, xpow, r, 0, N); in _arb_sin_sum_bs_powtab() 162 flint_bitcnt_t Qexp[1]; in arb_sin_cos_fmpz_div_2exp_bsplit() local 177 _arb_sin_sum_bs_powtab(T, Q, Qexp, x, r, N); in arb_sin_cos_fmpz_div_2exp_bsplit() 181 Qexp[0] += r; in arb_sin_cos_fmpz_div_2exp_bsplit() 184 if (Qexp[0] >= prec) in arb_sin_cos_fmpz_div_2exp_bsplit() [all …]
|
H A D | atan_sum_bs_powtab.c | 19 bsplit(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in bsplit() argument 32 *Qexp = 2 * r; in bsplit() 45 *Qexp = 4 * r; in bsplit() 59 bsplit(T, Q, Qexp, xexp, xpow, r, a, m); in bsplit() 72 *Qexp = *Qexp + *Q2exp; in bsplit() 78 _arb_atan_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in _arb_atan_sum_bs_powtab() argument 120 bsplit(T, Q, Qexp, xexp, xpow, r, 0, N); in _arb_atan_sum_bs_powtab()
|
H A D | exp_sum_bs_powtab.c | 127 bsplit(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in bsplit() argument 137 *Qexp = r + cc; in bsplit() 149 *Qexp = r + cc; in bsplit() 153 *Qexp += r + cc; in bsplit() 167 bsplit(T, Q, Qexp, xexp, xpow, r, a, m); in bsplit() 179 *Qexp = *Qexp + *Q2exp; in bsplit() 185 _arb_exp_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, in _arb_exp_sum_bs_powtab() argument 227 bsplit(T, Q, Qexp, xexp, xpow, r, 0, N); in _arb_exp_sum_bs_powtab()
|
H A D | exp_arf_bb.c | 48 flint_bitcnt_t Qexp[1]; in arb_exp_arf_bb() local 123 _arb_exp_sum_bs_powtab(T, Q, Qexp, u, r, N); in arb_exp_arf_bb() 126 if (*Qexp >= wp) in arb_exp_arf_bb() 128 fmpz_tdiv_q_2exp(T, T, *Qexp - wp); in arb_exp_arf_bb() 133 fmpz_mul_2exp(T, T, wp - *Qexp); in arb_exp_arf_bb()
|
H A D | atan_arf_bb.c | 156 flint_bitcnt_t Qexp[1]; in arb_atan_arf_bb() local 269 _arb_atan_sum_bs_powtab(P, Q, Qexp, u, r, N); in arb_atan_arf_bb() 273 *Qexp += r; in arb_atan_arf_bb() 276 if (*Qexp >= wp) in arb_atan_arf_bb() 278 fmpz_tdiv_q_2exp(P, P, *Qexp - wp); in arb_atan_arf_bb() 283 fmpz_mul_2exp(P, P, wp - *Qexp); in arb_atan_arf_bb()
|
/dports/math/arb/arb-2.21.1/arb/test/ |
H A D | t-atan_sum_bs_powtab.c | 29 flint_bitcnt_t Qexp, r; in main() local 41 _arb_atan_sum_bs_simple(T, Q, &Qexp, x, r, N); in main() 43 fmpq_div_2exp(S, S, Qexp); in main() 45 _arb_atan_sum_bs_powtab(T, Q, &Qexp, x, r, N); in main() 47 fmpq_div_2exp(V, V, Qexp); in main()
|
H A D | t-exp_sum_bs_powtab.c | 29 flint_bitcnt_t Qexp, r; in main() local 41 _arb_exp_sum_bs_simple(T, Q, &Qexp, x, r, N); in main() 43 fmpq_div_2exp(S, S, Qexp); in main() 45 _arb_exp_sum_bs_powtab(T, Q, &Qexp, x, r, N); in main() 47 fmpq_div_2exp(V, V, Qexp); in main()
|
/dports/science/jdftx/jdftx-1.6.0/jdftx/fluid/ |
H A D | PCM_internal.h | 256 …raint(const ScalarField& muPlus, const ScalarField& muMinus, const ScalarField& shape, double Qexp, 263 double mu0 = -(Qexp + Qdiff)/Qsum; 265 double mu0_Qsum = (Qexp + Qdiff)/(Qsum*Qsum); 280 double disc = sqrt(Qexp*Qexp - 4.*Qplus*Qminus); //discriminant for quadratic 282 if(Qexp<0) 283 { mu0 = log((disc-Qexp)/(2.*Qplus)); 284 mu0_Qplus = -2.*Qminus/(disc*(disc-Qexp)) - 1./Qplus; 285 mu0_Qminus = -2.*Qplus/(disc*(disc-Qexp)); 288 { mu0 = log(-2.*Qminus/(disc+Qexp)); 289 mu0_Qplus = 2.*Qminus/(disc*(disc+Qexp)); [all …]
|
H A D | NonlinearPCM.cpp | 191 double mu0 = 0., Qexp = 0., Adiel_Qexp = 0.; in operator ()() local 196 Qexp = integral(rhoExplicitTilde); in operator ()() 197 mu0 = screeningEval->neutralityConstraint(muPlus, muMinus, shape.back(), Qexp); in operator ()() 237 …screeningEval->neutralityConstraint(muPlus, muMinus, shape.back(), Qexp, &mu0_muPlus, &mu0_muMinus… in operator ()() 340 double Qexp = integral(rhoExplicitTilde); in dumpDensities() local 341 double mu0 = screeningEval->neutralityConstraint(muPlus, muMinus, shape.back(), Qexp); in dumpDensities()
|
/dports/math/flint2/flint-2.8.4/nmod_mpoly/ |
H A D | divides_heap_threaded.c | 849 ulong * Qexp = Q->exps; in _nmod_mpoly_divides_stripe1() local 906 &Qexp, &Q->exps_alloc, 1, Alen + 1); in _nmod_mpoly_divides_stripe1() 970 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 995 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1029 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1045 Q->exps = Qexp; in _nmod_mpoly_divides_stripe1() 1053 Q->exps = Qexp; in _nmod_mpoly_divides_stripe1() 1076 ulong * Qexp = Q->exps; in _nmod_mpoly_divides_stripe() local 1136 &Qexp, &Q->exps_alloc, N, Qlen + 1); in _nmod_mpoly_divides_stripe() 1296 Q->exps = Qexp; in _nmod_mpoly_divides_stripe() [all …]
|
H A D | mpolyn_divides_threaded.c | 920 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe1() local 982 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1, S->ctx); in _nmod_mpolyn_divides_stripe1() 984 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], mask); in _nmod_mpolyn_divides_stripe1() 1045 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1070 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1110 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1132 *Q_exp = Qexp; in _nmod_mpolyn_divides_stripe1() 1159 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe() local 1232 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N, S->ctx); in _nmod_mpolyn_divides_stripe() 1234 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _nmod_mpolyn_divides_stripe() [all …]
|
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/nmod_mpoly/ |
H A D | divides_heap_threaded.c | 855 ulong * Qexp = * Q_exp; in _nmod_mpoly_divides_stripe1() local 911 _nmod_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1); in _nmod_mpoly_divides_stripe1() 913 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], mask); in _nmod_mpoly_divides_stripe1() 974 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 999 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1033 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1053 *Q_exp = Qexp; in _nmod_mpoly_divides_stripe1() 1080 ulong * Qexp = * Q_exp; in _nmod_mpoly_divides_stripe() local 1158 _nmod_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N); in _nmod_mpoly_divides_stripe() 1228 Qexp + N*x->j, N); in _nmod_mpoly_divides_stripe() [all …]
|
H A D | mpolyn_divides_threaded.c | 934 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe1() local 996 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1, S->ctx); in _nmod_mpolyn_divides_stripe1() 998 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], mask); in _nmod_mpolyn_divides_stripe1() 1059 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1084 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1124 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1146 *Q_exp = Qexp; in _nmod_mpolyn_divides_stripe1() 1173 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe() local 1246 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N, S->ctx); in _nmod_mpolyn_divides_stripe() 1248 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _nmod_mpolyn_divides_stripe() [all …]
|
/dports/math/e-antic/flint2-ae7ec89/nmod_mpoly/ |
H A D | divides_heap_threaded.c | 855 ulong * Qexp = * Q_exp; in _nmod_mpoly_divides_stripe1() local 911 _nmod_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1); in _nmod_mpoly_divides_stripe1() 913 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], mask); in _nmod_mpoly_divides_stripe1() 974 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 999 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1033 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpoly_divides_stripe1() 1053 *Q_exp = Qexp; in _nmod_mpoly_divides_stripe1() 1080 ulong * Qexp = * Q_exp; in _nmod_mpoly_divides_stripe() local 1158 _nmod_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N); in _nmod_mpoly_divides_stripe() 1228 Qexp + N*x->j, N); in _nmod_mpoly_divides_stripe() [all …]
|
H A D | mpolyn_divides_threaded.c | 934 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe1() local 996 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1, S->ctx); in _nmod_mpolyn_divides_stripe1() 998 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], mask); in _nmod_mpolyn_divides_stripe1() 1059 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1084 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1124 texp = Bexp[x->i] + Qexp[x->j]; in _nmod_mpolyn_divides_stripe1() 1146 *Q_exp = Qexp; in _nmod_mpolyn_divides_stripe1() 1173 ulong * Qexp = * Q_exp; in _nmod_mpolyn_divides_stripe() local 1246 _nmod_mpolyn_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N, S->ctx); in _nmod_mpolyn_divides_stripe() 1248 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _nmod_mpolyn_divides_stripe() [all …]
|
/dports/math/e-antic/flint2-ae7ec89/fmpz_mpoly/ |
H A D | divides_heap_threaded.c | 933 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe1() local 1015 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1); in _fmpz_mpoly_divides_stripe1() 1107 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1132 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1221 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1244 *Q_exp = Qexp; in _fmpz_mpoly_divides_stripe1() 1271 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe() local 1379 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N); in _fmpz_mpoly_divides_stripe() 1382 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() 1385 lt_divides = mpoly_monomial_divides_mp(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() [all …]
|
H A D | mpolyu_divides_threaded.c | 737 ulong * Qexp = * Q_exp; in _fmpz_mpolyuu_divides_stripe() local 793 _fmpz_mpolyu_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, S->minor_bits, S->ctx); in _fmpz_mpolyuu_divides_stripe() 795 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], S->main_overflow_mask); in _fmpz_mpolyuu_divides_stripe() 873 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 898 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 950 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 972 *Q_exp = Qexp; in _fmpz_mpolyuu_divides_stripe()
|
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/fmpz_mpoly/ |
H A D | divides_heap_threaded.c | 933 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe1() local 1015 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1); in _fmpz_mpoly_divides_stripe1() 1107 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1132 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1221 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1244 *Q_exp = Qexp; in _fmpz_mpoly_divides_stripe1() 1271 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe() local 1379 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N); in _fmpz_mpoly_divides_stripe() 1382 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() 1385 lt_divides = mpoly_monomial_divides_mp(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() [all …]
|
H A D | mpolyu_divides_threaded.c | 737 ulong * Qexp = * Q_exp; in _fmpz_mpolyuu_divides_stripe() local 793 _fmpz_mpolyu_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, S->minor_bits, S->ctx); in _fmpz_mpolyuu_divides_stripe() 795 lt_divides = mpoly_monomial_divides1(Qexp + Qlen, exp, Bexp[0], S->main_overflow_mask); in _fmpz_mpolyuu_divides_stripe() 873 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 898 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 950 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpolyuu_divides_stripe() 972 *Q_exp = Qexp; in _fmpz_mpolyuu_divides_stripe()
|
/dports/math/flint2/flint-2.8.4/fmpz_mpoly/ |
H A D | divides_heap_threaded.c | 933 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe1() local 1015 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, 1); in _fmpz_mpoly_divides_stripe1() 1107 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1132 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1221 texp = Bexp[x->i] + Qexp[x->j]; in _fmpz_mpoly_divides_stripe1() 1244 *Q_exp = Qexp; in _fmpz_mpoly_divides_stripe1() 1271 ulong * Qexp = * Q_exp; in _fmpz_mpoly_divides_stripe() local 1359 _fmpz_mpoly_fit_length(&Qcoeff, &Qexp, &Qalloc, Qlen + 1, N); in _fmpz_mpoly_divides_stripe() 1367 lt_divides = mpoly_monomial_divides(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() 1374 lt_divides = mpoly_monomial_divides_mp(Qexp + N*Qlen, exp, in _fmpz_mpoly_divides_stripe() [all …]
|
/dports/math/arb/arb-2.21.1/ |
H A D | arb.h | 929 void _arb_exp_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, 932 void _arb_exp_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, 1004 void _arb_atan_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, 1007 void _arb_atan_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp,
|
/dports/math/arb/arb-2.21.1/doc/source/ |
H A D | arb.rst | 1697 .. function:: void _arb_exp_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, const fmpz_t x… 1699 .. function:: void _arb_exp_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, const fmpz_t x… 1701 Computes *T*, *Q* and *Qexp* such that 1702 `T / (Q 2^{\text{Qexp}}) = \sum_{k=1}^N (x/2^r)^k/k!` using binary splitting. 1714 .. function:: void _arb_atan_sum_bs_simple(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, const fmpz_t … 1716 .. function:: void _arb_atan_sum_bs_powtab(fmpz_t T, fmpz_t Q, flint_bitcnt_t * Qexp, const fmpz_t … 1718 Computes *T*, *Q* and *Qexp* such that 1719 `T / (Q 2^{\text{Qexp}}) = \sum_{k=1}^N (-1)^k (x/2^r)^{2k} / (2k+1)`
|