diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2006-12-19 22:32:09 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2006-12-19 22:32:09 +0000 |
commit | 718e0f2bfd39efeb0d6ed578300f4cf49c8394f1 (patch) | |
tree | bcb7044db3c683f924c66f4e67ba1f35bd138c34 /erfc.c | |
parent | e939a4de0d5f82064ab734b7fb753554f32c151e (diff) | |
download | mpfr-718e0f2bfd39efeb0d6ed578300f4cf49c8394f1.tar.gz |
Replaced mpfr_mul_2exp and mpfr_div_2exp by mpfr_mul_2ui and
mpfr_div_2ui respectively (when this makes sense, of course).
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4317 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'erfc.c')
-rw-r--r-- | erfc.c | 6 |
1 files changed, 3 insertions, 3 deletions
@@ -42,7 +42,7 @@ mpfr_erfc_asympt (mpfr_ptr y, mpfr_srcptr x) with a bound for err */ mpfr_mul (xx, x, x, GMP_RNDD); /* err <= 1 */ mpfr_ui_div (xx, 1, xx, GMP_RNDU); /* upper bound for 1/(2x^2), err <= 2 */ - mpfr_div_2exp (xx, xx, 1, GMP_RNDU); /* exact */ + mpfr_div_2ui (xx, xx, 1, GMP_RNDU); /* exact */ mpfr_set_ui (t, 1, GMP_RNDN); /* current term, exact */ mpfr_set (y, t, GMP_RNDN); /* current sum */ mpfr_set_ui (err, 0, GMP_RNDN); @@ -70,9 +70,9 @@ mpfr_erfc_asympt (mpfr_ptr y, mpfr_srcptr x) } /* the error on y is bounded by err*ulp(y) */ mpfr_mul (t, x, x, GMP_RNDU); /* rel. err <= 2^(1-p) */ - mpfr_div_2exp (err, err, 3, GMP_RNDU); /* err/8 */ + mpfr_div_2ui (err, err, 3, GMP_RNDU); /* err/8 */ mpfr_add (err, err, t, GMP_RNDU); /* err/8 + xx */ - mpfr_mul_2exp (err, err, 3, GMP_RNDU); /* err + 8*xx */ + mpfr_mul_2ui (err, err, 3, GMP_RNDU); /* err + 8*xx */ mpfr_exp (t, t, GMP_RNDU); /* err <= 1/2*ulp(t) + err(x*x)*t <= 1/2*ulp(t)+2*|x*x|*ulp(t) <= (2*|x*x|+1/2)*ulp(t) */ |