lib/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-07-29 11:41:21 +03:00
Code Activity

Update.

Browse Source 
* math/test-misc.c (main): More tests for frexp and some for
	fpclassify and isnormal.
This commit is contained in:
Ulrich Drepper
2000-12-03 23:38:18 +00:00
parent be2925fc8a
commit 0df757b10a
2 changed files with 63 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
ChangeLog
View File

@ -1,5 +1,8 @@
2000-12-03 Ulrich Drepper <drepper@redhat.com> 2000-12-03 Ulrich Drepper <drepper@redhat.com>
* math/test-misc.c (main): More tests for frexp and some for
fpclassify and isnormal.
* sysdeps/ieee754/flt-32/s_fpclassifyf.c (__fpclassifyf): Correct * sysdeps/ieee754/flt-32/s_fpclassifyf.c (__fpclassifyf): Correct
test for subnormal. test for subnormal.
Reported by Fred J. Tydeman <tydeman@tybor.com>. Reported by Fred J. Tydeman <tydeman@tybor.com>.

67
math/test-misc.c
View File

@ -42,8 +42,10 @@ main (void)
# if __GNUC__ >= 3 || __GNUC_MINOR__ >= 96 # if __GNUC__ >= 3 || __GNUC_MINOR__ >= 96
{ {
long double x = LDBL_MAX / ldexpl (1.0L, LDBL_MANT_DIG + 1); long double x;
long double m; long double m;
long double r;
int e;
int i; int i;
# if LDBL_MANT_DIG == 64 # if LDBL_MANT_DIG == 64
@ -52,12 +54,9 @@ main (void)
# error "Please adjust" # error "Please adjust"
# endif # endif
for (i = 0; i < LDBL_MANT_DIG + 1; ++i, x *= 2.0L) for (i = LDBL_MAX_EXP, x = LDBL_MAX; i >= LDBL_MIN_EXP; --i, x /= 2.0L)
{ {
long double r; printf ("2^%d: ", i);
int e;
printf ("2^%d: ", LDBL_MAX_EXP - (LDBL_MANT_DIG + 1) + i);
r = frexpl (x, &e); r = frexpl (x, &e);
if (r != m) if (r != m)
@ -66,7 +65,7 @@ main (void)
result = 1; result = 1;
continue; continue;
} }
if (e != LDBL_MAX_EXP - (LDBL_MANT_DIG + 1) + i) if (e != i)
{ {
printf ("exponent wrong %d (%.20Lg)\n", e, x); printf ("exponent wrong %d (%.20Lg)\n", e, x);
result = 1; result = 1;
@ -91,5 +90,59 @@ main (void)
} }
} }
if (fpclassify (FLT_MIN) != FP_NORMAL)
{
printf ("fpclassify (FLT_MIN) failed: %d\n", fpclassify (FLT_MIN));
result = 1;
}
if (fpclassify (nextafterf (FLT_MIN, FLT_MIN / 2.0f)) != FP_SUBNORMAL)
{
printf ("fpclassify (FLT_MIN-epsilon) failed: %d\n",
fpclassify (nextafterf (FLT_MIN, FLT_MIN / 2.0f)));
result = 1;
}
if (fpclassify (DBL_MIN) != FP_NORMAL)
{
printf ("fpclassify (DBL_MIN) failed: %d\n", fpclassify (DBL_MIN));
result = 1;
}
if (fpclassify (nextafter (DBL_MIN, DBL_MIN / 2.0)) != FP_SUBNORMAL)
{
printf ("fpclassify (DBL_MIN-epsilon) failed: %d\n",
fpclassify (nextafter (DBL_MIN, DBL_MIN / 2.0)));
result = 1;
}
#ifndef NO_LONG_DOUBLE
if (fpclassify (LDBL_MIN) != FP_NORMAL)
{
printf ("fpclassify (LDBL_MIN) failed: %d\n", fpclassify (LDBL_MIN));
result = 1;
}
if (fpclassify (nextafterl (LDBL_MIN, LDBL_MIN / 2.0)) != FP_SUBNORMAL)
{
printf ("fpclassify (LDBL_MIN-epsilon) failed: %d\n",
fpclassify (nextafterl (LDBL_MIN, LDBL_MIN / 2.0)));
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
return result; return result;
} }