diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2003-06-13 08:22:53 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2003-06-13 08:22:53 +0000 |
commit | 26e1a5b3a26dff669a0cf1f63880d97c2c195e0b (patch) | |
tree | 4be9251a084fe36ce39a06de75ce69aaf7b90c29 /erf.c | |
parent | 79e3cdfa0402607afe8dbaab26c120f31ecf5e3e (diff) | |
download | mpfr-26e1a5b3a26dff669a0cf1f63880d97c2c195e0b.tar.gz |
improved error analysis on erf
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@2332 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'erf.c')
-rw-r--r-- | erf.c | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -175,9 +175,9 @@ mpfr_erf_0 (mpfr_ptr res, mpfr_srcptr x, mp_rnd_t rnd_mode) if ((nuk < - (mp_exp_t) m) && ((double) k >= xf * xf)) break; - /* tauk <- 1/2 + tauk * 2^sigmak + 2^(2k+3+nuk) */ + /* tauk <- 1/2 + tauk * 2^sigmak + (1+8k)*2^nuk */ tauk = 0.5 + mul_2exp (tauk, sigmak) - + mul_2exp (1.0, 2 * k + 3 + nuk); + + mul_2exp (1.0 + 8.0 * (double) k, nuk); } mpfr_mul (s, x, s, GMP_RNDU); |