diff options
author | Paul Zimmermann <Paul.Zimmermann@inria.fr> | 2018-06-20 11:11:52 +0200 |
---|---|---|
committer | Paul Zimmermann <Paul.Zimmermann@inria.fr> | 2018-06-20 11:11:52 +0200 |
commit | 389c43dfbf14f5200ad88170d0e14f90df180066 (patch) | |
tree | 7bc70a78f6b67659fccdbbc90b0f7e2c49ffb53e | |
parent | 54b8d5f09a3cf6f5d88d616118996b20413c7acd (diff) | |
parent | e9438ddfc68d8952c02d2a1ed53d13a1d8960b3f (diff) | |
download | mpc-git-389c43dfbf14f5200ad88170d0e14f90df180066.tar.gz |
Merge branch 'master' of git+ssh://scm.gforge.inria.fr//gitroot/mpc/mpc
-rw-r--r-- | src/dot.c | 2 | ||||
-rw-r--r-- | tests/tdot.c | 12 |
2 files changed, 11 insertions, 3 deletions
@@ -53,7 +53,7 @@ mpc_dot (mpc_ptr res, const mpc_ptr *x, const mpc_ptr *y, mpfr_mul (z[i], mpc_realref (x[i]), mpc_realref (y[i]), MPFR_RNDZ); /* idem for z[n+i]: we allocate with prec_x_im + prec_y_max bits */ mpfr_init2 (z[n+i], prec_x_im + prec_y_max); - mpfr_set_prec (z[i], prec_x_im + prec_y_im); + mpfr_set_prec (z[n+i], prec_x_im + prec_y_im); mpfr_mul (z[n+i], mpc_imagref (x[i]), mpc_imagref (y[i]), MPFR_RNDZ); mpfr_neg (z[n+i], z[n+i], MPFR_RNDZ); } diff --git a/tests/tdot.c b/tests/tdot.c index b5a61d6..c6d1f39 100644 --- a/tests/tdot.c +++ b/tests/tdot.c @@ -38,21 +38,29 @@ check_special (void) /* dot product of empty vectors is 0 */ inex = mpc_dot (res, t, t, 0, MPC_RNDNN); MPC_ASSERT (inex == 0); - MPC_ASSERT (mpfr_cmp_ui (mpc_realref (res), 0) == 0); - MPC_ASSERT (mpfr_cmp_ui (mpc_imagref (res), 0) == 0); + MPC_ASSERT (mpfr_zero_p (mpc_realref (res))); + MPC_ASSERT (mpfr_zero_p (mpc_imagref (res))); + MPC_ASSERT (mpfr_signbit (mpc_realref (res)) == 0); + MPC_ASSERT (mpfr_signbit (mpc_imagref (res)) == 0); /* (1,2)*(1,2) = (-3,4) */ inex = mpc_dot (res, t, t, 1, MPC_RNDNN); MPC_ASSERT (inex == 0); + MPC_ASSERT (mpfr_regular_p (mpc_realref (res))); + MPC_ASSERT (mpfr_regular_p (mpc_imagref (res))); MPC_ASSERT (mpfr_cmp_si (mpc_realref (res), -3) == 0); MPC_ASSERT (mpfr_cmp_ui (mpc_imagref (res), 4) == 0); /* (1,2)*(1,2) + (2,3)*(2,3) = (-8,16) */ inex = mpc_dot (res, t, t, 2, MPC_RNDNN); MPC_ASSERT (inex == 0); + MPC_ASSERT (mpfr_regular_p (mpc_realref (res))); + MPC_ASSERT (mpfr_regular_p (mpc_imagref (res))); MPC_ASSERT (mpfr_cmp_si (mpc_realref (res), -8) == 0); MPC_ASSERT (mpfr_cmp_ui (mpc_imagref (res), 16) == 0); /* (1,2)*(1,2) + (2,3)*(2,3) + (3,4)*(3,4) = (-15,40) */ inex = mpc_dot (res, t, t, 3, MPC_RNDNN); MPC_ASSERT (inex == 0); + MPC_ASSERT (mpfr_regular_p (mpc_realref (res))); + MPC_ASSERT (mpfr_regular_p (mpc_imagref (res))); MPC_ASSERT (mpfr_cmp_si (mpc_realref (res), -15) == 0); MPC_ASSERT (mpfr_cmp_ui (mpc_imagref (res), 40) == 0); for (i = 0; i < 3; i++) |