hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-04 17:01:40 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

rel accuracy version of an internal function

Browse source
This commit is contained in:
fredrik 2022-01-13 11:12:26 +01:00
parent f6d1112d25
commit f7849aa1d8
1 changed files with 33 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

40
arb_hypgeom/gamma_upper_fmpq.c
View file

@ -44,7 +44,7 @@ mag_pow_fmpq_fast(mag_t res, const mag_t x, const fmpq_t e)
}
slong
_arb_hypgeom_gamma_upper_fmpq_inf_choose_N(mag_t err, const fmpq_t a, const arb_t z, const mag_t abs_tol)
_arb_hypgeom_gamma_upper_fmpq_inf_choose_N_1(mag_t err, const fmpq_t a, const arb_t z, int prefactor, const mag_t abs_tol)
{
slong N, aa, ab;
fmpz_t az1, az2;
@ -71,12 +71,19 @@ _arb_hypgeom_gamma_upper_fmpq_inf_choose_N(mag_t err, const fmpq_t a, const arb_
ab = fmpz_get_si(az2);
/* prefactor z^(a-1) * exp(-z) */
arb_get_mag_lower(t, z);
mag_expinv(t, t);
fmpq_sub_ui(a1, a, 1);
arb_get_mag(u, z);
mag_pow_fmpq_fast(u, u, a1);
mag_mul(err, t, u);
if (prefactor)
{
arb_get_mag_lower(t, z);
mag_expinv(t, t);
fmpq_sub_ui(a1, a, 1);
arb_get_mag(u, z);
mag_pow_fmpq_fast(u, u, a1);
mag_mul(err, t, u);
}
else
{
mag_one(err);
}
if (mag_is_inf(err))
{
@ -121,6 +128,24 @@ _arb_hypgeom_gamma_upper_fmpq_inf_choose_N(mag_t err, const fmpq_t a, const arb_
return N;
}
slong
_arb_hypgeom_gamma_upper_fmpq_inf_choose_N(mag_t err, const fmpq_t a, const arb_t z, const mag_t abs_tol)
{
return _arb_hypgeom_gamma_upper_fmpq_inf_choose_N_1(err, a, z, 1, abs_tol);
}
slong
_arb_hypgeom_gamma_upper_fmpq_inf_choose_N_rel(mag_t err, const fmpq_t a, const arb_t z, slong prec)
{
mag_t tol;
slong N;
mag_init(tol);
mag_set_ui_2exp_si(tol, 1, -prec);
N = _arb_hypgeom_gamma_upper_fmpq_inf_choose_N_1(err, a, z, 0, tol);
mag_clear(tol);
return N;
}
static void
upper_bsplit(arb_t M, arb_t S, arb_t Q, const fmpz_t ap, const fmpz_t aq, const arb_t z, slong na, slong nb, int cont, slong prec)
{
@ -282,6 +307,7 @@ _arb_hypgeom_gamma_lower_fmpq_0_choose_N(mag_t err, const fmpq_t a, const arb_t
return N;
}
/* todo: squaring opt; power table */
static void
lower_bsplit(arb_t M, arb_t S, arb_t Q, const fmpz_t ap, const fmpz_t aq, const arb_t z, slong na, slong nb, int cont, slong prec)