2000-01-20 09:27:38 +01:00
|
|
|
/* Miscellaneous tests which don't fit anywhere else.
|
2007-01-16 00:43:04 +01:00
|
|
|
Copyright (C) 2000, 2001, 2004, 2005, 2007 Free Software Foundation, Inc.
|
2000-01-20 09:27:38 +01:00
|
|
|
This file is part of the GNU C Library.
|
|
|
|
|
|
|
|
The GNU C Library is free software; you can redistribute it and/or
|
2001-07-06 06:58:11 +02:00
|
|
|
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.
|
2000-01-20 09:27:38 +01:00
|
|
|
|
|
|
|
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
|
2001-07-06 06:58:11 +02:00
|
|
|
Lesser General Public License for more details.
|
2000-01-20 09:27:38 +01:00
|
|
|
|
2001-07-06 06:58:11 +02:00
|
|
|
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. */
|
2000-01-20 09:27:38 +01:00
|
|
|
|
2000-12-05 07:49:26 +01:00
|
|
|
#include <fenv.h>
|
2001-01-27 20:09:14 +01:00
|
|
|
#include <float.h>
|
2000-12-21 05:17:49 +01:00
|
|
|
#include <ieee754.h>
|
2000-01-20 09:27:38 +01:00
|
|
|
#include <math.h>
|
|
|
|
#include <stdio.h>
|
2000-12-05 07:36:53 +01:00
|
|
|
#include <string.h>
|
2000-01-20 09:27:38 +01:00
|
|
|
|
|
|
|
|
|
|
|
int
|
|
|
|
main (void)
|
|
|
|
{
|
|
|
|
int result = 0;
|
|
|
|
|
2000-10-25 19:19:53 +02:00
|
|
|
#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;
|
|
|
|
}
|
|
|
|
}
|
2000-12-03 11:00:46 +01:00
|
|
|
|
|
|
|
{
|
2000-12-04 00:38:18 +01:00
|
|
|
long double x;
|
2000-12-03 11:00:46 +01:00
|
|
|
long double m;
|
2000-12-04 00:38:18 +01:00
|
|
|
long double r;
|
|
|
|
int e;
|
2000-12-03 11:00:46 +01:00
|
|
|
int i;
|
|
|
|
|
2007-01-16 00:43:04 +01:00
|
|
|
# if LDBL_MANT_DIG == 64
|
2000-12-03 11:00:46 +01:00
|
|
|
m = 0xf.fffffffffffffffp-4L;
|
2007-01-16 00:43:04 +01:00
|
|
|
# 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;
|
2007-01-16 00:43:04 +01:00
|
|
|
# elif LDBL_MANT_DIG == 113
|
2000-12-13 15:04:38 +01:00
|
|
|
m = 0x1.ffffffffffffffffffffffffffffp-1L;
|
2007-01-16 00:43:04 +01:00
|
|
|
# else
|
|
|
|
# error "Please adjust"
|
|
|
|
# endif
|
2000-12-03 11:00:46 +01:00
|
|
|
|
2000-12-04 00:38:18 +01:00
|
|
|
for (i = LDBL_MAX_EXP, x = LDBL_MAX; i >= LDBL_MIN_EXP; --i, x /= 2.0L)
|
2000-12-03 11:00:46 +01:00
|
|
|
{
|
2000-12-04 00:38:18 +01:00
|
|
|
printf ("2^%d: ", i);
|
2000-12-03 11:00:46 +01:00
|
|
|
|
|
|
|
r = frexpl (x, &e);
|
|
|
|
if (r != m)
|
|
|
|
{
|
|
|
|
printf ("mantissa incorrect: %.20La\n", r);
|
|
|
|
result = 1;
|
|
|
|
continue;
|
|
|
|
}
|
2000-12-04 00:38:18 +01:00
|
|
|
if (e != i)
|
2000-12-03 11:00:46 +01:00
|
|
|
{
|
|
|
|
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");
|
|
|
|
}
|
|
|
|
|
2000-12-03 11:00:46 +01:00
|
|
|
}
|
2000-12-05 07:36:53 +01:00
|
|
|
|
2007-01-16 00:43:04 +01:00
|
|
|
# if 0
|
2000-12-05 07:36:53 +01:00
|
|
|
{
|
|
|
|
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;
|
|
|
|
}
|
|
|
|
}
|
2007-01-16 00:43:04 +01:00
|
|
|
# endif
|
2000-10-25 19:19:53 +02:00
|
|
|
#endif
|
2000-01-20 09:27:38 +01:00
|
|
|
|
|
|
|
{
|
|
|
|
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;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2000-12-19 06:07:31 +01:00
|
|
|
{
|
2000-12-21 05:17:49 +01:00
|
|
|
union ieee754_float v1;
|
|
|
|
union ieee754_float v2;
|
2000-12-19 06:07:31 +01:00
|
|
|
float f;
|
|
|
|
|
2000-12-21 05:17:49 +01:00
|
|
|
v1.f = f = FLT_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.f = f = nextafterf (f, FLT_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
if (fpclassify (f) != FP_NORMAL)
|
|
|
|
{
|
2000-12-19 06:48:11 +01:00
|
|
|
printf ("fpclassify (FLT_MIN-epsilon+epsilon) failed: %d\n",
|
2000-12-19 06:07:31 +01:00
|
|
|
fpclassify (f));
|
|
|
|
result = 1;
|
|
|
|
}
|
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
|
|
|
|
|
|
|
v1.f = f = -FLT_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.f = f = nextafterf (f, -FLT_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
if (fpclassify (f) != FP_NORMAL)
|
|
|
|
{
|
2000-12-19 06:48:11 +01:00
|
|
|
printf ("fpclassify (-FLT_MIN-epsilon+epsilon) failed: %d\n",
|
2000-12-19 06:07:31 +01:00
|
|
|
fpclassify (f));
|
|
|
|
result = 1;
|
|
|
|
}
|
2000-12-19 06:48:11 +01:00
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
2001-06-16 06:34:51 +02:00
|
|
|
|
|
|
|
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-19 06:07:31 +01:00
|
|
|
}
|
2000-12-21 09:09:22 +01:00
|
|
|
|
2000-12-19 06:07:31 +01:00
|
|
|
{
|
2000-12-21 05:17:49 +01:00
|
|
|
union ieee754_double v1;
|
|
|
|
union ieee754_double v2;
|
2000-12-19 06:07:31 +01:00
|
|
|
double d;
|
|
|
|
|
2000-12-21 05:17:49 +01:00
|
|
|
v1.d = d = DBL_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.d = d = nextafter (d, DBL_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
if (fpclassify (d) != FP_NORMAL)
|
|
|
|
{
|
|
|
|
printf ("fpclassify (DBL_MIN-epsilon+epsilon) failed: %d\n",
|
|
|
|
fpclassify (d));
|
|
|
|
result = 1;
|
|
|
|
}
|
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
|
|
|
|
|
|
|
v1.d = d = -DBL_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.d = d = nextafter (d, -DBL_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
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
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
2001-06-16 06:34:51 +02:00
|
|
|
|
|
|
|
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-19 06:07:31 +01:00
|
|
|
}
|
2000-12-21 09:09:22 +01:00
|
|
|
|
2000-12-04 00:38:18 +01:00
|
|
|
#ifndef NO_LONG_DOUBLE
|
2000-12-19 06:07:31 +01:00
|
|
|
{
|
2000-12-21 05:17:49 +01:00
|
|
|
union ieee854_long_double v1;
|
|
|
|
union ieee854_long_double v2;
|
2000-12-19 06:07:31 +01:00
|
|
|
long double ld;
|
|
|
|
|
2000-12-21 05:17:49 +01:00
|
|
|
v1.d = ld = LDBL_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.d = ld = nextafterl (ld, LDBL_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
if (fpclassify (ld) != FP_NORMAL)
|
|
|
|
{
|
|
|
|
printf ("fpclassify (LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
|
|
|
|
fpclassify (ld), ld);
|
|
|
|
result = 1;
|
|
|
|
}
|
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
|
|
|
|
|
|
|
v1.d = ld = -LDBL_MIN;
|
2000-12-19 06:07:31 +01:00
|
|
|
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;
|
|
|
|
}
|
2000-12-21 05:17:49 +01:00
|
|
|
v2.d = ld = nextafterl (ld, -LDBL_MIN);
|
2000-12-19 06:07:31 +01:00
|
|
|
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
|
|
|
|
2000-12-21 05:17:49 +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;
|
|
|
|
}
|
2001-06-16 06:34:51 +02:00
|
|
|
|
|
|
|
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;
|
|
|
|
}
|
2000-12-19 06:07:31 +01:00
|
|
|
}
|
2000-12-04 00:38:18 +01:00
|
|
|
#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
|
|
|
|
|
2000-12-04 02:27:16 +01:00
|
|
|
#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;
|
|
|
|
}
|
|
|
|
}
|
2000-12-05 07:36:53 +01:00
|
|
|
|
|
|
|
/* 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;
|
|
|
|
}
|
|
|
|
}
|
2000-12-04 02:27:16 +01:00
|
|
|
#endif
|
|
|
|
|
2000-12-05 07:49:26 +01:00
|
|
|
#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
|
|
|
|
|
2007-01-16 00:43:04 +01:00
|
|
|
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
|
|
|
|
|
2007-05-03 14:15:04 +02:00
|
|
|
#if !defined NO_LONG_DOUBLE && LDBL_MANT_DIG >= DBL_MANT_DIG + 4
|
2007-05-03 18:36:49 +02:00
|
|
|
int oldmode = fegetround ();
|
|
|
|
int j;
|
|
|
|
for (j = 0; j < 4; j++)
|
2007-05-03 14:15:04 +02:00
|
|
|
{
|
2007-05-03 18:36:49 +02:00
|
|
|
int mode;
|
|
|
|
int i;
|
|
|
|
int k = 0;
|
|
|
|
const char *mstr;
|
|
|
|
switch (j)
|
2007-05-03 14:15:04 +02:00
|
|
|
{
|
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;
|
|
|
|
}
|
2007-05-03 14:15:04 +02:00
|
|
|
}
|
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;
|
|
|
|
}
|
2007-05-03 14:15:04 +02:00
|
|
|
#endif
|
|
|
|
|
2000-01-20 09:27:38 +01:00
|
|
|
return result;
|
|
|
|
}
|