diff options
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); |