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

fix issue #404: nan fp computation resulting in unreasonably large N in erf on some systems

Browse source
This commit is contained in:
fredrik 2022-04-02 13:35:38 +02:00
parent 076249eaee
commit 7de55e6b43
1 changed files with 2 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
arb_hypgeom/erf.c
View file

@ -97,6 +97,8 @@ arb_hypgeom_erf_1f1b(arb_t res, const arb_t z, slong prec)
dz = arf_get_d(arb_midref(z), ARF_RND_DOWN); dz = arf_get_d(arb_midref(z), ARF_RND_DOWN);
dz = fabs(dz); dz = fabs(dz);
u = -dz * dz + prec * LOG2 + log(dz); u = -dz * dz + prec * LOG2 + log(dz);
if (dz < 1.0)
u = FLINT_MAX(u, 1e-6);
u = u / d_lambertw(u / (EXP1 * dz * dz)); u = u / d_lambertw(u / (EXP1 * dz * dz));
N = u + 1; N = u + 1;
} }