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

increase working precision in dilog

Browse source
This commit is contained in:
Fredrik Johansson 2017-02-26 00:08:30 +01:00
parent 3b1e0a2e01
commit 882c0ec62d
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
acb_hypgeom/dilog.c
View file

@ -16,7 +16,7 @@ acb_hypgeom_dilog(acb_t res, const acb_t z, slong prec)
{
double a, b, best, mz, mz1, t, u;
int algorithm;
slong acc;
slong acc, inprec;
if (!acb_is_finite(z))
{
@ -34,6 +34,7 @@ acb_hypgeom_dilog(acb_t res, const acb_t z, slong prec)
acc = FLINT_MAX(acc, 0);
acc = FLINT_MIN(acc, prec);
prec = FLINT_MIN(prec, acc + 30);
inprec = prec;
/* first take care of exponents that may overflow doubles */
if (arf_cmpabs_2exp_si(arb_midref(acb_realref(z)), -20) <= 0 &&
@ -50,6 +51,8 @@ acb_hypgeom_dilog(acb_t res, const acb_t z, slong prec)
return;
}
prec = 1.005 * prec + 5;
a = arf_get_d(arb_midref(acb_realref(z)), ARF_RND_DOWN);
b = arf_get_d(arb_midref(acb_imagref(z)), ARF_RND_DOWN);
@ -149,5 +152,7 @@ acb_hypgeom_dilog(acb_t res, const acb_t z, slong prec)
acb_hypgeom_dilog_transform(res, z, algorithm, prec);
else /* (algorithm == 8) */
acb_hypgeom_dilog_bernoulli(res, z, prec);
acb_set_round(res, res, inprec);
}