1 /* Copyright (C) 1997, 1998, 1999 Free Software Foundation, Inc.
2 This file is part of the GNU C Library.
3 Contributed by Andreas Jaeger <aj@arthur.rhein-neckar.de>, 1997.
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Library General Public License as
7 published by the Free Software Foundation; either version 2 of the
8 License, or (at your option) any later version.
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Library General Public License for more details.
15 You should have received a copy of the GNU Library General Public
16 License along with the GNU C Library; see the file COPYING.LIB. If not,
17 write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
18 Boston, MA 02111-1307, USA. */
20 /* Part of testsuite for libm.
22 This file is processed by a perl script. The resulting file has to
23 be included by a master file that defines:
26 FUNC(function): converts general function name (like cos) to
27 name with correct suffix (e.g. cosl or cosf)
28 MATHCONST(x): like FUNC but for constants (e.g convert 0.0 to 0.0L)
29 FLOAT: floating point type to test
30 - TEST_MSG: informal message to be displayed
31 CHOOSE(Clongdouble,Cdouble,Cfloat,Cinlinelongdouble,Cinlinedouble,Cinlinefloat):
32 chooses one of the parameters as delta for testing
34 PRINTF_EXPR Floating point conversion specification to print a variable
35 of type FLOAT with printf. PRINTF_EXPR just contains
36 the specifier, not the percent and width arguments,
38 PRINTF_XEXPR Like PRINTF_EXPR, but print in hexadecimal format.
39 PRINTF_NEXPR Like PRINTF_EXPR, but print nice. */
41 /* This testsuite has currently tests for:
42 acos, acosh, asin, asinh, atan, atan2, atanh,
43 cbrt, ceil, copysign, cos, cosh, erf, erfc, exp, exp10, exp2, expm1,
44 fabs, fdim, floor, fma, fmax, fmin, fmod, fpclassify,
46 ilogb, isfinite, isinf, isnan, isnormal,
47 isless, islessequal, isgreater, isgreaterequal, islessgreater, isunordered,
49 ldexp, lgamma, log, log10, log1p, log2, logb,
50 modf, nearbyint, nextafter,
51 pow, remainder, remquo, rint, lrint, llrint,
52 round, lround, llround,
53 scalb, scalbn, scalbln, signbit, sin, sincos, sinh, sqrt, tan, tanh, tgamma, trunc,
56 and for the following complex math functions:
57 cabs, cacos, cacosh, carg, casin, casinh, catan, catanh,
58 ccos, ccosh, cexp, clog, cpow, cproj, csin, csinh, csqrt, ctan, ctanh.
60 At the moment the following functions aren't tested:
61 conj, cimag, creal, drem,
65 Parameter handling is primitive in the moment:
66 --verbose=[0..3] for different levels of output:
68 1: basic report on failed tests (default)
69 2: full report on all tests
70 -v for full output (equals --verbose=3)
71 -u for generation of an ULPs file
76 This suite tests some aspects of the correct implementation of
77 mathematical functions in libm. Some simple, specific parameters
78 are tested for correctness but there's no exhaustive
79 testing. Handling of specific inputs (e.g. infinity, not-a-number)
80 is also tested. Correct handling of exceptions is checked
81 against. These implemented tests should check all cases that are
82 specified in ISO C 9X.
84 Exception testing: At the moment only divide-by-zero and invalid
85 exceptions are tested. Overflow/underflow and inexact exceptions
86 aren't checked at the moment.
88 NaN values: There exist signalling and quiet NaNs. This implementation
89 only uses signalling NaN as parameter but does not differenciate
90 between the two kinds of NaNs as result.
92 Inline functions: Inlining functions should give an improvement in
93 speed - but not in precission. The inlined functions return
94 reasonable values for a reasonable range of input values. The
95 result is not necessarily correct for all values and exceptions are
96 not correctly raised in all cases. Problematic input and return
97 values are infinity, not-a-number and minus zero. This suite
98 therefore does not check these specific inputs and the exception
99 handling for inlined mathematical functions - just the "reasonable"
102 Beware: The tests might fail for any of the following reasons:
104 - Functions are wrong
105 - Floating Point Unit not working properly
106 - Compiler has errors
108 With e.g. gcc 2.7.2.2 the test for cexp fails because of a compiler error.
115 #include "libm-test-ulps.h"
126 /* Possible exceptions */
127 #define NO_EXCEPTION 0x0
128 #define INVALID_EXCEPTION 0x1
129 #define DIVIDE_BY_ZERO_EXCEPTION 0x2
130 /* The next flags signals that those exceptions are allowed but not required. */
131 #define INVALID_EXCEPTION_OK 0x4
132 #define DIVIDE_BY_ZERO_EXCEPTION_OK 0x8
133 #define EXCEPTIONS_OK INVALID_EXCEPTION_OK+DIVIDE_BY_ZERO_EXCEPTION_OK
134 /* Some special test flags, passed togther with exceptions. */
135 #define IGNORE_ZERO_INF_SIGN 0x10
137 /* Various constants (we must supply them precalculated for accuracy). */
138 #define M_PI_6l .52359877559829887308L
139 #define M_E2l 7.389056098930650227230L
140 #define M_E3l 20.08553692318766774093L
141 #define M_2_SQRT_PIl 3.5449077018110320545963L /* 2 sqrt (M_PIl) */
142 #define M_SQRT_PIl 1.77245385090551602729817L /* sqrt (M_PIl) */
143 #define M_LOG_SQRT_PIl 0.572364942924700087072L /* log(sqrt(M_PIl)) */
144 #define M_LOG_2_SQRT_PIl 1.265512123484645396489L /* log(2*sqrt(M_PIl)) */
145 #define M_PI_34l (M_PIl - M_PI_4l) /* 3*pi/4 */
146 #define M_PI_34_LOG10El (M_PIl - M_PI_4l) * M_LOG10El
147 #define M_PI2_LOG10El M_PI_2l * M_LOG10El
148 #define M_PI4_LOG10El M_PI_4l * M_LOG10El
149 #define M_PI_LOG10El M_PIl * M_LOG10El
151 static FILE *ulps_file; /* File to document difference. */
152 static int output_ulps; /* Should ulps printed? */
154 static int noErrors; /* number of errors */
155 static int noTests; /* number of tests (without testing exceptions) */
156 static int noExcTests; /* number of tests for exception flags */
157 static int noXFails; /* number of expected failures. */
158 static int noXPasses; /* number of unexpected passes. */
161 static int output_max_error; /* Should the maximal errors printed? */
162 static int output_points; /* Should the single function results printed? */
164 static FLOAT minus_zero, plus_zero;
165 static FLOAT plus_infty, minus_infty, nan_value;
167 static FLOAT max_error;
170 #define BUILD_COMPLEX(real, imag) \
171 ({ __complex__ FLOAT __retval; \
172 __real__ __retval = (real); \
173 __imag__ __retval = (imag); \
177 #define MANT_DIG CHOOSE ((LDBL_MANT_DIG-1), (DBL_MANT_DIG-1), (FLT_MANT_DIG-1), \
178 (LDBL_MANT_DIG-1), (DBL_MANT_DIG-1), (FLT_MANT_DIG-1))
181 init_max_error (void)
187 set_max_error (FLOAT current)
189 if (current > max_error)
194 /* Should the message print to screen? This depends on the verbose flag,
195 and the test status. */
197 print_screen (int ok, int xfail)
201 || ((verbose == 1) && (ok == xfail))))
207 /* Should the message print to screen? This depends on the verbose flag,
208 and the test status. */
210 print_screen_max_error (int ok, int xfail)
214 || ((verbose == 1) && (ok == xfail))))
219 /* Update statistic counters. */
221 update_stats (int ok, int xfail)
226 else if (!ok && xfail)
228 else if (!ok && !xfail)
233 print_ulps (const char *test_name, FLOAT ulp)
237 fprintf (ulps_file, "Test \"%s\":\n", test_name);
238 fprintf (ulps_file, "%s: % .4" PRINTF_NEXPR "\n",
239 CHOOSE("ldouble", "double", "float",
240 "ildouble", "idouble", "ifloat"), ulp);
245 print_function_ulps (const char *function_name, FLOAT ulp)
249 fprintf (ulps_file, "Function \"%s\":\n", function_name);
250 fprintf (ulps_file, "%s: % .4" PRINTF_NEXPR "\n",
251 CHOOSE("ldouble", "double", "float",
252 "ildouble", "idouble", "ifloat"), ulp);
258 /* Test if Floating-Point stack hasn't changed */
260 fpstack_test (const char *test_name)
263 static int old_stack;
266 asm ("fnstsw" : "=a" (sw));
272 printf ("FP-Stack wrong after test %s (%d, should be %d)\n",
273 test_name, sw, old_stack);
282 print_max_error (const char *func_name, FLOAT allowed, int xfail)
287 sprintf (str, "Maximal error of `%s':", func_name);
289 if (max_error <= allowed)
295 print_function_ulps (func_name, max_error);
298 if (print_screen_max_error (ok, xfail))
300 printf ("%s\n", str);
301 printf (" is : % .4" PRINTF_NEXPR " ulp\n", max_error);
302 printf (" accepted: % .4" PRINTF_NEXPR " ulp\n", allowed);
305 update_stats (ok, xfail);
309 /* Test whether a given exception was raised. */
311 test_single_exception (const char *test_name,
315 const char *flag_name)
319 if (exception & exc_flag)
321 if (fetestexcept (fe_flag))
323 if (print_screen (1, 0))
324 printf ("Pass: %s: Exception \"%s\" set\n", test_name, flag_name);
329 if (print_screen (0, 0))
330 printf ("Failure: %s: Exception \"%s\" not set\n",
331 test_name, flag_name);
336 if (fetestexcept (fe_flag))
339 if (print_screen (0, 0))
340 printf ("Failure: %s: Exception \"%s\" set\n",
341 test_name, flag_name);
345 if (print_screen (1, 0))
346 printf ("%s: Exception \"%s\" not set\n", test_name,
357 /* Test whether exceptions given by EXCEPTION are raised. Ignore thereby
358 allowed but not required exceptions.
361 test_exceptions (const char *test_name, int exception)
365 if ((exception & DIVIDE_BY_ZERO_EXCEPTION_OK) == 0)
366 test_single_exception (test_name, exception,
367 DIVIDE_BY_ZERO_EXCEPTION, FE_DIVBYZERO,
371 if ((exception & INVALID_EXCEPTION_OK) == 0)
372 test_single_exception (test_name, exception, INVALID_EXCEPTION, FE_INVALID,
373 "Invalid operation");
375 feclearexcept (FE_ALL_EXCEPT);
380 check_float (const char *test_name, FLOAT computed, FLOAT expected,
381 FLOAT max_ulp, int xfail, int exceptions)
388 test_exceptions (test_name, exceptions);
389 if (isnan (computed) && isnan (expected))
391 else if (isinf (computed) && isinf (expected))
393 /* Test for sign of infinities. */
394 if (((exceptions & IGNORE_ZERO_INF_SIGN) == 0)
395 && (isinf (computed) != isinf (expected)))
398 printf ("infinity has wrong sign.\n");
403 /* Don't calc ulp for NaNs or infinities. */
404 else if (isinf (computed) || isnan (computed) || isinf (expected) || isnan (expected))
408 diff = FUNC(fabs) (computed - expected);
409 /* ilogb (0) isn't allowed. */
411 ulp = diff / FUNC(ldexp) (1.0, - MANT_DIG);
413 ulp = diff / FUNC(ldexp) (1.0, FUNC(ilogb) (expected) - MANT_DIG);
416 if (((exceptions & IGNORE_ZERO_INF_SIGN) == 0)
417 && (computed == 0.0 && expected == 0.0
418 && signbit(computed) != signbit (expected)))
421 print_ulps (test_name, ulp);
423 else if (ulp <= max_ulp)
428 print_ulps (test_name, ulp);
432 if (print_screen (ok, xfail))
435 printf ("Failure: ");
436 printf ("Test: %s\n", test_name);
437 printf ("Result:\n");
438 printf (" is: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR "\n",
440 printf (" should be: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR "\n",
444 printf (" difference: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR
446 printf (" ulp : % .4" PRINTF_NEXPR "\n", ulp);
447 printf (" max.ulp : % .4" PRINTF_NEXPR "\n", max_ulp);
450 update_stats (ok, xfail);
452 fpstack_test (test_name);
457 check_complex (const char *test_name, FLOAT __complex__ computed,
458 FLOAT __complex__ expected,
459 FLOAT max_ulp, int xfail, int exception)
461 FLOAT part_comp, part_exp;
464 sprintf (str, "Real part of: %s", test_name);
465 part_comp = __real__ computed;
466 part_exp = __real__ expected;
467 check_float (str, part_comp, part_exp, max_ulp, xfail, exception);
469 sprintf (str, "Imaginary part of: %s", test_name);
470 part_comp = __imag__ computed;
471 part_exp = __imag__ expected;
472 /* Don't check again for exceptions, just pass through the
473 zero/inf sign test. */
474 check_float (str, part_comp, part_exp, max_ulp, xfail,
475 exception & IGNORE_ZERO_INF_SIGN);
479 /* Check that computed and expected values are equal (int values). */
481 check_int (const char *test_name, int computed, int expected, int max_ulp,
482 int xfail, int exceptions)
484 int diff = computed - expected;
487 test_exceptions (test_name, exceptions);
489 if (abs (diff) <= max_ulp)
493 print_ulps (test_name, diff);
495 if (print_screen (ok, xfail))
498 printf ("Failure: ");
499 printf ("Test: %s\n", test_name);
500 printf ("Result:\n");
501 printf (" is: %d\n", computed);
502 printf (" should be: %d\n", expected);
505 update_stats (ok, xfail);
506 fpstack_test (test_name);
510 /* Check that computed and expected values are equal (long int values). */
512 check_long (const char *test_name, long int computed, long int expected,
513 long int max_ulp, int xfail, int exceptions)
515 long int diff = computed - expected;
518 test_exceptions (test_name, exceptions);
520 if (labs (diff) <= max_ulp)
524 print_ulps (test_name, diff);
526 if (print_screen (ok, xfail))
529 printf ("Failure: ");
530 printf ("Test: %s\n", test_name);
531 printf ("Result:\n");
532 printf (" is: %ld\n", computed);
533 printf (" should be: %ld\n", expected);
536 update_stats (ok, xfail);
537 fpstack_test (test_name);
541 /* Check that computed value is true/false. */
543 check_bool (const char *test_name, int computed, int expected,
544 long int max_ulp, int xfail, int exceptions)
548 test_exceptions (test_name, exceptions);
550 if ((computed == 0) == (expected == 0))
553 if (print_screen (ok, xfail))
556 printf ("Failure: ");
557 printf ("Test: %s\n", test_name);
558 printf ("Result:\n");
559 printf (" is: %d\n", computed);
560 printf (" should be: %d\n", expected);
563 update_stats (ok, xfail);
564 fpstack_test (test_name);
568 /* check that computed and expected values are equal (long int values) */
570 check_longlong (const char *test_name, long long int computed,
571 long long int expected,
572 long long int max_ulp, int xfail,
575 long long int diff = computed - expected;
578 test_exceptions (test_name, exceptions);
580 if (llabs (diff) <= max_ulp)
584 print_ulps (test_name, diff);
586 if (print_screen (ok, xfail))
590 printf ("Test: %s\n", test_name);
591 printf ("Result:\n");
592 printf (" is: %lld\n", computed);
593 printf (" should be: %lld\n", expected);
596 update_stats (ok, xfail);
597 fpstack_test (test_name);
602 /* This is to prevent messages from the SVID libm emulation. */
604 matherr (struct exception *x __attribute__ ((unused)))
610 /****************************************************************************
611 Tests for single functions of libm.
612 Please keep them alphabetically sorted!
613 ****************************************************************************/
620 TEST_f_f (acos, plus_infty, nan_value, INVALID_EXCEPTION);
621 TEST_f_f (acos, minus_infty, nan_value, INVALID_EXCEPTION);
622 TEST_f_f (acos, nan_value, nan_value);
625 TEST_f_f (acos, 1.1, nan_value, INVALID_EXCEPTION);
626 TEST_f_f (acos, -1.1, nan_value, INVALID_EXCEPTION);
628 TEST_f_f (acos, 0, M_PI_2l);
629 TEST_f_f (acos, minus_zero, M_PI_2l);
630 TEST_f_f (acos, 1, 0);
631 TEST_f_f (acos, -1, M_PIl);
632 TEST_f_f (acos, 0.5, M_PI_6l*2.0);
633 TEST_f_f (acos, -0.5, M_PI_6l*4.0);
634 TEST_f_f (acos, 0.7, 0.7953988301841435554L);
644 TEST_f_f (acosh, plus_infty, plus_infty);
645 TEST_f_f (acosh, minus_infty, nan_value, INVALID_EXCEPTION);
648 TEST_f_f (acosh, -1.1, nan_value, INVALID_EXCEPTION);
650 TEST_f_f (acosh, 1, 0);
651 TEST_f_f (acosh, 7, 2.6339157938496334172L);
662 TEST_f_f (asin, plus_infty, nan_value, INVALID_EXCEPTION);
663 TEST_f_f (asin, minus_infty, nan_value, INVALID_EXCEPTION);
664 TEST_f_f (asin, nan_value, nan_value);
666 /* asin x == NaN plus invalid exception for |x| > 1. */
667 TEST_f_f (asin, 1.1, nan_value, INVALID_EXCEPTION);
668 TEST_f_f (asin, -1.1, nan_value, INVALID_EXCEPTION);
670 TEST_f_f (asin, 0, 0);
671 TEST_f_f (asin, minus_zero, minus_zero);
672 TEST_f_f (asin, 0.5, M_PI_6l);
673 TEST_f_f (asin, -0.5, -M_PI_6l);
674 TEST_f_f (asin, 1.0, M_PI_2l);
675 TEST_f_f (asin, -1.0, -M_PI_2l);
676 TEST_f_f (asin, 0.7, 0.7753974966107530637L);
686 TEST_f_f (asinh, 0, 0);
687 TEST_f_f (asinh, minus_zero, minus_zero);
688 TEST_f_f (asinh, plus_infty, plus_infty);
689 TEST_f_f (asinh, minus_infty, minus_infty);
690 TEST_f_f (asinh, 0.7, 0.652666566082355786L);
701 TEST_f_f (atan, 0, 0);
702 TEST_f_f (atan, minus_zero, minus_zero);
704 TEST_f_f (atan, plus_infty, M_PI_2l);
705 TEST_f_f (atan, minus_infty, -M_PI_2l);
707 TEST_f_f (atan, 1, M_PI_4l);
708 TEST_f_f (atan, -1, -M_PI_4l);
710 TEST_f_f (atan, 0.7, 0.6107259643892086165L);
724 TEST_f_f (atanh, 0, 0);
725 TEST_f_f (atanh, minus_zero, minus_zero);
727 TEST_f_f (atanh, 1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
728 TEST_f_f (atanh, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
730 /* atanh (x) == NaN plus invalid exception if |x| > 1. */
731 TEST_f_f (atanh, 1.1, nan_value, INVALID_EXCEPTION);
732 TEST_f_f (atanh, -1.1, nan_value, INVALID_EXCEPTION);
734 TEST_f_f (atanh, 0.7, 0.8673005276940531944L);
744 /* atan2 (0,x) == 0 for x > 0. */
745 TEST_ff_f (atan2, 0, 1, 0);
747 /* atan2 (-0,x) == -0 for x > 0. */
748 TEST_ff_f (atan2, minus_zero, 1, minus_zero);
750 TEST_ff_f (atan2, 0, 0, 0);
751 TEST_ff_f (atan2, minus_zero, 0, minus_zero);
753 /* atan2 (+0,x) == +pi for x < 0. */
754 TEST_ff_f (atan2, 0, -1, M_PIl);
756 /* atan2 (-0,x) == -pi for x < 0. */
757 TEST_ff_f (atan2, minus_zero, -1, -M_PIl);
759 TEST_ff_f (atan2, 0, minus_zero, M_PIl);
760 TEST_ff_f (atan2, minus_zero, minus_zero, -M_PIl);
762 /* atan2 (y,+0) == pi/2 for y > 0. */
763 TEST_ff_f (atan2, 1, 0, M_PI_2l);
765 /* atan2 (y,-0) == pi/2 for y > 0. */
766 TEST_ff_f (atan2, 1, minus_zero, M_PI_2l);
768 /* atan2 (y,+0) == -pi/2 for y < 0. */
769 TEST_ff_f (atan2, -1, 0, -M_PI_2l);
771 /* atan2 (y,-0) == -pi/2 for y < 0. */
772 TEST_ff_f (atan2, -1, minus_zero, -M_PI_2l);
774 /* atan2 (y,inf) == +0 for finite y > 0. */
775 TEST_ff_f (atan2, 1, plus_infty, 0);
777 /* atan2 (y,inf) == -0 for finite y < 0. */
778 TEST_ff_f (atan2, -1, plus_infty, minus_zero);
780 /* atan2(+inf, x) == pi/2 for finite x. */
781 TEST_ff_f (atan2, plus_infty, -1, M_PI_2l);
783 /* atan2(-inf, x) == -pi/2 for finite x. */
784 TEST_ff_f (atan2, minus_infty, 1, -M_PI_2l);
786 /* atan2 (y,-inf) == +pi for finite y > 0. */
787 TEST_ff_f (atan2, 1, minus_infty, M_PIl);
789 /* atan2 (y,-inf) == -pi for finite y < 0. */
790 TEST_ff_f (atan2, -1, minus_infty, -M_PIl);
792 TEST_ff_f (atan2, plus_infty, plus_infty, M_PI_4l);
793 TEST_ff_f (atan2, minus_infty, plus_infty, -M_PI_4l);
794 TEST_ff_f (atan2, plus_infty, minus_infty, M_PI_34l);
795 TEST_ff_f (atan2, minus_infty, minus_infty, -M_PI_34l);
797 TEST_ff_f (atan2, 0.7, 1, 0.6107259643892086165L);
798 TEST_ff_f (atan2, 0.4, 0.0003, 1.5700463269355215718L);
809 /* cabs (x + iy) is specified as hypot (x,y) */
811 /* cabs (+inf + i x) == +inf. */
812 TEST_c_f (cabs, plus_infty, 1.0, plus_infty);
813 /* cabs (-inf + i x) == +inf. */
814 TEST_c_f (cabs, minus_infty, 1.0, plus_infty);
816 TEST_c_f (cabs, minus_infty, nan_value, plus_infty);
817 TEST_c_f (cabs, minus_infty, nan_value, plus_infty);
819 TEST_c_f (cabs, nan_value, nan_value, nan_value);
821 /* cabs (x,y) == cabs (y,x). */
822 TEST_c_f (cabs, 0.7L, 12.4L, 12.41974234837422060118L);
823 /* cabs (x,y) == cabs (-x,y). */
824 TEST_c_f (cabs, -12.4L, 0.7L, 12.41974234837422060118L);
825 /* cabs (x,y) == cabs (-y,x). */
826 TEST_c_f (cabs, -0.7L, 12.4L, 12.41974234837422060118L);
827 /* cabs (x,y) == cabs (-x,-y). */
828 TEST_c_f (cabs, -12.4L, -0.7L, 12.41974234837422060118L);
829 /* cabs (x,y) == cabs (-y,-x). */
830 TEST_c_f (cabs, -0.7L, -12.4L, 12.41974234837422060118L);
831 /* cabs (x,0) == fabs (x). */
832 TEST_c_f (cabs, -0.7L, 0, 0.7L);
833 TEST_c_f (cabs, 0.7L, 0, 0.7L);
834 TEST_c_f (cabs, -1.0L, 0, 1.0L);
835 TEST_c_f (cabs, 1.0L, 0, 1.0L);
836 TEST_c_f (cabs, -5.7e7L, 0, 5.7e7L);
837 TEST_c_f (cabs, 5.7e7L, 0, 5.7e7L);
839 TEST_c_f (cabs, 0.7, 1.2, 1.3892443989449804508L);
850 TEST_c_c (cacos, 0, 0, M_PI_2l, minus_zero);
851 TEST_c_c (cacos, minus_zero, 0, M_PI_2l, minus_zero);
852 TEST_c_c (cacos, minus_zero, minus_zero, M_PI_2l, 0.0);
853 TEST_c_c (cacos, 0, minus_zero, M_PI_2l, 0.0);
855 TEST_c_c (cacos, minus_infty, plus_infty, M_PI_34l, minus_infty);
856 TEST_c_c (cacos, minus_infty, minus_infty, M_PI_34l, plus_infty);
858 TEST_c_c (cacos, plus_infty, plus_infty, M_PI_4l, minus_infty);
859 TEST_c_c (cacos, plus_infty, minus_infty, M_PI_4l, plus_infty);
861 TEST_c_c (cacos, -10.0, plus_infty, M_PI_2l, minus_infty);
862 TEST_c_c (cacos, -10.0, minus_infty, M_PI_2l, plus_infty);
863 TEST_c_c (cacos, 0, plus_infty, M_PI_2l, minus_infty);
864 TEST_c_c (cacos, 0, minus_infty, M_PI_2l, plus_infty);
865 TEST_c_c (cacos, 0.1, plus_infty, M_PI_2l, minus_infty);
866 TEST_c_c (cacos, 0.1, minus_infty, M_PI_2l, plus_infty);
868 TEST_c_c (cacos, minus_infty, 0, M_PIl, minus_infty);
869 TEST_c_c (cacos, minus_infty, minus_zero, M_PIl, plus_infty);
870 TEST_c_c (cacos, minus_infty, 100, M_PIl, minus_infty);
871 TEST_c_c (cacos, minus_infty, -100, M_PIl, plus_infty);
873 TEST_c_c (cacos, plus_infty, 0, 0.0, minus_infty);
874 TEST_c_c (cacos, plus_infty, minus_zero, 0.0, plus_infty);
875 TEST_c_c (cacos, plus_infty, 0.5, 0.0, minus_infty);
876 TEST_c_c (cacos, plus_infty, -0.5, 0.0, plus_infty);
878 TEST_c_c (cacos, plus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
879 TEST_c_c (cacos, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
881 TEST_c_c (cacos, 0, nan_value, M_PI_2l, nan_value);
882 TEST_c_c (cacos, minus_zero, nan_value, M_PI_2l, nan_value);
884 TEST_c_c (cacos, nan_value, plus_infty, nan_value, minus_infty);
885 TEST_c_c (cacos, nan_value, minus_infty, nan_value, plus_infty);
887 TEST_c_c (cacos, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
888 TEST_c_c (cacos, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
890 TEST_c_c (cacos, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
891 TEST_c_c (cacos, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
893 TEST_c_c (cacos, nan_value, nan_value, nan_value, nan_value);
895 TEST_c_c (cacos, 0.7, 1.2, 1.1351827477151551089L, -1.0927647857577371459L);
896 TEST_c_c (cacos, -2, -3, 2.1414491111159960199L, 1.9833870299165354323L);
908 TEST_c_c (cacosh, 0, 0, 0.0, M_PI_2l);
909 TEST_c_c (cacosh, minus_zero, 0, 0.0, M_PI_2l);
910 TEST_c_c (cacosh, 0, minus_zero, 0.0, -M_PI_2l);
911 TEST_c_c (cacosh, minus_zero, minus_zero, 0.0, -M_PI_2l);
912 TEST_c_c (cacosh, minus_infty, plus_infty, plus_infty, M_PI_34l);
913 TEST_c_c (cacosh, minus_infty, minus_infty, plus_infty, -M_PI_34l);
915 TEST_c_c (cacosh, plus_infty, plus_infty, plus_infty, M_PI_4l);
916 TEST_c_c (cacosh, plus_infty, minus_infty, plus_infty, -M_PI_4l);
918 TEST_c_c (cacosh, -10.0, plus_infty, plus_infty, M_PI_2l);
919 TEST_c_c (cacosh, -10.0, minus_infty, plus_infty, -M_PI_2l);
920 TEST_c_c (cacosh, 0, plus_infty, plus_infty, M_PI_2l);
921 TEST_c_c (cacosh, 0, minus_infty, plus_infty, -M_PI_2l);
922 TEST_c_c (cacosh, 0.1, plus_infty, plus_infty, M_PI_2l);
923 TEST_c_c (cacosh, 0.1, minus_infty, plus_infty, -M_PI_2l);
925 TEST_c_c (cacosh, minus_infty, 0, plus_infty, M_PIl);
926 TEST_c_c (cacosh, minus_infty, minus_zero, plus_infty, -M_PIl);
927 TEST_c_c (cacosh, minus_infty, 100, plus_infty, M_PIl);
928 TEST_c_c (cacosh, minus_infty, -100, plus_infty, -M_PIl);
930 TEST_c_c (cacosh, plus_infty, 0, plus_infty, 0.0);
931 TEST_c_c (cacosh, plus_infty, minus_zero, plus_infty, minus_zero);
932 TEST_c_c (cacosh, plus_infty, 0.5, plus_infty, 0.0);
933 TEST_c_c (cacosh, plus_infty, -0.5, plus_infty, minus_zero);
935 TEST_c_c (cacosh, plus_infty, nan_value, plus_infty, nan_value);
936 TEST_c_c (cacosh, minus_infty, nan_value, plus_infty, nan_value);
938 TEST_c_c (cacosh, 0, nan_value, nan_value, nan_value);
939 TEST_c_c (cacosh, minus_zero, nan_value, nan_value, nan_value);
941 TEST_c_c (cacosh, nan_value, plus_infty, plus_infty, nan_value);
942 TEST_c_c (cacosh, nan_value, minus_infty, plus_infty, nan_value);
944 TEST_c_c (cacosh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
945 TEST_c_c (cacosh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
947 TEST_c_c (cacosh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
948 TEST_c_c (cacosh, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
950 TEST_c_c (cacosh, nan_value, nan_value, nan_value, nan_value);
952 TEST_c_c (cacosh, 0.7, 1.2, 1.0927647857577371459L, 1.1351827477151551089L);
953 TEST_c_c (cacosh, -2, -3, -1.9833870299165354323L, 2.1414491111159960199L);
963 /* carg (x + iy) is specified as atan2 (y, x) */
965 /* carg (x + i 0) == 0 for x > 0. */
966 TEST_c_f (carg, 2.0, 0, 0);
967 /* carg (x - i 0) == -0 for x > 0. */
968 TEST_c_f (carg, 2.0, minus_zero, minus_zero);
970 TEST_c_f (carg, 0, 0, 0);
971 TEST_c_f (carg, 0, minus_zero, minus_zero);
973 /* carg (x + i 0) == +pi for x < 0. */
974 TEST_c_f (carg, -2.0, 0, M_PIl);
976 /* carg (x - i 0) == -pi for x < 0. */
977 TEST_c_f (carg, -2.0, minus_zero, -M_PIl);
979 TEST_c_f (carg, minus_zero, 0, M_PIl);
980 TEST_c_f (carg, minus_zero, minus_zero, -M_PIl);
982 /* carg (+0 + i y) == pi/2 for y > 0. */
983 TEST_c_f (carg, 0, 2.0, M_PI_2l);
985 /* carg (-0 + i y) == pi/2 for y > 0. */
986 TEST_c_f (carg, minus_zero, 2.0, M_PI_2l);
988 /* carg (+0 + i y) == -pi/2 for y < 0. */
989 TEST_c_f (carg, 0, -2.0, -M_PI_2l);
991 /* carg (-0 + i y) == -pi/2 for y < 0. */
992 TEST_c_f (carg, minus_zero, -2.0, -M_PI_2l);
994 /* carg (inf + i y) == +0 for finite y > 0. */
995 TEST_c_f (carg, plus_infty, 2.0, 0);
997 /* carg (inf + i y) == -0 for finite y < 0. */
998 TEST_c_f (carg, plus_infty, -2.0, minus_zero);
1000 /* carg(x + i inf) == pi/2 for finite x. */
1001 TEST_c_f (carg, 10.0, plus_infty, M_PI_2l);
1003 /* carg(x - i inf) == -pi/2 for finite x. */
1004 TEST_c_f (carg, 10.0, minus_infty, -M_PI_2l);
1006 /* carg (-inf + i y) == +pi for finite y > 0. */
1007 TEST_c_f (carg, minus_infty, 10.0, M_PIl);
1009 /* carg (-inf + i y) == -pi for finite y < 0. */
1010 TEST_c_f (carg, minus_infty, -10.0, -M_PIl);
1012 TEST_c_f (carg, plus_infty, plus_infty, M_PI_4l);
1014 TEST_c_f (carg, plus_infty, minus_infty, -M_PI_4l);
1016 TEST_c_f (carg, minus_infty, plus_infty, 3 * M_PI_4l);
1018 TEST_c_f (carg, minus_infty, minus_infty, -3 * M_PI_4l);
1028 TEST_c_c (casin, 0, 0, 0.0, 0.0);
1029 TEST_c_c (casin, minus_zero, 0, minus_zero, 0.0);
1030 TEST_c_c (casin, 0, minus_zero, 0.0, minus_zero);
1031 TEST_c_c (casin, minus_zero, minus_zero, minus_zero, minus_zero);
1033 TEST_c_c (casin, plus_infty, plus_infty, M_PI_4l, plus_infty);
1034 TEST_c_c (casin, plus_infty, minus_infty, M_PI_4l, minus_infty);
1035 TEST_c_c (casin, minus_infty, plus_infty, -M_PI_4l, plus_infty);
1036 TEST_c_c (casin, minus_infty, minus_infty, -M_PI_4l, minus_infty);
1038 TEST_c_c (casin, -10.0, plus_infty, minus_zero, plus_infty);
1039 TEST_c_c (casin, -10.0, minus_infty, minus_zero, minus_infty);
1040 TEST_c_c (casin, 0, plus_infty, 0.0, plus_infty);
1041 TEST_c_c (casin, 0, minus_infty, 0.0, minus_infty);
1042 TEST_c_c (casin, minus_zero, plus_infty, minus_zero, plus_infty);
1043 TEST_c_c (casin, minus_zero, minus_infty, minus_zero, minus_infty);
1044 TEST_c_c (casin, 0.1, plus_infty, 0.0, plus_infty);
1045 TEST_c_c (casin, 0.1, minus_infty, 0.0, minus_infty);
1047 TEST_c_c (casin, minus_infty, 0, -M_PI_2l, plus_infty);
1048 TEST_c_c (casin, minus_infty, minus_zero, -M_PI_2l, minus_infty);
1049 TEST_c_c (casin, minus_infty, 100, -M_PI_2l, plus_infty);
1050 TEST_c_c (casin, minus_infty, -100, -M_PI_2l, minus_infty);
1052 TEST_c_c (casin, plus_infty, 0, M_PI_2l, plus_infty);
1053 TEST_c_c (casin, plus_infty, minus_zero, M_PI_2l, minus_infty);
1054 TEST_c_c (casin, plus_infty, 0.5, M_PI_2l, plus_infty);
1055 TEST_c_c (casin, plus_infty, -0.5, M_PI_2l, minus_infty);
1057 TEST_c_c (casin, nan_value, plus_infty, nan_value, plus_infty);
1058 TEST_c_c (casin, nan_value, minus_infty, nan_value, minus_infty);
1060 TEST_c_c (casin, 0.0, nan_value, 0.0, nan_value);
1061 TEST_c_c (casin, minus_zero, nan_value, minus_zero, nan_value);
1063 TEST_c_c (casin, plus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
1064 TEST_c_c (casin, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
1066 TEST_c_c (casin, nan_value, 10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1067 TEST_c_c (casin, nan_value, -10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1069 TEST_c_c (casin, 0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1070 TEST_c_c (casin, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1072 TEST_c_c (casin, nan_value, nan_value, nan_value, nan_value);
1074 TEST_c_c (casin, 0.7, 1.2, 0.4356135790797415103L, 1.0927647857577371459L);
1075 TEST_c_c (casin, -2, -3, -0.5706527843210994007L, -1.9833870299165354323L);
1086 TEST_c_c (casinh, 0, 0, 0.0, 0.0);
1087 TEST_c_c (casinh, minus_zero, 0, minus_zero, 0);
1088 TEST_c_c (casinh, 0, minus_zero, 0.0, minus_zero);
1089 TEST_c_c (casinh, minus_zero, minus_zero, minus_zero, minus_zero);
1091 TEST_c_c (casinh, plus_infty, plus_infty, plus_infty, M_PI_4l);
1092 TEST_c_c (casinh, plus_infty, minus_infty, plus_infty, -M_PI_4l);
1093 TEST_c_c (casinh, minus_infty, plus_infty, minus_infty, M_PI_4l);
1094 TEST_c_c (casinh, minus_infty, minus_infty, minus_infty, -M_PI_4l);
1096 TEST_c_c (casinh, -10.0, plus_infty, minus_infty, M_PI_2l);
1097 TEST_c_c (casinh, -10.0, minus_infty, minus_infty, -M_PI_2l);
1098 TEST_c_c (casinh, 0, plus_infty, plus_infty, M_PI_2l);
1099 TEST_c_c (casinh, 0, minus_infty, plus_infty, -M_PI_2l);
1100 TEST_c_c (casinh, minus_zero, plus_infty, minus_infty, M_PI_2l);
1101 TEST_c_c (casinh, minus_zero, minus_infty, minus_infty, -M_PI_2l);
1102 TEST_c_c (casinh, 0.1, plus_infty, plus_infty, M_PI_2l);
1103 TEST_c_c (casinh, 0.1, minus_infty, plus_infty, -M_PI_2l);
1105 TEST_c_c (casinh, minus_infty, 0, minus_infty, 0.0);
1106 TEST_c_c (casinh, minus_infty, minus_zero, minus_infty, minus_zero);
1107 TEST_c_c (casinh, minus_infty, 100, minus_infty, 0.0);
1108 TEST_c_c (casinh, minus_infty, -100, minus_infty, minus_zero);
1110 TEST_c_c (casinh, plus_infty, 0, plus_infty, 0.0);
1111 TEST_c_c (casinh, plus_infty, minus_zero, plus_infty, minus_zero);
1112 TEST_c_c (casinh, plus_infty, 0.5, plus_infty, 0.0);
1113 TEST_c_c (casinh, plus_infty, -0.5, plus_infty, minus_zero);
1115 TEST_c_c (casinh, plus_infty, nan_value, plus_infty, nan_value);
1116 TEST_c_c (casinh, minus_infty, nan_value, minus_infty, nan_value);
1118 TEST_c_c (casinh, nan_value, 0, nan_value, 0.0);
1119 TEST_c_c (casinh, nan_value, minus_zero, nan_value, minus_zero);
1121 TEST_c_c (casinh, nan_value, plus_infty, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
1122 TEST_c_c (casinh, nan_value, minus_infty, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
1124 TEST_c_c (casinh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1125 TEST_c_c (casinh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1127 TEST_c_c (casinh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
1128 TEST_c_c (casinh, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1130 TEST_c_c (casinh, nan_value, nan_value, nan_value, nan_value);
1132 TEST_c_c (casinh, 0.7, 1.2, 0.9786545955936738768L, 0.9113541895315601156L);
1133 TEST_c_c (casinh, -2, -3, -1.9686379257930962917L, -0.9646585044076027920L);
1144 TEST_c_c (catan, 0, 0, 0, 0);
1145 TEST_c_c (catan, minus_zero, 0, minus_zero, 0);
1146 TEST_c_c (catan, 0, minus_zero, 0, minus_zero);
1147 TEST_c_c (catan, minus_zero, minus_zero, minus_zero, minus_zero);
1149 TEST_c_c (catan, plus_infty, plus_infty, M_PI_2l, 0);
1150 TEST_c_c (catan, plus_infty, minus_infty, M_PI_2l, minus_zero);
1151 TEST_c_c (catan, minus_infty, plus_infty, -M_PI_2l, 0);
1152 TEST_c_c (catan, minus_infty, minus_infty, -M_PI_2l, minus_zero);
1155 TEST_c_c (catan, plus_infty, -10.0, M_PI_2l, minus_zero);
1156 TEST_c_c (catan, minus_infty, -10.0, -M_PI_2l, minus_zero);
1157 TEST_c_c (catan, plus_infty, minus_zero, M_PI_2l, minus_zero);
1158 TEST_c_c (catan, minus_infty, minus_zero, -M_PI_2l, minus_zero);
1159 TEST_c_c (catan, plus_infty, 0.0, M_PI_2l, 0);
1160 TEST_c_c (catan, minus_infty, 0.0, -M_PI_2l, 0);
1161 TEST_c_c (catan, plus_infty, 0.1, M_PI_2l, 0);
1162 TEST_c_c (catan, minus_infty, 0.1, -M_PI_2l, 0);
1164 TEST_c_c (catan, 0.0, minus_infty, M_PI_2l, minus_zero);
1165 TEST_c_c (catan, minus_zero, minus_infty, -M_PI_2l, minus_zero);
1166 TEST_c_c (catan, 100.0, minus_infty, M_PI_2l, minus_zero);
1167 TEST_c_c (catan, -100.0, minus_infty, -M_PI_2l, minus_zero);
1169 TEST_c_c (catan, 0.0, plus_infty, M_PI_2l, 0);
1170 TEST_c_c (catan, minus_zero, plus_infty, -M_PI_2l, 0);
1171 TEST_c_c (catan, 0.5, plus_infty, M_PI_2l, 0);
1172 TEST_c_c (catan, -0.5, plus_infty, -M_PI_2l, 0);
1174 TEST_c_c (catan, nan_value, 0.0, nan_value, 0);
1175 TEST_c_c (catan, nan_value, minus_zero, nan_value, minus_zero);
1177 TEST_c_c (catan, nan_value, plus_infty, nan_value, 0);
1178 TEST_c_c (catan, nan_value, minus_infty, nan_value, minus_zero);
1180 TEST_c_c (catan, 0.0, nan_value, nan_value, nan_value);
1181 TEST_c_c (catan, minus_zero, nan_value, nan_value, nan_value);
1183 TEST_c_c (catan, plus_infty, nan_value, M_PI_2l, 0, IGNORE_ZERO_INF_SIGN);
1184 TEST_c_c (catan, minus_infty, nan_value, -M_PI_2l, 0, IGNORE_ZERO_INF_SIGN);
1186 TEST_c_c (catan, nan_value, 10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1187 TEST_c_c (catan, nan_value, -10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1189 TEST_c_c (catan, 0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1190 TEST_c_c (catan, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1192 TEST_c_c (catan, nan_value, nan_value, nan_value, nan_value);
1194 TEST_c_c (catan, 0.7, 1.2, 1.0785743834118921877L, 0.5770573776534306764L);
1196 TEST_c_c (catan, -2, -3, -1.4099210495965755225L, -0.2290726829685387662L);
1206 TEST_c_c (catanh, 0, 0, 0.0, 0.0);
1207 TEST_c_c (catanh, minus_zero, 0, minus_zero, 0.0);
1208 TEST_c_c (catanh, 0, minus_zero, 0.0, minus_zero);
1209 TEST_c_c (catanh, minus_zero, minus_zero, minus_zero, minus_zero);
1211 TEST_c_c (catanh, plus_infty, plus_infty, 0.0, M_PI_2l);
1212 TEST_c_c (catanh, plus_infty, minus_infty, 0.0, -M_PI_2l);
1213 TEST_c_c (catanh, minus_infty, plus_infty, minus_zero, M_PI_2l);
1214 TEST_c_c (catanh, minus_infty, minus_infty, minus_zero, -M_PI_2l);
1216 TEST_c_c (catanh, -10.0, plus_infty, minus_zero, M_PI_2l);
1217 TEST_c_c (catanh, -10.0, minus_infty, minus_zero, -M_PI_2l);
1218 TEST_c_c (catanh, minus_zero, plus_infty, minus_zero, M_PI_2l);
1219 TEST_c_c (catanh, minus_zero, minus_infty, minus_zero, -M_PI_2l);
1220 TEST_c_c (catanh, 0, plus_infty, 0.0, M_PI_2l);
1221 TEST_c_c (catanh, 0, minus_infty, 0.0, -M_PI_2l);
1222 TEST_c_c (catanh, 0.1, plus_infty, 0.0, M_PI_2l);
1223 TEST_c_c (catanh, 0.1, minus_infty, 0.0, -M_PI_2l);
1225 TEST_c_c (catanh, minus_infty, 0, minus_zero, M_PI_2l);
1226 TEST_c_c (catanh, minus_infty, minus_zero, minus_zero, -M_PI_2l);
1227 TEST_c_c (catanh, minus_infty, 100, minus_zero, M_PI_2l);
1228 TEST_c_c (catanh, minus_infty, -100, minus_zero, -M_PI_2l);
1230 TEST_c_c (catanh, plus_infty, 0, 0.0, M_PI_2l);
1231 TEST_c_c (catanh, plus_infty, minus_zero, 0.0, -M_PI_2l);
1232 TEST_c_c (catanh, plus_infty, 0.5, 0.0, M_PI_2l);
1233 TEST_c_c (catanh, plus_infty, -0.5, 0.0, -M_PI_2l);
1235 TEST_c_c (catanh, 0, nan_value, 0.0, nan_value);
1236 TEST_c_c (catanh, minus_zero, nan_value, minus_zero, nan_value);
1238 TEST_c_c (catanh, plus_infty, nan_value, 0.0, nan_value);
1239 TEST_c_c (catanh, minus_infty, nan_value, minus_zero, nan_value);
1241 TEST_c_c (catanh, nan_value, 0, nan_value, nan_value);
1242 TEST_c_c (catanh, nan_value, minus_zero, nan_value, nan_value);
1244 TEST_c_c (catanh, nan_value, plus_infty, 0.0, M_PI_2l, IGNORE_ZERO_INF_SIGN);
1245 TEST_c_c (catanh, nan_value, minus_infty, 0.0, -M_PI_2l, IGNORE_ZERO_INF_SIGN);
1247 TEST_c_c (catanh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1248 TEST_c_c (catanh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1250 TEST_c_c (catanh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
1251 TEST_c_c (catanh, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
1253 TEST_c_c (catanh, nan_value, nan_value, nan_value, nan_value);
1255 TEST_c_c (catanh, 0.7, 1.2, 0.2600749516525135959L, 0.9702403077950989849L);
1256 TEST_c_c (catanh, -2, -3, -0.1469466662255297520L, -1.3389725222944935611L);
1266 TEST_f_f (cbrt, 0.0, 0.0);
1267 TEST_f_f (cbrt, minus_zero, minus_zero);
1269 TEST_f_f (cbrt, plus_infty, plus_infty);
1270 TEST_f_f (cbrt, minus_infty, minus_infty);
1271 TEST_f_f (cbrt, nan_value, nan_value);
1273 TEST_f_f (cbrt, -0.001, -0.1);
1274 TEST_f_f (cbrt, 8, 2);
1275 TEST_f_f (cbrt, -27.0, -3.0);
1276 TEST_f_f (cbrt, 0.970299, 0.99);
1277 TEST_f_f (cbrt, 0.7, 0.8879040017426007084L);
1288 TEST_c_c (ccos, 0.0, 0.0, 1.0, minus_zero);
1289 TEST_c_c (ccos, minus_zero, 0.0, 1.0, 0.0);
1290 TEST_c_c (ccos, 0.0, minus_zero, 1.0, 0.0);
1291 TEST_c_c (ccos, minus_zero, minus_zero, 1.0, minus_zero);
1293 TEST_c_c (ccos, plus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1294 TEST_c_c (ccos, plus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1295 TEST_c_c (ccos, minus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1296 TEST_c_c (ccos, minus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1298 TEST_c_c (ccos, 0.0, plus_infty, plus_infty, minus_zero);
1299 TEST_c_c (ccos, 0.0, minus_infty, plus_infty, 0.0);
1300 TEST_c_c (ccos, minus_zero, plus_infty, plus_infty, 0.0);
1301 TEST_c_c (ccos, minus_zero, minus_infty, plus_infty, minus_zero);
1303 TEST_c_c (ccos, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1304 TEST_c_c (ccos, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1305 TEST_c_c (ccos, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1306 TEST_c_c (ccos, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1308 TEST_c_c (ccos, 4.625, plus_infty, minus_infty, plus_infty);
1309 TEST_c_c (ccos, 4.625, minus_infty, minus_infty, minus_infty);
1310 TEST_c_c (ccos, -4.625, plus_infty, minus_infty, minus_infty);
1311 TEST_c_c (ccos, -4.625, minus_infty, minus_infty, plus_infty);
1313 TEST_c_c (ccos, plus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
1314 TEST_c_c (ccos, plus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
1315 TEST_c_c (ccos, minus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
1316 TEST_c_c (ccos, minus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
1318 TEST_c_c (ccos, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1319 TEST_c_c (ccos, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1321 TEST_c_c (ccos, nan_value, plus_infty, plus_infty, nan_value);
1322 TEST_c_c (ccos, nan_value, minus_infty, plus_infty, nan_value);
1324 TEST_c_c (ccos, nan_value, 9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1325 TEST_c_c (ccos, nan_value, -9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1327 TEST_c_c (ccos, 0.0, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1328 TEST_c_c (ccos, minus_zero, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1330 TEST_c_c (ccos, 10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1331 TEST_c_c (ccos, -10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1333 TEST_c_c (ccos, plus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1334 TEST_c_c (ccos, minus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1336 TEST_c_c (ccos, nan_value, nan_value, nan_value, nan_value);
1338 TEST_c_c (ccos, 0.7, 1.2, 1.3848657645312111080L, -0.97242170335830028619L);
1340 TEST_c_c (ccos, -2, -3, -4.1896256909688072301L, -9.1092278937553365979L);
1352 TEST_c_c (ccosh, 0.0, 0.0, 1.0, 0.0);
1353 TEST_c_c (ccosh, minus_zero, 0.0, 1.0, minus_zero);
1354 TEST_c_c (ccosh, 0.0, minus_zero, 1.0, minus_zero);
1355 TEST_c_c (ccosh, minus_zero, minus_zero, 1.0, 0.0);
1357 TEST_c_c (ccosh, 0.0, plus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1358 TEST_c_c (ccosh, minus_zero, plus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1359 TEST_c_c (ccosh, 0.0, minus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1360 TEST_c_c (ccosh, minus_zero, minus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1362 TEST_c_c (ccosh, plus_infty, 0.0, plus_infty, 0.0);
1363 TEST_c_c (ccosh, minus_infty, 0.0, plus_infty, minus_zero);
1364 TEST_c_c (ccosh, plus_infty, minus_zero, plus_infty, minus_zero);
1365 TEST_c_c (ccosh, minus_infty, minus_zero, plus_infty, 0.0);
1367 TEST_c_c (ccosh, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1368 TEST_c_c (ccosh, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1369 TEST_c_c (ccosh, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1370 TEST_c_c (ccosh, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
1372 TEST_c_c (ccosh, plus_infty, 4.625, minus_infty, minus_infty);
1373 TEST_c_c (ccosh, minus_infty, 4.625, minus_infty, plus_infty);
1374 TEST_c_c (ccosh, plus_infty, -4.625, minus_infty, plus_infty);
1375 TEST_c_c (ccosh, minus_infty, -4.625, minus_infty, minus_infty);
1377 TEST_c_c (ccosh, 6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1378 TEST_c_c (ccosh, -6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1379 TEST_c_c (ccosh, 6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1380 TEST_c_c (ccosh, -6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1382 TEST_c_c (ccosh, 0.0, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1383 TEST_c_c (ccosh, minus_zero, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1385 TEST_c_c (ccosh, plus_infty, nan_value, plus_infty, nan_value);
1386 TEST_c_c (ccosh, minus_infty, nan_value, plus_infty, nan_value);
1388 TEST_c_c (ccosh, 9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1389 TEST_c_c (ccosh, -9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1391 TEST_c_c (ccosh, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1392 TEST_c_c (ccosh, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1394 TEST_c_c (ccosh, nan_value, 10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1395 TEST_c_c (ccosh, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1397 TEST_c_c (ccosh, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
1398 TEST_c_c (ccosh, nan_value, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
1400 TEST_c_c (ccosh, nan_value, nan_value, nan_value, nan_value);
1402 TEST_c_c (ccosh, 0.7, 1.2, 0.4548202223691477654L, 0.7070296600921537682L);
1404 TEST_c_c (ccosh, -2, -3, -3.7245455049153225654L, 0.5118225699873846088L);
1415 TEST_f_f (ceil, 0.0, 0.0);
1416 TEST_f_f (ceil, minus_zero, minus_zero);
1417 TEST_f_f (ceil, plus_infty, plus_infty);
1418 TEST_f_f (ceil, minus_infty, minus_infty);
1420 TEST_f_f (ceil, M_PIl, 4.0);
1421 TEST_f_f (ceil, -M_PIl, -3.0);
1432 TEST_c_c (cexp, plus_zero, plus_zero, 1, 0.0);
1433 TEST_c_c (cexp, minus_zero, plus_zero, 1, 0.0);
1434 TEST_c_c (cexp, plus_zero, minus_zero, 1, minus_zero);
1435 TEST_c_c (cexp, minus_zero, minus_zero, 1, minus_zero);
1437 TEST_c_c (cexp, plus_infty, plus_zero, plus_infty, 0.0);
1438 TEST_c_c (cexp, plus_infty, minus_zero, plus_infty, minus_zero);
1440 TEST_c_c (cexp, minus_infty, plus_zero, 0.0, 0.0);
1441 TEST_c_c (cexp, minus_infty, minus_zero, 0.0, minus_zero);
1443 TEST_c_c (cexp, 0.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1444 TEST_c_c (cexp, minus_zero, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1446 TEST_c_c (cexp, 0.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1447 TEST_c_c (cexp, minus_zero, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1449 TEST_c_c (cexp, 100.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1450 TEST_c_c (cexp, -100.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1452 TEST_c_c (cexp, 100.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1453 TEST_c_c (cexp, -100.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1455 TEST_c_c (cexp, minus_infty, 2.0, minus_zero, 0.0);
1456 TEST_c_c (cexp, minus_infty, 4.0, minus_zero, minus_zero);
1457 TEST_c_c (cexp, plus_infty, 2.0, minus_infty, plus_infty);
1458 TEST_c_c (cexp, plus_infty, 4.0, minus_infty, minus_infty);
1460 TEST_c_c (cexp, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1461 TEST_c_c (cexp, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1463 TEST_c_c (cexp, minus_infty, plus_infty, 0.0, 0.0, IGNORE_ZERO_INF_SIGN);
1464 TEST_c_c (cexp, minus_infty, minus_infty, 0.0, minus_zero, IGNORE_ZERO_INF_SIGN);
1466 TEST_c_c (cexp, minus_infty, nan_value, 0, 0, IGNORE_ZERO_INF_SIGN);
1468 TEST_c_c (cexp, plus_infty, nan_value, plus_infty, nan_value);
1470 TEST_c_c (cexp, nan_value, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1471 TEST_c_c (cexp, nan_value, 1.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1473 TEST_c_c (cexp, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
1474 TEST_c_c (cexp, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1475 TEST_c_c (cexp, 1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1476 TEST_c_c (cexp, nan_value, nan_value, nan_value, nan_value);
1478 TEST_c_c (cexp, 0.7, 1.2, 0.7296989091503236012L, 1.8768962328348102821L);
1479 TEST_c_c (cexp, -2.0, -3.0, -0.1339809149295426134L, -0.0190985162611351964L);
1489 TEST_c_c (clog, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
1490 TEST_c_c (clog, minus_zero, minus_zero, minus_infty, -M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
1492 TEST_c_c (clog, 0, 0, minus_infty, 0.0, DIVIDE_BY_ZERO_EXCEPTION);
1493 TEST_c_c (clog, 0, minus_zero, minus_infty, minus_zero, DIVIDE_BY_ZERO_EXCEPTION);
1495 TEST_c_c (clog, minus_infty, plus_infty, plus_infty, M_PI_34l);
1496 TEST_c_c (clog, minus_infty, minus_infty, plus_infty, -M_PI_34l);
1498 TEST_c_c (clog, plus_infty, plus_infty, plus_infty, M_PI_4l);
1499 TEST_c_c (clog, plus_infty, minus_infty, plus_infty, -M_PI_4l);
1501 TEST_c_c (clog, 0, plus_infty, plus_infty, M_PI_2l);
1502 TEST_c_c (clog, 3, plus_infty, plus_infty, M_PI_2l);
1503 TEST_c_c (clog, minus_zero, plus_infty, plus_infty, M_PI_2l);
1504 TEST_c_c (clog, -3, plus_infty, plus_infty, M_PI_2l);
1505 TEST_c_c (clog, 0, minus_infty, plus_infty, -M_PI_2l);
1506 TEST_c_c (clog, 3, minus_infty, plus_infty, -M_PI_2l);
1507 TEST_c_c (clog, minus_zero, minus_infty, plus_infty, -M_PI_2l);
1508 TEST_c_c (clog, -3, minus_infty, plus_infty, -M_PI_2l);
1510 TEST_c_c (clog, minus_infty, 0, plus_infty, M_PIl);
1511 TEST_c_c (clog, minus_infty, 1, plus_infty, M_PIl);
1512 TEST_c_c (clog, minus_infty, minus_zero, plus_infty, -M_PIl);
1513 TEST_c_c (clog, minus_infty, -1, plus_infty, -M_PIl);
1515 TEST_c_c (clog, plus_infty, 0, plus_infty, 0.0);
1516 TEST_c_c (clog, plus_infty, 1, plus_infty, 0.0);
1517 TEST_c_c (clog, plus_infty, minus_zero, plus_infty, minus_zero);
1518 TEST_c_c (clog, plus_infty, -1, plus_infty, minus_zero);
1520 TEST_c_c (clog, plus_infty, nan_value, plus_infty, nan_value);
1521 TEST_c_c (clog, minus_infty, nan_value, plus_infty, nan_value);
1523 TEST_c_c (clog, nan_value, plus_infty, plus_infty, nan_value);
1524 TEST_c_c (clog, nan_value, minus_infty, plus_infty, nan_value);
1526 TEST_c_c (clog, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1527 TEST_c_c (clog, 3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1528 TEST_c_c (clog, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1529 TEST_c_c (clog, -3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1531 TEST_c_c (clog, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1532 TEST_c_c (clog, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1533 TEST_c_c (clog, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
1534 TEST_c_c (clog, nan_value, -5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1536 TEST_c_c (clog, nan_value, nan_value, nan_value, nan_value);
1537 TEST_c_c (clog, -2, -3, 1.2824746787307683680L, -2.1587989303424641704L);
1548 TEST_c_c (clog10, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
1549 TEST_c_c (clog10, minus_zero, minus_zero, minus_infty, -M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
1551 TEST_c_c (clog10, 0, 0, minus_infty, 0.0, DIVIDE_BY_ZERO_EXCEPTION);
1552 TEST_c_c (clog10, 0, minus_zero, minus_infty, minus_zero, DIVIDE_BY_ZERO_EXCEPTION);
1554 TEST_c_c (clog10, minus_infty, plus_infty, plus_infty, M_PI_34_LOG10El);
1556 TEST_c_c (clog10, plus_infty, plus_infty, plus_infty, M_PI4_LOG10El);
1557 TEST_c_c (clog10, plus_infty, minus_infty, plus_infty, -M_PI4_LOG10El);
1559 TEST_c_c (clog10, 0, plus_infty, plus_infty, M_PI2_LOG10El);
1560 TEST_c_c (clog10, 3, plus_infty, plus_infty, M_PI2_LOG10El);
1561 TEST_c_c (clog10, minus_zero, plus_infty, plus_infty, M_PI2_LOG10El);
1562 TEST_c_c (clog10, -3, plus_infty, plus_infty, M_PI2_LOG10El);
1563 TEST_c_c (clog10, 0, minus_infty, plus_infty, -M_PI2_LOG10El);
1564 TEST_c_c (clog10, 3, minus_infty, plus_infty, -M_PI2_LOG10El);
1565 TEST_c_c (clog10, minus_zero, minus_infty, plus_infty, -M_PI2_LOG10El);
1566 TEST_c_c (clog10, -3, minus_infty, plus_infty, -M_PI2_LOG10El);
1568 TEST_c_c (clog10, minus_infty, 0, plus_infty, M_PI_LOG10El);
1569 TEST_c_c (clog10, minus_infty, 1, plus_infty, M_PI_LOG10El);
1570 TEST_c_c (clog10, minus_infty, minus_zero, plus_infty, -M_PI_LOG10El);
1571 TEST_c_c (clog10, minus_infty, -1, plus_infty, -M_PI_LOG10El);
1573 TEST_c_c (clog10, plus_infty, 0, plus_infty, 0.0);
1574 TEST_c_c (clog10, plus_infty, 1, plus_infty, 0.0);
1575 TEST_c_c (clog10, plus_infty, minus_zero, plus_infty, minus_zero);
1576 TEST_c_c (clog10, plus_infty, -1, plus_infty, minus_zero);
1578 TEST_c_c (clog10, plus_infty, nan_value, plus_infty, nan_value);
1579 TEST_c_c (clog10, minus_infty, nan_value, plus_infty, nan_value);
1581 TEST_c_c (clog10, nan_value, plus_infty, plus_infty, nan_value);
1582 TEST_c_c (clog10, nan_value, minus_infty, plus_infty, nan_value);
1584 TEST_c_c (clog10, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1585 TEST_c_c (clog10, 3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1586 TEST_c_c (clog10, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1587 TEST_c_c (clog10, -3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1589 TEST_c_c (clog10, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1590 TEST_c_c (clog10, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1591 TEST_c_c (clog10, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
1592 TEST_c_c (clog10, nan_value, -5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1594 TEST_c_c (clog10, nan_value, nan_value, nan_value, nan_value);
1596 TEST_c_c (clog10, 0.7, 1.2, 0.1427786545038868803L, 0.4528483579352493248L);
1597 TEST_c_c (clog10, -2, -3, 0.5569716761534183846L, -0.9375544629863747085L);
1603 copysign_test (void)
1607 TEST_ff_f (copysign, 0, 4, 0);
1608 TEST_ff_f (copysign, 0, -4, minus_zero);
1609 TEST_ff_f (copysign, minus_zero, 4, 0);
1610 TEST_ff_f (copysign, minus_zero, -4, minus_zero);
1612 TEST_ff_f (copysign, plus_infty, 0, plus_infty);
1613 TEST_ff_f (copysign, plus_infty, minus_zero, minus_infty);
1614 TEST_ff_f (copysign, minus_infty, 0, plus_infty);
1615 TEST_ff_f (copysign, minus_infty, minus_zero, minus_infty);
1617 TEST_ff_f (copysign, 0, plus_infty, 0);
1618 TEST_ff_f (copysign, 0, minus_zero, minus_zero);
1619 TEST_ff_f (copysign, minus_zero, plus_infty, 0);
1620 TEST_ff_f (copysign, minus_zero, minus_zero, minus_zero);
1622 /* XXX More correctly we would have to check the sign of the NaN. */
1623 TEST_ff_f (copysign, nan_value, 0, nan_value);
1624 TEST_ff_f (copysign, nan_value, minus_zero, nan_value);
1625 TEST_ff_f (copysign, -nan_value, 0, nan_value);
1626 TEST_ff_f (copysign, -nan_value, minus_zero, nan_value);
1636 TEST_f_f (cos, 0, 1);
1637 TEST_f_f (cos, minus_zero, 1);
1638 TEST_f_f (cos, plus_infty, nan_value, INVALID_EXCEPTION);
1639 TEST_f_f (cos, minus_infty, nan_value, INVALID_EXCEPTION);
1641 TEST_f_f (cos, M_PI_6l * 2.0, 0.5);
1642 TEST_f_f (cos, M_PI_6l * 4.0, -0.5);
1643 TEST_f_f (cos, M_PI_2l, 0);
1645 TEST_f_f (cos, 0.7, 0.7648421872844884262L);
1654 TEST_f_f (cosh, 0, 1);
1655 TEST_f_f (cosh, minus_zero, 1);
1657 TEST_f_f (cosh, plus_infty, plus_infty);
1658 TEST_f_f (cosh, minus_infty, plus_infty);
1660 TEST_f_f (cosh, 0.7, 1.255169005630943018L);
1670 TEST_cc_c (cpow, 1, 0, 0, 0, 1.0, 0.0);
1671 TEST_cc_c (cpow, 2, 0, 10, 0, 1024.0, 0.0);
1673 TEST_cc_c (cpow, M_El, 0, 0, 2 * M_PIl, 1.0, 0.0);
1674 TEST_cc_c (cpow, 2, 3, 4, 0, -119.0, -120.0);
1683 TEST_c_c (cproj, 0.0, 0.0, 0.0, 0.0);
1684 TEST_c_c (cproj, minus_zero, minus_zero, minus_zero, minus_zero);
1685 TEST_c_c (cproj, 0.0, minus_zero, 0.0, minus_zero);
1686 TEST_c_c (cproj, minus_zero, 0.0, minus_zero, 0.0);
1688 TEST_c_c (cproj, nan_value, nan_value, nan_value, nan_value);
1690 TEST_c_c (cproj, plus_infty, plus_infty, plus_infty, 0.0);
1691 TEST_c_c (cproj, plus_infty, minus_infty, plus_infty, minus_zero);
1692 TEST_c_c (cproj, minus_infty, plus_infty, plus_infty, 0.0);
1693 TEST_c_c (cproj, minus_infty, minus_infty, plus_infty, minus_zero);
1695 TEST_c_c (cproj, 1.0, 0.0, 1.0, 0.0);
1696 TEST_c_c (cproj, 2.0, 3.0, 0.28571428571428571429L, .42857142857142857143L);
1708 TEST_c_c (csin, 0.0, 0.0, 0.0, 0.0);
1709 TEST_c_c (csin, minus_zero, 0.0, minus_zero, 0.0);
1710 TEST_c_c (csin, 0.0, minus_zero, 0, minus_zero);
1711 TEST_c_c (csin, minus_zero, minus_zero, minus_zero, minus_zero);
1713 TEST_c_c (csin, 0.0, plus_infty, 0.0, plus_infty);
1714 TEST_c_c (csin, minus_zero, plus_infty, minus_zero, plus_infty);
1715 TEST_c_c (csin, 0.0, minus_infty, 0.0, minus_infty);
1716 TEST_c_c (csin, minus_zero, minus_infty, minus_zero, minus_infty);
1718 TEST_c_c (csin, plus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1719 TEST_c_c (csin, minus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1720 TEST_c_c (csin, plus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1721 TEST_c_c (csin, minus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1723 TEST_c_c (csin, plus_infty, plus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1724 TEST_c_c (csin, minus_infty, plus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1725 TEST_c_c (csin, plus_infty, minus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1726 TEST_c_c (csin, minus_infty, minus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1728 TEST_c_c (csin, plus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
1729 TEST_c_c (csin, plus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
1730 TEST_c_c (csin, minus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
1731 TEST_c_c (csin, minus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
1733 TEST_c_c (csin, 4.625, plus_infty, minus_infty, minus_infty);
1734 TEST_c_c (csin, 4.625, minus_infty, minus_infty, plus_infty);
1735 TEST_c_c (csin, -4.625, plus_infty, plus_infty, minus_infty);
1736 TEST_c_c (csin, -4.625, minus_infty, plus_infty, plus_infty);
1738 TEST_c_c (csin, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1739 TEST_c_c (csin, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
1741 TEST_c_c (csin, nan_value, plus_infty, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
1742 TEST_c_c (csin, nan_value, minus_infty, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
1744 TEST_c_c (csin, nan_value, 9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1745 TEST_c_c (csin, nan_value, -9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1747 TEST_c_c (csin, 0.0, nan_value, 0.0, nan_value);
1748 TEST_c_c (csin, minus_zero, nan_value, minus_zero, nan_value);
1750 TEST_c_c (csin, 10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1751 TEST_c_c (csin, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1753 TEST_c_c (csin, plus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1754 TEST_c_c (csin, minus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1756 TEST_c_c (csin, nan_value, nan_value, nan_value, nan_value);
1758 TEST_c_c (csin, 0.7, 1.2, 1.1664563419657581376L, 1.1544997246948547371L);
1760 TEST_c_c (csin, -2, -3, -9.1544991469114295734L, 4.1689069599665643507L);
1772 TEST_c_c (csinh, 0.0, 0.0, 0.0, 0.0);
1773 TEST_c_c (csinh, minus_zero, 0.0, minus_zero, 0.0);
1774 TEST_c_c (csinh, 0.0, minus_zero, 0.0, minus_zero);
1775 TEST_c_c (csinh, minus_zero, minus_zero, minus_zero, minus_zero);
1777 TEST_c_c (csinh, 0.0, plus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1778 TEST_c_c (csinh, minus_zero, plus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1779 TEST_c_c (csinh, 0.0, minus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1780 TEST_c_c (csinh, minus_zero, minus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1782 TEST_c_c (csinh, plus_infty, 0.0, plus_infty, 0.0);
1783 TEST_c_c (csinh, minus_infty, 0.0, minus_infty, 0.0);
1784 TEST_c_c (csinh, plus_infty, minus_zero, plus_infty, minus_zero);
1785 TEST_c_c (csinh, minus_infty, minus_zero, minus_infty, minus_zero);
1787 TEST_c_c (csinh, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1788 TEST_c_c (csinh, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1789 TEST_c_c (csinh, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1790 TEST_c_c (csinh, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
1792 TEST_c_c (csinh, plus_infty, 4.625, minus_infty, minus_infty);
1793 TEST_c_c (csinh, minus_infty, 4.625, plus_infty, minus_infty);
1794 TEST_c_c (csinh, plus_infty, -4.625, minus_infty, plus_infty);
1795 TEST_c_c (csinh, minus_infty, -4.625, plus_infty, plus_infty);
1797 TEST_c_c (csinh, 6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1798 TEST_c_c (csinh, -6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1799 TEST_c_c (csinh, 6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1800 TEST_c_c (csinh, -6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1802 TEST_c_c (csinh, 0.0, nan_value, 0.0, nan_value, IGNORE_ZERO_INF_SIGN);
1803 TEST_c_c (csinh, minus_zero, nan_value, 0.0, nan_value, IGNORE_ZERO_INF_SIGN);
1805 TEST_c_c (csinh, plus_infty, nan_value, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
1806 TEST_c_c (csinh, minus_infty, nan_value, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
1808 TEST_c_c (csinh, 9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION);
1809 TEST_c_c (csinh, -9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION);
1811 TEST_c_c (csinh, nan_value, 0.0, nan_value, 0.0);
1812 TEST_c_c (csinh, nan_value, minus_zero, nan_value, minus_zero);
1814 TEST_c_c (csinh, nan_value, 10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1815 TEST_c_c (csinh, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1817 TEST_c_c (csinh, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
1818 TEST_c_c (csinh, nan_value, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
1820 TEST_c_c (csinh, nan_value, nan_value, nan_value, nan_value);
1822 TEST_c_c (csinh, 0.7, 1.2, 0.27487868678117583582L, 1.1698665727426565139L);
1823 TEST_c_c (csinh, -2, -3, 3.5905645899857799520L, -0.5309210862485198052L);
1833 TEST_c_c (csqrt, 0, 0, 0.0, 0.0);
1834 TEST_c_c (csqrt, 0, minus_zero, 0, minus_zero);
1835 TEST_c_c (csqrt, minus_zero, 0, 0.0, 0.0);
1836 TEST_c_c (csqrt, minus_zero, minus_zero, 0.0, minus_zero);
1838 TEST_c_c (csqrt, minus_infty, 0, 0.0, plus_infty);
1839 TEST_c_c (csqrt, minus_infty, 6, 0.0, plus_infty);
1840 TEST_c_c (csqrt, minus_infty, minus_zero, 0.0, minus_infty);
1841 TEST_c_c (csqrt, minus_infty, -6, 0.0, minus_infty);
1843 TEST_c_c (csqrt, plus_infty, 0, plus_infty, 0.0);
1844 TEST_c_c (csqrt, plus_infty, 6, plus_infty, 0.0);
1845 TEST_c_c (csqrt, plus_infty, minus_zero, plus_infty, minus_zero);
1846 TEST_c_c (csqrt, plus_infty, -6, plus_infty, minus_zero);
1848 TEST_c_c (csqrt, 0, plus_infty, plus_infty, plus_infty);
1849 TEST_c_c (csqrt, 4, plus_infty, plus_infty, plus_infty);
1850 TEST_c_c (csqrt, plus_infty, plus_infty, plus_infty, plus_infty);
1851 TEST_c_c (csqrt, minus_zero, plus_infty, plus_infty, plus_infty);
1852 TEST_c_c (csqrt, -4, plus_infty, plus_infty, plus_infty);
1853 TEST_c_c (csqrt, minus_infty, plus_infty, plus_infty, plus_infty);
1854 TEST_c_c (csqrt, 0, minus_infty, plus_infty, minus_infty);
1855 TEST_c_c (csqrt, 4, minus_infty, plus_infty, minus_infty);
1856 TEST_c_c (csqrt, plus_infty, minus_infty, plus_infty, minus_infty);
1857 TEST_c_c (csqrt, minus_zero, minus_infty, plus_infty, minus_infty);
1858 TEST_c_c (csqrt, -4, minus_infty, plus_infty, minus_infty);
1859 TEST_c_c (csqrt, minus_infty, minus_infty, plus_infty, minus_infty);
1861 TEST_c_c (csqrt, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
1863 TEST_c_c (csqrt, plus_infty, nan_value, plus_infty, nan_value);
1865 TEST_c_c (csqrt, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1866 TEST_c_c (csqrt, 1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1867 TEST_c_c (csqrt, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1868 TEST_c_c (csqrt, -1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1870 TEST_c_c (csqrt, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1871 TEST_c_c (csqrt, nan_value, 8, nan_value, nan_value, INVALID_EXCEPTION_OK);
1872 TEST_c_c (csqrt, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
1873 TEST_c_c (csqrt, nan_value, -8, nan_value, nan_value, INVALID_EXCEPTION_OK);
1875 TEST_c_c (csqrt, nan_value, nan_value, nan_value, nan_value);
1877 TEST_c_c (csqrt, 16.0, -30.0, 5.0, -3.0);
1878 TEST_c_c (csqrt, -1, 0, 0.0, 1.0);
1879 TEST_c_c (csqrt, 0, 2, 1.0, 1.0);
1880 TEST_c_c (csqrt, 119, 120, 12.0, 5.0);
1881 TEST_c_c (csqrt, 0.7, 1.2, 1.0220676100300264507L, 0.5870453129635652115L);
1882 TEST_c_c (csqrt, -2, -3, 0.8959774761298381247L, -1.6741492280355400404L);
1883 TEST_c_c (csqrt, -2, 3, 0.8959774761298381247L, 1.6741492280355400404L);
1893 TEST_c_c (ctan, 0, 0, 0.0, 0.0);
1894 TEST_c_c (ctan, 0, minus_zero, 0.0, minus_zero);
1895 TEST_c_c (ctan, minus_zero, 0, minus_zero, 0.0);
1896 TEST_c_c (ctan, minus_zero, minus_zero, minus_zero, minus_zero);
1898 TEST_c_c (ctan, 0, plus_infty, 0.0, 1.0);
1899 TEST_c_c (ctan, 1, plus_infty, 0.0, 1.0);
1900 TEST_c_c (ctan, minus_zero, plus_infty, minus_zero, 1.0);
1901 TEST_c_c (ctan, -1, plus_infty, minus_zero, 1.0);
1903 TEST_c_c (ctan, 0, minus_infty, 0.0, -1.0);
1904 TEST_c_c (ctan, 1, minus_infty, 0.0, -1.0);
1905 TEST_c_c (ctan, minus_zero, minus_infty, minus_zero, -1.0);
1906 TEST_c_c (ctan, -1, minus_infty, minus_zero, -1.0);
1908 TEST_c_c (ctan, plus_infty, 0, nan_value, nan_value, INVALID_EXCEPTION);
1909 TEST_c_c (ctan, plus_infty, 2, nan_value, nan_value, INVALID_EXCEPTION);
1910 TEST_c_c (ctan, minus_infty, 0, nan_value, nan_value, INVALID_EXCEPTION);
1911 TEST_c_c (ctan, minus_infty, 2, nan_value, nan_value, INVALID_EXCEPTION);
1912 TEST_c_c (ctan, plus_infty, minus_zero, nan_value, nan_value, INVALID_EXCEPTION);
1913 TEST_c_c (ctan, plus_infty, -2, nan_value, nan_value, INVALID_EXCEPTION);
1914 TEST_c_c (ctan, minus_infty, minus_zero, nan_value, nan_value, INVALID_EXCEPTION);
1915 TEST_c_c (ctan, minus_infty, -2, nan_value, nan_value, INVALID_EXCEPTION);
1917 TEST_c_c (ctan, nan_value, plus_infty, 0.0, 1.0, IGNORE_ZERO_INF_SIGN);
1918 TEST_c_c (ctan, nan_value, minus_infty, 0.0, -1.0, IGNORE_ZERO_INF_SIGN);
1920 TEST_c_c (ctan, 0, nan_value, 0.0, nan_value);
1921 TEST_c_c (ctan, minus_zero, nan_value, minus_zero, nan_value);
1923 TEST_c_c (ctan, 0.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1924 TEST_c_c (ctan, -4.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1926 TEST_c_c (ctan, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
1927 TEST_c_c (ctan, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1928 TEST_c_c (ctan, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
1929 TEST_c_c (ctan, nan_value, -0.25, nan_value, nan_value, INVALID_EXCEPTION_OK);
1931 TEST_c_c (ctan, nan_value, nan_value, nan_value, nan_value);
1933 TEST_c_c (ctan, 0.7, 1.2, 0.1720734197630349001L, 0.9544807059989405538L);
1934 TEST_c_c (ctan, -2, -3, 0.0037640256415042482L, -1.0032386273536098014L);
1945 TEST_c_c (ctanh, 0, 0, 0.0, 0.0);
1946 TEST_c_c (ctanh, 0, minus_zero, 0.0, minus_zero);
1947 TEST_c_c (ctanh, minus_zero, 0, minus_zero, 0.0);
1948 TEST_c_c (ctanh, minus_zero, minus_zero, minus_zero, minus_zero);
1950 TEST_c_c (ctanh, plus_infty, 0, 1.0, 0.0);
1951 TEST_c_c (ctanh, plus_infty, 1, 1.0, 0.0);
1952 TEST_c_c (ctanh, plus_infty, minus_zero, 1.0, minus_zero);
1953 TEST_c_c (ctanh, plus_infty, -1, 1.0, minus_zero);
1954 TEST_c_c (ctanh, minus_infty, 0, -1.0, 0.0);
1955 TEST_c_c (ctanh, minus_infty, 1, -1.0, 0.0);
1956 TEST_c_c (ctanh, minus_infty, minus_zero, -1.0, minus_zero);
1957 TEST_c_c (ctanh, minus_infty, -1, -1.0, minus_zero);
1959 TEST_c_c (ctanh, 0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1960 TEST_c_c (ctanh, 2, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1961 TEST_c_c (ctanh, 0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1962 TEST_c_c (ctanh, 2, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1963 TEST_c_c (ctanh, minus_zero, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1964 TEST_c_c (ctanh, -2, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1965 TEST_c_c (ctanh, minus_zero, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1966 TEST_c_c (ctanh, -2, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
1968 TEST_c_c (ctanh, plus_infty, nan_value, 1.0, 0.0, IGNORE_ZERO_INF_SIGN);
1969 TEST_c_c (ctanh, minus_infty, nan_value, -1.0, 0.0, IGNORE_ZERO_INF_SIGN);
1971 TEST_c_c (ctanh, nan_value, 0, nan_value, 0.0);
1972 TEST_c_c (ctanh, nan_value, minus_zero, nan_value, minus_zero);
1974 TEST_c_c (ctanh, nan_value, 0.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1975 TEST_c_c (ctanh, nan_value, -4.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
1977 TEST_c_c (ctanh, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1978 TEST_c_c (ctanh, 5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1979 TEST_c_c (ctanh, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1980 TEST_c_c (ctanh, -0.25, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
1982 TEST_c_c (ctanh, nan_value, nan_value, nan_value, nan_value);
1984 TEST_c_c (ctanh, 0, M_PI_4l, 0.0, 1.0);
1986 TEST_c_c (ctanh, 0.7, 1.2, 1.3472197399061191630L, 0.4778641038326365540L);
1987 TEST_c_c (ctanh, -2, -3, -0.9653858790221331242L, 0.0098843750383224937L);
1997 if (errno == ENOSYS)
1998 /* Function not implemented. */
2003 TEST_f_f (erf, 0, 0);
2004 TEST_f_f (erf, minus_zero, minus_zero);
2005 TEST_f_f (erf, plus_infty, 1);
2006 TEST_f_f (erf, minus_infty, -1);
2008 TEST_f_f (erf, 0.7, 0.67780119383741847297L);
2010 TEST_f_f (erf, 1.2, 0.91031397822963538024L);
2011 TEST_f_f (erf, 2.0, 0.99532226501895273416L);
2012 TEST_f_f (erf, 4.1, 0.99999999329997234592L);
2013 TEST_f_f (erf, 27, 1.0L);
2024 if (errno == ENOSYS)
2025 /* Function not implemented. */
2030 TEST_f_f (erfc, plus_infty, 0.0);
2031 TEST_f_f (erfc, minus_infty, 2.0);
2032 TEST_f_f (erfc, 0.0, 1.0);
2033 TEST_f_f (erfc, minus_zero, 1.0);
2035 TEST_f_f (erfc, 0.7, 0.32219880616258152702L);
2037 TEST_f_f (erfc, 1.2, 0.089686021770364619762L);
2038 TEST_f_f (erfc, 2.0, 0.0046777349810472658379L);
2039 TEST_f_f (erfc, 4.1, 0.67000276540848983727e-8L);
2040 TEST_f_f (erfc, 9, 0.41370317465138102381e-36L);
2050 TEST_f_f (exp, 0, 1);
2051 TEST_f_f (exp, minus_zero, 1);
2053 TEST_f_f (exp, plus_infty, plus_infty);
2054 TEST_f_f (exp, minus_infty, 0);
2055 TEST_f_f (exp, 1, M_El);
2057 TEST_f_f (exp, 2, M_E2l);
2058 TEST_f_f (exp, 3, M_E3l);
2059 TEST_f_f (exp, 0.7, 2.0137527074704765216L);
2070 if (errno == ENOSYS)
2071 /* Function not implemented. */
2076 TEST_f_f (exp10, 0, 1);
2077 TEST_f_f (exp10, minus_zero, 1);
2079 TEST_f_f (exp10, plus_infty, plus_infty);
2080 TEST_f_f (exp10, minus_infty, 0);
2081 TEST_f_f (exp10, 3, 1000);
2082 TEST_f_f (exp10, -1, 0.1);
2083 TEST_f_f (exp10, 1e6, plus_infty);
2084 TEST_f_f (exp10, -1e6, 0);
2085 TEST_f_f (exp10, 0.7, 5.0118723362727228500L);
2095 if (errno == ENOSYS)
2096 /* Function not implemented. */
2101 TEST_f_f (exp2, 0, 1);
2102 TEST_f_f (exp2, minus_zero, 1);
2103 TEST_f_f (exp2, plus_infty, plus_infty);
2104 TEST_f_f (exp2, minus_infty, 0);
2105 TEST_f_f (exp2, 10, 1024);
2106 TEST_f_f (exp2, -1, 0.5);
2107 TEST_f_f (exp2, 1e6, plus_infty);
2108 TEST_f_f (exp2, -1e6, 0);
2109 TEST_f_f (exp2, 0.7, 1.6245047927124710452L);
2119 TEST_f_f (expm1, 0, 0);
2120 TEST_f_f (expm1, minus_zero, minus_zero);
2122 TEST_f_f (expm1, plus_infty, plus_infty);
2123 TEST_f_f (expm1, minus_infty, -1);
2125 TEST_f_f (expm1, 1, M_El - 1.0);
2126 TEST_f_f (expm1, 0.7, 1.0137527074704765216L);
2136 TEST_f_f (fabs, 0, 0);
2137 TEST_f_f (fabs, minus_zero, 0);
2139 TEST_f_f (fabs, plus_infty, plus_infty);
2140 TEST_f_f (fabs, minus_infty, plus_infty);
2142 TEST_f_f (fabs, 38.0, 38.0);
2143 TEST_f_f (fabs, -M_El, M_El);
2153 TEST_ff_f (fdim, 0, 0, 0);
2154 TEST_ff_f (fdim, 9, 0, 9);
2155 TEST_ff_f (fdim, 0, 9, 0);
2156 TEST_ff_f (fdim, -9, 0, 0);
2157 TEST_ff_f (fdim, 0, -9, 9);
2159 TEST_ff_f (fdim, plus_infty, 9, plus_infty);
2160 TEST_ff_f (fdim, plus_infty, -9, plus_infty);
2161 TEST_ff_f (fdim, minus_infty, 9, 0);
2162 TEST_ff_f (fdim, minus_infty, -9, 0);
2163 TEST_ff_f (fdim, 9, minus_infty, plus_infty);
2164 TEST_ff_f (fdim, -9, minus_infty, plus_infty);
2165 TEST_ff_f (fdim, 9, plus_infty, 0);
2166 TEST_ff_f (fdim, -9, plus_infty, 0);
2168 TEST_ff_f (fdim, 0, nan_value, nan_value);
2169 TEST_ff_f (fdim, 9, nan_value, nan_value);
2170 TEST_ff_f (fdim, -9, nan_value, nan_value);
2171 TEST_ff_f (fdim, nan_value, 9, nan_value);
2172 TEST_ff_f (fdim, nan_value, -9, nan_value);
2173 TEST_ff_f (fdim, plus_infty, nan_value, nan_value);
2174 TEST_ff_f (fdim, minus_infty, nan_value, nan_value);
2175 TEST_ff_f (fdim, nan_value, plus_infty, nan_value);
2176 TEST_ff_f (fdim, nan_value, minus_infty, nan_value);
2177 TEST_ff_f (fdim, nan_value, nan_value, nan_value);
2187 TEST_f_f (floor, 0.0, 0.0);
2188 TEST_f_f (floor, minus_zero, minus_zero);
2189 TEST_f_f (floor, plus_infty, plus_infty);
2190 TEST_f_f (floor, minus_infty, minus_infty);
2192 TEST_f_f (floor, M_PIl, 3.0);
2193 TEST_f_f (floor, -M_PIl, -4.0);
2203 TEST_fff_f (fma, 1.0, 2.0, 3.0, 5.0);
2204 TEST_fff_f (fma, nan_value, 2.0, 3.0, nan_value);
2205 TEST_fff_f (fma, 1.0, nan_value, 3.0, nan_value);
2206 TEST_fff_f (fma, 1.0, 2.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
2207 TEST_fff_f (fma, plus_infty, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
2208 TEST_fff_f (fma, minus_infty, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
2209 TEST_fff_f (fma, 0.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
2210 TEST_fff_f (fma, 0.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
2211 TEST_fff_f (fma, plus_infty, 0.0, 1.0, nan_value, INVALID_EXCEPTION);
2212 TEST_fff_f (fma, minus_infty, 0.0, 1.0, nan_value, INVALID_EXCEPTION);
2213 TEST_fff_f (fma, 0.0, plus_infty, 1.0, nan_value, INVALID_EXCEPTION);
2214 TEST_fff_f (fma, 0.0, minus_infty, 1.0, nan_value, INVALID_EXCEPTION);
2216 TEST_fff_f (fma, plus_infty, plus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
2217 TEST_fff_f (fma, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
2218 TEST_fff_f (fma, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
2219 TEST_fff_f (fma, minus_infty, minus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
2230 TEST_ff_f (fmax, 0, 0, 0);
2231 TEST_ff_f (fmax, minus_zero, minus_zero, minus_zero);
2232 TEST_ff_f (fmax, 9, 0, 9);
2233 TEST_ff_f (fmax, 0, 9, 9);
2234 TEST_ff_f (fmax, -9, 0, 0);
2235 TEST_ff_f (fmax, 0, -9, 0);
2237 TEST_ff_f (fmax, plus_infty, 9, plus_infty);
2238 TEST_ff_f (fmax, 0, plus_infty, plus_infty);
2239 TEST_ff_f (fmax, -9, plus_infty, plus_infty);
2240 TEST_ff_f (fmax, plus_infty, -9, plus_infty);
2242 TEST_ff_f (fmax, minus_infty, 9, 9);
2243 TEST_ff_f (fmax, minus_infty, -9, -9);
2244 TEST_ff_f (fmax, 9, minus_infty, 9);
2245 TEST_ff_f (fmax, -9, minus_infty, -9);
2247 TEST_ff_f (fmax, 0, nan_value, 0);
2248 TEST_ff_f (fmax, 9, nan_value, 9);
2249 TEST_ff_f (fmax, -9, nan_value, -9);
2250 TEST_ff_f (fmax, nan_value, 0, 0);
2251 TEST_ff_f (fmax, nan_value, 9, 9);
2252 TEST_ff_f (fmax, nan_value, -9, -9);
2253 TEST_ff_f (fmax, plus_infty, nan_value, plus_infty);
2254 TEST_ff_f (fmax, minus_infty, nan_value, minus_infty);
2255 TEST_ff_f (fmax, nan_value, plus_infty, plus_infty);
2256 TEST_ff_f (fmax, nan_value, minus_infty, minus_infty);
2257 TEST_ff_f (fmax, nan_value, nan_value, nan_value);
2268 TEST_ff_f (fmin, 0, 0, 0);
2269 TEST_ff_f (fmin, minus_zero, minus_zero, minus_zero);
2270 TEST_ff_f (fmin, 9, 0, 0);
2271 TEST_ff_f (fmin, 0, 9, 0);
2272 TEST_ff_f (fmin, -9, 0, -9);
2273 TEST_ff_f (fmin, 0, -9, -9);
2275 TEST_ff_f (fmin, plus_infty, 9, 9);
2276 TEST_ff_f (fmin, 9, plus_infty, 9);
2277 TEST_ff_f (fmin, plus_infty, -9, -9);
2278 TEST_ff_f (fmin, -9, plus_infty, -9);
2279 TEST_ff_f (fmin, minus_infty, 9, minus_infty);
2280 TEST_ff_f (fmin, minus_infty, -9, minus_infty);
2281 TEST_ff_f (fmin, 9, minus_infty, minus_infty);
2282 TEST_ff_f (fmin, -9, minus_infty, minus_infty);
2284 TEST_ff_f (fmin, 0, nan_value, 0);
2285 TEST_ff_f (fmin, 9, nan_value, 9);
2286 TEST_ff_f (fmin, -9, nan_value, -9);
2287 TEST_ff_f (fmin, nan_value, 0, 0);
2288 TEST_ff_f (fmin, nan_value, 9, 9);
2289 TEST_ff_f (fmin, nan_value, -9, -9);
2290 TEST_ff_f (fmin, plus_infty, nan_value, plus_infty);
2291 TEST_ff_f (fmin, minus_infty, nan_value, minus_infty);
2292 TEST_ff_f (fmin, nan_value, plus_infty, plus_infty);
2293 TEST_ff_f (fmin, nan_value, minus_infty, minus_infty);
2294 TEST_ff_f (fmin, nan_value, nan_value, nan_value);
2306 /* fmod (+0, y) == +0 for y != 0. */
2307 TEST_ff_f (fmod, 0, 3, 0);
2309 /* fmod (-0, y) == -0 for y != 0. */
2310 TEST_ff_f (fmod, minus_zero, 3, minus_zero);
2312 /* fmod (+inf, y) == NaN plus invalid exception. */
2313 TEST_ff_f (fmod, plus_infty, 3, nan_value, INVALID_EXCEPTION);
2314 /* fmod (-inf, y) == NaN plus invalid exception. */
2315 TEST_ff_f (fmod, minus_infty, 3, nan_value, INVALID_EXCEPTION);
2316 /* fmod (x, +0) == NaN plus invalid exception. */
2317 TEST_ff_f (fmod, 3, 0, nan_value, INVALID_EXCEPTION);
2318 /* fmod (x, -0) == NaN plus invalid exception. */
2319 TEST_ff_f (fmod, 3, minus_zero, nan_value, INVALID_EXCEPTION);
2321 /* fmod (x, +inf) == x for x not infinite. */
2322 TEST_ff_f (fmod, 3.0, plus_infty, 3.0);
2323 /* fmod (x, -inf) == x for x not infinite. */
2324 TEST_ff_f (fmod, 3.0, minus_infty, 3.0);
2326 TEST_ff_f (fmod, 6.5, 2.3, 1.9);
2327 TEST_ff_f (fmod, -6.5, 2.3, -1.9);
2328 TEST_ff_f (fmod, 6.5, -2.3, 1.9);
2329 TEST_ff_f (fmod, -6.5, -2.3, -1.9);
2335 fpclassify_test (void)
2339 TEST_f_i (fpclassify, nan_value, FP_NAN);
2340 TEST_f_i (fpclassify, plus_infty, FP_INFINITE);
2341 TEST_f_i (fpclassify, minus_infty, FP_INFINITE);
2342 TEST_f_i (fpclassify, plus_zero, FP_ZERO);
2343 TEST_f_i (fpclassify, minus_zero, FP_ZERO);
2344 TEST_f_i (fpclassify, 1000, FP_NORMAL);
2357 TEST_fI_f1 (frexp, plus_infty, plus_infty, IGNORE);
2358 TEST_fI_f1 (frexp, minus_infty, minus_infty, IGNORE);
2359 TEST_fI_f1 (frexp, nan_value, nan_value, IGNORE);
2361 TEST_fI_f1 (frexp, 0.0, 0.0, 0.0);
2362 TEST_fI_f1 (frexp, minus_zero, minus_zero, 0.0);
2364 TEST_fI_f1 (frexp, 12.8, 0.8, 4);
2365 TEST_fI_f1 (frexp, -27.34, -0.854375, 5);
2377 if (errno == ENOSYS)
2378 /* Function not implemented. */
2380 feclearexcept (FE_ALL_EXCEPT);
2384 TEST_f_f (gamma, plus_infty, plus_infty);
2385 TEST_f_f (gamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2386 TEST_f_f (gamma, -3, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2387 TEST_f_f (gamma, minus_infty, plus_infty);
2389 TEST_f_f1 (gamma, 1, 0, 1);
2390 TEST_f_f1 (gamma, 3, M_LN2l, 1);
2392 TEST_f_f1 (gamma, 0.5, M_LOG_SQRT_PIl, 1);
2393 TEST_f_f1 (gamma, -0.5, M_LOG_2_SQRT_PIl, -1);
2403 TEST_ff_f (hypot, plus_infty, 1, plus_infty, IGNORE_ZERO_INF_SIGN);
2404 TEST_ff_f (hypot, minus_infty, 1, plus_infty, IGNORE_ZERO_INF_SIGN);
2406 TEST_ff_f (hypot, plus_infty, nan_value, plus_infty);
2407 TEST_ff_f (hypot, minus_infty, nan_value, plus_infty);
2408 TEST_ff_f (hypot, nan_value, plus_infty, plus_infty);
2409 TEST_ff_f (hypot, nan_value, minus_infty, plus_infty);
2411 TEST_ff_f (hypot, nan_value, nan_value, nan_value);
2413 /* hypot (x,y) == hypot (+-x, +-y) */
2414 TEST_ff_f (hypot, 0.7, 12.4, 12.41974234837422060118L);
2415 TEST_ff_f (hypot, -0.7, 12.4, 12.41974234837422060118L);
2416 TEST_ff_f (hypot, 0.7, -12.4, 12.41974234837422060118L);
2417 TEST_ff_f (hypot, -0.7, -12.4, 12.41974234837422060118L);
2418 TEST_ff_f (hypot, 12.4, 0.7, 12.41974234837422060118L);
2419 TEST_ff_f (hypot, -12.4, 0.7, 12.41974234837422060118L);
2420 TEST_ff_f (hypot, 12.4, -0.7, 12.41974234837422060118L);
2421 TEST_ff_f (hypot, -12.4, -0.7, 12.41974234837422060118L);
2423 /* hypot (x,0) == fabs (x) */
2424 TEST_ff_f (hypot, 0.7, 0, 0.7);
2425 TEST_ff_f (hypot, -0.7, 0, 0.7);
2426 TEST_ff_f (hypot, -5.7e7, 0, 5.7e7);
2428 TEST_ff_f (hypot, 0.7, 1.2, 1.3892443989449804508L);
2439 TEST_f_i (ilogb, 1, 0);
2440 TEST_f_i (ilogb, M_El, 1);
2441 TEST_f_i (ilogb, 1024, 10);
2442 TEST_f_i (ilogb, -2000, 10);
2444 /* XXX We have a problem here: the standard does not tell us whether
2445 exceptions are allowed/required. ignore them for now. */
2447 TEST_f_i (ilogb, 0.0, FP_ILOGB0, EXCEPTIONS_OK);
2448 TEST_f_i (ilogb, nan_value, FP_ILOGBNAN, EXCEPTIONS_OK);
2454 isfinite_test (void)
2458 TEST_f_b (isfinite, 0, 1);
2459 TEST_f_b (isfinite, minus_zero, 1);
2460 TEST_f_b (isfinite, 10, 1);
2461 TEST_f_b (isfinite, plus_infty, 0);
2462 TEST_f_b (isfinite, minus_infty, 0);
2463 TEST_f_b (isfinite, nan_value, 0);
2469 isnormal_test (void)
2473 TEST_f_b (isnormal, 0, 0);
2474 TEST_f_b (isnormal, minus_zero, 0);
2475 TEST_f_b (isnormal, 10, 1);
2476 TEST_f_b (isnormal, plus_infty, 0);
2477 TEST_f_b (isnormal, minus_infty, 0);
2478 TEST_f_b (isnormal, nan_value, 0);
2488 if (errno == ENOSYS)
2489 /* Function not implemented. */
2494 /* j0 is the Bessel function of the first kind of order 0 */
2495 TEST_f_f (j0, nan_value, nan_value);
2496 TEST_f_f (j0, plus_infty, 0);
2497 TEST_f_f (j0, -1.0, 0.76519768655796655145);
2498 TEST_f_f (j0, 0.0, 1.0);
2499 TEST_f_f (j0, 0.1, 0.99750156206604003228);
2500 TEST_f_f (j0, 0.7, 0.88120088860740528084);
2501 TEST_f_f (j0, 1.0, 0.76519768655796655145);
2502 TEST_f_f (j0, 1.5, 0.51182767173591812875);
2503 TEST_f_f (j0, 2.0, 0.22389077914123566805);
2504 TEST_f_f (j0, 8.0, 0.17165080713755390609);
2505 TEST_f_f (j0, 10.0, -0.24593576445134833520);
2516 if (errno == ENOSYS)
2517 /* Function not implemented. */
2520 /* j1 is the Bessel function of the first kind of order 1 */
2524 TEST_f_f (j1, nan_value, nan_value);
2525 TEST_f_f (j1, plus_infty, 0);
2527 TEST_f_f (j1, -1.0, -0.44005058574493351596);
2528 TEST_f_f (j1, 0.0, 0.0);
2529 TEST_f_f (j1, 0.1, 0.049937526036241997556);
2530 TEST_f_f (j1, 0.7, 0.32899574154005894785);
2531 TEST_f_f (j1, 1.0, 0.44005058574493351596);
2532 TEST_f_f (j1, 1.5, 0.55793650791009964199);
2533 TEST_f_f (j1, 2.0, 0.57672480775687338720);
2534 TEST_f_f (j1, 8.0, 0.23463634685391462438);
2535 TEST_f_f (j1, 10.0, 0.043472746168861436670);
2545 if (errno == ENOSYS)
2546 /* Function not implemented. */
2549 /* jn is the Bessel function of the first kind of order n. */
2552 /* jn (0, x) == j0 (x) */
2553 TEST_ff_f (jn, 0, nan_value, nan_value);
2554 TEST_ff_f (jn, 0, plus_infty, 0);
2555 TEST_ff_f (jn, 0, -1.0, 0.76519768655796655145);
2556 TEST_ff_f (jn, 0, 0.0, 1.0);
2557 TEST_ff_f (jn, 0, 0.1, 0.99750156206604003228);
2558 TEST_ff_f (jn, 0, 0.7, 0.88120088860740528084);
2559 TEST_ff_f (jn, 0, 1.0, 0.76519768655796655145);
2560 TEST_ff_f (jn, 0, 1.5, 0.51182767173591812875);
2561 TEST_ff_f (jn, 0, 2.0, 0.22389077914123566805);
2562 TEST_ff_f (jn, 0, 8.0, 0.17165080713755390609);
2563 TEST_ff_f (jn, 0, 10.0, -0.24593576445134833520);
2565 /* jn (1, x) == j1 (x) */
2566 TEST_ff_f (jn, 1, nan_value, nan_value);
2567 TEST_ff_f (jn, 1, plus_infty, 0);
2569 TEST_ff_f (jn, 1, -1.0, -0.44005058574493351596);
2570 TEST_ff_f (jn, 1, 0.0, 0.0);
2571 TEST_ff_f (jn, 1, 0.1, 0.049937526036241997556);
2572 TEST_ff_f (jn, 1, 0.7, 0.32899574154005894785);
2573 TEST_ff_f (jn, 1, 1.0, 0.44005058574493351596);
2574 TEST_ff_f (jn, 1, 1.5, 0.55793650791009964199);
2575 TEST_ff_f (jn, 1, 2.0, 0.57672480775687338720);
2576 TEST_ff_f (jn, 1, 8.0, 0.23463634685391462438);
2577 TEST_ff_f (jn, 1, 10.0, 0.043472746168861436670);
2580 TEST_ff_f (jn, 3, nan_value, nan_value);
2581 TEST_ff_f (jn, 3, plus_infty, 0);
2583 TEST_ff_f (jn, 3, -1.0, -0.019563353982668405919);
2584 TEST_ff_f (jn, 3, 0.0, 0.0);
2585 TEST_ff_f (jn, 3, 0.1, 0.000020820315754756261429);
2586 TEST_ff_f (jn, 3, 0.7, 0.0069296548267508408077);
2587 TEST_ff_f (jn, 3, 1.0, 0.019563353982668405919);
2588 TEST_ff_f (jn, 3, 2.0, 0.12894324947440205110);
2589 TEST_ff_f (jn, 3, 10.0, 0.058379379305186812343);
2592 TEST_ff_f (jn, 10, nan_value, nan_value);
2593 TEST_ff_f (jn, 10, plus_infty, 0);
2595 TEST_ff_f (jn, 10, -1.0, 0.26306151236874532070e-9);
2596 TEST_ff_f (jn, 10, 0.0, 0.0);
2597 TEST_ff_f (jn, 10, 0.1, 0.26905328954342155795e-19);
2598 TEST_ff_f (jn, 10, 0.7, 0.75175911502153953928e-11);
2599 TEST_ff_f (jn, 10, 1.0, 0.26306151236874532070e-9);
2600 TEST_ff_f (jn, 10, 2.0, 0.25153862827167367096e-6);
2601 TEST_ff_f (jn, 10, 10.0, 0.20748610663335885770);
2610 TEST_ff_f (ldexp, 0, 0, 0);
2611 TEST_ff_f (ldexp, minus_zero, 0, minus_zero);
2613 TEST_ff_f (ldexp, plus_infty, 1, plus_infty);
2614 TEST_ff_f (ldexp, minus_infty, 1, minus_infty);
2615 TEST_ff_f (ldexp, nan_value, 1, nan_value);
2617 TEST_ff_f (ldexp, 0.8L, 4, 12.8L);
2618 TEST_ff_f (ldexp, -0.854375L, 5, -27.34L);
2620 /* ldexp (x, 0) == x. */
2621 TEST_ff_f (ldexp, 1.0L, 0L, 1.0L);
2629 if (errno == ENOSYS)
2630 /* Function not implemented. */
2632 feclearexcept (FE_ALL_EXCEPT);
2636 TEST_f_f (lgamma, plus_infty, plus_infty);
2637 TEST_f_f (lgamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2639 /* lgamma (x) == +inf plus divide by zero exception for integer x <= 0. */
2640 TEST_f_f (lgamma, -3, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2641 TEST_f_f (lgamma, minus_infty, plus_infty);
2643 TEST_f_f1 (lgamma, 1, 0, 1);
2645 TEST_f_f1 (lgamma, 3, M_LN2l, 1);
2647 TEST_f_f1 (lgamma, 0.5, M_LOG_SQRT_PIl, 1);
2648 TEST_f_f1 (lgamma, -0.5, M_LOG_2_SQRT_PIl, -1);
2649 TEST_f_f1 (lgamma, 0.7, 0.26086724653166651439L, 1);
2650 TEST_f_f1 (lgamma, 1.2, -0.853740900033158497197e-1L, 1);
2658 /* XXX this test is incomplete. We need to have a way to specifiy
2659 the rounding method and test the critical cases. So far, only
2660 unproblematic numbers are tested. */
2664 TEST_f_l (lrint, 0.0, 0);
2665 TEST_f_l (lrint, minus_zero, 0);
2666 TEST_f_l (lrint, 0.2, 0);
2667 TEST_f_l (lrint, -0.2, 0);
2669 TEST_f_l (lrint, 1.4, 1);
2670 TEST_f_l (lrint, -1.4, -1);
2672 TEST_f_l (lrint, 8388600.3, 8388600);
2673 TEST_f_l (lrint, -8388600.3, -8388600);
2681 /* XXX this test is incomplete. We need to have a way to specifiy
2682 the rounding method and test the critical cases. So far, only
2683 unproblematic numbers are tested. */
2687 TEST_f_L (llrint, 0.0, 0);
2688 TEST_f_L (llrint, minus_zero, 0);
2689 TEST_f_L (llrint, 0.2, 0);
2690 TEST_f_L (llrint, -0.2, 0);
2692 TEST_f_L (llrint, 1.4, 1);
2693 TEST_f_L (llrint, -1.4, -1);
2695 TEST_f_L (llrint, 8388600.3, 8388600);
2696 TEST_f_L (llrint, -8388600.3, -8388600);
2698 /* Test boundary conditions. */
2700 TEST_f_L (llrint, 2097151.0,2097151LL);
2702 TEST_f_L (llrint, 8388608.0, 8388608LL);
2704 TEST_f_L (llrint, 16777216.0, 16777216LL);
2706 TEST_f_L (llrint, 2199023255552.0, 2199023255552LL);
2708 TEST_f_L (llrint, 4398046511104.0, 4398046511104LL);
2709 /* 0x10000000000000 */
2710 TEST_f_L (llrint, 4503599627370496.0, 4503599627370496LL);
2711 /* 0x10000080000000 */
2712 TEST_f_L (llrint, 4503601774854144.0, 4503601774854144LL);
2713 /* 0x20000000000000 */
2714 TEST_f_L (llrint, 9007199254740992.0, 9007199254740992LL);
2715 /* 0x80000000000000 */
2716 TEST_f_L (llrint, 36028797018963968.0, 36028797018963968LL);
2717 /* 0x100000000000000 */
2718 TEST_f_L (llrint, 72057594037927936.0, 72057594037927936LL);
2728 TEST_f_f (log, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2729 TEST_f_f (log, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2731 TEST_f_f (log, 1, 0);
2733 TEST_f_f (log, -1, nan_value, INVALID_EXCEPTION);
2734 TEST_f_f (log, plus_infty, plus_infty);
2736 TEST_f_f (log, M_El, 1);
2737 TEST_f_f (log, 1.0 / M_El, -1);
2738 TEST_f_f (log, 2, M_LN2l);
2739 TEST_f_f (log, 10, M_LN10l);
2740 TEST_f_f (log, 0.7, -0.35667494393873237891L);
2751 TEST_f_f (log10, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2752 TEST_f_f (log10, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2754 TEST_f_f (log10, 1, 0);
2756 /* log10 (x) == NaN plus invalid exception if x < 0. */
2757 TEST_f_f (log10, -1, nan_value, INVALID_EXCEPTION);
2759 TEST_f_f (log10, plus_infty, plus_infty);
2761 TEST_f_f (log10, 0.1L, -1);
2762 TEST_f_f (log10, 10.0, 1);
2763 TEST_f_f (log10, 100.0, 2);
2764 TEST_f_f (log10, 10000.0, 4);
2765 TEST_f_f (log10, M_El, M_LOG10El);
2766 TEST_f_f (log10, 0.7, -0.15490195998574316929L);
2777 TEST_f_f (log1p, 0, 0);
2778 TEST_f_f (log1p, minus_zero, minus_zero);
2780 TEST_f_f (log1p, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2781 TEST_f_f (log1p, -2, nan_value, INVALID_EXCEPTION);
2783 TEST_f_f (log1p, plus_infty, plus_infty);
2785 TEST_f_f (log1p, M_El - 1.0, 1);
2787 TEST_f_f (log1p, -0.3, -0.35667494393873237891L);
2798 TEST_f_f (log2, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2799 TEST_f_f (log2, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2801 TEST_f_f (log2, 1, 0);
2803 TEST_f_f (log2, -1, nan_value, INVALID_EXCEPTION);
2805 TEST_f_f (log2, plus_infty, plus_infty);
2807 TEST_f_f (log2, M_El, M_LOG2El);
2808 TEST_f_f (log2, 2.0, 1);
2809 TEST_f_f (log2, 16.0, 4);
2810 TEST_f_f (log2, 256.0, 8);
2811 TEST_f_f (log2, 0.7, -0.51457317282975824043L);
2822 TEST_f_f (logb, plus_infty, plus_infty);
2823 TEST_f_f (logb, minus_infty, plus_infty);
2825 TEST_f_f (logb, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2827 TEST_f_f (logb, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
2829 TEST_f_f (logb, 1, 0);
2830 TEST_f_f (logb, M_El, 1);
2831 TEST_f_f (logb, 1024, 10);
2832 TEST_f_f (logb, -2000, 10);
2842 TEST_f_l (lround, 0, 0);
2843 TEST_f_l (lround, minus_zero, 0);
2844 TEST_f_l (lround, 0.2, 0.0);
2845 TEST_f_l (lround, -0.2, 0);
2846 TEST_f_l (lround, 0.5, 1);
2847 TEST_f_l (lround, -0.5, -1);
2848 TEST_f_l (lround, 0.8, 1);
2849 TEST_f_l (lround, -0.8, -1);
2850 TEST_f_l (lround, 1.5, 2);
2851 TEST_f_l (lround, -1.5, -2);
2852 TEST_f_l (lround, 22514.5, 22515);
2853 TEST_f_l (lround, -22514.5, -22515);
2855 TEST_f_l (lround, 2097152.5, 2097153);
2856 TEST_f_l (lround, -2097152.5, -2097153);
2867 TEST_f_L (llround, 0, 0);
2868 TEST_f_L (llround, minus_zero, 0);
2869 TEST_f_L (llround, 0.2, 0.0);
2870 TEST_f_L (llround, -0.2, 0);
2871 TEST_f_L (llround, 0.5, 1);
2872 TEST_f_L (llround, -0.5, -1);
2873 TEST_f_L (llround, 0.8, 1);
2874 TEST_f_L (llround, -0.8, -1);
2875 TEST_f_L (llround, 1.5, 2);
2876 TEST_f_L (llround, -1.5, -2);
2877 TEST_f_L (llround, 22514.5, 22515);
2878 TEST_f_L (llround, -22514.5, -22515);
2880 TEST_f_L (llround, 2097152.5, 2097153);
2881 TEST_f_L (llround, -2097152.5, -2097153);
2882 TEST_f_L (llround, 34359738368.5, 34359738369ll);
2883 TEST_f_L (llround, -34359738368.5, -34359738369ll);
2886 /* Test boundary conditions. */
2888 TEST_f_L (llround, 2097151.0, 2097151LL);
2890 TEST_f_L (llround, 8388608.0, 8388608LL);
2892 TEST_f_L (llround, 16777216.0, 16777216LL);
2894 TEST_f_L (llround, 2199023255552.0, 2199023255552LL);
2896 TEST_f_L (llround, 4398046511104.0, 4398046511104LL);
2897 /* 0x10000000000000 */
2898 TEST_f_L (llround, 4503599627370496.0, 4503599627370496LL);
2899 /* 0x10000080000000 */
2900 TEST_f_L (llrint, 4503601774854144.0, 4503601774854144LL);
2901 /* 0x20000000000000 */
2902 TEST_f_L (llround, 9007199254740992.0, 9007199254740992LL);
2903 /* 0x80000000000000 */
2904 TEST_f_L (llround, 36028797018963968.0, 36028797018963968LL);
2905 /* 0x100000000000000 */
2906 TEST_f_L (llround, 72057594037927936.0, 72057594037927936LL);
2918 TEST_fF_f1 (modf, plus_infty, 0, plus_infty);
2919 TEST_fF_f1 (modf, minus_infty, minus_zero, minus_infty);
2920 TEST_fF_f1 (modf, nan_value, nan_value, nan_value);
2921 TEST_fF_f1 (modf, 0, 0, 0);
2922 TEST_fF_f1 (modf, 1.5, 0.5, 1);
2923 TEST_fF_f1 (modf, 2.5, 0.5, 2);
2924 TEST_fF_f1 (modf, -2.5, -0.5, -2);
2925 TEST_fF_f1 (modf, 20, 0, 20);
2926 TEST_fF_f1 (modf, 21, 0, 21);
2927 TEST_fF_f1 (modf, 89.5, 0.5, 89);
2934 nearbyint_test (void)
2938 TEST_f_f (nearbyint, 0.0, 0.0);
2939 TEST_f_f (nearbyint, minus_zero, minus_zero);
2940 TEST_f_f (nearbyint, plus_infty, plus_infty);
2941 TEST_f_f (nearbyint, minus_infty, minus_infty);
2947 nextafter_test (void)
2952 TEST_ff_f (nextafter, 0, 0, 0);
2953 TEST_ff_f (nextafter, minus_zero, 0, 0);
2954 TEST_ff_f (nextafter, 0, minus_zero, minus_zero);
2955 TEST_ff_f (nextafter, minus_zero, minus_zero, minus_zero);
2957 TEST_ff_f (nextafter, 9, 9, 9);
2958 TEST_ff_f (nextafter, -9, -9, -9);
2959 TEST_ff_f (nextafter, plus_infty, plus_infty, plus_infty);
2960 TEST_ff_f (nextafter, minus_infty, minus_infty, minus_infty);
2962 TEST_ff_f (nextafter, nan_value, 1.1, nan_value);
2963 TEST_ff_f (nextafter, 1.1, nan_value, nan_value);
2964 TEST_ff_f (nextafter, nan_value, nan_value, nan_value);
2966 /* XXX We need the hexadecimal FP number representation here for further
2974 nexttoward_test (void)
2977 TEST_ff_f (nexttoward, 0, 0, 0);
2978 TEST_ff_f (nexttoward, minus_zero, 0, 0);
2979 TEST_ff_f (nexttoward, 0, minus_zero, minus_zero);
2980 TEST_ff_f (nexttoward, minus_zero, minus_zero, minus_zero);
2982 TEST_ff_f (nexttoward, 9, 9, 9);
2983 TEST_ff_f (nexttoward, -9, -9, -9);
2984 TEST_ff_f (nexttoward, plus_infty, plus_infty, plus_infty);
2985 TEST_ff_f (nexttoward, minus_infty, minus_infty, minus_infty);
2987 TEST_ff_f (nexttoward, nan_value, 1.1, nan_value);
2988 TEST_ff_f (nexttoward, 1.1, nan_value, nan_value);
2989 TEST_ff_f (nexttoward, nan_value, nan_value, nan_value);
2991 /* XXX We need the hexadecimal FP number representation here for further
3004 TEST_ff_f (pow, 0, 0, 1);
3005 TEST_ff_f (pow, 0, minus_zero, 1);
3006 TEST_ff_f (pow, minus_zero, 0, 1);
3007 TEST_ff_f (pow, minus_zero, minus_zero, 1);
3009 TEST_ff_f (pow, 10, 0, 1);
3010 TEST_ff_f (pow, 10, minus_zero, 1);
3011 TEST_ff_f (pow, -10, 0, 1);
3012 TEST_ff_f (pow, -10, minus_zero, 1);
3014 TEST_ff_f (pow, nan_value, 0, 1);
3015 TEST_ff_f (pow, nan_value, minus_zero, 1);
3018 TEST_ff_f (pow, 1.1, plus_infty, plus_infty);
3019 TEST_ff_f (pow, plus_infty, plus_infty, plus_infty);
3020 TEST_ff_f (pow, -1.1, plus_infty, plus_infty);
3021 TEST_ff_f (pow, minus_infty, plus_infty, plus_infty);
3023 TEST_ff_f (pow, 0.9L, plus_infty, 0);
3024 TEST_ff_f (pow, 1e-7L, plus_infty, 0);
3025 TEST_ff_f (pow, -0.9L, plus_infty, 0);
3026 TEST_ff_f (pow, -1e-7L, plus_infty, 0);
3028 TEST_ff_f (pow, 1.1, minus_infty, 0);
3029 TEST_ff_f (pow, plus_infty, minus_infty, 0);
3030 TEST_ff_f (pow, -1.1, minus_infty, 0);
3031 TEST_ff_f (pow, minus_infty, minus_infty, 0);
3033 TEST_ff_f (pow, 0.9L, minus_infty, plus_infty);
3034 TEST_ff_f (pow, 1e-7L, minus_infty, plus_infty);
3035 TEST_ff_f (pow, -0.9L, minus_infty, plus_infty);
3036 TEST_ff_f (pow, -1e-7L, minus_infty, plus_infty);
3038 TEST_ff_f (pow, plus_infty, 1e-7L, plus_infty);
3039 TEST_ff_f (pow, plus_infty, 1, plus_infty);
3040 TEST_ff_f (pow, plus_infty, 1e7L, plus_infty);
3042 TEST_ff_f (pow, plus_infty, -1e-7L, 0);
3043 TEST_ff_f (pow, plus_infty, -1, 0);
3044 TEST_ff_f (pow, plus_infty, -1e7L, 0);
3046 TEST_ff_f (pow, minus_infty, 1, minus_infty);
3047 TEST_ff_f (pow, minus_infty, 11, minus_infty);
3048 TEST_ff_f (pow, minus_infty, 1001, minus_infty);
3050 TEST_ff_f (pow, minus_infty, 2, plus_infty);
3051 TEST_ff_f (pow, minus_infty, 12, plus_infty);
3052 TEST_ff_f (pow, minus_infty, 1002, plus_infty);
3053 TEST_ff_f (pow, minus_infty, 0.1, plus_infty);
3054 TEST_ff_f (pow, minus_infty, 1.1, plus_infty);
3055 TEST_ff_f (pow, minus_infty, 11.1, plus_infty);
3056 TEST_ff_f (pow, minus_infty, 1001.1, plus_infty);
3058 TEST_ff_f (pow, minus_infty, -1, minus_zero);
3059 TEST_ff_f (pow, minus_infty, -11, minus_zero);
3060 TEST_ff_f (pow, minus_infty, -1001, minus_zero);
3062 TEST_ff_f (pow, minus_infty, -2, 0);
3063 TEST_ff_f (pow, minus_infty, -12, 0);
3064 TEST_ff_f (pow, minus_infty, -1002, 0);
3065 TEST_ff_f (pow, minus_infty, -0.1, 0);
3066 TEST_ff_f (pow, minus_infty, -1.1, 0);
3067 TEST_ff_f (pow, minus_infty, -11.1, 0);
3068 TEST_ff_f (pow, minus_infty, -1001.1, 0);
3070 TEST_ff_f (pow, nan_value, nan_value, nan_value);
3071 TEST_ff_f (pow, 0, nan_value, nan_value);
3072 TEST_ff_f (pow, 1, nan_value, nan_value);
3073 TEST_ff_f (pow, -1, nan_value, nan_value);
3074 TEST_ff_f (pow, nan_value, 1, nan_value);
3075 TEST_ff_f (pow, nan_value, -1, nan_value);
3077 /* pow (x, NaN) == NaN. */
3078 TEST_ff_f (pow, 3.0, nan_value, nan_value);
3080 TEST_ff_f (pow, 1, plus_infty, nan_value, INVALID_EXCEPTION);
3081 TEST_ff_f (pow, -1, plus_infty, nan_value, INVALID_EXCEPTION);
3082 TEST_ff_f (pow, 1, minus_infty, nan_value, INVALID_EXCEPTION);
3083 TEST_ff_f (pow, -1, minus_infty, nan_value, INVALID_EXCEPTION);
3085 TEST_ff_f (pow, -0.1, 1.1, nan_value, INVALID_EXCEPTION);
3086 TEST_ff_f (pow, -0.1, -1.1, nan_value, INVALID_EXCEPTION);
3087 TEST_ff_f (pow, -10.1, 1.1, nan_value, INVALID_EXCEPTION);
3088 TEST_ff_f (pow, -10.1, -1.1, nan_value, INVALID_EXCEPTION);
3090 TEST_ff_f (pow, 0, -1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3091 TEST_ff_f (pow, 0, -11, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3092 TEST_ff_f (pow, minus_zero, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3093 TEST_ff_f (pow, minus_zero, -11, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3095 TEST_ff_f (pow, 0, -2, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3096 TEST_ff_f (pow, 0, -11.1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3097 TEST_ff_f (pow, minus_zero, -2, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3098 TEST_ff_f (pow, minus_zero, -11.1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
3101 TEST_ff_f (pow, 0, 1, 0);
3102 TEST_ff_f (pow, 0, 11, 0);
3104 TEST_ff_f (pow, minus_zero, 1, minus_zero);
3105 TEST_ff_f (pow, minus_zero, 11, minus_zero);
3108 TEST_ff_f (pow, 0, 2, 0);
3109 TEST_ff_f (pow, 0, 11.1, 0);
3112 TEST_ff_f (pow, minus_zero, 2, 0);
3113 TEST_ff_f (pow, minus_zero, 11.1, 0);
3115 /* pow (x, +inf) == +inf for |x| > 1. */
3116 TEST_ff_f (pow, 1.5, plus_infty, plus_infty);
3118 /* pow (x, +inf) == +0 for |x| < 1. */
3119 TEST_ff_f (pow, 0.5, plus_infty, 0.0);
3121 /* pow (x, -inf) == +0 for |x| > 1. */
3122 TEST_ff_f (pow, 1.5, minus_infty, 0.0);
3124 /* pow (x, -inf) == +inf for |x| < 1. */
3125 TEST_ff_f (pow, 0.5, minus_infty, plus_infty);
3127 /* pow (+inf, y) == +inf for y > 0. */
3128 TEST_ff_f (pow, plus_infty, 2, plus_infty);
3130 /* pow (+inf, y) == +0 for y < 0. */
3131 TEST_ff_f (pow, plus_infty, -1, 0.0);
3133 /* pow (-inf, y) == -inf for y an odd integer > 0. */
3134 TEST_ff_f (pow, minus_infty, 27, minus_infty);
3136 /* pow (-inf, y) == +inf for y > 0 and not an odd integer. */
3137 TEST_ff_f (pow, minus_infty, 28, plus_infty);
3139 /* pow (-inf, y) == -0 for y an odd integer < 0. */
3140 TEST_ff_f (pow, minus_infty, -3, minus_zero);
3141 /* pow (-inf, y) == +0 for y < 0 and not an odd integer. */
3142 TEST_ff_f (pow, minus_infty, -2.0, 0.0);
3144 /* pow (+0, y) == +0 for y an odd integer > 0. */
3145 TEST_ff_f (pow, 0.0, 27, 0.0);
3147 /* pow (-0, y) == -0 for y an odd integer > 0. */
3148 TEST_ff_f (pow, minus_zero, 27, minus_zero);
3150 /* pow (+0, y) == +0 for y > 0 and not an odd integer. */
3151 TEST_ff_f (pow, 0.0, 4, 0.0);
3153 /* pow (-0, y) == +0 for y > 0 and not an odd integer. */
3154 TEST_ff_f (pow, minus_zero, 4, 0.0);
3156 TEST_ff_f (pow, 0.7, 1.2, 0.65180494056638638188L);
3159 TEST_ff_f (pow, -7.49321e+133, -9.80818e+16, 0);
3166 remainder_test (void)
3171 TEST_ff_f (remainder, 1, 0, nan_value, INVALID_EXCEPTION);
3172 TEST_ff_f (remainder, 1, minus_zero, nan_value, INVALID_EXCEPTION);
3173 TEST_ff_f (remainder, plus_infty, 1, nan_value, INVALID_EXCEPTION);
3174 TEST_ff_f (remainder, minus_infty, 1, nan_value, INVALID_EXCEPTION);
3176 TEST_ff_f (remainder, 1.625, 1.0, -0.375);
3177 TEST_ff_f (remainder, -1.625, 1.0, 0.375);
3178 TEST_ff_f (remainder, 1.625, -1.0, -0.375);
3179 TEST_ff_f (remainder, -1.625, -1.0, 0.375);
3180 TEST_ff_f (remainder, 5.0, 2.0, 1.0);
3181 TEST_ff_f (remainder, 3.0, 2.0, -1.0);
3194 TEST_ffI_f1 (remquo, 1, 0, nan_value, IGNORE, INVALID_EXCEPTION);
3195 TEST_ffI_f1 (remquo, 1, minus_zero, nan_value, IGNORE, INVALID_EXCEPTION);
3196 TEST_ffI_f1 (remquo, plus_infty, 1, nan_value, IGNORE, INVALID_EXCEPTION);
3197 TEST_ffI_f1 (remquo, minus_infty, 1, nan_value, IGNORE, INVALID_EXCEPTION);
3199 TEST_ffI_f1 (remquo, 1.625, 1.0, -0.375, 2);
3200 TEST_ffI_f1 (remquo, -1.625, 1.0, 0.375, -2);
3201 TEST_ffI_f1 (remquo, 1.625, -1.0, -0.375, -2);
3202 TEST_ffI_f1 (remquo, -1.625, -1.0, 0.375, 2);
3204 TEST_ffI_f1 (remquo, 5, 2, 1, 2);
3205 TEST_ffI_f1 (remquo, 3, 2, -1, 2);
3215 TEST_f_f (rint, 0.0, 0.0);
3216 TEST_f_f (rint, minus_zero, minus_zero);
3217 TEST_f_f (rint, plus_infty, plus_infty);
3218 TEST_f_f (rint, minus_infty, minus_infty);
3228 TEST_f_f (round, 0, 0);
3229 TEST_f_f (round, minus_zero, minus_zero);
3230 TEST_f_f (round, 0.2, 0.0);
3231 TEST_f_f (round, -0.2, minus_zero);
3232 TEST_f_f (round, 0.5, 1.0);
3233 TEST_f_f (round, -0.5, -1.0);
3234 TEST_f_f (round, 0.8, 1.0);
3235 TEST_f_f (round, -0.8, -1.0);
3236 TEST_f_f (round, 1.5, 2.0);
3237 TEST_f_f (round, -1.5, -2.0);
3238 TEST_f_f (round, 2097152.5, 2097153);
3239 TEST_f_f (round, -2097152.5, -2097153);
3251 TEST_ff_f (scalb, 2.0, 0.5, nan_value);
3252 TEST_ff_f (scalb, 3.0, -2.5, nan_value);
3254 TEST_ff_f (scalb, 0, nan_value, nan_value);
3255 TEST_ff_f (scalb, 1, nan_value, nan_value);
3257 TEST_ff_f (scalb, 1, 0, 1);
3258 TEST_ff_f (scalb, -1, 0, -1);
3260 TEST_ff_f (scalb, 0, plus_infty, nan_value, INVALID_EXCEPTION);
3261 TEST_ff_f (scalb, minus_zero, plus_infty, nan_value, INVALID_EXCEPTION);
3263 TEST_ff_f (scalb, 0, 2, 0);
3264 TEST_ff_f (scalb, minus_zero, -4, minus_zero);
3265 TEST_ff_f (scalb, 0, 0, 0);
3266 TEST_ff_f (scalb, minus_zero, 0, minus_zero);
3267 TEST_ff_f (scalb, 0, -1, 0);
3268 TEST_ff_f (scalb, minus_zero, -10, minus_zero);
3269 TEST_ff_f (scalb, 0, minus_infty, 0);
3270 TEST_ff_f (scalb, minus_zero, minus_infty, minus_zero);
3272 TEST_ff_f (scalb, plus_infty, -1, plus_infty);
3273 TEST_ff_f (scalb, minus_infty, -10, minus_infty);
3274 TEST_ff_f (scalb, plus_infty, 0, plus_infty);
3275 TEST_ff_f (scalb, minus_infty, 0, minus_infty);
3276 TEST_ff_f (scalb, plus_infty, 2, plus_infty);
3277 TEST_ff_f (scalb, minus_infty, 100, minus_infty);
3279 TEST_ff_f (scalb, 0.1, minus_infty, 0.0);
3280 TEST_ff_f (scalb, -0.1, minus_infty, minus_zero);
3282 TEST_ff_f (scalb, 1, plus_infty, plus_infty);
3283 TEST_ff_f (scalb, -1, plus_infty, minus_infty);
3284 TEST_ff_f (scalb, plus_infty, plus_infty, plus_infty);
3285 TEST_ff_f (scalb, minus_infty, plus_infty, minus_infty);
3287 TEST_ff_f (scalb, plus_infty, minus_infty, nan_value);
3288 TEST_ff_f (scalb, minus_infty, minus_infty, nan_value);
3290 TEST_ff_f (scalb, nan_value, 1, nan_value);
3291 TEST_ff_f (scalb, 1, nan_value, nan_value);
3292 TEST_ff_f (scalb, nan_value, 0, nan_value);
3293 TEST_ff_f (scalb, 0, nan_value, nan_value);
3294 TEST_ff_f (scalb, nan_value, plus_infty, nan_value);
3295 TEST_ff_f (scalb, plus_infty, nan_value, nan_value);
3296 TEST_ff_f (scalb, nan_value, nan_value, nan_value);
3298 TEST_ff_f (scalb, 0.8L, 4, 12.8L);
3299 TEST_ff_f (scalb, -0.854375L, 5, -27.34L);
3311 TEST_fi_f (scalbn, 0, 0, 0);
3312 TEST_fi_f (scalbn, minus_zero, 0, minus_zero);
3314 TEST_fi_f (scalbn, plus_infty, 1, plus_infty);
3315 TEST_fi_f (scalbn, minus_infty, 1, minus_infty);
3316 TEST_fi_f (scalbn, nan_value, 1, nan_value);
3318 TEST_fi_f (scalbn, 0.8L, 4, 12.8L);
3319 TEST_fi_f (scalbn, -0.854375L, 5, -27.34L);
3321 TEST_fi_f (scalbn, 1, 0L, 1);
3332 TEST_fl_f (scalbln, 0, 0, 0);
3333 TEST_fl_f (scalbln, minus_zero, 0, minus_zero);
3335 TEST_fl_f (scalbln, plus_infty, 1, plus_infty);
3336 TEST_fl_f (scalbln, minus_infty, 1, minus_infty);
3337 TEST_fl_f (scalbln, nan_value, 1, nan_value);
3339 TEST_fl_f (scalbln, 0.8L, 4, 12.8L);
3340 TEST_fl_f (scalbln, -0.854375L, 5, -27.34L);
3342 TEST_fl_f (scalbln, 1, 0L, 1);
3353 TEST_f_b (signbit, 0, 0);
3354 TEST_f_b (signbit, minus_zero, 1);
3355 TEST_f_b (signbit, plus_infty, 0);
3356 TEST_f_b (signbit, minus_infty, 1);
3358 /* signbit (x) != 0 for x < 0. */
3359 TEST_f_b (signbit, -1, 1);
3360 /* signbit (x) == 0 for x >= 0. */
3361 TEST_f_b (signbit, 1, 0);
3371 TEST_f_f (sin, 0, 0);
3372 TEST_f_f (sin, minus_zero, minus_zero);
3373 TEST_f_f (sin, plus_infty, nan_value, INVALID_EXCEPTION);
3374 TEST_f_f (sin, minus_infty, nan_value, INVALID_EXCEPTION);
3376 TEST_f_f (sin, M_PI_6l, 0.5);
3377 TEST_f_f (sin, -M_PI_6l, -0.5);
3378 TEST_f_f (sin, M_PI_2l, 1);
3379 TEST_f_f (sin, -M_PI_2l, -1);
3380 TEST_f_f (sin, 0.7, 0.64421768723769105367L);
3389 FLOAT sin_res, cos_res;
3393 /* sincos is treated differently because it returns void. */
3394 TEST_extra (sincos, 0, 0, 1);
3396 TEST_extra (sincos, minus_zero, minus_zero, 1);
3397 TEST_extra (sincos, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
3398 TEST_extra (sincos, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
3400 TEST_extra (sincos, M_PI_2l, 1, 0);
3401 TEST_extra (sincos, M_PI_6l, 0.5, 0.866025403784438646764L);
3402 TEST_extra (sincos, M_PI_6l*2.0, 0.866025403784438646764L, 0.5);
3403 TEST_extra (sincos, 0.7, 0.64421768723769105367L, 0.76484218728448842626L);
3412 TEST_f_f (sinh, 0, 0);
3413 TEST_f_f (sinh, minus_zero, minus_zero);
3415 TEST_f_f (sinh, plus_infty, plus_infty);
3416 TEST_f_f (sinh, minus_infty, minus_infty);
3418 TEST_f_f (sinh, 0.7, 0.75858370183953350346L);
3428 TEST_f_f (sqrt, 0, 0);
3429 TEST_f_f (sqrt, nan_value, nan_value);
3430 TEST_f_f (sqrt, plus_infty, plus_infty);
3432 TEST_f_f (sqrt, minus_zero, minus_zero);
3434 /* sqrt (x) == NaN plus invalid exception for x < 0. */
3435 TEST_f_f (sqrt, -1, nan_value, INVALID_EXCEPTION);
3436 TEST_f_f (sqrt, minus_infty, nan_value, INVALID_EXCEPTION);
3438 TEST_f_f (sqrt, 2209, 47);
3439 TEST_f_f (sqrt, 4, 2);
3440 TEST_f_f (sqrt, 2, M_SQRT2l);
3441 TEST_f_f (sqrt, 0.25, 0.5);
3442 TEST_f_f (sqrt, 6642.25, 81.5);
3443 TEST_f_f (sqrt, 15239.9025, 123.45);
3444 TEST_f_f (sqrt, 0.7, 0.83666002653407554798L);
3454 TEST_f_f (tan, 0, 0);
3455 TEST_f_f (tan, minus_zero, minus_zero);
3456 TEST_f_f (tan, plus_infty, nan_value, INVALID_EXCEPTION);
3457 TEST_f_f (tan, minus_infty, nan_value, INVALID_EXCEPTION);
3459 TEST_f_f (tan, M_PI_4l, 1);
3460 TEST_f_f (tan, 0.7, 0.84228838046307944813L);
3470 TEST_f_f (tanh, 0, 0);
3471 TEST_f_f (tanh, minus_zero, minus_zero);
3473 TEST_f_f (tanh, plus_infty, 1);
3474 TEST_f_f (tanh, minus_infty, -1);
3476 TEST_f_f (tanh, 0.7, 0.60436777711716349631L);
3486 if (errno == ENOSYS)
3487 /* Function not implemented. */
3489 feclearexcept (FE_ALL_EXCEPT);
3493 TEST_f_f (tgamma, plus_infty, plus_infty);
3494 TEST_f_f (tgamma, 0, nan_value, INVALID_EXCEPTION);
3495 TEST_f_f (tgamma, minus_zero, nan_value, INVALID_EXCEPTION);
3496 /* tgamma (x) == NaN plus invalid exception for integer x <= 0. */
3497 TEST_f_f (tgamma, -2, nan_value, INVALID_EXCEPTION);
3498 TEST_f_f (tgamma, minus_infty, nan_value, INVALID_EXCEPTION);
3500 TEST_f_f (tgamma, 0.5, M_SQRT_PIl);
3501 TEST_f_f (tgamma, -0.5, -M_2_SQRT_PIl);
3503 TEST_f_f (tgamma, 1, 1);
3504 TEST_f_f (tgamma, 4, 6);
3506 TEST_f_f (tgamma, 0.7, 1.29805533264755778568L);
3507 TEST_f_f (tgamma, 1.2, 0.91816874239976061064L);
3517 TEST_f_f (trunc, plus_infty, plus_infty);
3518 TEST_f_f (trunc, minus_infty, minus_infty);
3519 TEST_f_f (trunc, nan_value, nan_value);
3521 TEST_f_f (trunc, 0, 0);
3522 TEST_f_f (trunc, minus_zero, minus_zero);
3523 TEST_f_f (trunc, 0.625, 0);
3524 TEST_f_f (trunc, -0.625, minus_zero);
3525 TEST_f_f (trunc, 1, 1);
3526 TEST_f_f (trunc, -1, -1);
3527 TEST_f_f (trunc, 1.625, 1);
3528 TEST_f_f (trunc, -1.625, -1);
3530 TEST_f_f (trunc, 1048580.625L, 1048580L);
3531 TEST_f_f (trunc, -1048580.625L, -1048580L);
3533 TEST_f_f (trunc, 8388610.125L, 8388610.0L);
3534 TEST_f_f (trunc, -8388610.125L, -8388610.0L);
3536 TEST_f_f (trunc, 4294967296.625L, 4294967296.0L);
3537 TEST_f_f (trunc, -4294967296.625L, -4294967296.0L);
3548 if (errno == ENOSYS)
3549 /* Function not implemented. */
3552 /* y0 is the Bessel function of the second kind of order 0 */
3555 TEST_f_f (y0, -1.0, minus_infty);
3556 TEST_f_f (y0, 0.0, minus_infty);
3557 TEST_f_f (y0, nan_value, nan_value);
3558 TEST_f_f (y0, plus_infty, 0);
3560 TEST_f_f (y0, 0.1, -1.5342386513503668441);
3561 TEST_f_f (y0, 0.7, -0.19066492933739506743);
3562 TEST_f_f (y0, 1.0, 0.088256964215676957983);
3563 TEST_f_f (y0, 1.5, 0.38244892379775884396);
3564 TEST_f_f (y0, 2.0, 0.51037567264974511960);
3565 TEST_f_f (y0, 8.0, 0.22352148938756622053);
3566 TEST_f_f (y0, 10.0, 0.055671167283599391424);
3577 if (errno == ENOSYS)
3578 /* Function not implemented. */
3581 /* y1 is the Bessel function of the second kind of order 1 */
3584 TEST_f_f (y1, -1.0, minus_infty);
3585 TEST_f_f (y1, 0.0, minus_infty);
3586 TEST_f_f (y1, plus_infty, 0);
3587 TEST_f_f (y1, nan_value, nan_value);
3589 TEST_f_f (y1, 0.1, -6.4589510947020269877);
3590 TEST_f_f (y1, 0.7, -1.1032498719076333697);
3591 TEST_f_f (y1, 1.0, -0.78121282130028871655);
3592 TEST_f_f (y1, 1.5, -0.41230862697391129595);
3593 TEST_f_f (y1, 2.0, -0.10703243154093754689);
3594 TEST_f_f (y1, 8.0, -0.15806046173124749426);
3595 TEST_f_f (y1, 10.0, 0.24901542420695388392);
3605 if (errno == ENOSYS)
3606 /* Function not implemented. */
3609 /* yn is the Bessel function of the second kind of order n */
3612 /* yn (0, x) == y0 (x) */
3613 TEST_ff_f (yn, 0, -1.0, minus_infty);
3614 TEST_ff_f (yn, 0, 0.0, minus_infty);
3615 TEST_ff_f (yn, 0, nan_value, nan_value);
3616 TEST_ff_f (yn, 0, plus_infty, 0);
3618 TEST_ff_f (yn, 0, 0.1, -1.5342386513503668441);
3619 TEST_ff_f (yn, 0, 0.7, -0.19066492933739506743);
3620 TEST_ff_f (yn, 0, 1.0, 0.088256964215676957983);
3621 TEST_ff_f (yn, 0, 1.5, 0.38244892379775884396);
3622 TEST_ff_f (yn, 0, 2.0, 0.51037567264974511960);
3623 TEST_ff_f (yn, 0, 8.0, 0.22352148938756622053);
3624 TEST_ff_f (yn, 0, 10.0, 0.055671167283599391424);
3626 /* yn (1, x) == y1 (x) */
3627 TEST_ff_f (yn, 1, -1.0, minus_infty);
3628 TEST_ff_f (yn, 1, 0.0, minus_infty);
3629 TEST_ff_f (yn, 1, plus_infty, 0);
3630 TEST_ff_f (yn, 1, nan_value, nan_value);
3632 TEST_ff_f (yn, 1, 0.1, -6.4589510947020269877);
3633 TEST_ff_f (yn, 1, 0.7, -1.1032498719076333697);
3634 TEST_ff_f (yn, 1, 1.0, -0.78121282130028871655);
3635 TEST_ff_f (yn, 1, 1.5, -0.41230862697391129595);
3636 TEST_ff_f (yn, 1, 2.0, -0.10703243154093754689);
3637 TEST_ff_f (yn, 1, 8.0, -0.15806046173124749426);
3638 TEST_ff_f (yn, 1, 10.0, 0.24901542420695388392);
3641 TEST_ff_f (yn, 3, plus_infty, 0);
3642 TEST_ff_f (yn, 3, nan_value, nan_value);
3644 TEST_ff_f (yn, 3, 0.1, -5099.3323786129048894);
3645 TEST_ff_f (yn, 3, 0.7, -15.819479052819633505);
3646 TEST_ff_f (yn, 3, 1.0, -5.8215176059647288478);
3647 TEST_ff_f (yn, 3, 2.0, -1.1277837768404277861);
3648 TEST_ff_f (yn, 3, 10.0, -0.25136265718383732978);
3651 TEST_ff_f (yn, 10, plus_infty, 0);
3652 TEST_ff_f (yn, 10, nan_value, nan_value);
3654 TEST_ff_f (yn, 10, 0.1, -0.11831335132045197885e19);
3655 TEST_ff_f (yn, 10, 0.7, -0.42447194260703866924e10);
3656 TEST_ff_f (yn, 10, 1.0, -0.12161801427868918929e9);
3657 TEST_ff_f (yn, 10, 2.0, -129184.54220803928264);
3658 TEST_ff_f (yn, 10, 10.0, -0.35981415218340272205);
3669 fpstack_test ("start *init*");
3671 nan_value = plus_zero / plus_zero; /* Suppress GCC warning */
3673 minus_zero = FUNC(copysign) (0.0, -1.0);
3674 plus_infty = CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF,
3675 HUGE_VALL, HUGE_VAL, HUGE_VALF);
3676 minus_infty = CHOOSE (-HUGE_VALL, -HUGE_VAL, -HUGE_VALF,
3677 -HUGE_VALL, -HUGE_VAL, -HUGE_VALF);
3683 (void) &minus_infty;
3685 /* Clear all exceptions. From now on we must not get random exceptions. */
3686 feclearexcept (FE_ALL_EXCEPT);
3688 /* Test to make sure we start correctly. */
3689 fpstack_test ("end *init*");
3692 /* Definitions of arguments for argp functions. */
3693 static const struct argp_option options[] =
3695 { "verbose", 'v', "NUMBER", 0, "Level of verbosity (0..3)"},
3696 { "ulps-file", 'u', NULL, 0, "Output ulps to file ULPs"},
3697 { "no-max-error", 'f', NULL, 0,
3698 "Don't output maximal errors of functions"},
3699 { "no-points", 'p', NULL, 0,
3700 "Don't output results of functions invocations"},
3701 { NULL, 0, NULL, 0, NULL }
3704 /* Short description of program. */
3705 static const char doc[] = "Math test suite: " TEST_MSG ;
3707 /* Prototype for option handler. */
3708 static error_t parse_opt (int key, char *arg, struct argp_state *state);
3710 /* Data structure to communicate with argp functions. */
3711 static struct argp argp =
3713 options, parse_opt, NULL, doc,
3717 /* Handle program arguments. */
3719 parse_opt (int key, char *arg, struct argp_state *state)
3727 output_max_error = 0;
3734 verbose = (unsigned int) strtoul (optarg, NULL, 0);
3739 return ARGP_ERR_UNKNOWN;
3745 /* function to check our ulp calculation. */
3752 /* This gives one ulp. */
3753 u = FUNC(nextafter) (10, 20);
3754 check_equal (10.0, u, 1, &diff, &ulp);
3755 printf ("One ulp: % .4" PRINTF_NEXPR "\n", ulp);
3757 /* This gives one more ulp. */
3758 u = FUNC(nextafter) (u, 20);
3759 check_equal (10.0, u, 2, &diff, &ulp);
3760 printf ("two ulp: % .4" PRINTF_NEXPR "\n", ulp);
3762 /* And now calculate 100 ulp. */
3763 for (i = 2; i < 100; i++)
3764 u = FUNC(nextafter) (u, 20);
3765 check_equal (10.0, u, 100, &diff, &ulp);
3766 printf ("100 ulp: % .4" PRINTF_NEXPR "\n", ulp);
3771 main (int argc, char **argv)
3778 output_max_error = 1;
3781 /* Parse and process arguments. */
3782 argp_parse (&argp, argc, argv, 0, &remaining, NULL);
3784 if (remaining != argc)
3786 fprintf (stderr, "wrong number of arguments");
3787 argp_help (&argp, stdout, ARGP_HELP_SEE, program_invocation_short_name);
3788 exit (EXIT_FAILURE);
3793 ulps_file = fopen ("ULPs", "a");
3794 if (ulps_file == NULL)
3796 perror ("can't open file `ULPs' for writing: ");
3809 /* Keep the tests a wee bit ordered (according to ISO 9X). */
3810 /* Classification macros: */
3816 /* Trigonometric functions: */
3826 /* Hyperbolic functions: */
3834 /* Exponential and logarithmic functions: */
3852 /* Power and absolute value functions: */
3859 /* Error and gamma functions: */
3866 /* Nearest integer functions: */
3878 /* Remainder functions: */
3883 /* Manipulation functions: */
3888 /* maximum, minimum and positive difference functions */
3893 /* Multiply and add: */
3896 /* Complex functions: */
3918 /* Bessel functions: */
3929 printf ("\nTest suite completed:\n");
3930 printf (" %d test cases plus %d tests for exception flags executed.\n",
3931 noTests, noExcTests);
3933 printf (" %d expected failures occured.\n", noXFails);
3935 printf (" %d unexpected passes occured.\n", noXPasses);
3938 printf (" %d errors occured.\n", noErrors);
3941 printf (" All tests passed successfully.\n");