Home
last modified time | relevance | path

Searched refs:Qexp (Results 1 – 24 of 24) sorted by relevance

/dports/math/arb/arb-2.21.1/arb/
H A Dexp_sum_bs_simple.c15 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 Datan_sum_bs_simple.c15 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 Dsin_cos_arf_bb.c18 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 Datan_sum_bs_powtab.c19 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 Dexp_sum_bs_powtab.c127 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 Dexp_arf_bb.c48 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 Datan_arf_bb.c156 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 Dt-atan_sum_bs_powtab.c29 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 Dt-exp_sum_bs_powtab.c29 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 DPCM_internal.h256 …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 DNonlinearPCM.cpp191 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 Ddivides_heap_threaded.c849 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 Dmpolyn_divides_threaded.c920 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 Ddivides_heap_threaded.c855 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 Dmpolyn_divides_threaded.c934 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 Ddivides_heap_threaded.c855 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 Dmpolyn_divides_threaded.c934 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 Ddivides_heap_threaded.c933 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 Dmpolyu_divides_threaded.c737 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 Ddivides_heap_threaded.c933 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 Dmpolyu_divides_threaded.c737 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 Ddivides_heap_threaded.c933 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 Darb.h929 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 Darb.rst1697 .. 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)`