diff options
author | Philippe Theveny <philippe.theveny@laposte.net> | 2008-05-20 15:56:15 +0000 |
---|---|---|
committer | Philippe Theveny <philippe.theveny@laposte.net> | 2008-05-20 15:56:15 +0000 |
commit | 8d4114fc6c70eb8591ee8e78fd5027f2e045b0ab (patch) | |
tree | bf86f7e7aeb766c08f6f2814a9aaeb8013e2a319 /tests/tcos.c | |
parent | 9dc5b73974ba3d973b3fcd70a3352986c493fee7 (diff) | |
download | mpc-git-8d4114fc6c70eb8591ee8e78fd5027f2e045b0ab.tar.gz |
Fix wrong sign.
git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/mpc/trunk@134 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'tests/tcos.c')
-rw-r--r-- | tests/tcos.c | 37 |
1 files changed, 19 insertions, 18 deletions
diff --git a/tests/tcos.c b/tests/tcos.c index 9d47a69..16ce677 100644 --- a/tests/tcos.c +++ b/tests/tcos.c @@ -294,8 +294,8 @@ special () static void pure_real_argument () { - /* cos(x -i*0) = cos(x) -i*0*sin(x) */ - /* sin(x +i*0) = cos(x) +i*0*sin(x) */ + /* cos(x -i*0) = cos(x) +i*0*sin(x) */ + /* cos(x +i*0) = cos(x) -i*0*sin(x) */ mpc_t u; mpc_t z; mpc_t cos_z; @@ -304,59 +304,60 @@ pure_real_argument () mpc_init (z); mpc_init (cos_z); - /* cos(1 +i*0) = cos(1) +i*0*sin(1) */ + /* cos(1 +i*0) = cos(1) -i*0*sin(1) */ mpc_set_ui_ui (z, 1, 0, MPC_RNDNN); mpfr_cos (MPC_RE (u), MPC_RE (z), GMP_RNDN); mpfr_set_ui (MPC_IM (u), 0, GMP_RNDN); + mpc_conj (u, u, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(1 -i*0) = cos(1) -i*0*sin(1) */ + /* cos(1 -i*0) = cos(1) +i*0*sin(1) */ mpc_conj (z, z, MPC_RNDNN); mpc_conj (u, u, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(-1 +i*0) = cos(-1) +i*0*sin(-1) */ + /* cos(-1 +i*0) = cos(-1) -i*0*sin(-1) */ mpc_neg (z, z, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(-1 -i*0) = cos(-1) -i*0*sin(-1) */ + /* cos(-1 -i*0) = cos(-1) +i*0*sin(-1) */ mpc_conj (z, z, MPC_RNDNN); mpc_conj (u, u, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(2 +i*0) = cos(2) +i*0*sin(2) */ + /* cos(2 +i*0) = cos(2) -i*0*sin(2) */ mpc_set_ui_ui (z, 2, 0, MPC_RNDNN); mpfr_cos (MPC_RE (u), MPC_RE (z), GMP_RNDN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(2 -i*0) = cos(2) -i*0*sin(2) */ + /* cos(2 -i*0) = cos(2) +i*0*sin(2) */ mpc_conj (z, z, MPC_RNDNN); mpc_conj (u, u, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(-2 +i*0) = cos(-2) +i*0*sin(-2) */ + /* cos(-2 +i*0) = cos(-2) -i*0*sin(-2) */ mpc_neg (z, z, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); - /* cos(-2 -i*0) = cos(-2) -i*0*sin(-2) */ + /* cos(-2 -i*0) = cos(-2) +i*0*sin(-2) */ mpc_conj (z, z, MPC_RNDNN); mpc_conj (u, u, MPC_RNDNN); mpc_cos (cos_z, z, MPC_RNDNN); - if (mpc_cmp (cos_z, u) != 0 || mpfr_signbit (MPC_IM (cos_z))) + if (mpc_cmp (cos_z, u) != 0 || !mpfr_signbit (MPC_IM (cos_z))) test_failed (z, cos_z, u); mpc_clear (cos_z); |