diff options
author | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2010-07-30 15:35:43 +0000 |
---|---|---|
committer | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2010-07-30 15:35:43 +0000 |
commit | 77a6d6e0f8dd1b3cabe40f7a3b90cc2423a4b6d8 (patch) | |
tree | c613656dd08c8a0aaa9d463f23cf1f98ecc0623d /src/mul.c | |
parent | f5e44f7b814899e0c8a01cf4803e99c37de9fecb (diff) | |
download | mpc-77a6d6e0f8dd1b3cabe40f7a3b90cc2423a4b6d8.tar.gz |
code cleanup:
removed macro MPFR_PREC from mpc-impl.h
replaced MPFR_PREC (MPC_RE (op)) by MPC_PREC_RE (op) and
MPFR_PREC (MPC_IM (op)) by MPC_PREC_IM (op)
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@814 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'src/mul.c')
-rw-r--r-- | src/mul.c | 17 |
1 files changed, 8 insertions, 9 deletions
@@ -1,6 +1,6 @@ /* mpc_mul -- Multiply two complex numbers -Copyright (C) 2002, 2004, 2005, 2008, 2009 Andreas Enge, Paul Zimmermann, Philippe Th\'eveny +Copyright (C) 2002, 2004, 2005, 2008, 2009, 2010 Andreas Enge, Paul Zimmermann, Philippe Th\'eveny This file is part of the MPC Library. @@ -222,7 +222,7 @@ mul_pure_imaginary (mpc_ptr a, mpc_srcptr u, mpfr_srcptr y, mpc_rnd_t rnd, mpc_t result; if (overlap) - mpc_init3 (result, MPFR_PREC (MPC_RE (a)), MPFR_PREC (MPC_IM (a))); + mpc_init3 (result, MPC_PREC_RE (a), MPC_PREC_IM (a)); else result [0] = a [0]; @@ -259,7 +259,7 @@ mpc_mul_naive (mpc_ptr a, mpc_srcptr b, mpc_srcptr c, mpc_rnd_t rnd) if (overlap) { - mpfr_init2 (t, MPFR_PREC(MPC_RE(a))); + mpfr_init2 (t, MPC_PREC_RE(a)); inex_re = mpfr_sub (t, u, v, MPC_RND_RE(rnd)); } else @@ -309,8 +309,7 @@ mpc_mul_karatsuba (mpc_ptr rop, mpc_srcptr op1, mpc_srcptr op2, mpc_rnd_t rnd) overlap = (rop == op1) || (rop == op2); if (overlap) - mpc_init3 (result, MPFR_PREC (MPC_RE (rop)), - MPFR_PREC (MPC_IM (rop))); + mpc_init3 (result, MPC_PREC_RE (rop), MPC_PREC_IM (rop)); else result [0] = rop [0]; @@ -342,12 +341,12 @@ mpc_mul_karatsuba (mpc_ptr rop, mpc_srcptr op1, mpc_srcptr op2, mpc_rnd_t rnd) /* find the precision and rounding mode for the new real part */ if (mul_i % 2) { - prec_re = MPFR_PREC(MPC_IM(rop)); + prec_re = MPC_PREC_IM(rop); rnd_re = MPC_RND_IM(rnd); } else /* mul_i = 0 or 2 */ { - prec_re = MPFR_PREC(MPC_RE(rop)); + prec_re = MPC_PREC_RE(rop); rnd_re = MPC_RND_RE(rnd); } @@ -358,8 +357,8 @@ mpc_mul_karatsuba (mpc_ptr rop, mpc_srcptr op1, mpc_srcptr op2, mpc_rnd_t rnd) prec = MPC_MAX_PREC(rop); mpfr_init2 (u, 2); - mpfr_init2 (v, prec_v = MPFR_PREC(a) + MPFR_PREC(d)); - mpfr_init2 (w, prec_w = MPFR_PREC(b) + MPFR_PREC(c)); + mpfr_init2 (v, prec_v = mpfr_get_prec (a) + mpfr_get_prec (d)); + mpfr_init2 (w, prec_w = mpfr_get_prec (b) + mpfr_get_prec (c)); mpfr_init2 (x, 2); mpfr_mul (v, a, d, GMP_RNDN); /* exact */ |