[div.c] fixed bug in div(x,x,x) when real part x is zero

git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@657 211d60ee-9f03-0410-a15a-8952a2c7a4e4

1 files changed, 2 insertions, 1 deletions

diff --git a/src/div.c b/src/div.c
index b788ecd..1e346ef 100644
--- a/src/div.c
+++ b/src/div.c
@@ -206,6 +206,7 @@ mpc_div (mpc_ptr a, mpc_srcptr b, mpc_srcptr c, mpc_rnd_t rnd)
    {
       /* (re_b+i*im_b)/(i*c) = im_b/c - i * (re_b/c) */
       int overlap = (a == b) || (a == c);
+      int imag_b = mpfr_zero_p (MPC_RE (b));
       mpfr_t cloc;
 
       if (overlap)
@@ -226,7 +227,7 @@ mpc_div (mpc_ptr a, mpc_srcptr b, mpc_srcptr c, mpc_rnd_t rnd)
       if (mpfr_zero_p (MPC_RE (a)))
          mpfr_setsign (MPC_RE (a), MPC_RE (a), (brs != crs && bis != cis),
             GMP_RNDN); /* exact */
-      if (mpfr_zero_p (MPC_RE (b)))
+      if (imag_b)
          mpfr_setsign (MPC_IM (a), MPC_IM (a), (bis != crs && brs == cis),
             GMP_RNDN);