summaryrefslogtreecommitdiff
path: root/erf.c
diff options
context:
space:
mode:
authorzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2003-06-13 08:22:53 +0000
committerzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2003-06-13 08:22:53 +0000
commit26e1a5b3a26dff669a0cf1f63880d97c2c195e0b (patch)
tree4be9251a084fe36ce39a06de75ce69aaf7b90c29 /erf.c
parent79e3cdfa0402607afe8dbaab26c120f31ecf5e3e (diff)
downloadmpfr-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.c4
1 files changed, 2 insertions, 2 deletions
diff --git a/erf.c b/erf.c
index f435f8274..b13ef47ec 100644
--- a/erf.c
+++ b/erf.c
@@ -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);