summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorzimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4>2009-11-20 20:00:50 +0000
committerzimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4>2009-11-20 20:00:50 +0000
commite50f474738de94f86c172f4e8fe56c4610bbcca2 (patch)
treed929974e91999c3336f64ca540986aa445efb518
parentc1bce31b9e8704cac4bc857ab13512f45c8e5bc8 (diff)
downloadmpc-e50f474738de94f86c172f4e8fe56c4610bbcca2.tar.gz
[atan.c] fixed problem with variable y which was initialized twice
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@725 211d60ee-9f03-0410-a15a-8952a2c7a4e4
-rw-r--r--src/atan.c5
1 files changed, 2 insertions, 3 deletions
diff --git a/src/atan.c b/src/atan.c
index f99be1c..8f6ab79 100644
--- a/src/atan.c
+++ b/src/atan.c
@@ -198,7 +198,7 @@ mpc_atan (mpc_ptr rop, mpc_srcptr op, mpc_rnd_t rnd)
mp_exp_t op_im_exp;
mp_rnd_t rnd1, rnd2;
- mpfr_inits (a, b, x, y, (mpfr_ptr)0);
+ mpfr_inits (a, b, x, y, (mpfr_ptr) 0);
/* real part: Re(arctan(x+i*y)) = [arctan2(x,1-y) - arctan2(-x,1+y)]/2 */
minus_op_re[0] = MPC_RE (op)[0];
@@ -286,7 +286,6 @@ mpc_atan (mpc_ptr rop, mpc_srcptr op, mpc_rnd_t rnd)
/* Imaginary part
Im(atan(x+I*y)) = 1/4 * [log(x^2+(1+y)^2) - log (x^2 +(1-y)^2)] */
prec = mpfr_get_prec (MPC_IM (rop)); /* result precision */
- mpfr_init (y);
/* a = o(1+y) error(a) < 1 ulp(a)
b = o(a^2) error(b) < 5 ulp(b)
@@ -352,7 +351,7 @@ mpc_atan (mpc_ptr rop, mpc_srcptr op, mpc_rnd_t rnd)
inex = mpc_set_fr_fr (rop, x, y, rnd);
- mpfr_clears (a, b, x, y, (mpfr_ptr)0);
+ mpfr_clears (a, b, x, y, (mpfr_ptr) 0);
return inex;
}
}