math_self_complex: Add test case

diff --git a/math/t_self_complex.c b/math/t_self_complex.c
index 4399ec1..7fed4ac 100644 (file)
--- a/math/t_self_complex.c
+++ b/math/t_self_complex.c
@@ -2,8 +2,8 @@
 
 #include <atf-c.h>
 #include <math.h>
-#include <stdio.h>
 
+#include "config.h"
 #include "subr_atf.h"
 #include "subr_fpcmp.h"
 #include "subr_random.h"
@@ -44,10 +44,49 @@ ATF_TC_BODY(test_self_complex1, tc)
        }
 }
 
+ATF_TC(test_self_complex2);
+ATF_TC_HEAD(test_self_complex2, tc)
+{
+       atf_tc_set_md_var(tc,
+           "descr",
+           "Check that our random generators work fine "
+           "for complex numbers");
+}
+ATF_TC_BODY(test_self_complex2, tc)
+{
+       float complex fcx;
+       double complex dcx;
+       long double complex ldcx;
+       long i, N;
+
+       N = atf_tc_get_config_var_as_long(tc, "iterations");
+       ATF_REQUIRE(N > 0);
+
+       ATF_FOR_LOOP(i, N, i++) {
+               /* float */
+               fcx = random_float_complex(FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(crealf(fcx)) == FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(cimagf(fcx)) == FP_NORMAL);
+
+               /* double */
+               dcx = random_double_complex(FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(creal(dcx)) == FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(cimag(dcx)) == FP_NORMAL);
+
+               /* long double */
+#if defined(HAVE_CREALL) && defined(HAVE_CIMAGL)
+               ldcx = random_long_double_complex(FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(creall(ldcx)) == FP_NORMAL);
+               ATF_PASS_OR_BREAK(fpclassify(cimagl(ldcx)) == FP_NORMAL);
+#endif
+       }
+}
+
 /* Add test cases to test program */
 ATF_TP_ADD_TCS(tp)
 {
        ATF_TP_ADD_TC(tp, test_self_complex1);
+       ATF_TP_ADD_TC(tp, test_self_complex2);
 
        return atf_no_error();
 }