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 (not elegant) for fragile argument reduction in incomplete elliptic integrals

Browse source
This commit is contained in:
fredrik 2018-10-24 22:36:30 +02:00
parent 512ada9674
commit eb60d7ac4e
3 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
acb_elliptic/e_inc.c
View file

@ -142,7 +142,7 @@ acb_elliptic_e_inc(acb_t res, const acb_t phi, const acb_t m, int times_pi, slon
acb_zero(z);
arf_set_mag(arb_midref(acb_realref(z)), arb_radref(d));
mag_zero(arb_radref(d));
arb_sub(d, d, acb_realref(z), prec);
arb_sub(d, d, acb_realref(z), 2 * prec + 100); /* meant to be exact */
arb_floor(d, d, prec);
/* w = 2 E(m) */

2
acb_elliptic/f.c
View file

@ -127,7 +127,7 @@ acb_elliptic_f(acb_t res, const acb_t phi, const acb_t m, int times_pi, slong pr
acb_zero(z);
arf_set_mag(arb_midref(acb_realref(z)), arb_radref(d));
mag_zero(arb_radref(d));
arb_sub(d, d, acb_realref(z), prec);
arb_sub(d, d, acb_realref(z), 2 * prec + 100); /* meant to be exact */
arb_floor(d, d, prec);
/* w = 2 K(m) */

2
acb_elliptic/pi.c
View file

@ -182,7 +182,7 @@ acb_elliptic_pi_inc(acb_t res, const acb_t n, const acb_t phi, const acb_t m, in
acb_zero(z);
arf_set_mag(arb_midref(acb_realref(z)), arb_radref(d));
mag_zero(arb_radref(d));
arb_sub(d, d, acb_realref(z), prec);
arb_sub(d, d, acb_realref(z), 2 * prec + 100); /* meant to be exact */
arb_floor(d, d, prec);
/* w = 2 Pi(n, m) */