diff options
author | zimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2011-09-08 10:18:04 +0000 |
---|---|---|
committer | zimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2011-09-08 10:18:04 +0000 |
commit | d622cf4115f155cfaab532ca6929dd1f647b2b24 (patch) | |
tree | f6940b368240c2867a625c884ed969ff8e0e878a /tests/tpow_fr.c | |
parent | e30b69480b4585d565926fe5ae82ae4437a41686 (diff) | |
download | mpc-d622cf4115f155cfaab532ca6929dd1f647b2b24.tar.gz |
[tests/tpow_fr.c] added test case exposing reuse bug (I don't know how to enter
such a test in pow_fr.dat)
[src/pow.c] fixed corresponding bug
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1088 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'tests/tpow_fr.c')
-rw-r--r-- | tests/tpow_fr.c | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/tests/tpow_fr.c b/tests/tpow_fr.c index bccc093..d9968a8 100644 --- a/tests/tpow_fr.c +++ b/tests/tpow_fr.c @@ -20,12 +20,40 @@ along with this program. If not, see http://www.gnu.org/licenses/ . #include "mpc-tests.h" +static void +test_reuse (void) +{ + mpc_t z; + mpfr_t y; + int inex; + + mpfr_init2 (y, 2); + mpc_init2 (z, 2); + mpc_set_si_si (z, 0, -1, MPC_RNDNN); + mpfr_neg (mpc_realref (z), mpc_realref (z), GMP_RNDN); + mpc_div_2exp (z, z, 4, MPC_RNDNN); + mpfr_set_ui (y, 512, GMP_RNDN); + inex = mpc_pow_fr (z, z, y, MPC_RNDNN); + if (MPC_INEX_RE(inex) != 0 || MPC_INEX_IM(inex) != 0 || + mpfr_cmp_ui_2exp (MPC_RE(z), 1, -2048) != 0 || + mpfr_cmp_ui (MPC_IM(z), 0) != 0 || mpfr_signbit (MPC_IM(z)) == 0) + { + printf ("Error in test_reuse, wrong ternary value or output\n"); + printf ("inex=(%d %d)\n", MPC_INEX_RE(inex), MPC_INEX_IM(inex)); + printf ("z="); mpc_out_str (stdout, 2, 0, z, MPC_RNDNN); printf ("\n"); + exit (1); + } + mpfr_clear (y); + mpc_clear (z); +} + int main (void) { DECL_FUNC (CCF, f, mpc_pow_fr); test_start (); + test_reuse (); data_check (f, "pow_fr.dat"); tgeneric (f, 2, 1024, 7, 10); |