Home
last modified time | relevance | path

Searched refs:term_prec (Results 1 – 14 of 14) sorted by relevance

/dports/math/arb/arb-2.21.1/acb_hypgeom/
H A Dgamma_stirling_sum_improved.c31 slong n, term_mag, term_prec; in acb_hypgeom_gamma_stirling_sum_improved() local
232 term_prec = FLINT_MAX(term_prec, 10); in acb_hypgeom_gamma_stirling_sum_improved()
303 term_prec = prec + term_mag; in acb_hypgeom_gamma_stirling_sum_improved()
304 term_prec = FLINT_MIN(term_prec, prec); in acb_hypgeom_gamma_stirling_sum_improved()
305 term_prec = FLINT_MAX(term_prec, 10); in acb_hypgeom_gamma_stirling_sum_improved()
337 term_prec = FLINT_MIN(term_prec, prec); in acb_hypgeom_gamma_stirling_sum_improved()
338 term_prec = FLINT_MAX(term_prec, 10); in acb_hypgeom_gamma_stirling_sum_improved()
363 term_prec = FLINT_MIN(term_prec, prec); in acb_hypgeom_gamma_stirling_sum_improved()
364 term_prec = FLINT_MAX(term_prec, 10); in acb_hypgeom_gamma_stirling_sum_improved()
388 term_prec = FLINT_MIN(term_prec, prec); in acb_hypgeom_gamma_stirling_sum_improved()
[all …]
H A Dgamma_stirling_sum_horner.c22 slong n, term_mag, term_prec; in acb_hypgeom_gamma_stirling_sum_horner() local
50 term_prec = prec + term_mag; in acb_hypgeom_gamma_stirling_sum_horner()
51 term_prec = FLINT_MIN(term_prec, prec); in acb_hypgeom_gamma_stirling_sum_horner()
52 term_prec = FLINT_MAX(term_prec, 10); in acb_hypgeom_gamma_stirling_sum_horner()
54 if (prec - term_prec > 200) in acb_hypgeom_gamma_stirling_sum_horner()
56 acb_set_round(t, w, term_prec); in acb_hypgeom_gamma_stirling_sum_horner()
57 acb_mul(s, s, t, term_prec); in acb_hypgeom_gamma_stirling_sum_horner()
60 acb_mul(s, s, w, term_prec); in acb_hypgeom_gamma_stirling_sum_horner()
62 arb_gamma_stirling_coeff(b, n, 0, term_prec); in acb_hypgeom_gamma_stirling_sum_horner()
63 acb_add(s, s, b, term_prec); in acb_hypgeom_gamma_stirling_sum_horner()
H A Dgamma_taylor.c16 evaluate_rect(acb_t res, const short * term_prec, slong len, const acb_t x, slong prec) in evaluate_rect() argument
57 acb_set_round(t, xs + m, term_prec[n1]); in evaluate_rect()
58 acb_mul(res, res, t, term_prec[n1]); in evaluate_rect()
59 acb_add(res, res, s, term_prec[n1]); in evaluate_rect()
61 acb_mul(res, res, xs + m, term_prec[n1]); in evaluate_rect()
62 acb_add(res, res, s, term_prec[n1]); in evaluate_rect()
130 short term_prec[ARB_HYPGEOM_GAMMA_TAB_NUM]; in acb_hypgeom_gamma_taylor() local
203 term_prec[0] = wp; in acb_hypgeom_gamma_taylor()
216 term_prec[i] = FLINT_MIN(FLINT_MAX(wp + tail_bound, 2), wp); in acb_hypgeom_gamma_taylor()
218 if (term_prec[i] > arb_hypgeom_gamma_coeffs[i].nlimbs * FLINT_BITS) in acb_hypgeom_gamma_taylor()
[all …]
/dports/math/arb/arb-2.21.1/arb_hypgeom/
H A Dgamma_stirling_sum_improved.c20 slong n, term_mag, term_prec; in arb_hypgeom_gamma_stirling_sum_improved() local
221 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_gamma_stirling_sum_improved()
292 term_prec = prec + term_mag; in arb_hypgeom_gamma_stirling_sum_improved()
293 term_prec = FLINT_MIN(term_prec, prec); in arb_hypgeom_gamma_stirling_sum_improved()
294 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_gamma_stirling_sum_improved()
326 term_prec = FLINT_MIN(term_prec, prec); in arb_hypgeom_gamma_stirling_sum_improved()
327 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_gamma_stirling_sum_improved()
352 term_prec = FLINT_MIN(term_prec, prec); in arb_hypgeom_gamma_stirling_sum_improved()
353 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_gamma_stirling_sum_improved()
377 term_prec = FLINT_MIN(term_prec, prec); in arb_hypgeom_gamma_stirling_sum_improved()
[all …]
H A Dgamma_stirling_sum_horner.c21 slong n, term_mag, term_prec; in arb_hypgeom_gamma_stirling_sum_horner() local
49 term_prec = prec + term_mag; in arb_hypgeom_gamma_stirling_sum_horner()
50 term_prec = FLINT_MIN(term_prec, prec); in arb_hypgeom_gamma_stirling_sum_horner()
51 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_gamma_stirling_sum_horner()
53 if (prec - term_prec > 200) in arb_hypgeom_gamma_stirling_sum_horner()
55 arb_set_round(t, w, term_prec); in arb_hypgeom_gamma_stirling_sum_horner()
56 arb_mul(s, s, t, term_prec); in arb_hypgeom_gamma_stirling_sum_horner()
59 arb_mul(s, s, w, term_prec); in arb_hypgeom_gamma_stirling_sum_horner()
61 arb_gamma_stirling_coeff(b, n, 0, term_prec); in arb_hypgeom_gamma_stirling_sum_horner()
62 arb_add(s, s, b, term_prec); in arb_hypgeom_gamma_stirling_sum_horner()
H A Dcentral_bin_ui.c29 slong j, k, term_prec, wp; in arb_hypgeom_central_bin_ui_asymp() local
65 term_prec = wp + term_mag; in arb_hypgeom_central_bin_ui_asymp()
66 term_prec = FLINT_MIN(term_prec, wp); in arb_hypgeom_central_bin_ui_asymp()
67 term_prec = FLINT_MAX(term_prec, 10); in arb_hypgeom_central_bin_ui_asymp()
69 arb_gamma_stirling_coeff(t, j + 1, 0, term_prec); in arb_hypgeom_central_bin_ui_asymp()
72 arb_sub(t, u, t, term_prec); in arb_hypgeom_central_bin_ui_asymp()
H A Dgamma_taylor.c251 short term_prec[ARB_HYPGEOM_GAMMA_TAB_NUM]; in arb_hypgeom_gamma_taylor() local
530 term_prec[0] = wp; in arb_hypgeom_gamma_taylor()
543 term_prec[i] = FLINT_MIN(FLINT_MAX(wp + tail_bound, 2), wp); in arb_hypgeom_gamma_taylor()
560 for (i = n - 1; i >= 1 && term_prec[i] <= 53; i--) in arb_hypgeom_gamma_taylor()
587 prec, du, log2u, i, term_prec[i]); in arb_hypgeom_gamma_taylor()
591 if (term_prec[i] < wp - 128) in arb_hypgeom_gamma_taylor()
593 arf_set_round(v, u, term_prec[i], ARF_RND_DOWN); in arb_hypgeom_gamma_taylor()
594 arf_mul(s, s, v, term_prec[i], ARF_RND_DOWN); in arb_hypgeom_gamma_taylor()
595 arf_add(s, s, c, term_prec[i], ARF_RND_DOWN); in arb_hypgeom_gamma_taylor()
599 arf_mul(s, s, u, term_prec[i], ARF_RND_DOWN); in arb_hypgeom_gamma_taylor()
[all …]
/dports/math/arb/arb-2.21.1/acb/
H A Dgamma.c32 slong k, term_prec; in acb_gamma_stirling_eval() local
57 term_prec = prec + term_mag; in acb_gamma_stirling_eval()
58 term_prec = FLINT_MIN(term_prec, prec); in acb_gamma_stirling_eval()
59 term_prec = FLINT_MAX(term_prec, 10); in acb_gamma_stirling_eval()
61 arb_gamma_stirling_coeff(b, k, digamma, term_prec); in acb_gamma_stirling_eval()
65 acb_set_round(t, zinv2, term_prec); in acb_gamma_stirling_eval()
66 acb_mul(s, s, t, term_prec); in acb_gamma_stirling_eval()
69 acb_mul(s, s, zinv2, term_prec); in acb_gamma_stirling_eval()
71 arb_add(acb_realref(s), acb_realref(s), b, term_prec); in acb_gamma_stirling_eval()
/dports/math/arb/arb-2.21.1/acb_modular/
H A Dtheta_const_sum_rs.c87 slong k, term_prec, i, e, eprev; in acb_modular_theta_const_sum_rs() local
138 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in acb_modular_theta_const_sum_rs()
139 … _acb_modular_mul(qpow + k, tmp1, tmp2, qpow + tab[k], qpow + k - tab[k], term_prec, prec); in acb_modular_theta_const_sum_rs()
145 term_prec = prec; in acb_modular_theta_const_sum_rs()
153 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in acb_modular_theta_const_sum_rs()
159 _acb_modular_mul(theta2, tmp1, tmp2, theta2, qpow + m2, term_prec, prec); in acb_modular_theta_const_sum_rs()
170 term_prec = prec; in acb_modular_theta_const_sum_rs()
178 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in acb_modular_theta_const_sum_rs()
184 _acb_modular_mul(theta3, tmp1, tmp2, theta3, qpow + m3, term_prec, prec); in acb_modular_theta_const_sum_rs()
187 _acb_modular_mul(theta4, tmp1, tmp2, theta4, qpow + m3, term_prec, prec); in acb_modular_theta_const_sum_rs()
H A Deta_sum.c65 slong e, e1, e2, k, k1, k2, num, term_prec; in _acb_modular_eta_sum_basecase() local
115 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in _acb_modular_eta_sum_basecase()
126 _acb_modular_mul(qpow + k, tmp1, tmp2, qpow + k1, qpow + k2, term_prec, prec); in _acb_modular_eta_sum_basecase()
130 _acb_modular_mul(qpow + k, tmp1, tmp2, qpow + k1, qpow + k1, term_prec, prec); in _acb_modular_eta_sum_basecase()
131 _acb_modular_mul(qpow + k, tmp1, tmp2, qpow + k, qpow + k2, term_prec, prec); in _acb_modular_eta_sum_basecase()
158 slong k, term_prec, i, e, eprev; in _acb_modular_eta_sum_rs() local
195 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in _acb_modular_eta_sum_rs()
196 … _acb_modular_mul(qpow + k, tmp1, tmp2, qpow + tab[k], qpow + k - tab[k], term_prec, prec); in _acb_modular_eta_sum_rs()
202 term_prec = prec; in _acb_modular_eta_sum_rs()
210 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in _acb_modular_eta_sum_rs()
[all …]
H A Dtheta_sum.c24 slong e, e1, e2, k, k1, k2, r, n, N, WN, term_prec; in acb_modular_theta_sum() local
231 qpow + k1, qpow + k2, term_prec, prec); in acb_modular_theta_sum()
236 qpow + k1, qpow + k1, term_prec, prec); in acb_modular_theta_sum()
238 qpow + k, qpow + k2, term_prec, prec); in acb_modular_theta_sum()
279 acb_add(tmp1, wpow + n, vpow + n, term_prec); in acb_modular_theta_sum()
280 acb_mul(tmp1, qpow + k, tmp1, term_prec); in acb_modular_theta_sum()
285 acb_sub(tmp2, wpow + n, vpow + n, term_prec); in acb_modular_theta_sum()
286 acb_mul(tmp2, qpow + k, tmp2, term_prec); in acb_modular_theta_sum()
296 acb_mul_ui(term, term, 2 * n, term_prec); in acb_modular_theta_sum()
321 acb_mul(tmp1, qpow + k, tmp1, term_prec); in acb_modular_theta_sum()
[all …]
H A Dtheta_const_sum_basecase.c19 slong k, term_prec; in acb_modular_theta_const_sum_basecase() local
141 term_prec = FLINT_MIN(FLINT_MAX(prec + log2term_approx + 16.0, 16.0), prec); in acb_modular_theta_const_sum_basecase()
142 _acb_modular_mul(qpow + k, t1, t2, qpow + tab[k], qpow + k - tab[k], term_prec, prec); in acb_modular_theta_const_sum_basecase()
/dports/math/arb/arb-2.21.1/arb/
H A Dgamma.c188 slong k, term_prec; in arb_gamma_stirling_eval() local
214 term_prec = prec + term_mag; in arb_gamma_stirling_eval()
215 term_prec = FLINT_MIN(term_prec, prec); in arb_gamma_stirling_eval()
216 term_prec = FLINT_MAX(term_prec, 10); in arb_gamma_stirling_eval()
220 arb_set_round(t, zinv2, term_prec); in arb_gamma_stirling_eval()
221 arb_mul(s, s, t, term_prec); in arb_gamma_stirling_eval()
224 arb_mul(s, s, zinv2, term_prec); in arb_gamma_stirling_eval()
226 arb_gamma_stirling_coeff(b, k, digamma, term_prec); in arb_gamma_stirling_eval()
227 arb_add(s, s, b, term_prec); in arb_gamma_stirling_eval()
H A Ddot.c237 slong shift_bits, shift_limbs, term_prec; in _arb_dot_addmul_generic() local
256 term_prec = sn * FLINT_BITS - shift; in _arb_dot_addmul_generic()
261 term_prec += FLINT_BITS; in _arb_dot_addmul_generic()
262 nn = (term_prec + FLINT_BITS - 1) / FLINT_BITS; in _arb_dot_addmul_generic()
272 if (term_prec >= MUL_MPFR_MIN_LIMBS * FLINT_BITS && in _arb_dot_addmul_generic()
274 xn * FLINT_BITS > 0.9 * term_prec && in _arb_dot_addmul_generic()
275 yn * FLINT_BITS > 0.9 * term_prec) in _arb_dot_addmul_generic()
363 slong shift_bits, shift_limbs, term_prec; in _arb_dot_add_generic() local
371 term_prec = sn * FLINT_BITS - shift; in _arb_dot_add_generic()
372 term_prec += FLINT_BITS; in _arb_dot_add_generic()
[all …]