diff options
author | thevenyp <thevenyp@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2008-04-21 15:29:37 +0000 |
---|---|---|
committer | thevenyp <thevenyp@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2008-04-21 15:29:37 +0000 |
commit | ab49659dc541afdf1431a66346e29b547d2f6bbd (patch) | |
tree | f2fa5889eb7a4592f706638b9739ae3ca0a02af8 /tests/tsin.c | |
parent | c81321dced4bb67b00ab1257a9a8c393379103a5 (diff) | |
download | mpc-ab49659dc541afdf1431a66346e29b547d2f6bbd.tar.gz |
Add special values processing in mpc_sin as in the C99 standard
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@115 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'tests/tsin.c')
-rw-r--r-- | tests/tsin.c | 228 |
1 files changed, 227 insertions, 1 deletions
diff --git a/tests/tsin.c b/tests/tsin.c index ff3ce84..22098c8 100644 --- a/tests/tsin.c +++ b/tests/tsin.c @@ -28,6 +28,231 @@ MA 02111-1307, USA. */ #define TEST_FUNCTION mpc_sin #include "tgeneric.c" +static void +test_failed (mpc_t op, mpc_t get, mpc_t expected) +{ + printf ("mpc_sin(op) failed with "); + OUT (op); + printf (" "); + OUT (get); + OUT (expected); + exit (1); +} + +/* check special values as defined in C99 standard */ +static void +special () +{ + mpc_t z; + mpc_t s; + mpc_t c99; + + mpc_init (z); + mpc_init (s); + mpc_init (c99); + + mpfr_set_nan (MPC_RE (z)); + /* sin(NaN +i*NaN) = NaN +i*NaN */ + mpfr_set_nan (MPC_IM (z)); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(NaN +i*infinity) = NaN +/-i*infinity */ + mpfr_set_inf (MPC_IM (z), +1); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(NaN -i*infinity) = NaN -/+i*infinity */ + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(NaN +i*y) = NaN +i*NaN where 0<|y|<infinity */ + mpfr_set_ui (MPC_IM (z), 1, GMP_RNDN); + mpfr_set_nan (MPC_RE (c99)); + mpfr_set_nan (MPC_IM (c99)); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(NaN +i*0) = NaN +/-i*0 */ + mpfr_set_ui (MPC_IM (z), 0, GMP_RNDN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, z); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(+0 +i*NaN) = +0 +i*NaN */ + mpc_mul_i (z, z, +1, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_zero_p (MPC_RE (s)) || mpfr_signbit (MPC_RE (s)) != 0 + || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(-0 +i*NaN) = -0 +i*NaN */ + mpc_neg (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_zero_p (MPC_RE (s)) || mpfr_signbit (MPC_RE (s)) == 0 + || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, z); + + /* sin(x +i*NaN) = NaN +i*NaN where x!=0 */ + mpfr_set_inf (MPC_RE (z), +1); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_neg (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpfr_set_ui (MPC_RE (z), +1, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_neg (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) && !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + + mpfr_set_inf (MPC_RE(z), -1); + /* sin(-infinity -i*infinity) = NaN -/+i*infinity */ + mpfr_set_inf (MPC_IM(z), -1); + mpfr_set_nan (MPC_RE(c99)); + mpfr_set_inf (MPC_IM(c99), -1); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(-infinity +i*infinity) = NaN +/-i*infinity */ + mpc_conj (z, z, MPC_RNDNN); + mpc_conj (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(+infinity -i*infinity) = NaN -/+i*infinity */ + mpc_neg (z, z, MPC_RNDNN); + mpc_neg (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(+infinity +i*infinity) = NaN +/-i*infinity */ + mpc_conj (z, z, MPC_RNDNN); + mpc_conj (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_inf_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(x +/-i*infinity) = infinity*(sin x -i*cos x) where 0<|x|<infinity */ + mpfr_set_ui (MPC_RE (z), 1, GMP_RNDN); + mpfr_set_inf (MPC_RE (c99), +1); + mpfr_set_inf (MPC_IM (c99), +1); + mpc_sin (s, z, MPC_RNDUU); + if (mpc_cmp (s, c99) != 0) + test_failed (z, s, c99); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDUU); + if (mpc_cmp (s, c99) != 0) + test_failed (z, s, c99); + mpfr_neg (MPC_RE (z), MPC_RE (z), GMP_RNDN); + mpfr_neg (MPC_RE (c99), MPC_RE (c99), GMP_RNDN); + mpc_sin (s, z, MPC_RNDUU); + if (mpc_cmp (s, c99) != 0) + test_failed (z, s, c99); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDUU); + if (mpc_cmp (s, c99) != 0) + test_failed (z, s, c99); + + /* sin(+0 -i*infinity) = +0 -i*infinity */ + mpfr_set_ui (MPC_RE (z), 0, GMP_RNDN); + mpc_sin (s, z, MPC_RNDNN); + if (mpc_cmp (s, z) != 0) + test_failed (z, s, z); + + /* sin(+0 +i*infinity) = +0 +i*infinity */ + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (mpc_cmp (s, z) != 0) + test_failed (z, s, z); + + /* sin(-0 -i*infinity) = -0 -i*infinity */ + mpc_neg (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (mpc_cmp (s, z) != 0) + test_failed (z, s, z); + + /* sin(-0 +i*infinity) = -0 +i*infinity */ + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (mpc_cmp (s, z) != 0) + test_failed (z, s, z); + + /* sin(-infinity +i*0) = NaN +/-i*0 */ + mpfr_set_inf (MPC_RE (z), -1); + mpfr_set_ui (MPC_IM (z), 0, GMP_RNDN); + mpfr_set_nan (MPC_RE (c99)); + mpfr_set_ui (MPC_IM (c99), 0, GMP_RNDN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(-infinity -i*0) = NaN -/+i*0 */ + mpc_conj (z, z, MPC_RNDNN); + mpc_conj (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(+infinity +i*0) = NaN +/-i*0 */ + mpc_neg (z, z, MPC_RNDNN); + mpc_neg (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(+infinity -i*0) = NaN -/+i*0 */ + mpc_conj (z, z, MPC_RNDNN); + mpc_conj (c99, c99, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_zero_p (MPC_IM (s))) + test_failed (z, s, c99); + + /* sin(+/-infinity +i*y) = NaN +i*NaN where 0<|y|<infinity*/ + mpfr_set_ui (MPC_IM (z), 1, GMP_RNDN); + mpfr_set_nan (MPC_IM (c99)); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_neg (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + mpc_conj (z, z, MPC_RNDNN); + mpc_sin (s, z, MPC_RNDNN); + if (!mpfr_nan_p (MPC_RE (s)) || !mpfr_nan_p (MPC_IM (s))) + test_failed (z, s, c99); + + mpc_clear (c99); + mpc_clear (s); + mpc_clear (z); +} + int main() { @@ -39,7 +264,7 @@ main() mpc_init (z); mpfr_init (g); - for (prec = 2; prec <= 1000; prec++) + for (prec = 2; prec <= 1000; prec+=2) { mpc_set_prec (x, prec); mpc_set_prec (z, prec); @@ -65,6 +290,7 @@ main() } } + special (); tgeneric (); mpc_clear (x); |