glibc/math/test-misc.c

1144 lines
28 KiB
C
Raw Normal View History

/* Miscellaneous tests which don't fit anywhere else.
2001-01-27 19:09:14 +00:00
Copyright (C) 2000, 2001 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 19:09:14 +00: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;
}
}
# if __GNUC__ >= 3 || __GNUC_MINOR__ >= 96
{
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 == 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 12:55:46 +00: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");
}
}
# endif
#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 05:48:11 +00: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 05:48:11 +00:00
printf ("fpclassify (-FLT_MIN-epsilon+epsilon) failed: %d\n",
fpclassify (f));
result = 1;
}
2000-12-19 05:48:11 +00: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 05:48:11 +00: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 08:09:22 +00: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 08:09:22 +00: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 05:48:11 +00: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 05:48:11 +00: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 08:09:22 +00: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 08:09:22 +00: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 05:48:11 +00: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 05:48:11 +00: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 08:09:22 +00: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
return result;
}