glibc/math/test-misc.c

1320 lines
33 KiB
C
Raw Normal View History

/* Miscellaneous tests which don't fit anywhere else.
Copyright (C) 2000, 2001, 2004, 2005, 2007 Free Software Foundation, Inc.
This file is part of the GNU C Library.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 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
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, write to the Free
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
02111-1307 USA. */
#include <fenv.h>
2001-01-27 20:09:14 +01:00
#include <float.h>
#include <ieee754.h>
#include <math.h>
#include <stdio.h>
#include <string.h>
int
main (void)
{
int result = 0;
#ifndef NO_LONG_DOUBLE
{
long double x = 0x100000001ll + (long double) 0.5;
long double q;
long double r;
r = modfl (x, &q);
if (q != (long double) 0x100000001ll || r != 0.5)
{
printf ("modfl (%Lg, ...) failed\n", x);
result = 1;
}
}
{
long double x;
long double m;
long double r;
int e;
int i;
# if LDBL_MANT_DIG == 64
m = 0xf.fffffffffffffffp-4L;
# elif LDBL_MANT_DIG == 106
Update. 2004-05-03 Jakub Jelinek <jakub@redhat.com> * sysdeps/unix/bsd/bits/fcntl.h (F_SETOWN, F_GETOWN): Define if __USE_BSD or __USE_UNIX98. * sysdeps/unix/bsd/ultrix4/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/bsd/bsd4.4/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/bsd/sun/sunos4/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/common/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/aix/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/irix4/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/alpha/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/s390/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/cris/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. *sysdeps/unix/sysv/linux/powerpc/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/sparc/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/mips/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/sh/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/i386/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/m68k/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/ia64/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/arm/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/hppa/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/unix/sysv/linux/x86_64/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/generic/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * sysdeps/mach/hurd/bits/fcntl.h (F_SETOWN, F_GETOWN): Likewise. * io/sys/stat.h (S_ISSOCK, S_IFSOCK): Likewise. 2004-05-01 Jakub Jelinek <jakub@redhat.com> * posix/Versions (libc): Remove __libc_wait, __libc_waitpid, __libc_pause, __libc_nanosleep, __libc_pread, __libc_pread64, __libc_pwrite64, __waitid and __pselect @@GLIBC_PRIVATE. * stdlib/Versions (libc): Remove __on_exit@@GLIBC_PRIVATE. * sysdeps/unix/sysv/linux/Versions (libc): Remove __libc_sigaction@@GLIBC_PRIVATE. * sysdeps/unix/sysv/linux/x86_64/Versions (libc): Remove __modify_ldt@@GLIBC_PRIVATE. * socket/Versions (libc): Remove __libc_accept, __libc_send, __libc_recvfrom, __libc_recvmsg, __libc_sendmsg, __libc_recv, __libc_sendto and __libc_connect @@GLIBC_PRIVATE. * stdio-common/Versions (libc): Remove _itoa_upper_digits@@GLIBC_PRIVATE. * resolv/Versions (libresolv): Remove __ns_samename@@GLIBC_PRIVATE. * misc/Versions (libc): Remove __libc_fsync, __libc_msync, __libc_readv and __libc_writev @@GLIBC_PRIVATE. * termios/Versions (libc): Remove __libc_tcdrain@@GLIBC_PRIVATE. * io/Versions (libc): Remove __libc_open, __libc_close, __libc_read, __libc_write, __libc_lseek, __libc_fcntl, __libc_open64 and __libc_lseek64 @@GLIBC_PRIVATE. 2004-04-30 Jakub Jelinek <jakub@redhat.com> * elf/dl-load.c (open_verify): Move e_phentsize check after e_type check. 2004-04-29 Steven Munroe <sjmunroe@us.ibm.com> * sysdeps/powerpc/powerpc64/fpu/s_ceil.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_ceilf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_floor.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_floorf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_llrint.c: Removed. * sysdeps/powerpc/powerpc64/fpu/s_llrint.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_llrintf.c: Removed. * sysdeps/powerpc/powerpc64/fpu/s_llrintf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_llround.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_llroundf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_lround.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_lroundf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_rint.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_rintf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_round.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_roundf.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_trunc.S: New file. * sysdeps/powerpc/powerpc64/fpu/s_truncf.S: New file. * math/test-misc.c [LDBL_MANT_DIG == 106](main): Correct LDBL_MAX mantissa for AIX long double format. * misc/qefgcvt.c [LDBL_MANT_DIG == 106] (NDIGIT_MAX): Define for AIX long double format. * misc/qefgcvt_r.c [LDBL_MANT_DIG == 106] (NDIGIT_MAX): Likewise. * stdlib/fpioconst.c [!__NO_LONG_DOUBLE_MATH && __LDBL_MAX_EXP__>1024] (_fpioconst_pow10): AIX long double format has same exponent range as double. * stdlib/fpioconst.h [!__NO_LONG_DOUBLE_MATH && __LDBL_MAX_EXP__>1024] (LDBL_MAX_10_EXP_LOG): AIX long double format has same exponent range as double. 2004-04-23 Art Haas <ahaas@airmail.net> * sysdeps/unix/sysv/linux/kernel-features.h: Add 'defined'.
2004-05-03 23:40:12 +02:00
/* This has to match the mantissa of LDBL_MAX which actually does have a
missing bit in the middle. */
m = 0x1.fffffffffffff7ffffffffffff8p-1L;
# elif LDBL_MANT_DIG == 113
m = 0x1.ffffffffffffffffffffffffffffp-1L;
# else
# error "Please adjust"
# endif
for (i = LDBL_MAX_EXP, x = LDBL_MAX; i >= LDBL_MIN_EXP; --i, x /= 2.0L)
{
printf ("2^%d: ", i);
r = frexpl (x, &e);
if (r != m)
{
printf ("mantissa incorrect: %.20La\n", r);
result = 1;
continue;
}
if (e != i)
{
printf ("exponent wrong %d (%.20Lg)\n", e, x);
result = 1;
continue;
}
puts ("ok");
}
Update. 2001-06-04 Bruno Haible <haible@clisp.cons.org> * iconv/loop.c (UNICODE_TAG_HANDLER): New macro. * iconv/gconv_simple.c (__gconv_transform_internal_ascii): Invoke UNICODE_TAG_HANDLER. (__gconv_transform_internal_ucs2): Likewise. (__gconv_transform_internal_ucs2reverse): Likewise. * iconvdata/8bit-gap.c (BODY for TO_LOOP): Invoke UNICODE_TAG_HANDLER. * iconvdata/8bit-generic.c (BODY for TO_LOOP): Likewise. * iconvdata/ansi_x3.110.c (BODY for TO_LOOP): Likewise. * iconvdata/big5.c (BODY for TO_LOOP): Likewise. * iconvdata/big5hkscs.c (BODY for TO_LOOP): Likewise. * iconvdata/cp1255.c (BODY for TO_LOOP): Likewise. * iconvdata/cp1258.c (BODY for TO_LOOP): Likewise. * iconvdata/euc-cn.c (BODY for TO_LOOP): Likewise. * iconvdata/euc-jp.c (BODY for TO_LOOP): Likewise. * iconvdata/euc-kr.c (BODY for TO_LOOP): Likewise. * iconvdata/euc-tw.c (BODY for TO_LOOP): Likewise. * iconvdata/gbk.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm930.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm932.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm933.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm935.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm937.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm939.c (BODY for TO_LOOP): Likewise. * iconvdata/ibm943.c (BODY for TO_LOOP): Likewise. * iconvdata/iso646.c (BODY for TO_LOOP): Likewise. * iconvdata/iso8859-1.c (BODY for TO_LOOP): Likewise. * iconvdata/iso_6937.c (BODY for TO_LOOP): Likewise. * iconvdata/iso_6937-2.c (BODY for TO_LOOP): Likewise. * iconvdata/iso-2022-cn.c (BODY for TO_LOOP): Likewise. * iconvdata/iso-2022-cn-ext.c (BODY for TO_LOOP): Likewise. * iconvdata/iso-2022-kr.c (BODY for TO_LOOP): Likewise. * iconvdata/johab.c (BODY for TO_LOOP): Likewise. * iconvdata/sjis.c (BODY for TO_LOOP): Likewise. * iconvdata/t.61.c (BODY for TO_LOOP): Likewise. * iconvdata/uhc.c (BODY for TO_LOOP): Likewise. * iconvdata/unicode.c (BODY for TO_LOOP): Likewise. * iconvdata/iso-2022-jp.c (TAG_none, TAG_language, TAG_language_j, TAG_language_ja, TAG_language_k, TAG_language_ko, TAG_language_z, TAG_language_zh, CURRENT_TAG_MASK): New enum values. (EMIT_SHIFT_TO_INIT): Don't emit an escape sequence if ASCII_set is already selected but set2 or tag are set. (conversion): New enum type. (cvlist_t): New type. (CVLIST, CVLIST_FIRST, CVLIST_REST): New macros. (conversion_lists): New array. (BODY for TO_LOOP): Keep track of Unicode 3.1 language tag. If "ja", prefer conversion to Japanese character sets. If "zh", prefer conversion to GB2312. If "ko", prefer conversion to KSC5601. Small optimizations. (INIT_PARAMS): Add tag. (UPDATE_PARAMS): Add tag. 2001-06-04 Bruno Haible <haible@clisp.cons.org> * locale/programs/locfile.c (write_locale_data): Before creat(), unlink the file, to avoid crashing the processes that mmap it. Change a double slash to a single slash. Free fname in case of error return. 2001-06-02 Jakub Jelinek <jakub@redhat.com> * sysdeps/i386/fpu/s_frexpl.S (__frexpl): Mostly revert 2000-12-03 changes, do the special handling for denormal numbers, not for normalized numbers (patch by <trevin@xmission.com>). * math/test-misc.c (main): Test frexpl with denormal arguments. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * math/libm-test.inc (llround_test): Add two new llround tests. * sysdeps/ieee754/ldbl-96/s_llroundl.c (__llroundl): Don't allow overflow when rounding away from zero. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * math/Makefile (libm-calls): Add e_log2, w_log2, remove s_log2. * math/math_private.h (__ieee754_log2, __ieee754_log2f, __ieee754_log2l): New prototypes. * sysdeps/generic/w_log2.c: New file. * sysdeps/generic/w_log2f.c: New file. * sysdeps/generic/w_log2l.c: New file. * sysdeps/generic/s_log2l.c: Move... * sysdeps/generic/e_log2l.c: ...to here. Rename to __ieee754_log2l. * sysdeps/ieee754/k_standard.c (__kernel_standard): Handle log2(0) and log2(x < 0). * sysdeps/i386/fpu/s_log2.S: Move... * sysdeps/i386/fpu/e_log2.S: ...to here. Rename to __ieee754_log2. * sysdeps/i386/fpu/s_log2f.S: Move... * sysdeps/i386/fpu/e_log2f.S: ...to here. Rename to __ieee754_log2f. * sysdeps/i386/fpu/s_log2l.S: Move... * sysdeps/i386/fpu/e_log2l.S: ...to here. Rename to __ieee754_log2l. * sysdeps/m68k/fpu/s_log2.S: Move... * sysdeps/m68k/fpu/e_log2.S: ...to here. Rename to __ieee754_log2. * sysdeps/m68k/fpu/s_log2f.S: Move... * sysdeps/m68k/fpu/e_log2f.S: ...to here. Rename to __ieee754_log2f. * sysdeps/m68k/fpu/s_log2l.S: Move... * sysdeps/m68k/fpu/e_log2l.S: ...to here. Rename to __ieee754_log2l. * sysdeps/ieee754/dbl-64/s_log2.c: Move... * sysdeps/ieee754/dbl-64/e_log2.c: ...to here. Rename to __ieee754_log2. * sysdeps/ieee754/flt-32/s_log2f.c: Move... * sysdeps/ieee754/flt-32/e_log2f.c: ...to here. Rename to __ieee754_log2f. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * sysdeps/generic/w_exp2.c (u_threshold): Lower threshold so that even arguments which result in denormalized exp2 are accepted. (__exp2): Arguments equal to u_threshold already result into underflow. * sysdeps/generic/w_exp2f.c (u_threshold, __exp2f): Likewise. * sysdeps/generic/w_exp2l.c (u_threshold, __exp2l): Likewise. * sysdeps/ieee754/dbl-64/e_exp2.c (__ieee754_exp2): Lomark was too low, with corrected lowmark use greaterequal, not greater. * sysdeps/ieee754/flt-32/e_exp2f.c (__ieee754_exp2f): Likewise. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * math/libm-test.inc (ilogb_test): Test that ilogb(+-Inf) == INT_MAX. * sysdeps/i386/fpu/s_ilogb.S (__ilogb): Return INT_MAX for +-Inf. * sysdeps/i386/fpu/s_ilogbf.S (__ilogbf): Likewise. * sysdeps/i386/fpu/s_ilogbl.S (__ilogbl): Likewise. * sysdeps/ieee754/dbl-64/s_ilogb.c (__ilogb): Likewise. * sysdeps/ieee754/flt-32/s_ilogbf.c (__ilogbf): Likewise. * sysdeps/ieee754/ldbl-128/s_ilogbl.c (__ilogbl): Likewise. * sysdeps/ieee754/ldbl-96/s_ilogbl.c (__ilogbl): Likewise. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * sysdeps/generic/w_coshl.c (__coshl): Test if finite argument gave non-finite result instead of using constant in generic version. * sysdeps/generic/w_coshf.c (__coshf): Likewise. * sysdeps/generic/w_cosh.c (__cosh): Likewise. * sysdeps/generic/w_exp10.c (o_threshold, u_threshold): Remove. (__exp10): Test if finite argument gave non-finite result. * sysdeps/generic/w_exp10f.c (o_threshold, u_threshold, __exp10f): Likewise. * sysdeps/generic/w_exp10l.c (o_threshold, u_threshold, __exp10l): Likewise. 2001-06-04 Jakub Jelinek <jakub@redhat.com> * sysdeps/ieee754/ldbl-96/e_coshl.c (__ieee754_coshl): Fix overflow threshold constant (log(LDBL_MAX)+M_LN2l). 2001-05-29 Bruno Haible <haible@clisp.cons.org> * locale/programs/ld-ctype.c (idx_table): New struct type. (idx_table_init, idx_table_get, idx_table_add): New functions. (MAX_CHARNAMES_IDX): Remove macro. (locale_ctype_t): Change type of charnames_idx field. (ctype_startup): Change initialization of charnames_idx field. (find_idx): Use idx_table_get and idx_table_add for speed. * locale/programs/charmap.c (charmap_new_char): Fix ucs4 value computation of characters in a range. 2001-05-29 Bruno Haible <haible@clisp.cons.org> * iconvdata/gb18030.c (__fourbyte_to_ucs1): Add mappings for <U03F4>, <U03F5>. (__ucs_to_gb18030_tab1): Likewise. (BODY for FROM_LOOP): Add mapping for <U00010000>..<U0010FFFF>. (BODY for TO_LOOP): Likewise. * iconvdata/tst-table-charmap.sh: Update for charmaps containing <U00xxxxxx> syntax. * iconvdata/tst-table-from.c (bmp_only): New variable. (utf8_decode): If bmp_only, don't return characters outside Unicode plane 0. (main): When testing UTF-8 or GB18030, set bmp_only to 1. Don't print a conversion line if utf8_decode returns NULL. * iconvdata/tst-table-to.c (main): When testing encodings other than UTF-8 and GB18030, loop upto U+30000 instead of U+10000. Use UTF-8 instead of UCS-2 as input. * iconvdata/tst-table.sh: For GB18030, use only the part < 0x10000 of the charmap. 2001-05-29 Bruno Haible <haible@clisp.cons.org> * iconvdata/cns11643l1.c: Update to Unicode 3.1. (__cns11643l1_to_ucs4_tab): Regenerated. (__cns11643l1_from_ucs4_tab12): Regenerated. * iconvdata/cns11643.c: Update to Unicode 3.1. (__cns11643l14_to_ucs4_tab): Remove array. (__cns11643l3_to_ucs4_tab, __cns11643l4_to_ucs4_tab, __cns11643l5_to_ucs4_tab, __cns11643l6_to_ucs4_tab, __cns11643l7_to_ucs4_tab, __cns11643l15_to_ucs4_tab): New arrays. (__cns11643_from_ucs4p0_tab): Renamed from __cns11643_from_ucs4_tab. (__cns11643_from_ucs4p2_tab): New array. * iconvdata/cns11643.h (__cns11643l14_to_ucs4_tab): Remove declaration. (__cns11643l3_to_ucs4_tab, __cns11643l4_to_ucs4_tab, __cns11643l5_to_ucs4_tab, __cns11643l6_to_ucs4_tab, __cns11643l7_to_ucs4_tab, __cns11643l15_to_ucs4_tab): New declarations. (cns11643_to_ucs4): Treat planes 3, 4, 5, 6, 7, 15 instead of 14. (__cns11643_from_ucs4_tab): Remove declaration. (__cns11643_from_ucs4p0_tab, __cns11643_from_ucs4p2_tab): New declarations. (ucs4_to_cns11643): Update for new arrays. Treat U+3400..U+4DFF and U+20000..U+2A6D6. * iconvdata/cns11643l2.h (__cns11643_from_ucs4_tab): Remove declaration. (__cns11643_from_ucs4p0_tab): New declaration. (ucs4_to_cns11643l2): Update for new arrays. * iconvdata/iso-2022-cn-ext.c (BODY for FROM_LOOP): Handle planes 3 to 7. (BODY for TO_LOOP): Handle planes 3 to 7, instead of plane 14. * iconvdata/EUC-TW.irreversible: New file. * iconvdata/tst-table.sh: Use it. * iconvdata/Makefile (distribute): Add CP1255.irreversible, CP1258.irreversible, EUC-TW.irreversible. 2001-05-29 Bruno Haible <haible@clisp.cons.org> * locale/C-translit.h.in: Add transliterations for new Unicode 3.1 mathematical symbols.
2001-06-06 14:55:46 +02:00
for (i = LDBL_MIN_EXP, x = LDBL_MIN; i >= LDBL_MIN_EXP - LDBL_MANT_DIG + 1;
--i, x /= 2.0L)
{
printf ("2^%d: ", i);
r = frexpl (x, &e);
if (r != 0.5L)
{
printf ("mantissa incorrect: %.20La\n", r);
result = 1;
continue;
}
if (e != i)
{
printf ("exponent wrong %d (%.20Lg)\n", e, x);
result = 1;
continue;
}
puts ("ok");
}
}
# if 0
{
int e;
long double r = frexpl (LDBL_MIN * LDBL_EPSILON, &e);
if (r != 0.5)
{
printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): mantissa wrong: %Lg\n",
r);
result = 1;
}
else if (e != -16444)
{
printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): exponent wrong: %d\n",
e);
result = 1;
}
}
# endif
#endif
{
double x = 0x100000001ll + (double) 0.5;
double q;
double r;
r = modf (x, &q);
if (q != (double) 0x100000001ll || r != 0.5)
{
printf ("modf (%g, ...) failed\n", x);
result = 1;
}
}
{
union ieee754_float v1;
union ieee754_float v2;
float f;
v1.f = f = FLT_MIN;
if (fpclassify (f) != FP_NORMAL)
{
printf ("fpclassify (FLT_MIN) failed: %d\n", fpclassify (f));
result = 1;
}
f = nextafterf (f, FLT_MIN / 2.0f);
if (fpclassify (f) != FP_SUBNORMAL)
{
printf ("fpclassify (FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
result = 1;
}
v2.f = f = nextafterf (f, FLT_MIN);
if (fpclassify (f) != FP_NORMAL)
{
2000-12-19 06:48:11 +01:00
printf ("fpclassify (FLT_MIN-epsilon+epsilon) failed: %d\n",
fpclassify (f));
result = 1;
}
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("FLT_MIN: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("FLT_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("FLT_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.f = f = -FLT_MIN;
if (fpclassify (f) != FP_NORMAL)
{
printf ("fpclassify (-FLT_MIN) failed: %d\n", fpclassify (f));
result = 1;
}
f = nextafterf (f, -FLT_MIN / 2.0f);
if (fpclassify (f) != FP_SUBNORMAL)
{
printf ("fpclassify (-FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
result = 1;
}
v2.f = f = nextafterf (f, -FLT_MIN);
if (fpclassify (f) != FP_NORMAL)
{
2000-12-19 06:48:11 +01:00
printf ("fpclassify (-FLT_MIN-epsilon+epsilon) failed: %d\n",
fpclassify (f));
result = 1;
}
2000-12-19 06:48:11 +01:00
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("-FLT_MIN: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-FLT_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-FLT_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
2000-12-19 06:48:11 +01:00
f = FLT_MAX;
if (fpclassify (f) != FP_NORMAL)
{
printf ("fpclassify (FLT_MAX) failed: %d\n", fpclassify (f));
result = 1;
}
f = nextafterf (f, INFINITY);
if (fpclassify (f) != FP_INFINITE)
{
printf ("fpclassify (FLT_MAX+epsilon) failed: %d\n", fpclassify (f));
result = 1;
}
f = -FLT_MAX;
if (fpclassify (f) != FP_NORMAL)
{
printf ("fpclassify (-FLT_MAX) failed: %d\n", fpclassify (f));
result = 1;
}
f = nextafterf (f, -INFINITY);
if (fpclassify (f) != FP_INFINITE)
{
printf ("fpclassify (-FLT_MAX-epsilon) failed: %d\n", fpclassify (f));
result = 1;
}
2000-12-21 09:09:22 +01:00
v1.f = f = 0.0625;
f = nextafterf (f, 0.0);
v2.f = f = nextafterf (f, 1.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("0.0625f down: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625f down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625f down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.f = f = 0.0625;
f = nextafterf (f, 1.0);
v2.f = f = nextafterf (f, 0.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("0.0625f up: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625f up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625f up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.f = f = -0.0625;
f = nextafterf (f, 0.0);
v2.f = f = nextafterf (f, -1.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("-0.0625f up: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625f up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625f up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.f = f = -0.0625;
f = nextafterf (f, -1.0);
v2.f = f = nextafterf (f, 0.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("-0.0625f down: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625f down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625f down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.f = f = 0.0f;
f = nextafterf (f, 1.0);
v2.f = nextafterf (f, -1.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("0.0f up: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0f up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (0 != v2.ieee.negative)
{
printf ("0.0f up: negative differs: 0 vs %d\n",
v2.ieee.negative);
result = 1;
}
v1.f = f = 0.0f;
f = nextafterf (f, -1.0);
v2.f = nextafterf (f, 1.0);
if (v1.ieee.mantissa != v2.ieee.mantissa)
{
printf ("0.0f down: mantissa differs: %8x vs %8x\n",
v1.ieee.mantissa, v2.ieee.mantissa);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0f down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (1 != v2.ieee.negative)
{
printf ("0.0f down: negative differs: 1 vs %d\n",
v2.ieee.negative);
result = 1;
}
if (nextafterf (0.0f, INFINITY) != nextafterf (0.0f, 1.0f)
|| nextafterf (-0.0f, INFINITY) != nextafterf (-0.0f, 1.0f)
|| nextafterf (0.0f, -INFINITY) != nextafterf (0.0f, -1.0f)
|| nextafterf (-0.0f, -INFINITY) != nextafterf (-0.0f, -1.0f))
{
printf ("nextafterf (+-0, +-Inf) != nextafterf (+-0, +-1)\n");
result = 1;
}
if (nexttowardf (0.0f, INFINITY) != nexttowardf (0.0f, 1.0f)
|| nexttowardf (-0.0f, INFINITY) != nexttowardf (-0.0f, 1.0f)
|| nexttowardf (0.0f, -INFINITY) != nexttowardf (0.0f, -1.0f)
|| nexttowardf (-0.0f, -INFINITY) != nexttowardf (-0.0f, -1.0f))
{
printf ("nexttowardf (+-0, +-Inf) != nexttowardf (+-0, +-1)\n");
result = 1;
}
}
2000-12-21 09:09:22 +01:00
{
union ieee754_double v1;
union ieee754_double v2;
double d;
v1.d = d = DBL_MIN;
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (DBL_MIN) failed: %d\n", fpclassify (d));
result = 1;
}
d = nextafter (d, DBL_MIN / 2.0);
if (fpclassify (d) != FP_SUBNORMAL)
{
printf ("fpclassify (DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
result = 1;
}
v2.d = d = nextafter (d, DBL_MIN);
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (DBL_MIN-epsilon+epsilon) failed: %d\n",
fpclassify (d));
result = 1;
}
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("DBL_MIN: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("DBL_MIN: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("DBL_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("DBL_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = d = -DBL_MIN;
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (-DBL_MIN) failed: %d\n", fpclassify (d));
result = 1;
}
d = nextafter (d, -DBL_MIN / 2.0);
if (fpclassify (d) != FP_SUBNORMAL)
{
printf ("fpclassify (-DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
result = 1;
}
v2.d = d = nextafter (d, -DBL_MIN);
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (-DBL_MIN-epsilon+epsilon) failed: %d\n",
fpclassify (d));
result = 1;
}
2000-12-19 06:48:11 +01:00
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-DBL_MIN: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-DBL_MIN: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-DBL_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-DBL_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
2000-12-19 06:48:11 +01:00
d = DBL_MAX;
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (DBL_MAX) failed: %d\n", fpclassify (d));
result = 1;
}
d = nextafter (d, INFINITY);
if (fpclassify (d) != FP_INFINITE)
{
printf ("fpclassify (DBL_MAX+epsilon) failed: %d\n", fpclassify (d));
result = 1;
}
d = -DBL_MAX;
if (fpclassify (d) != FP_NORMAL)
{
printf ("fpclassify (-DBL_MAX) failed: %d\n", fpclassify (d));
result = 1;
}
d = nextafter (d, -INFINITY);
if (fpclassify (d) != FP_INFINITE)
{
printf ("fpclassify (-DBL_MAX-epsilon) failed: %d\n", fpclassify (d));
result = 1;
}
2000-12-21 09:09:22 +01:00
v1.d = d = 0.0625;
d = nextafter (d, 0.0);
v2.d = d = nextafter (d, 1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0625 down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0625 down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625 down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625 down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = d = 0.0625;
d = nextafter (d, 1.0);
v2.d = d = nextafter (d, 0.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0625 up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0625 up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625 up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625 up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = d = -0.0625;
d = nextafter (d, 0.0);
v2.d = d = nextafter (d, -1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-0.0625 up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-0.0625 up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625 up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625 up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = d = -0.0625;
d = nextafter (d, -1.0);
v2.d = d = nextafter (d, 0.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-0.0625 down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-0.0625 down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625 down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625 down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = d = 0.0;
d = nextafter (d, 1.0);
v2.d = nextafter (d, -1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0 up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0 up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0 up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (0 != v2.ieee.negative)
{
printf ("0.0 up: negative differs: 0 vs %d\n",
v2.ieee.negative);
result = 1;
}
v1.d = d = 0.0;
d = nextafter (d, -1.0);
v2.d = nextafter (d, 1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0 down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0 down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0 down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (1 != v2.ieee.negative)
{
printf ("0.0 down: negative differs: 1 vs %d\n",
v2.ieee.negative);
result = 1;
}
if (nextafter (0.0, INFINITY) != nextafter (0.0, 1.0)
|| nextafter (-0.0, INFINITY) != nextafter (-0.0, 1.0)
|| nextafter (0.0, -INFINITY) != nextafter (0.0, -1.0)
|| nextafter (-0.0, -INFINITY) != nextafter (-0.0, -1.0))
{
printf ("nextafter (+-0, +-Inf) != nextafter (+-0, +-1)\n");
result = 1;
}
if (nexttoward (0.0, INFINITY) != nexttoward (0.0, 1.0)
|| nexttoward (-0.0, INFINITY) != nexttoward (-0.0, 1.0)
|| nexttoward (0.0, -INFINITY) != nexttoward (0.0, -1.0)
|| nexttoward (-0.0, -INFINITY) != nexttoward (-0.0, -1.0))
{
printf ("nexttoward (+-0, +-Inf) != nexttoward (+-0, +-1)\n");
result = 1;
}
}
2000-12-21 09:09:22 +01:00
#ifndef NO_LONG_DOUBLE
{
union ieee854_long_double v1;
union ieee854_long_double v2;
long double ld;
v1.d = ld = LDBL_MIN;
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (LDBL_MIN) failed: %d\n", fpclassify (ld));
result = 1;
}
ld = nextafterl (ld, LDBL_MIN / 2.0);
if (fpclassify (ld) != FP_SUBNORMAL)
{
printf ("fpclassify (LDBL_MIN-epsilon) failed: %d (%La)\n",
fpclassify (ld), ld);
result = 1;
}
v2.d = ld = nextafterl (ld, LDBL_MIN);
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
fpclassify (ld), ld);
result = 1;
}
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("LDBL_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("LDBL_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = ld = -LDBL_MIN;
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (-LDBL_MIN) failed: %d\n", fpclassify (ld));
result = 1;
}
ld = nextafterl (ld, -LDBL_MIN / 2.0);
if (fpclassify (ld) != FP_SUBNORMAL)
{
printf ("fpclassify (-LDBL_MIN-epsilon) failed: %d (%La)\n",
fpclassify (ld), ld);
result = 1;
}
v2.d = ld = nextafterl (ld, -LDBL_MIN);
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (-LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
fpclassify (ld), ld);
result = 1;
}
2000-12-19 06:48:11 +01:00
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-LDBL_MIN: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-LDBL_MIN: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
2000-12-19 06:48:11 +01:00
ld = LDBL_MAX;
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (LDBL_MAX) failed: %d\n", fpclassify (ld));
result = 1;
}
ld = nextafterl (ld, INFINITY);
if (fpclassify (ld) != FP_INFINITE)
{
printf ("fpclassify (LDBL_MAX+epsilon) failed: %d\n", fpclassify (ld));
result = 1;
}
ld = -LDBL_MAX;
if (fpclassify (ld) != FP_NORMAL)
{
printf ("fpclassify (-LDBL_MAX) failed: %d\n", fpclassify (ld));
result = 1;
}
ld = nextafterl (ld, -INFINITY);
if (fpclassify (ld) != FP_INFINITE)
{
printf ("fpclassify (-LDBL_MAX-epsilon) failed: %d\n",
fpclassify (ld));
result = 1;
}
2000-12-21 09:09:22 +01:00
v1.d = ld = 0.0625;
ld = nextafterl (ld, 0.0);
v2.d = ld = nextafterl (ld, 1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0625L down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0625L down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625L down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625L down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = ld = 0.0625;
ld = nextafterl (ld, 1.0);
v2.d = ld = nextafterl (ld, 0.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0625L up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0625L up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0625L up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("0.0625L up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = ld = -0.0625;
ld = nextafterl (ld, 0.0);
v2.d = ld = nextafterl (ld, -1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-0.0625L up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-0.0625L up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625L up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625L up: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = ld = -0.0625;
ld = nextafterl (ld, -1.0);
v2.d = ld = nextafterl (ld, 0.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("-0.0625L down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("-0.0625L down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("-0.0625L down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (v1.ieee.negative != v2.ieee.negative)
{
printf ("-0.0625L down: negative differs: %d vs %d\n",
v1.ieee.negative, v2.ieee.negative);
result = 1;
}
v1.d = ld = 0.0;
ld = nextafterl (ld, 1.0);
v2.d = nextafterl (ld, -1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0L up: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0L up: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0L up: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (0 != v2.ieee.negative)
{
printf ("0.0L up: negative differs: 0 vs %d\n",
v2.ieee.negative);
result = 1;
}
v1.d = ld = 0.0;
ld = nextafterl (ld, -1.0);
v2.d = nextafterl (ld, 1.0);
if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
{
printf ("0.0L down: mantissa0 differs: %8x vs %8x\n",
v1.ieee.mantissa0, v2.ieee.mantissa0);
result = 1;
}
if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
{
printf ("0.0L down: mantissa1 differs: %8x vs %8x\n",
v1.ieee.mantissa1, v2.ieee.mantissa1);
result = 1;
}
if (v1.ieee.exponent != v2.ieee.exponent)
{
printf ("0.0L down: exponent differs: %4x vs %4x\n",
v1.ieee.exponent, v2.ieee.exponent);
result = 1;
}
if (1 != v2.ieee.negative)
{
printf ("0.0L down: negative differs: 1 vs %d\n",
v2.ieee.negative);
result = 1;
}
if (nextafterl (0.0, INFINITY) != nextafterl (0.0, 1.0)
|| nextafterl (-0.0, INFINITY) != nextafterl (-0.0, 1.0)
|| nextafterl (0.0, -INFINITY) != nextafterl (0.0, -1.0)
|| nextafterl (-0.0, -INFINITY) != nextafterl (-0.0, -1.0))
{
printf ("nextafterl (+-0, +-Inf) != nextafterl (+-0, +-1)\n");
result = 1;
}
if (nexttowardl (0.0L, INFINITY) != nexttowardl (0.0L, 1.0L)
|| nexttowardl (-0.0L, INFINITY) != nexttowardl (-0.0L, 1.0L)
|| nexttowardl (0.0L, -INFINITY) != nexttowardl (0.0L, -1.0L)
|| nexttowardl (-0.0L, -INFINITY) != nexttowardl (-0.0L, -1.0L))
{
printf ("nexttowardl (+-0, +-Inf) != nexttowardl (+-0, +-1)\n");
result = 1;
}
}
#endif
if (! isnormal (FLT_MIN))
{
puts ("isnormal (FLT_MIN) failed");
result = 1;
}
if (! isnormal (DBL_MIN))
{
puts ("isnormal (DBL_MIN) failed");
result = 1;
}
#ifndef NO_LONG_DOUBLE
if (! isnormal (LDBL_MIN))
{
puts ("isnormal (LDBL_MIN) failed");
result = 1;
}
#endif
#ifdef __i386__
/* This is a test for the strange long doubles in x86 FPUs. */
{
union
{
char b[10];
long double d;
} u =
{ .b = { 0, 0, 0, 0, 0, 0, 0, 0x80, 0, 0 } };
if (fpclassify (u.d) != FP_NORMAL)
{
printf ("fpclassify (0x00008000000000000000) failed: %d (%Lg)\n",
fpclassify (u.d), u.d);
result = 1;
}
}
/* Special NaNs in x86 long double. Test for scalbl. */
{
union
{
char b[10];
long double d;
} u =
{ .b = { 0, 1, 0, 0, 0, 0, 0, 0xc0, 0xff, 0x7f } };
long double r;
r = scalbl (u.d, 0.0);
if (!isnan (r))
{
puts ("scalbl(NaN, 0) does not return NaN");
result = 1;
}
else if (memcmp (&r, &u.d, sizeof (double)) != 0)
{
puts ("scalbl(NaN, 0) does not return the same NaN");
result = 1;
}
}
#endif
#ifndef NO_LONG_DOUBLE
{
long double r;
feclearexcept (FE_ALL_EXCEPT);
r = scalbl (LDBL_MIN, 2147483647);
if (! isinf (r))
{
puts ("scalbl (LDBL_MIN, 2147483647) does not return Inf");
result = 1;
}
else if (signbit (r) != 0)
{
puts ("scalbl (LDBL_MIN, 2147483647) returns -Inf");
result = 1;
}
else if (fetestexcept (FE_UNDERFLOW))
{
puts ("scalbl(NaN, 0) raises underflow exception");
result = 1;
}
feclearexcept (FE_ALL_EXCEPT);
r = scalbl (LDBL_MAX, -2147483647);
if (r != 0.0)
{
puts ("scalbl (LDBL_MAX, -2147483647) does not return 0");
result = 1;
}
else if (signbit (r) != 0)
{
puts ("scalbl (LDBL_MAX, -2147483647) returns -Inf");
result = 1;
}
else if (fetestexcept (FE_OVERFLOW))
{
puts ("scalbl(NaN, 0) raises overflow exception");
result = 1;
}
}
#endif
[BZ #779] 2005-03-10 Jakub Jelinek <jakub@redhat.com> * math/test-misc.c (main): Add some more tests. 2005-03-17 Jakub Jelinek <jakub@redhat.com> * posix/regcomp.c (re_compile_fastmap_iter): Fix check for failed __wcrtomb. Check return values of other __wcrtomb calls. * posix/regex_internal.c (build_wcs_buffer, re_string_skip_chars): Change mbclen type to size_t. (build_wcs_upper_buffer): Change mbclen and mbcdlen type to size_t. Handle mb chars whose upper case doesn't have multibyte representation in locale's charset. 2005-03-15 Jakub Jelinek <jakub@redhat.com> * malloc/malloc.c (_int_icalloc, _int_icomalloc, iALLOc, public_iCALLOc, public_iCALLOc, public_iCOMALLOc): Protect with #ifndef _LIBC. [BZ #779] * malloc/malloc.c (public_mTRIm): Initialize malloc if not yet initialized. 2005-03-10 Jakub Jelinek <jakub@redhat.com> * misc/sys/cdefs.h (__always_inline): Define. * posix/bits/unistd.h (read, pread, pread64, readlink, getcwd, getwd): Use __always_inline instead of __inline. * socket/bits/socket2.h (recv, recvfrom): Likewise. * libio/bits/stdio2.h (gets, fgets, fgets_unlocked): Likewise. * string/bits/string3.h (__memcpy_ichk, __memmove_ichk, __mempcpy_ichk, __memset_ichk, __strcpy_ichk, __stpcpy_ichk, __strncpy_ichk, __strcat_ichk, __strncat_ichk): Use __always_inline instead of __inline__ __attribute__ ((__always_inline__)). 2005-03-09 Jakub Jelinek <jakub@redhat.com> * debug/tst-chk1.c: Include sys/socket.h and sys/un.h. (do_test): Add new tests for recv, recvfrom, getcwd, getwd and readlink. Add some more tests for read, pread, pread64, fgets and fgets_unlocked. * posix/bits/unistd.h (read, pread, pread64, readlink, getcwd, getwd): Change macros into extern inline functions. (__read_alias, __pread_alias, __pread64_alias, __readlink_alias, __getcwd_alias, __getwd_alias): New prototypes. * socket/bits/socket2.h (recv, recvfrom): Change macros into extern inline functions. (__recv_alias, __recvfrom_alias): New prototypes. * libio/bits/stdio2.h (gets, fgets, fgets_unlocked): Change macros into extern inline functions. (__gets_alias, __fgets_alias, __fgets_unlocked_alias): New prototypes. * debug/pread_chk.c (__pread_chk): Fix order of arguments passed to __pread. * debug/pread64_chk.c (__pread64_chk): Fix order of arguments passed to __pread64.
2005-03-19 01:28:51 +01:00
/* The tests here are very similar to tests earlier in this file,
the important difference is just that there are no intervening
union variables that cause some GCC versions to hide possible
bugs in nextafter* implementation. */
if (nextafterf (nextafterf (FLT_MIN, FLT_MIN / 2.0), FLT_MIN) != FLT_MIN)
{
puts ("nextafterf FLT_MIN test failed");
result = 1;
}
if (nextafterf (nextafterf (-FLT_MIN, -FLT_MIN / 2.0), -FLT_MIN)
!= -FLT_MIN)
{
puts ("nextafterf -FLT_MIN test failed");
result = 1;
}
if (nextafter (nextafter (DBL_MIN, DBL_MIN / 2.0), DBL_MIN) != DBL_MIN)
{
puts ("nextafter DBL_MIN test failed");
result = 1;
}
if (nextafter (nextafter (-DBL_MIN, -DBL_MIN / 2.0), -DBL_MIN) != -DBL_MIN)
{
puts ("nextafter -DBL_MIN test failed");
result = 1;
}
#ifndef NO_LONG_DOUBLE
if (nextafterl (nextafterl (LDBL_MIN, LDBL_MIN / 2.0), LDBL_MIN)
!= LDBL_MIN)
{
puts ("nextafterl LDBL_MIN test failed");
result = 1;
}
if (nextafterl (nextafterl (-LDBL_MIN, -LDBL_MIN / 2.0), -LDBL_MIN)
!= -LDBL_MIN)
{
puts ("nextafterl -LDBL_MIN test failed");
result = 1;
}
#endif
volatile float f1 = FLT_MAX;
volatile float f2 = FLT_MAX / 2;
(void) &f1;
(void) &f2;
feclearexcept (FE_ALL_EXCEPT);
f2 += f1;
int fe = fetestexcept (FE_ALL_EXCEPT);
if (fe != (FE_OVERFLOW | FE_INEXACT))
{
printf ("float overflow test failed: %x\n", fe);
result = 1;
}
volatile double d1 = DBL_MAX;
volatile double d2 = DBL_MAX / 2;
(void) &d1;
(void) &d2;
feclearexcept (FE_ALL_EXCEPT);
d2 += d1;
fe = fetestexcept (FE_ALL_EXCEPT);
if (fe != (FE_OVERFLOW | FE_INEXACT))
{
printf ("double overflow test failed: %x\n", fe);
result = 1;
}
#ifndef NO_LONG_DOUBLE
volatile long double ld1 = LDBL_MAX;
volatile long double ld2 = LDBL_MAX / 2;
(void) &ld1;
(void) &ld2;
feclearexcept (FE_ALL_EXCEPT);
ld2 += ld1;
fe = fetestexcept (FE_ALL_EXCEPT);
if (fe != (FE_OVERFLOW | FE_INEXACT))
{
printf ("long double overflow test failed: %x\n", fe);
result = 1;
}
#endif
#if !defined NO_LONG_DOUBLE && LDBL_MANT_DIG == 113
volatile long double ld3 = 0x1.0000000000010000000100000001p+1;
volatile long double ld4 = 0x1.0000000000000000000000000001p+1;
(void) &ld3;
(void) &ld4;
ld3 -= ld4;
if (ld3 != 0x1.0p-47)
{
printf ("long double subtraction test failed %.28La\n", ld3);
result = 1;
}
#endif
#if !defined NO_LONG_DOUBLE && LDBL_MANT_DIG >= DBL_MANT_DIG + 4
* soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. * soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. 2007-04-14 Uros Bizjak <ubizjak@gmail.com> * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. 2007-04-16 Uros Bizjak <ubizjak@gmail.com> Jakub Jelinek <jakub@redhat.com> * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. 2007-05-03 Jakub Jelinek <jakub@redhat.com>
2007-05-03 18:36:49 +02:00
int oldmode = fegetround ();
int j;
for (j = 0; j < 4; j++)
{
* soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. * soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. 2007-04-14 Uros Bizjak <ubizjak@gmail.com> * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. 2007-04-16 Uros Bizjak <ubizjak@gmail.com> Jakub Jelinek <jakub@redhat.com> * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. 2007-05-03 Jakub Jelinek <jakub@redhat.com>
2007-05-03 18:36:49 +02:00
int mode;
int i;
int k = 0;
const char *mstr;
switch (j)
{
* soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. * soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. 2007-04-14 Uros Bizjak <ubizjak@gmail.com> * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. 2007-04-16 Uros Bizjak <ubizjak@gmail.com> Jakub Jelinek <jakub@redhat.com> * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. 2007-05-03 Jakub Jelinek <jakub@redhat.com>
2007-05-03 18:36:49 +02:00
#ifdef FE_TONEAREST
case 0:
mode = FE_TONEAREST;
mstr = "nearest";
k = 8;
break;
#endif
#ifdef FE_DOWNWARD
case 1:
mode = FE_DOWNWARD;
mstr = "-inf";
break;
#endif
#ifdef FE_UPWARD
case 2:
mode = FE_UPWARD;
mstr = "+inf";
k = 15;
break;
#endif
#ifdef FE_TOWARDZERO
case 3:
mode = FE_TOWARDZERO;
mstr = "0";
break;
#endif
default:
continue;
}
volatile long double ld5 = nextafter (0.0, 1.0) / 16.0L;
volatile double d5;
(void) &ld5;
for (i = 0; i <= 32; i++)
{
if (fesetround (mode))
{
printf ("failed to set rounding mode to %s\n", mstr);
result = 1;
break;
}
d5 = ld5 * i;
(void) &d5;
fesetround (oldmode);
if (d5 != ((j == 0 && i == 8) ? 0 : (i + k) / 16)
* nextafter (0.0, 1.0))
{
printf ("%La incorrectly rounded to %s as %a\n",
ld5 * i, mstr, d5);
result = 1;
}
}
* soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. * soft-fp/op-common.h (FP_TRUNC): Replace raising of FP_EX_INEXACT with setting the sticky bit. * math/test-misc.c (main): Add more truncation tests. 2007-04-14 Uros Bizjak <ubizjak@gmail.com> * soft-fp/floatunsidf.c (__floatunsidf): Use DFtype instead of double in the function declaration. * soft-fp/floatundidf.c (__floatundidf): Use DFtype instead of double in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/floatunsisf.c (__floatunsisf): Use SFtype instead of float in the function declaration. * soft-fp/extended.h (FP_UNPACK_RAW_E): Do not increase X##_e for denormal operands. Do not generate FP_EX_DENORM exception. (FP_UNPACK_RAW_EP): Ditto. (FP_UNPACK_SEMIRAW_E): Use FP_UNPACK_RAW_E instead of undefined _FP_UNPACK_RAW_E. (FP_UNPACK_SEMIRAW_EP): Use FP_UNPACK_RAW_EP instead of undefined _FP_UNPACK_RAW_EP. (FP_PACK_SEMIRAW_E): Use FP_PACK_RAW_E instead of undefined _FP_PACK_RAW_E. (FP_PACK_SEMIRAW_EP): Use FP_PACK_RAW_EP instead of undefined _FP_PACK_RAW_EP. * op-2.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_2. * op-4.h (_FP_FRAC_COPY_2_2): Define as alias to _FP_FRAC_COPY_4. 2007-04-16 Uros Bizjak <ubizjak@gmail.com> Jakub Jelinek <jakub@redhat.com> * soft-fp/op-common.h (FP_EXTEND): Do not abort when _FP_EXPBIAS_##dfs == _FP_EXPBIAS_##sfs. Handle denormals for this case. * soft-fp/op-common.h (FP_TRUNC): Ditto. 2007-05-03 Jakub Jelinek <jakub@redhat.com>
2007-05-03 18:36:49 +02:00
}
volatile long double ld7 = nextafterl (0.0L, 1.0L);
volatile double d7;
(void) &ld7;
fesetround (FE_UPWARD);
d7 = ld7;
(void) &d7;
fesetround (oldmode);
if (d7 != nextafter (0.0, 1.0))
{
printf ("%La incorrectly rounded upward to %a\n", ld7, d7);
result = 1;
}
#endif
return result;
}