diff options
Diffstat (limited to 'tests/tacos.c')
| -rw-r--r-- | tests/tacos.c | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/tests/tacos.c b/tests/tacos.c index 6ad1cc6..2f764b1 100644 --- a/tests/tacos.c +++ b/tests/tacos.c @@ -28,11 +28,46 @@ along with this program. If not, see http://www.gnu.org/licenses/ . #include "data_check.tpl" #include "tgeneric.tpl" +/* test with reduced exponent range */ +static void +test_reduced (void) +{ + mpfr_exp_t emin = mpfr_get_emin (); + mpfr_exp_t emax = mpfr_get_emax (); + mpc_t x, z, zr; + + mpfr_set_emin (-148); + mpfr_set_emax (128); + mpc_init2 (x, 24); + mpc_init2 (z, 24); + mpc_init2 (zr, 24); + mpfr_set_flt (mpc_realref (x), -0x0.01f28c10p0f, MPFR_RNDN); + mpfr_set_flt (mpc_imagref (x), -0x6.79cdd0p-68f, MPFR_RNDN); + mpc_acos (z, x, MPC_RNDNN); + /* square root is 0.00038652126908433 - 1.37634353022868e12*I */ + mpfr_set_flt (mpc_realref (zr), 0x1.941242p0f, MPFR_RNDN); + mpfr_set_flt (mpc_imagref (zr), 0x6.79da18p-68f, MPFR_RNDN); + if (mpc_cmp (z, zr)) + { + printf ("Incorrect acos with reduced exponent range:\n"); + mpfr_printf ("Expected (%Re,%Re)\n", mpc_realref (zr), mpc_imagref (zr)); + mpfr_printf ("Got (%Re,%Re)\n", mpc_realref (z), mpc_imagref (z)); + exit (1); + } + mpc_clear (x); + mpc_clear (z); + mpc_clear (zr); + mpfr_set_emin (emin); + mpfr_set_emax (emax); +} + int main (void) { test_start (); + test_reduced (); + data_check_template ("acos.dsc", "acos.dat"); tgeneric_template ("acos.dsc", 2, 512, 7, 7); |