/dports/math/arb/arb-2.21.1/acb_hypgeom/ |
H A D | gamma_stirling_sum_improved.c | 31 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 D | gamma_stirling_sum_horner.c | 22 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 D | gamma_taylor.c | 16 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 D | gamma_stirling_sum_improved.c | 20 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 D | gamma_stirling_sum_horner.c | 21 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 D | central_bin_ui.c | 29 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 D | gamma_taylor.c | 251 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 D | gamma.c | 32 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 D | theta_const_sum_rs.c | 87 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 D | eta_sum.c | 65 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 D | theta_sum.c | 24 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 D | theta_const_sum_basecase.c | 19 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 D | gamma.c | 188 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 D | dot.c | 237 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 …]
|