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

fix bug in fpwrap

Browse source
This commit is contained in:
fredrik 2021-12-10 13:36:41 +01:00
parent 2aa02b9171
commit 75df786325
2 changed files with 21 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
arb_fpwrap/fpwrap.c
View file

@ -590,7 +590,7 @@ int arb_fpwrap_double_4_int(double * res, arb_func_4_int func, double x1, double
arb_set_d(arb_x1, x1);
arb_set_d(arb_x2, x2);
arb_set_d(arb_x3, x3);
arb_set_d(arb_x3, x4);
arb_set_d(arb_x4, x4);
if (!arb_is_finite(arb_x1) || !arb_is_finite(arb_x2) || !arb_is_finite(arb_x3) || !arb_is_finite(arb_x4))
{

20
arb_fpwrap/test/t-fpwrap.c
View file

@ -148,6 +148,26 @@ int main()
}
}
{
double y;
arb_fpwrap_double_hypgeom_2f1(&y, 0.1, 0.2, 0.3, 0.4, 0, FPWRAP_CORRECT_ROUNDING);
if (fabs(y - 1.0341794015503748492) > 1e-16)
{
flint_printf("FAIL: value 6\n\n");
flint_abort();
}
arb_fpwrap_double_hypgeom_2f1(&y, 0.1, 0.2, 0.3, 0.4, 1, FPWRAP_CORRECT_ROUNDING);
if (fabs(y - 0.34569799520143110351) > 1e-16)
{
flint_printf("FAIL: value 6\n\n");
flint_abort();
}
}
{
double x, y, z;
complex_double cx, cy, cz, ctau;