diff options
author | Ulrich Drepper <drepper@redhat.com> | 1997-03-16 20:28:07 +0000 |
---|---|---|
committer | Ulrich Drepper <drepper@redhat.com> | 1997-03-16 20:28:07 +0000 |
commit | 0d8733c4fc3695c7987548d10c344445f3eb552d (patch) | |
tree | 4a2fe41c465963794a6f9d4d1717494972336e14 /math | |
parent | 1ed4222e5724ce863a95c15eb530dfab7ef0451c (diff) | |
download | glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.zip glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.tar.gz glibc-0d8733c4fc3695c7987548d10c344445f3eb552d.tar.bz2 |
Update.cvs/libc-ud-970316cvs/libc-970324cvs/libc-970323cvs/libc-970322cvs/libc-970321cvs/libc-970320cvs/libc-970319cvs/libc-970318cvs/libc-970317cvs/before-thomas-posix1996cvs/thomas-posix1996
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.c | 1520 | ||||
-rw-r--r-- | math/math.h | 4 |
2 files changed, 1522 insertions, 2 deletions
diff --git a/math/libm-test.c b/math/libm-test.c new file mode 100644 index 0000000..98e3cba --- /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 22274f0..89c1cfa 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) \ |