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

Airy functions: fix precision loss for huge input and high precision

Browse source
This commit is contained in:
fredrik 2018-05-01 11:02:49 +02:00
parent e0c823ab52
commit 08345b400c
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
acb_hypgeom/airy.c
View file

@ -227,7 +227,7 @@ acb_hypgeom_airy(acb_t ai, acb_t aip, acb_t bi, acb_t bip, const acb_t z, slong
zmag = fmpz_get_d(ARF_EXPREF(re));
else
zmag = fmpz_get_d(ARF_EXPREF(im));
zmag = (zmag + 1) * (1.0 / LOG2);
zmag = 3.0 * zmag + 1;
n = wp / (-zmag) + 1;
}
@ -246,7 +246,7 @@ acb_hypgeom_airy(acb_t ai, acb_t aip, acb_t bi, acb_t bip, const acb_t z, slong
{
x = fmpz_get_d(ARF_EXPREF(re));
y = fmpz_get_d(ARF_EXPREF(im));
zmag = (FLINT_MAX(x, y) - 2) * (1.0 / LOG2);
zmag = (FLINT_MAX(x, y) - 2) * LOG2;
n = asymp_pick_terms(wp, zmag);
n = FLINT_MAX(n, 1);
}