aboutsummaryrefslogtreecommitdiff
path: root/math
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>1997-03-16 20:28:07 +0000
committerUlrich Drepper <drepper@redhat.com>1997-03-16 20:28:07 +0000
commit0d8733c4fc3695c7987548d10c344445f3eb552d (patch)
tree4a2fe41c465963794a6f9d4d1717494972336e14 /math
parent1ed4222e5724ce863a95c15eb530dfab7ef0451c (diff)
downloadglibc-0d8733c4fc3695c7987548d10c344445f3eb552d.tar
glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.tar.gz
glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.tar.bz2
glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.zip
1997-03-16 18:43 Ulrich Drepper <drepper@cygnus.com> * manual/filesys.texi: Add documentation for scandir and alphasort. * math/math.c (fpclassify): Correct stupid typos. * math/libm-test.c: New file. libm test suite by Andreas Jaeger. * nss/nss_files/files-hosts.c: Add gethostbyname2 imlementation. * posix/Makefile (routines): Add bsd-getpgrp. * posix/bsd-getpgrp.c: New file. * posix/unistd.h [__FAVOR_BSD]: Define macro getpgrp which maps calls to __bsd_getpgrp. * sysdeps/generic/getpgrp.c: De-ANSI-declfy. * sysdeps/i386/huge_val.h: New file. ix87 specific infinity values. * sysdeps/m68k/huge_val.h: New file. m68k specific infinity values. * sysdeps/generic/huge_val.h: Remove definition of long double definition. Make it the same as the double definition. * sysdeps/libm-i387/e_acos.S: Fix bug in FPU stack handling. * sysdeps/libm-i387/e_acosf.S: Likewise. * sysdeps/libm-i387/e_acosl.S: Likewise. * sysdeps/libm-i387/e_asin.S: Likewise. * sysdeps/libm-i387/e_asinf.S: Likewise. * sysdeps/libm-i387/e_asinl.S: Likewise. * sysdeps/libm-i387/e_exp.S: Likewise. * sysdeps/libm-i387/e_expf.S: Likewise. * sysdeps/libm-i387/e_expl.S: Likewise. * sysdeps/libm-i387/e_scalbn.S: Likewise. * sysdeps/libm-i387/e_scalbnf.S: Likewise. * sysdeps/libm-i387/e_scalbnl.S: Likewise. * sysdeps/libm-i387/e_log.S: Optimize branch code. * sysdeps/libm-i387/e_logf.S: Likewise. * sysdeps/libm-i387/e_logl.S: Likewise. * sysdeps/libm-i387/e_log10.S: Likewise. * sysdeps/libm-i387/e_log10f.S: Likewise. * sysdeps/libm-i387/e_log10l.S: Likewise. * sysdeps/libm-i387/e_pow.S: Major rewrite to handle special cases. * sysdeps/libm-i387/e_powf.S: Likewise. * sysdeps/libm-i387/e_powl.S: Likewise. * sysdeps/libm-i387/e_expm1.S: Change return value for -inf argument to -1.0. * sysdeps/libm-i387/e_expm1f.S: Likewise. * sysdeps/libm-i387/e_expm1l.S: Likewise. * sysdeps/libm-i387/e_isinfl.c: Return -1 for -inf. * sysdeps/libm-i387/e_logbl.S: Correct return value. Discard first stack element after fxtract. * sysdeps/libm-ieee754/e_atan2l.c: New file. `long double' implementation for atan2 function. * sysdeps/libm-ieee754/k_standard.c: Return NAN for libm not in _SVID_ mode when acos, asin, atan2, log, log10 is called with argument out of range. Add new error case for pow(+0,neg). * sysdeps/libm-ieee754/s_fpclassifyf.c: Correct recognition of NaN and +-inf. * sysdeps/libm-ieee754/s_fpclassifyl.c: Mask out explicit leading digit in stupid 80 bit formats. * sysdeps/libm-ieee754/s_isinf.c: Rewrite to return -1 for -inf. * sysdeps/libm-ieee754/s_isinff.c: Likewise. * sysdeps/libm-ieee754/s_isinfl.c: Likewise. * sysdeps/libm-ieee754/s_scalbnl.c (huge, tiny): Adapt values for long double type. * sysdeps/libm-ieee754/w_atan2.c: Do not raise exception expect when in SVID mode. * sysdeps/libm-ieee754/w_atan2f.c: Likewise. * sysdeps/libm-ieee754/w_atan2l.c: Likewise. * sysdeps/libm-ieee754/w_pow.c: Distinguish error cases for x is +0 or -0. * sysdeps/posix/isfdtype.c: Add cast to prevent warning. * sysdeps/stub/fcntlbits.h: Update copyright. * sysdeps/unix/bsd/fcntlbits.h: Likewise. * sysdeps/unix/bsd/bsd4.4/fcntlbits.h: Likewise. * sysdeps/unix/bsd/sun/sunos4/fcntlbits.h: Likewise. * sysdeps/unix/bsd/ultrix4/fcntlbits.h: Likewise. * sysdeps/unix/common/fcntlbits.h: Likewise. * sysdeps/unix/sysv/fcntlbits.h: Likewise. Define O_FSYNC as alias of O_SYNC. Add BSD compatibility macros FAPPEND, FFSYNC, FNONBLOCK, and FNDELAY. * sysdeps/unix/sysv/irix4/fcntlbits.h: Likewise. * sysdeps/unix/readdir_r.c: Don't copy whole `struct dirent' record, only reclen bytes. * sysdeps/unix/sysv/linux/fcntlbits.h [__USE_GNU]: Add O_READ, O_WRITE and O_NORW. * sysdeps/unix/sysv/linux/alpha/fcntlbits.h: Likewise. * sysdeps/unix/sysv/linux/init-first.h: Add copyright. * sysdeps/unix/sysv/linux/fxstat.c: New file. Rewrite kernel-level struct stat to user-level form. * sysdeps/unix/sysv/linux/lxstat: New file. * sysdeps/unix/sysv/linux/xstat: New file. * sysdeps/unix/sysv/linux/kernel_stat.h: Define struct stat used in kernel. * sysdeps/unix/sysv/linux/statbuf.h (struct stat): Change definition to use prescribed types for elements. (_STAT_VER): Change to value 3. * sysdeps/unix/sysv/linux/alph/statbuf.h: Likewise. * sysdeps/unix/sysv/linux/Dist: Add kernel_stat.h. * sysdeps/unix/sysv/linux/alpha/Dist: Likewise. * time/Makefile: Correct dependencies for test-tz. 1997-03-16 14:59 Philip Blundell <phil@london.uk.eu.org> * resolv/netdb.h: Add prototypes for gai_strerror and getnameinfo (needed for IPv6 basic sockets API). 1997-03-16 15:02 a sun <asun@zoology.washington.edu> * sysdeps/unix/sysv/linux/net/if_ppp.h: Don't use incompatible kernel header. * sysdeps/unix/sysv/linux/net/ppp_defs.h: Likewise. 1997-03-14 17:15 Ulrich Drepper <drepper@cygnus.com> * db/hash/hash_bigkey.c (__big_delete): Don't call __free_ovflpage without testing for last_bfp to be NULL. Reported by fabsoft@fabserver1.zarm.uni-bremen.de. 1997-03-13 11:42 Jim Meyering <meyering@asic.sc.ti.com> * time/mktime.c (TIME_T_MIN): Work around a bug in Cray C 5.0.3.0. 1997-03-14 04:00 Kurt Garloff <garloff@kg1.ping.de> * sysdeps/unix/sysv/linux/fcntlbits.h (O_FSYNC): Make alias for O_SYNC. (FASYNC): Move to __USE_BSD section. Create new macro O_ASYNC. 1997-03-14 02:50 Ulrich Drepper <drepper@cygnus.com> * nis/nss_nis/nis-hosts.c (_nss_nis_gethostbyname2_r): New functions. Compare result for correct address type. (_nss_nis_gethostbyname_r): Use _nss_nis_gethostbyname2_r. Reported by Mirko Streckenbach <mirko@marian.hil.de>. 1997-02-17 01:40 Zlatko Calusic <zcalusic@srce.hr> * time/strptime.c (recursive): Return rp to caller. (strptime_internal): First check for long names, then abbreviated (month & weekday). 1997-03-10 19:44 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * Makeconfig: Remove useless definitions of ASFLAGS-%. * config.make.in (ASFLAGS-.so): Remove. * configure.in: Don't substitute ASFLAGS_SO. * sysdeps/sparc/configure.in: Remove file. * sysdeps/sparc/Makefile (ASFLAGS-.so): Define. 1997-03-11 17:00 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * time/strptime.c (strptime_internal) [case 'Y']: Always subtract 1900 from year, regardless of century. 1997-03-12 05:43 Ulrich Drepper <drepper@cygnus.com> * stdlib/strtod.c (_tens_in_limb) [BITS_PER_MP_LIMB > 32]: Make all numbers unsigned to make buggy gccs happy. Patch by Bryan W. Headley <bheadley@interaccess.com>. * sysdeps/unix/sysv/linux/netinet/ip.h: Add backward-compatibility definitions. Patch by a sun <asun@zoology.washington.edu>. Pretty print header. * Makerules (build-shlib): Also create symlink if library is versioned. based on a patch by H.J. Lu <hjl@gnu.ai.mit.edu>. Remove special rule to libc.so symlink. 1997-03-11 20:16 Andreas Jaeger <aj@arthur.pfalz.de> * manual/math.texi (Domain and Range Errors): Change descriptions according to recent changes for ISO C 9X. 1997-03-11 22:39 Ulrich Drepper <drepper@cygnus.com> * sysdeps/libm-ieee754/k_standard.c (__kernel_standard): Correct return values for acos, asin, and atan2. Reported by Andreas Jaeger <aj@arthur.pfalz.de>. 1997-03-10 18:16 Thorsten Kukuk <kukuk@vt.uni-paderborn.de> * ypclnt.c (__yp_bind): Fix possible buffer overflow. 1997-03-10 18:06 Bernd Schmidt <crux@Pool.Informatik.RWTH-Aachen.DE> * dirent/alphasort.c (alphasort): Interpret arguments as pointers to pointers to directory entries so that alphasort really can be used as argument for scandir. 1997-03-09 23:33 Andreas Jaeger <aj@arthur.pfalz.de> * string/strdup.c: Declare memcpy if !(_LIBC || STDC_HEADERS) instead of strcpy. 1997-03-10 03:34 Ulrich Drepper <drepper@cygnus.com> * catgets/catgets.c (catopen): Always add NLSPATH to search path for catalogs, not only if the envvar NLSPATH is not available. Reported by Andries.Brouwer@cwi.nl. 1997-03-10 02:46 Ulrich Drepper <drepper@cygnus.com> * Makeconfig (localtime-file): Don't define using installation directory. (inst_localtime-file): New variable. * time/Makefile (installed-localtime-file): Use inst_localtime-file. Reported by Edward Seidl <seidl@janed.com>. 1997-03-10 02:31 H.J. Lu <hjl@gnu.ai.mit.edu> * time/Makefile: Add source files to dependencies for test data. 1997-03-09 22:53 Thorsten Kukuk <kukuk@weber.uni-paderborn.de> * nis/nss_nis/nis-ethers.c: Don't ignore return value of yp_all. * nis/nss_nis/nis-proto.c: Likewise. * nis/nss_nis/nis-rpc.c: Likewise. * nis/nss_nis/nis-service.c: Likewise. 1997-03-08 14:37 Miguel de Icaza <miguel@nuclecu.unam.mx> * sysdeps/sparc/dl-machine.h (elf_machine_rela): Upgrade to versioning; Added missing R_SPARC_WDISP30 handling. (RTLD_START): Implement it. * sysdeps/unix/sysv/linux/sparc/brk.c: Fix. * sysdeps/unix/sysv/linux/sparc/start.c: Startup code for Linux/SPARC. 1997-03-02 18:06 Miguel de Icaza <miguel@nuclecu.unam.mx> * sysdeps/sparc/dl-machine.h (RTLD_START): Make arg as expected by the dynamic linker instead of having a new conditional define. Thanks to Richard Henderson for pointing this out. * elf/rtld.c: Remove usage of ELF_ADJUST_ARG. 1997-03-20 20:44 Thomas Bushnell, n/BSG <thomas@gnu.ai.mit.edu> * sysdeps/mach/hurd/euidaccess.c: Define as __euidaccess and make euidaccess weak alias. 1997-03-07 10:30 Thomas Bushnell, n/BSG <thomas@gnu.ai.mit.edu> * stdio-common/printf_fphex.c (MIN): New macro. * sysdeps/generic/netinet/in.h: Include <sys/types.h>. * sysdeps/generic/sys/mman.h (msync): Mention third arg. * sysdeps/generic/netinet/in.h: Add definitions for IPv6 basic API. (See change by Philip Blundell on Feb 16, 1997.) 1997-03-05 10:40 Thomas Bushnell, n/BSG <thomas@gnu.ai.mit.edu> * hurd/hurd.h (vpprintf): Include <stdarg.h>. New declaration. * hurd/set-host.c (_hurd_set_host_config): Cast second arg to __file_name_split. * mach/mach_error.c (mach_error_string_int): Give full prototype. * mach/errstring.c (mach_error_string_int): Likewise. * mach/error_compat.c (__mach_error_map_compat): Likewise. * hurd/vpprintf.c (pwrite, vpprintf): Likewise. * stdio/vasprintf.c (vasprintf): Likewise. * mach/mach/mach_traps.h: Include <mach/kern_return.h>. * mach/spin-solid.c: Include <mach/mach_traps.h>. * mach/spin-solid.c (__spin_lock_solid): Provide arg to swtch_pri. * mach/mach_init.c: Include <mach/mig_support.h>. * mach/mach_error.h (mach_error_string, mach_error, mach_error_type): Always provide prototypes. * mach/mach/error.h (mach_error_fn_t): Comment out declaration; it appears to be entirely unused dead code. * stdio/stdio.h (freopen): Fix spelling error. 1997-03-02 13:38 Miles Bader <miles@gnu.ai.mit.edu> * string/argz.h (__need_error_t): New macro, before including <errno.h> [!__const] (__const): New macro. [!__error_t_defined] (error_t): New typedef. * sysdeps/generic/socketbits.h: Add PF_FILE as synonym for PF_LOCAL * sysdeps/unix/sysv/linux/socketbits.h: Likewise.
Diffstat (limited to 'math')
-rw-r--r--math/libm-test.c1520
-rw-r--r--math/math.h4
2 files changed, 1522 insertions, 2 deletions
diff --git a/math/libm-test.c b/math/libm-test.c
new file mode 100644
index 0000000000..98e3cbad5d
--- /dev/null
+++ b/math/libm-test.c
@@ -0,0 +1,1520 @@
+/* Copyright (C) 1997 Free Software Foundation, Inc.
+ This file is part of the GNU C Library.
+ Contributed by Andreas Jaeger <aj@arthur.pfalz.de>, 1997.
+
+ The GNU C Library is free software; you can redistribute it and/or
+ modify it under the terms of the GNU Library General Public License as
+ published by the Free Software Foundation; either version 2 of the
+ License, or (at your option) any later version.
+
+ The GNU C Library is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ Library General Public License for more details.
+
+ You should have received a copy of the GNU Library General Public
+ License along with the GNU C Library; see the file COPYING.LIB. If not,
+ write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+ Boston, MA 02111-1307, USA. */
+
+
+/*
+ Part of testsuite for libm.
+
+ This file has to be included by a master file that defines:
+
+ Makros:
+ FUNC(function): converts general function name (like cos) to
+ name with correct suffix (e.g. cosl or cosf)
+ MATHCONST(x): like FUNC but for constants (e.g convert 0.0 to 0.0L)
+ MATHTYPE: floating point type to test
+ TEST_MSG: informal message to be displayed
+ CHOOSE(Clongdouble,Cdouble,Cfloat):
+ chooses one of the parameters as epsilon for testing
+ equality
+ PRINTF_EXPR Floating point conversion specification to print a variable
+ of type MATHTYPE with printf. PRINTF_EXPR just contains
+ the specifier, not the percent and width arguments,
+ e.g. "f"
+ */
+
+/* This program isn't finished yet.
+ It has tests for acos, acosh, asin, asinh, atan, atan2, atanh,
+ cbrt, ceil, cos, cosh, exp, exp2, expm1, fabs, floor, fpclassify,
+ frexp, ldexp,
+ log, log10, log1p, log2, logb,
+ pow, sin, sinh, tan, tanh, fabs, hypot.
+ Tests for the other libm-functions will come later.
+
+ The routines using random variables are still under construction. I don't
+ like it the way it's working now and will change it.
+
+ Exception handling has not been implemented so far so don't get fooled
+ that these tests pass.
+
+ Parameter handling is primitive in the moment:
+ --verbose=[0..3] for different levels of output:
+ 0: only error count
+ 1: basic report on failed tests
+ 2: full report on failed tests
+ 3: full report on failed and passed tests (default)
+ -v for full output (equals --verbose=3)
+ -s,--silent outputs only the error count (equals --verbose=0)
+ */
+
+/* Define if the following ISO C 9X functions are implemented: exp2,
+ log2. */
+#undef ISO_9X_IMPLEMENTED
+
+#define _GNU_SOURCE
+
+#include <math.h>
+#include <float.h>
+
+#include <errno.h>
+#include <stdlib.h>
+#include <stdio.h>
+#include <time.h>
+#include <getopt.h>
+
+/* TEST_EXCEPTION: tests if an exception as occured */
+/* for the moment: does nothing */
+/* Possible exceptions */
+#define NO_EXCEPTION 0x0
+#define INVALID_EXCEPTION 0x1
+#define DIVIDE_BY_ZERO_EXCEPTION 0x2
+
+#define PRINT 1
+#define NO_PRINT 0
+
+#define TEST_EXCEPTION(test) do {} while (0);
+/* As long as no exception code is available prevent warnings. */
+#define UNUSED __attribute__ ((unused))
+
+static int noErrors;
+
+static int verbose = 3;
+static MATHTYPE minus_zero, plus_zero;
+static MATHTYPE plus_infty, minus_infty, nan_value;
+
+typedef MATHTYPE (*mathfunc) (MATHTYPE);
+
+
+#define ISINF(x) \
+(sizeof (x) == sizeof (float) ? \
+ isinff (x) \
+ : sizeof (x) == sizeof (double) ? \
+ isinf (x) : isinfl (x))
+
+
+ /*
+ Test if Floating-Point stack hasn't changed
+ */
+static void
+fpstack_test (const char *test_name)
+{
+#ifdef i386
+ static int old_stack;
+ int sw;
+asm ("fnstsw":"=a" (sw));
+ sw >>= 11;
+ sw &= 7;
+ if (sw != old_stack)
+ {
+ printf ("FP-Stack wrong after test %s\n", test_name);
+ if (verbose > 2)
+ printf ("=======> stack = %d\n", sw);
+ ++noErrors;
+ old_stack = sw;
+ }
+#endif
+}
+
+
+/*
+ Call to an external function so that floating point registers
+ get moved to memory
+ */
+static void
+this_does_nothing (void)
+{
+ clock_t dummy;
+
+ dummy = clock ();
+}
+
+/*
+ Get a random value x with min_value < x < max_value
+ and min_value, max_value finite,
+ max_value and min_value shouldn't be too close together
+ */
+static MATHTYPE
+random_value (MATHTYPE min_value, MATHTYPE max_value)
+{
+ int r;
+ MATHTYPE x;
+
+ r = rand ();
+
+ x = (max_value - min_value) / RAND_MAX * (MATHTYPE) r + min_value;
+
+ if ((x <= min_value) || (x >= max_value) || !isfinite (x))
+ x = (max_value - min_value) / 2 + min_value;
+
+ return x;
+}
+
+/* Get a random value x with x > min_value. */
+static MATHTYPE
+random_greater (MATHTYPE min_value)
+{
+ return random_value (min_value, 1e6); /* CHOOSE (LDBL_MAX, DBL_MAX, FLT_MAX) */
+}
+
+/* Get a random value x with x < max_value. */
+static MATHTYPE
+random_less (MATHTYPE max_value)
+{
+ return random_value (-1e6, max_value);
+}
+
+
+/* Test if two floating point numbers are equal. */
+static int
+check_equal (MATHTYPE computed, MATHTYPE supplied, MATHTYPE eps, MATHTYPE * diff)
+{
+ /* Both plus Infinity or both minus infinity. */
+ if (ISINF (computed) && (ISINF (computed) == ISINF (supplied)))
+ return 1;
+
+ if (isnan (computed) && isnan (supplied)) /* isnan works for all types */
+ return 1;
+
+ *diff = FUNC(fabs) (computed - supplied);
+
+ if (*diff <= eps || signbit (computed) != signbit (supplied))
+ return 1;
+
+ return 0;
+}
+
+
+static void
+output_result_bool (const char *test_name, int result)
+{
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ noErrors++;
+ }
+
+ fpstack_test (test_name);
+}
+
+
+static void
+output_isvalue (const char *test_name, int result,
+ MATHTYPE value)
+{
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1)
+ printf (" Value: %.20" PRINTF_EXPR "\n", value);
+ noErrors++;
+ }
+
+ fpstack_test (test_name);
+}
+
+
+static void
+output_isvalue_ext (const char *test_name, int result,
+ MATHTYPE value, MATHTYPE parameter)
+{
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1)
+ {
+ printf (" Value: %.20" PRINTF_EXPR "\n", value);
+ printf (" Parameter: %.20" PRINTF_EXPR "\n", parameter);
+ }
+ noErrors++;
+ }
+
+ fpstack_test (test_name);
+}
+
+
+static void
+output_result (const char *test_name, int result,
+ MATHTYPE computed, MATHTYPE expected,
+ MATHTYPE difference,
+ int print_values, int print_diff)
+{
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1 && print_values)
+ {
+ printf ("Result:\n");
+ printf (" is: %.20" PRINTF_EXPR "\n", computed);
+ printf (" should be: %.20" PRINTF_EXPR "\n", expected);
+ if (print_diff)
+ printf (" difference: %.20" PRINTF_EXPR "\n", difference);
+ }
+ noErrors++;
+ }
+
+ fpstack_test (test_name);
+}
+
+
+static void
+output_result_ext (const char *test_name, int result,
+ MATHTYPE computed, MATHTYPE expected,
+ MATHTYPE difference,
+ MATHTYPE parameter,
+ int print_values, int print_diff)
+{
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1 && print_values)
+ {
+ printf ("Result:\n");
+ printf (" is: %.20" PRINTF_EXPR "\n", computed);
+ printf (" should be: %.20" PRINTF_EXPR "\n", expected);
+ if (print_diff)
+ printf (" difference: %.20" PRINTF_EXPR "\n", difference);
+ printf ("Parameter: %.20" PRINTF_EXPR "\n", parameter);
+ }
+ noErrors++;
+ }
+
+ fpstack_test (test_name);
+}
+
+
+static void
+check (const char *test_name, MATHTYPE computed, MATHTYPE expected)
+{
+ MATHTYPE diff;
+ int result;
+
+ result = check_equal (computed, expected, 0, &diff);
+ output_result (test_name, result,
+ computed, expected, diff, PRINT, PRINT);
+}
+
+
+static void
+check_ext (const char *test_name, MATHTYPE computed, MATHTYPE expected,
+ MATHTYPE parameter)
+{
+ MATHTYPE diff;
+ int result;
+
+ result = check_equal (computed, expected, 0, &diff);
+ output_result_ext (test_name, result,
+ computed, expected, diff, parameter, PRINT, PRINT);
+}
+
+
+static void
+check_eps (const char *test_name, MATHTYPE computed, MATHTYPE expected,
+ MATHTYPE epsilon)
+{
+ MATHTYPE diff;
+ int result;
+
+ result = check_equal (computed, expected, epsilon, &diff);
+ output_result (test_name, result,
+ computed, expected, diff, PRINT, PRINT);
+}
+
+
+static void
+check_bool (const char *test_name, int computed)
+{
+ output_result_bool (test_name, computed);
+}
+
+
+static void
+check_isnan (const char *test_name, MATHTYPE computed)
+{
+ output_isvalue (test_name, isnan (computed), computed);
+}
+
+
+static void
+check_isnan_exc (const char *test_name, MATHTYPE computed,
+ short exception UNUSED)
+{
+ output_isvalue (test_name, isnan (computed), computed);
+}
+
+
+static void
+check_isnan_ext (const char *test_name, MATHTYPE computed,
+ MATHTYPE parameter)
+{
+ output_isvalue_ext (test_name, isnan (computed), computed, parameter);
+}
+
+
+/* Tests if computed is +Inf */
+static void
+check_isinfp (const char *test_name, MATHTYPE computed)
+{
+ output_isvalue (test_name, (ISINF (computed) == +1), computed);
+}
+
+
+static void
+check_isinfp_ext (const char *test_name, MATHTYPE computed,
+ MATHTYPE parameter)
+{
+ output_isvalue_ext (test_name, (ISINF (computed) == +1), computed, parameter);
+}
+
+
+/* Tests if computed is +Inf */
+static void
+check_isinfp_exc (const char *test_name, MATHTYPE computed,
+ int exception UNUSED)
+{
+ output_isvalue (test_name, (ISINF (computed) == +1), computed);
+}
+
+/* Tests if computed is -Inf */
+static void
+check_isinfn (const char *test_name, MATHTYPE computed)
+{
+ output_isvalue (test_name, (ISINF (computed) == -1), computed);
+}
+
+
+static void
+check_isinfn_ext (const char *test_name, MATHTYPE computed,
+ MATHTYPE parameter)
+{
+ output_isvalue_ext (test_name, (ISINF (computed) == -1), computed, parameter);
+}
+
+
+/* Tests if computed is -Inf */
+static void
+check_isinfn_exc (const char *test_name, MATHTYPE computed,
+ int exception UNUSED)
+{
+ output_isvalue (test_name, (ISINF (computed) == -1), computed);
+}
+
+
+/****************************************************************************
+ Test for single functions of libm
+****************************************************************************/
+
+static void
+acos_test (void)
+{
+ MATHTYPE x;
+
+ check ("acos (1) == 0", FUNC(acos) (1), 0);
+
+ x = random_greater (1);
+ check_isnan_exc ("acos (x) == NaN + invalid exception for |x| > 1",
+ FUNC(acos) (x),
+ INVALID_EXCEPTION);
+}
+
+static void
+acosh_test (void)
+{
+ MATHTYPE x;
+
+ check ("acosh(1) == 0", FUNC(acosh) (1), 0);
+ check_isinfp ("acosh(+inf) == +inf", FUNC(acosh) (plus_infty));
+
+ x = random_less (1);
+ check_isnan_exc ("acosh(x) == NaN plus invalid exception if x < 1",
+ FUNC(acosh) (x), INVALID_EXCEPTION);
+}
+
+
+static void
+asin_test (void)
+{
+ MATHTYPE x;
+ check ("asin (0) == 0", FUNC(asin) (0), 0);
+
+ x = random_greater (1);
+ check_isnan_exc ("asin x == NaN + invalid exception for |x| > 1",
+ FUNC(asin) (x),
+ INVALID_EXCEPTION);
+}
+
+static void
+asinh_test (void)
+{
+
+ check ("asinh(+0) == +0", FUNC(asinh) (0), 0);
+ check ("asinh(-0) == -0", FUNC(asinh) (minus_zero), minus_zero);
+}
+
+
+static void
+atan_test (void)
+{
+ check ("atan (0) == 0", FUNC(atan) (0), 0);
+ check ("atan (-0) == -0", FUNC(atan) (minus_zero), minus_zero);
+
+ check ("atan (+inf) == pi/2", FUNC(atan) (plus_infty), M_PI_2);
+ check ("atan (-inf) == -pi/2", FUNC(atan) (minus_infty), -M_PI_2);
+
+}
+
+static void
+atan2_test (void)
+{
+ MATHTYPE x;
+
+ x = random_greater (0);
+ check ("atan2 (0,x) == 0 for x > 0",
+ FUNC(atan2) (0, x), 0);
+ x = random_greater (0);
+ check ("atan2 (-0,x) == -0 for x > 0",
+ FUNC(atan2) (minus_zero, x), minus_zero);
+
+ check ("atan2 (+0,+0) == +0", FUNC(atan2) (0, 0), 0);
+ check ("atan2 (-0,+0) == -0", FUNC(atan2) (minus_zero, 0), minus_zero);
+
+ x = -random_greater (0);
+ check ("atan2 (+0,x) == +pi for x < 0", FUNC(atan2) (0, x), M_PI);
+
+ x = -random_greater (0);
+ check ("atan2 (-0,x) == -pi for x < 0", FUNC(atan2) (minus_zero, x), -M_PI);
+
+ check ("atan2 (+0,-0) == +pi", FUNC(atan2) (0, minus_zero), M_PI);
+ check ("atan2 (-0,-0) == -pi", FUNC(atan2) (minus_zero, minus_zero), -M_PI);
+
+ x = random_greater (0);
+ check ("atan2 (y,+0) == pi/2 for y > 0", FUNC(atan2) (x, 0), M_PI_2);
+
+ x = random_greater (0);
+ check ("atan2 (y,-0) == pi/2 for y > 0", FUNC(atan2) (x, minus_zero), M_PI_2);
+
+ x = random_greater (0);
+ check ("atan2 (y,-inf) == +pi for finite y > 0",
+ FUNC(atan2) (x, minus_infty), M_PI);
+
+ x = -random_greater (0);
+ check ("atan2 (y,-inf) == -pi for finite y < 0",
+ FUNC(atan2) (x, minus_infty), -M_PI);
+
+ check ("atan2 (+inf,+inf) == +pi/4",
+ FUNC(atan2) (plus_infty, plus_infty), M_PI_4);
+
+ check ("atan2 (-inf,+inf) == -pi/4",
+ FUNC(atan2) (minus_infty, plus_infty), -M_PI_4);
+
+ check ("atan2 (+inf,-inf) == +3*pi/4",
+ FUNC(atan2) (plus_infty, minus_infty), 3 * M_PI_4);
+
+ check ("atan2 (-inf,-inf) == -3*pi/4",
+ FUNC(atan2) (minus_infty, minus_infty), -3 * M_PI_4);
+}
+
+
+static void
+atanh_test (void)
+{
+
+ check ("atanh(+0) == +0", FUNC(atanh) (0), 0);
+ check ("atanh(-0) == -0", FUNC(atanh) (minus_zero), minus_zero);
+ check_isinfp_exc ("atanh(+1) == +inf plus divide-by-zero exception",
+ FUNC(atanh) (1), DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("atanh(-1) == -inf plus divide-by-zero exception",
+ FUNC(atanh) (-1), DIVIDE_BY_ZERO_EXCEPTION);
+}
+
+
+static void
+cbrt_test (void)
+{
+ check ("cbrt (+0) == +0", FUNC(cbrt) (0.0), 0.0);
+ check ("cbrt (-0) == -0", FUNC(cbrt) (minus_zero), minus_zero);
+
+ check ("cbrt (8) == 2", FUNC(cbrt) (8), 2);
+ check ("cbrt (-27) == -3", FUNC(cbrt) (-27.0), -3.0);
+}
+
+
+static void
+ceil_test (void)
+{
+ check ("ceil (+0) == +0", FUNC(ceil) (0.0), 0.0);
+ check ("ceil (-0) == -0", FUNC(ceil) (minus_zero), minus_zero);
+ check_isinfp ("ceil (+inf) == +inf", FUNC(ceil) (plus_infty));
+ check_isinfn ("ceil (-inf) == -inf", FUNC(ceil) (minus_infty));
+
+ check ("ceil (pi) == 4", FUNC(ceil) (M_PI), 4.0);
+ check ("ceil (-pi) == -3", FUNC(ceil) (-M_PI), 3.0);
+}
+
+
+static void
+cos_test (void)
+{
+
+ check ("cos (+0) == 1", FUNC(cos) (0), 1);
+ check ("cos (-0) == 1", FUNC(cos) (minus_zero), 1);
+ check_isnan_exc ("cos (+inf) == NaN plus invalid exception",
+ FUNC(cos) (plus_infty),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("cos (-inf) == NaN plus invalid exception",
+ FUNC(cos) (minus_infty),
+ INVALID_EXCEPTION);
+
+ check_eps ("cos (pi/3) == 0.5", FUNC(cos) (M_PI / 3.0),
+ 0.5, CHOOSE (0, 1e-15L, 1e-7L));
+ check_eps ("cos (pi/2) == 0.5", FUNC(cos) (M_PI_2),
+ 0, CHOOSE (1e-19L, 1e-16L, 1e-7L));
+
+}
+
+static void
+cosh_test (void)
+{
+ check ("cosh (+0) == 1", FUNC(cosh) (0), 1);
+ check ("cosh (-0) == 1", FUNC(cosh) (minus_zero), 1);
+
+ check_isinfp ("cosh (+inf) == +inf", FUNC(cosh) (plus_infty));
+ check_isinfp ("cosh (-inf) == +inf", FUNC(cosh) (minus_infty));
+}
+
+
+static void
+exp_test (void)
+{
+ check ("exp (+0) == 1", FUNC(exp) (0), 1);
+ check ("exp (-0) == 1", FUNC(exp) (minus_zero), 1);
+
+ check_isinfp ("exp (+inf) == +inf", FUNC(exp) (plus_infty));
+ check ("exp (-inf) == 0", FUNC(exp) (minus_infty), 0);
+
+ check ("exp (1) == e", FUNC(exp) (1), M_E);
+}
+
+
+#ifdef ISO_9X_IMPLEMENTED
+static void
+exp2_test (void)
+{
+ check ("exp2 (+0) == 1", FUNC(exp2) (0), 1);
+ check ("exp2 (-0) == 1", FUNC(exp2) (minus_zero), 1);
+
+ check_isinfp ("exp2 (+inf) == +inf", FUNC(exp2) (plus_infty));
+ check ("exp2 (-inf) == 0", FUNC(exp2) (minus_infty), 0);
+ check ("exp2 (10) == 1024", FUNC(exp2) (10), 1024);
+}
+#endif
+
+
+static void
+expm1_test (void)
+{
+ check ("expm1 (+0) == 0", FUNC(expm1) (0), 0);
+ check ("expm1 (-0) == -0", FUNC(expm1) (minus_zero), minus_zero);
+
+ check_isinfp ("expm1 (+inf) == +inf", FUNC(expm1) (plus_infty));
+ check ("expm1 (-inf) == -1", FUNC(expm1) (minus_infty), -1);
+
+ check ("expm1 (1) == e-1", FUNC(expm1) (1), M_E - 1.0);
+}
+
+
+
+
+static void
+check_frexp (const char *test_name, MATHTYPE computed, MATHTYPE expected,
+ int comp_int, int exp_int)
+{
+ MATHTYPE diff;
+ int result;
+
+ result = (check_equal (computed, expected, 0, &diff)
+ && (comp_int == exp_int));
+
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1)
+ {
+ printf ("Result:\n");
+ printf (" is: %.20" PRINTF_EXPR " *2^%d\n", computed, comp_int);
+ printf (" should be: %.20" PRINTF_EXPR " *2^%d\n", expected, exp_int);
+ printf (" difference: %.20" PRINTF_EXPR "\n", diff);
+ }
+ noErrors++;
+ }
+ fpstack_test (test_name);
+ output_result (test_name, result,
+ computed, expected, diff, PRINT, PRINT);
+}
+
+
+static void
+frexp_test (void)
+{
+ int x_int;
+ MATHTYPE result;
+
+ result = FUNC(frexp) (plus_infty, &x_int);
+ check_isinfp ("frexp (+inf, expr) == +inf", result);
+
+ result = FUNC(frexp) (minus_infty, &x_int);
+ check_isinfn ("frexp (-inf, expr) == -inf", result);
+
+ result = FUNC(frexp) (nan_value, &x_int);
+ check_isnan ("frexp (Nan, expr) == NaN", result);
+
+ result = FUNC(frexp) (0, &x_int);
+ check_frexp ("frexp: +0 == 0 * 2^0", result, 0, x_int, 0);
+
+ result = FUNC(frexp) (minus_zero, &x_int);
+ check_frexp ("frexp: -0 == -0 * 2^0", result, minus_zero, x_int, 0);
+
+ result = FUNC(frexp) (12.8L, &x_int);
+ check_frexp ("frexp: 12.8 == 0.8 * 2^4", result, 0.8L, x_int, 4);
+
+ result = FUNC(frexp) (-27.34L, &x_int);
+ check_frexp ("frexp: -27.34 == -0.854375 * 2^5", result, -0.854375L, x_int, 5);
+
+}
+
+
+#if __GLIBC__ < 2 || (__GLIBC__ == 2 && __GLIBC_MINOR__ < 1)
+/* All floating-point numbers can be put in one of these categories. */
+enum
+{
+ FP_NAN,
+#define FP_NAN FP_NAN
+ FP_INFINITE,
+#define FP_INFINITE FP_INFINITE
+ FP_ZERO,
+#define FP_ZERO FP_ZERO
+ FP_SUBNORMAL,
+#define FP_SUBNORMAL FP_SUBNORMAL
+ FP_NORMAL
+#define FP_NORMAL FP_NORMAL
+};
+#endif
+
+
+static void
+fpclassify_test (void)
+{
+ MATHTYPE x;
+
+ /* fpclassify is a macro, don't give it constants as parameter */
+ check_bool ("fpclassify (NaN) == FP_NAN", fpclassify (nan_value) == FP_NAN);
+ check_bool ("fpclassify (+inf) == FP_INFINITE",
+ fpclassify (plus_infty) == FP_INFINITE);
+ check_bool ("fpclassify (-inf) == FP_INFINITE",
+ fpclassify (minus_infty) == FP_INFINITE);
+ check_bool ("fpclassify (+0) == FP_ZERO",
+ fpclassify (plus_zero) == FP_ZERO);
+ check_bool ("fpclassify (-0) == FP_ZERO",
+ fpclassify (minus_zero) == FP_ZERO);
+
+ x = 1000.0;
+ check_bool ("fpclassify (1000) == FP_NORMAL",
+ fpclassify (x) == FP_NORMAL);
+}
+
+
+static void
+ldexp_test (void)
+{
+ check ("ldexp (0, 0) == 0", FUNC(ldexp) (0, 0), 0);
+
+ check_isinfp ("ldexp (+inf, 1) == +inf", FUNC(ldexp) (plus_infty, 1));
+ check_isinfn ("ldexp (-inf, 1) == -inf", FUNC(ldexp) (minus_infty, 1));
+ check_isnan ("ldexp (NaN, 1) == NaN", FUNC(ldexp) (nan_value, 1));
+
+ check ("ldexp (0.8, 4) == 12.8", FUNC(ldexp) (0.8L, 4), 12.8L);
+ check ("ldexp (-0.854375, 5) == -27.34", FUNC(ldexp) (-0.854375L, 5), -27.34L);
+}
+
+
+static void
+log_test (void)
+{
+ check_isinfn_exc ("log (+0) == -inf", FUNC(log) (0),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("log (-0) == -inf", FUNC(log) (minus_zero),
+ DIVIDE_BY_ZERO_EXCEPTION);
+
+ check ("log (1) == 0", FUNC(log) (1), 0);
+
+ check_isnan_exc ("log (x) == NaN plus divide-by-zero exception if x < 0",
+ FUNC(log) (-1), INVALID_EXCEPTION);
+ check_isinfp ("log (+inf) == +inf", FUNC(log) (plus_infty));
+
+ check_eps ("log (e) == 1", FUNC(log) (M_E), 1, CHOOSE (0, 0, 9e-8L));
+ check ("log (1/e) == -1", FUNC(log) (1.0 / M_E), 1);
+ check ("log (2) == M_LN2", FUNC(log) (2), M_LN2);
+ check ("log (10) == M_LN10", FUNC(log) (10), M_LN10);
+}
+
+
+static void
+log10_test (void)
+{
+ check_isinfn_exc ("log10 (+0) == -inf", FUNC(log10) (0),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("log10 (-0) == -inf", FUNC(log10) (minus_zero),
+ DIVIDE_BY_ZERO_EXCEPTION);
+
+ check ("log10 (1) == +0", FUNC(log10) (1), 0);
+
+ check_isnan_exc ("log10 (x) == NaN plus divide-by-zero exception if x < 0",
+ FUNC(log10) (-1), INVALID_EXCEPTION);
+
+ check_isinfp ("log10 (+inf) == +inf", FUNC(log10) (plus_infty));
+
+ check ("log10 (0.1) == -1", FUNC(log10) (0.1L), -1);
+ check ("log10 (10) == 1", FUNC(log10) (10.0), 1);
+ check ("log10 (100) == 2", FUNC(log10) (100.0), 2);
+ check ("log10 (10000) == 4", FUNC(log10) (10000.0), 4);
+ check_eps ("log10 (e) == M_LOG10E", FUNC(log10) (M_E), M_LOG10E,
+ CHOOSE (0, 0, 9e-8));
+}
+
+
+static void
+log1p_test (void)
+{
+ check ("log1p (+0) == +0", FUNC(log1p) (0), 0);
+ check ("log1p (-0) == -0", FUNC(log1p) (minus_zero), minus_zero);
+
+ check_isinfn_exc ("log1p (-1) == -inf", FUNC(log1p) (-1),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isnan_exc ("log1p (x) == NaN plus divide-by-zero exception if x < -1",
+ FUNC(log1p) (-2), INVALID_EXCEPTION);
+
+ check_isinfp ("log1p (+inf) == +inf", FUNC(log1p) (plus_infty));
+
+ check ("log1p (e-1) == 1", FUNC(log1p) (M_E - 1.0), 1);
+
+}
+
+
+#ifdef ISO_9X_IMPLEMENTED
+static void
+log2_test (void)
+{
+ check_isinfn_exc ("log2 (+0) == -inf", FUNC(log2) (0),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("log2 (-0) == -inf", FUNC(log2) (minus_zero),
+ DIVIDE_BY_ZERO_EXCEPTION);
+
+ check ("log2 (1) == +0", FUNC(log2) (1), 0);
+
+ check_isnan_exc ("log2 (x) == NaN plus divide-by-zero exception if x < 0",
+ FUNC(log2) (-1), INVALID_EXCEPTION);
+
+ check_isinfp ("log2 (+inf) == +inf", FUNC(log2) (plus_infty));
+
+ check ("log2 (e) == M_LOG2E", FUNC(log2) (M_E), M_LOG2E);
+ check ("log2 (2) == 1", FUNC(log2) (2.0), 1);
+ check ("log2 (16) == 4", FUNC(log2) (16.0), 4);
+ check ("log2 (256) == 8", FUNC(log2) (256.0), 8);
+
+}
+#endif
+
+
+static void
+logb_test (void)
+{
+ check_isinfp ("logb (+inf) == +inf", FUNC(logb) (plus_infty));
+ check_isinfp ("logb (-inf) == +inf", FUNC(logb) (minus_infty));
+
+ check_isinfn_exc ("logb (+0) == -inf plus divide-by-zero exception",
+ FUNC(logb) (0), DIVIDE_BY_ZERO_EXCEPTION);
+
+ check_isinfn_exc ("logb (-0) == -inf plus divide-by-zero exception",
+ FUNC(logb) (minus_zero), DIVIDE_BY_ZERO_EXCEPTION);
+
+ check ("logb (1) == 0", FUNC(logb) (1), 0);
+ check ("logb (e) == 1", FUNC(logb) (M_E), 1);
+ check ("logb (1024) == 10", FUNC(logb) (1024), 10);
+ check ("logb (-2000) == 10", FUNC(logb) (-2000), 10);
+
+}
+
+
+static void
+sin_test (void)
+{
+ check ("sin (+0) == +0", FUNC(sin) (0), 0);
+ check ("sin (-0) == -0", FUNC(sin) (minus_zero), minus_zero);
+ check_isnan_exc ("sin (+inf) == NaN plus invalid exception",
+ FUNC(sin) (plus_infty),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("sin (-inf) == NaN plus invalid exception",
+ FUNC(sin) (minus_infty),
+ INVALID_EXCEPTION);
+
+ check ("sin (pi/6) == 0.5", FUNC(sin) (M_PI / 6.0), 0.5);
+ check ("sin (pi/2) == 1", FUNC(sin) (M_PI_2), 1);
+}
+
+
+static void
+sinh_test (void)
+{
+ check ("sinh (+0) == +0", FUNC(sinh) (0), 0);
+ check ("sinh (-0) == -0", FUNC(sinh) (minus_zero), minus_zero);
+
+ check_isinfp ("sinh (+inf) == +inf", FUNC(sinh) (plus_infty));
+ check_isinfn ("sinh (-inf) == -inf", FUNC(sinh) (minus_infty));
+}
+
+
+static void
+tan_test (void)
+{
+ check ("tan (+0) == -0", FUNC(tan) (0), 0);
+ check ("tan (-0) == -0", FUNC(tan) (minus_zero), minus_zero);
+ check_isnan_exc ("tan (+inf) == NaN plus invalid exception",
+ FUNC(tan) (plus_infty), INVALID_EXCEPTION);
+ check_isnan_exc ("tan (-inf) == NaN plus invalid exception",
+ FUNC(tan) (minus_infty), INVALID_EXCEPTION);
+
+ check_eps ("tan (pi/4) == 1", FUNC(tan) (M_PI_4), 1, CHOOSE (0, 1e-15L, 0));
+}
+
+
+static void
+tanh_test (void)
+{
+ check ("tanh (+0) == +0", FUNC(tanh) (0), 0);
+ check ("tanh (-0) == -0", FUNC(tanh) (minus_zero), minus_zero);
+
+ check ("tanh (+inf) == +1", FUNC(tanh) (plus_infty), 1);
+ check ("tanh (-inf) == -1", FUNC(tanh) (minus_infty), -1);
+}
+
+
+static void
+fabs_test (void)
+{
+ check ("fabs (+0) == +0", FUNC(fabs) (0), 0);
+ check ("fabs (-0) == +0", FUNC(fabs) (minus_zero), 0);
+
+ check_isinfp ("fabs (+inf) == +inf", FUNC(fabs) (plus_infty));
+ check_isinfp ("fabs (-inf) == +inf", FUNC(fabs) (minus_infty));
+
+ check ("fabs (+38) == 38", FUNC(fabs) (38.0), 38.0);
+ check ("fabs (-e) == e", FUNC(fabs) (-M_E), M_E);
+}
+
+
+static void
+floor_test (void)
+{
+ check ("floor (+0) == +0", FUNC(floor) (0.0), 0.0);
+ check ("floor (-0) == -0", FUNC(floor) (minus_zero), minus_zero);
+ check_isinfp ("floor (+inf) == +inf", FUNC(floor) (plus_infty));
+ check_isinfn ("floor (-inf) == -inf", FUNC(floor) (minus_infty));
+
+ check ("floor (pi) == 3", FUNC(floor) (M_PI), 3.0);
+ check ("floor (-pi) == -4", FUNC(floor) (-M_PI), 4.0);
+}
+
+
+static void
+hypot_report (const char *test_name, MATHTYPE computed, MATHTYPE expected)
+{
+ MATHTYPE diff;
+ int result;
+
+ result = check_equal (computed, expected, 0, &diff);
+
+ if (result)
+ {
+ if (verbose > 2)
+ printf ("Pass: %s\n", test_name);
+ }
+ else
+ {
+ if (verbose)
+ printf ("Fail: %s\n", test_name);
+ if (verbose > 1)
+ {
+ printf ("Result:\n");
+ printf (" is: %.20" PRINTF_EXPR, computed);
+ printf (" should be: %.20" PRINTF_EXPR, expected);
+ printf (" difference: %.20" PRINTF_EXPR "\n", diff);
+ }
+ noErrors++;
+ }
+ fpstack_test (test_name);
+ output_result (test_name, result,
+ computed, expected, diff, PRINT, PRINT);
+}
+
+
+static void
+hypot_test (void)
+{
+ MATHTYPE a = FUNC(hypot) (12.4L, 0.7L);
+
+ hypot_report ("hypot (x,y) == hypot (y,x)", FUNC(hypot) (0.7L, 12.4L), a);
+ hypot_report ("hypot (x,y) == hypot (-x,y)", FUNC(hypot) (-12.4L, 0.7L), a);
+ hypot_report ("hypot (x,y) == hypot (-y,x)", FUNC(hypot) (-0.7L, 12.4L), a);
+ hypot_report ("hypot (x,y) == hypot (-x,-y)", FUNC(hypot) (-12.4L, -0.7L), a);
+ hypot_report ("hypot (x,y) == hypot (-y,-x)", FUNC(hypot) (-0.7L, -12.4L), a);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (-0.7L, 0), 0.7L);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (0.7L, 0), 0.7L);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (-1.0L, 0), 1.0L);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (-1.0L, 0), 1.0L);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (-5.7e7L, 0), 5.7e7L);
+ check ("hypot (x,0) == fabs (x)", FUNC(hypot) (-5.7e7L, 0), 5.7e7L);
+}
+
+
+static void
+pow_test (void)
+{
+ MATHTYPE x;
+
+ check ("pow (+0, +0) == 1", FUNC(pow) (0, 0), 1);
+ check ("pow (+0, -0) == 1", FUNC(pow) (0, minus_zero), 1);
+ check ("pow (-0, +0) == 1", FUNC(pow) (minus_zero, 0), 1);
+ check ("pow (-0, -0) == 1", FUNC(pow) (minus_zero, minus_zero), 1);
+
+ check ("pow (+10, +0) == 1", FUNC(pow) (10, 0), 1);
+ check ("pow (+10, -0) == 1", FUNC(pow) (10, minus_zero), 1);
+ check ("pow (-10, +0) == 1", FUNC(pow) (-10, 0), 1);
+ check ("pow (-10, -0) == 1", FUNC(pow) (-10, minus_zero), 1);
+
+ check ("pow (NaN, +0) == 1", FUNC(pow) (nan_value, 0), 1);
+ check ("pow (NaN, -0) == 1", FUNC(pow) (nan_value, minus_zero), 1);
+
+ check_isinfp ("pow (+1.1, +inf) == +inf", FUNC(pow) (1.1, plus_infty));
+ check_isinfp ("pow (+inf, +inf) == +inf", FUNC(pow) (plus_infty, plus_infty));
+ check_isinfp ("pow (-1.1, +inf) == +inf", FUNC(pow) (-1.1, plus_infty));
+ check_isinfp ("pow (-inf, +inf) == +inf", FUNC(pow) (minus_infty, plus_infty));
+
+ check ("pow (0.9, +inf) == +0", FUNC(pow) (0.9L, plus_infty), 0);
+ check ("pow (1e-7, +inf) == +0", FUNC(pow) (1e-7L, plus_infty), 0);
+ check ("pow (-0.9, +inf) == +0", FUNC(pow) (-0.9L, plus_infty), 0);
+ check ("pow (-1e-7, +inf) == +0", FUNC(pow) (-1e-7L, plus_infty), 0);
+
+ check ("pow (+1.1, -inf) == +inf", FUNC(pow) (1.1, minus_infty), 0);
+ check ("pow (+inf, -inf) == +inf", FUNC(pow) (plus_infty, minus_infty), 0);
+ check ("pow (-1.1, -inf) == +inf", FUNC(pow) (-1.1, minus_infty), 0);
+ check ("pow (-inf, -inf) == +inf", FUNC(pow) (minus_infty, minus_infty), 0);
+
+ check_isinfp ("pow (0.9, -inf) == +0", FUNC(pow) (0.9L, minus_infty));
+ check_isinfp ("pow (1e-7, -inf) == +0", FUNC(pow) (1e-7L, minus_infty));
+ check_isinfp ("pow (-0.9, -inf) == +0", FUNC(pow) (-0.9L, minus_infty));
+ check_isinfp ("pow (-1e-7, -inf) == +0", FUNC(pow) (-1e-7L, minus_infty));
+
+ check_isinfp ("pow (+inf, 1e-7) == +inf", FUNC(pow) (plus_infty, 1e-7L));
+ check_isinfp ("pow (+inf, 1) == +inf", FUNC(pow) (plus_infty, 1));
+ check_isinfp ("pow (+inf, 1e7) == +inf", FUNC(pow) (plus_infty, 1e7L));
+
+ check ("pow (+inf, -1e-7) == 0", FUNC(pow) (plus_infty, -1e-7L), 0);
+ check ("pow (+inf, -1) == 0", FUNC(pow) (plus_infty, -1), 0);
+ check ("pow (+inf, -1e7) == 0", FUNC(pow) (plus_infty, -1e7L), 0);
+
+ check_isinfn ("pow (-inf, 1) == -inf", FUNC(pow) (minus_infty, 1));
+ check_isinfn ("pow (-inf, 11) == -inf", FUNC(pow) (minus_infty, 11));
+ check_isinfn ("pow (-inf, 1001) == -inf", FUNC(pow) (minus_infty, 1001));
+
+ check_isinfp ("pow (-inf, 2) == +inf", FUNC(pow) (minus_infty, 2));
+ check_isinfp ("pow (-inf, 12) == +inf", FUNC(pow) (minus_infty, 12));
+ check_isinfp ("pow (-inf, 1002) == +inf", FUNC(pow) (minus_infty, 1002));
+ check_isinfp ("pow (-inf, 0.1) == +inf", FUNC(pow) (minus_infty, 0.1));
+ check_isinfp ("pow (-inf, 1.1) == +inf", FUNC(pow) (minus_infty, 1.1));
+ check_isinfp ("pow (-inf, 11.1) == +inf", FUNC(pow) (minus_infty, 11.1));
+ check_isinfp ("pow (-inf, 1001.1) == +inf", FUNC(pow) (minus_infty, 1001.1));
+
+ check ("pow (-inf, -1) == -0", FUNC(pow) (-minus_infty, -1), minus_zero);
+ check ("pow (-inf, -11) == -0", FUNC(pow) (-minus_infty, -11), minus_zero);
+ check ("pow (-inf, -1001) == -0", FUNC(pow) (-minus_infty, -1001), minus_zero);
+
+ check ("pow (-inf, -2) == +0", FUNC(pow) (minus_infty, -2), 0);
+ check ("pow (-inf, -12) == +0", FUNC(pow) (minus_infty, -12), 0);
+ check ("pow (-inf, -1002) == +0", FUNC(pow) (minus_infty, -1002), 0);
+ check ("pow (-inf, -0.1) == +0", FUNC(pow) (minus_infty, -0.1), 0);
+ check ("pow (-inf, -1.1) == +0", FUNC(pow) (minus_infty, -1.1), 0);
+ check ("pow (-inf, -11.1) == +0", FUNC(pow) (minus_infty, -11.1), 0);
+ check ("pow (-inf, -1001.1) == +0", FUNC(pow) (minus_infty, -1001.1), 0);
+
+ check_isnan ("pow (NaN, NaN) == NaN", FUNC(pow) (nan_value, nan_value));
+ check_isnan ("pow (0, NaN) == NaN", FUNC(pow) (0, nan_value));
+ check_isnan ("pow (1, NaN) == NaN", FUNC(pow) (1, nan_value));
+ check_isnan ("pow (-1, NaN) == NaN", FUNC(pow) (-1, nan_value));
+ check_isnan ("pow (NaN, 1) == NaN", FUNC(pow) (nan_value, 1));
+ check_isnan ("pow (NaN, -1) == NaN", FUNC(pow) (nan_value, -1));
+
+ x = random_greater (0.0);
+ check_isnan_ext ("pow (x, NaN) == NaN", FUNC(pow) (x, nan_value), x);
+
+ check_isnan_exc ("pow (+1, +inf) == NaN", FUNC(pow) (1, plus_infty),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (-1, +inf) == NaN", FUNC(pow) (1, plus_infty),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (+1, -inf) == NaN", FUNC(pow) (1, plus_infty),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (-1, -inf) == NaN", FUNC(pow) (1, plus_infty),
+ INVALID_EXCEPTION);
+
+ check_isnan_exc ("pow (-0.1, 1.1) == NaN", FUNC(pow) (-0.1, 1.1),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (-0.1, -1.1) == NaN", FUNC(pow) (-0.1, -1.1),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (-10.1, 1.1) == NaN", FUNC(pow) (-10.1, 1.1),
+ INVALID_EXCEPTION);
+ check_isnan_exc ("pow (-10.1, -1.1) == NaN", FUNC(pow) (-10.1, -1.1),
+ INVALID_EXCEPTION);
+
+ check_isinfp_exc ("pow (+0, -1) == +inf", FUNC(pow) (0, -1),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfp_exc ("pow (+0, -11) == +inf", FUNC(pow) (0, -11),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("pow (-0, -1) == -inf", FUNC(pow) (minus_zero, -1),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfn_exc ("pow (-0, -11) == -inf", FUNC(pow) (minus_zero, -11),
+ DIVIDE_BY_ZERO_EXCEPTION);
+
+ check_isinfp_exc ("pow (+0, -2) == +inf", FUNC(pow) (0, -2),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfp_exc ("pow (+0, -11.1) == +inf", FUNC(pow) (0, -11.1),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfp_exc ("pow (-0, -2) == +inf", FUNC(pow) (minus_zero, -2),
+ DIVIDE_BY_ZERO_EXCEPTION);
+ check_isinfp_exc ("pow (-0, -11.1) == +inf", FUNC(pow) (minus_zero, -11.1),
+ DIVIDE_BY_ZERO_EXCEPTION);
+
+ check ("pow (+0, 1) == +0", FUNC(pow) (0, 1), 0);
+ check ("pow (+0, 11) == +0", FUNC(pow) (0, 11), 0);
+ check ("pow (-0, 1) == -0", FUNC(pow) (minus_zero, 1), minus_zero);
+ check ("pow (-0, 11) == -0", FUNC(pow) (minus_zero, 11), minus_zero);
+
+ check ("pow (+0, 2) == +0", FUNC(pow) (0, 2), 0);
+ check ("pow (+0, 11.1) == +0", FUNC(pow) (0, 11.1), 0);
+ check ("pow (-0, 2) == +0", FUNC(pow) (minus_zero, 2), 0);
+ check ("pow (-0, 11.1) == +0", FUNC(pow) (minus_zero, 11.1), 0);
+
+ x = random_greater (0.0);
+ check_isinfp_ext ("pow (x, +inf) == +inf for |x| > 1",
+ FUNC(pow) (x, plus_infty), x);
+
+ x = random_value (-1.0, 1.0);
+ check_ext ("pow (x, +inf) == +0 for |x| < 1",
+ FUNC(pow) (x, plus_infty), 0.0, x);
+
+ x = random_greater (0.0);
+ check_ext ("pow (x, -inf) == +0 for |x| > 1",
+ FUNC(pow) (x, minus_infty), 0.0, x);
+
+ x = random_value (-1.0, 1.0);
+ check_isinfp_ext ("pow (x, -inf) == +inf for |x| < 1",
+ FUNC(pow) (x, minus_infty), x);
+
+ x = random_greater (0.0);
+ check_isinfp_ext ("pow (+inf, y) == +inf for y > 0",
+ FUNC(pow) (plus_infty, x), x);
+
+ x = random_less (0.0);
+ check_ext ("pow (+inf, y) == +0 for y < 0",
+ FUNC(pow) (plus_infty, x), 0.0, x);
+
+ x = (rand () % 1000000) * 2.0 + 1; /* Get random odd integer > 0 */
+ check_isinfn_ext ("pow (-inf, y) == -inf for y an odd integer > 0",
+ FUNC(pow) (minus_infty, x), x);
+
+ x = ((rand () % 1000000) + 1) * 2.0; /* Get random even integer > 1 */
+ check_isinfp_ext ("pow (-inf, y) == +inf for y > 0 and not an odd integer",
+ FUNC(pow) (minus_infty, x), x);
+
+ x = -((rand () % 1000000) * 2.0 + 1); /* Get random odd integer < 0 */
+ check_ext ("pow (-inf, y) == -0 for y an odd integer < 0",
+ FUNC(pow) (minus_infty, x), minus_zero, x);
+
+ x = ((rand () % 1000000) + 1) * -2.0; /* Get random even integer < 0 */
+ check_ext ("pow (-inf, y) == -0 for y < 0 and not an odd integer",
+ FUNC(pow) (minus_infty, x), minus_zero, x);
+
+ x = (rand () % 1000000) * 2.0 + 1; /* Get random odd integer > 0 */
+ check_ext ("pow (+0, y) == +0 for y an odd integer > 0",
+ FUNC(pow) (0.0, x), 0.0, x);
+ x = (rand () % 1000000) * 2.0 + 1; /* Get random odd integer > 0 */
+ check_ext ("pow (-0, y) == -0 for y an odd integer > 0",
+ FUNC(pow) (minus_zero, x), minus_zero, x);
+
+ x = ((rand () % 1000000) + 1) * 2.0; /* Get random even integer > 1 */
+ check_ext ("pow (+0, y) == +0 for y > 0 and not an odd integer",
+ FUNC(pow) (0.0, x), 0.0, x);
+
+ x = ((rand () % 1000000) + 1) * 2.0; /* Get random even integer > 1 */
+ check_ext ("pow (-0, y) == +0 for y > 0 and not an odd integer",
+ FUNC(pow) (minus_zero, x), 0.0, x);
+}
+
+
+static void
+inverse_func_pair_test (const char *test_name,
+ mathfunc f1, mathfunc inverse,
+ MATHTYPE x, MATHTYPE epsilon)
+{
+ MATHTYPE a, b, difference;
+ int result;
+
+ a = f1 (x);
+ this_does_nothing ();
+ b = inverse (a);
+ this_does_nothing ();
+
+ result = check_equal (b, x, epsilon, &difference);
+ output_result (test_name, result,
+ b, x, difference, PRINT, PRINT);
+}
+
+
+static void
+inverse_functions (void)
+{
+ inverse_func_pair_test ("(asin(sin(x)) == x",
+ FUNC(sin), FUNC(asin), 1.0, CHOOSE (0, 0, 1e-7L));
+ inverse_func_pair_test ("(sin(asin(x)) == x",
+ FUNC(asin), FUNC(sin), 1.0, 0.0);
+
+ inverse_func_pair_test ("(acos(cos(x)) == x",
+ FUNC(cos), FUNC(acos), 1.0, CHOOSE (0, 1e-15L, 0));
+ inverse_func_pair_test ("(cos(acos(x)) == x",
+ FUNC(acos), FUNC(cos), 1.0, 0.0);
+ inverse_func_pair_test ("(atan(tan(x)) == x",
+ FUNC(tan), FUNC(atan), 1.0, 0.0);
+ inverse_func_pair_test ("(tan(atan(x)) == x",
+ FUNC(atan), FUNC(tan), 1.0, CHOOSE (0, 1e-15L, 0));
+
+ inverse_func_pair_test ("(asinh(sinh(x)) == x",
+ FUNC(sinh), FUNC(asinh), 1.0, CHOOSE (1e-18L, 0, 0));
+ inverse_func_pair_test ("(sinh(asinh(x)) == x",
+ FUNC(asinh), FUNC(sinh), 1.0, 0.0);
+
+ inverse_func_pair_test ("(acosh(cosh(x)) == x",
+ FUNC(cosh), FUNC(acosh), 1.0, CHOOSE (1e-18L, 1e-15L, 0));
+ inverse_func_pair_test ("(cosh(acosh(x)) == x",
+ FUNC(acosh), FUNC(cosh), 1.0, 0.0);
+
+ inverse_func_pair_test ("(atanh(tanh(x)) == x",
+ FUNC(tanh), FUNC(atanh), 1.0, CHOOSE (0, 1e-15L, 0));
+ inverse_func_pair_test ("(tanh(atanh(x)) == x",
+ FUNC(atanh), FUNC(tanh), 1.0, 0.0);
+
+}
+
+/* Test sin and cos with the identity: sin(x)^2 + cos(x)^2 = 1. */
+static void
+identities1_test (MATHTYPE x, MATHTYPE epsilon)
+{
+ MATHTYPE res1, res2, res3, diff;
+ int result;
+
+ res1 = FUNC(sin) (x);
+ this_does_nothing ();
+ res2 = FUNC(cos) (x);
+ this_does_nothing ();
+ res3 = res1 * res1 + res2 * res2;
+ this_does_nothing ();
+
+ result = check_equal (res3, 1.0, epsilon, &diff);
+ output_result_ext ("sin^2 + cos^2 == 1", result,
+ res3, 1.0, diff, x, PRINT, PRINT);
+}
+
+
+/* Test sin, cos, tan with the following relation: tan = sin/cos. */
+static void
+identities2_test (MATHTYPE x, MATHTYPE epsilon)
+{
+ MATHTYPE res1, res2, res3, res4, diff;
+ int result;
+
+ res1 = FUNC(sin) (x);
+ this_does_nothing ();
+ res2 = FUNC(cos) (x);
+ this_does_nothing ();
+ res3 = FUNC(tan) (x);
+ this_does_nothing ();
+ res4 = res1 / res2;
+ this_does_nothing ();
+
+ result = check_equal (res4, res3, epsilon, &diff);
+ output_result_ext ("sin/cos == tan", result,
+ res4, res3, diff, x, PRINT, PRINT);
+}
+
+
+/* Test cosh and sinh with the identity cosh^2 - sinh^2 = 1. */
+static void
+identities3_test (MATHTYPE x, MATHTYPE epsilon)
+{
+ MATHTYPE res1, res2, res3, diff;
+ int result;
+
+ res1 = FUNC(sinh) (x);
+ this_does_nothing ();
+ res2 = FUNC(cosh) (x);
+ this_does_nothing ();
+ res3 = res2 * res2 - res1 * res1;
+ this_does_nothing ();
+
+ result = check_equal (res3, 1.0, epsilon, &diff);
+ output_result_ext ("cosh^2 - sinh^2 == 1", result,
+ res3, 1.0, diff, x, PRINT, PRINT);
+}
+
+
+static void
+identities (void)
+{
+ identities1_test (0.2L, CHOOSE (1e-19L, 0, 0));
+ identities1_test (0.9L, 0);
+ identities1_test (0, 0);
+ identities1_test (-1, CHOOSE (0, 0, 1e-7));
+
+ identities2_test (0.2L, CHOOSE (0, 1e-16, 0));
+ identities2_test (0.9L, CHOOSE (0, 1e-15, 0));
+ identities2_test (0, 0);
+ identities2_test (-1, CHOOSE (1e-18L, 1e-15, 0));
+
+ identities3_test (0.2L, CHOOSE (0, 0, 1e-7));
+ identities3_test (0.9L, CHOOSE (1e-18L, 1e-15, 1e-6));
+ identities3_test (0, CHOOSE (0, 0, 1e-6));
+ identities3_test (-1, CHOOSE (1e-18L, 0, 1e-6));
+}
+
+
+/*
+ Let's test that basic arithmetic is working
+ tests: Infinity and NaN
+ */
+static void
+basic_tests (void)
+{
+ /* variables are declared volatile to forbid some compiler
+ optimizations */
+ volatile MATHTYPE Inf_var, NaN_var, zero_var, one_var;
+ MATHTYPE x1, x2;
+
+ zero_var = 0.0;
+ one_var = 1.0;
+ NaN_var = nan_value;
+ Inf_var = one_var / zero_var;
+
+ this_does_nothing ();
+
+ check_isinfp ("isinf (1/0) == +1", Inf_var);
+ check_isinfn ("isinf (-1/0) == -1", -Inf_var);
+ check_bool ("!isinf (1)", !(FUNC(isinf) (one_var)));
+ check_bool ("!isinf (0/0)", !(FUNC(isinf) (NaN_var)));
+
+ check_isnan ("isnan (0/0)", NaN_var);
+ check_isnan ("isnan (-(0/0))", -NaN_var);
+ check_bool ("!isnan (1)", !(FUNC(isnan) (one_var)));
+ check_bool ("!isnan (0/0)", !(FUNC(isnan) (Inf_var)));
+
+ check_bool ("inf == inf", Inf_var == Inf_var);
+ check_bool ("-inf == -inf", -Inf_var == -Inf_var);
+ check_bool ("inf != -inf", Inf_var != -Inf_var);
+ check_bool ("NaN != NaN", NaN_var != NaN_var);
+
+ /*
+ the same tests but this time with NAN from <nan.h>
+ NAN is a double const
+ */
+ check_bool ("isnan (NAN)", isnan (NAN));
+ check_bool ("isnan (-NAN)", isnan (-NAN));
+ check_bool ("!isinf (NAN)", !(isinf (NAN)));
+ check_bool ("!isinf (-NAN)", !(isinf (-NAN)));
+ check_bool ("NAN != NAN", NAN != NAN);
+
+ /* test if EPSILON is ok */
+ x1 = MATHCONST (1.0);
+ x2 = x1 + CHOOSE (LDBL_EPSILON, DBL_EPSILON, FLT_EPSILON);
+ check_bool ("1 != 1+EPSILON", x1 != x2);
+
+ x1 = MATHCONST (1.0);
+ x2 = x1 - CHOOSE (LDBL_EPSILON, DBL_EPSILON, FLT_EPSILON);
+ check_bool ("1 != 1-EPSILON", x1 != x2);
+
+ /* test if HUGE_VALx is ok */
+ x1 = CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF);
+ check_bool ("isinf (HUGE_VALx) == +1", ISINF (x1) == +1);
+ x1 = -CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF);
+ check_bool ("isinf (-HUGE_VALx) == -1", ISINF (x1) == -1);
+
+}
+
+
+static void
+initialize (void)
+{
+ plus_zero = 0.0;
+ nan_value = plus_zero / plus_zero; /* Suppress GCC warning */
+
+ minus_zero = copysign (0.0, -1.0);
+ plus_infty = CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF);
+ minus_infty = -CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF);
+
+ /* Test to make sure we start correctly. */
+ fpstack_test ("*init*");
+}
+
+
+static struct option long_options[] =
+{
+ {"verbose", optional_argument, NULL, 'v'},
+ {"silent", no_argument, NULL, 's'},
+ {0, 0, 0, 0}
+};
+
+
+static void
+parse_options (int argc, char *argv[])
+{
+ int c;
+ int option_index;
+
+ verbose = 1;
+
+ while (1)
+ {
+ c = getopt_long (argc, argv, "vs",
+ long_options, &option_index);
+
+ /* Detect the end of the options. */
+ if (c == -1)
+ break;
+
+ switch (c)
+ {
+ case 'v':
+ if (optarg)
+ verbose = (unsigned int) strtoul (optarg, NULL, 0);
+ else
+ verbose = 3;
+ break;
+ case 's':
+ verbose = 0;
+ default:
+ break;
+ }
+ }
+}
+
+
+int
+main (int argc, char *argv[])
+{
+ parse_options (argc, argv);
+
+ initialize ();
+ printf (TEST_MSG);
+ basic_tests ();
+
+ acos_test ();
+ acosh_test ();
+ asin_test ();
+ asinh_test ();
+ atan_test ();
+ atanh_test ();
+ atan2_test ();
+ cbrt_test ();
+ ceil_test ();
+ cos_test ();
+ cosh_test ();
+ exp_test ();
+#ifdef ISO_9X_IMPLEMENTED
+ exp2_test ();
+#endif
+ expm1_test ();
+ frexp_test ();
+ ldexp_test ();
+ log_test ();
+ log10_test ();
+ log1p_test ();
+#ifdef ISO_9X_IMPLEMENTED
+ log2_test ();
+#endif
+ logb_test ();
+ sin_test ();
+ sinh_test ();
+ tan_test ();
+ tanh_test ();
+ fabs_test ();
+ floor_test ();
+ fpclassify_test ();
+ hypot_test ();
+ pow_test ();
+
+ identities ();
+ inverse_functions ();
+
+ if (noErrors)
+ {
+ printf ("\n%d errors occured.\n", noErrors);
+ exit (1);
+ }
+ printf ("\n All tests passed sucessfully.\n");
+ exit (0);
+}
diff --git a/math/math.h b/math/math.h
index 22274f023e..89c1cfa9bb 100644
--- a/math/math.h
+++ b/math/math.h
@@ -128,9 +128,9 @@ enum
/* Return number of classification appropriate for X. */
#define fpclassify(x) \
(sizeof (x) == sizeof (float) ? \
- __finitef (x) \
+ __fpclassifyf (x) \
: sizeof (x) == sizeof (double) ? \
- __finite (x) : __finitel (x))
+ __fpclassify (x) : __fpclassifyl (x))
/* Return nonzero value if sign of X is negative. */
#define signbit(x) \