Mail Archives: djgpp-workers/2013/10/19/07:44:13
| X-Authentication-Warning: | delorie.com: mail set sender to djgpp-workers-bounces using -f
|
| X-Recipient: | djgpp-workers AT delorie DOT com
|
| Message-ID: | <52626FC1.9010706@gmx.de>
|
| Date: | Sat, 19 Oct 2013 13:40:49 +0200
|
| From: | Juan Manuel Guerrero <juan DOT guerrero AT gmx DOT de>
|
| User-Agent: | Mozilla/5.0 (X11; Linux x86_64; rv:16.0) Gecko/20121025 Thunderbird/16.0.2
|
| MIME-Version: | 1.0
|
| To: | djgpp-workers AT delorie DOT com
|
| Subject: | Implementation of the [l]lround[f|l] family of functions
|
| X-Provags-ID: | V03:K0:ayy+zvO+x5rM7D3uGEh62YJC4Td5h1rIGN96/4nnFbVp3KKaeCa
|
| rbRPXOewW9so7yy7SorILTk2DupEOxoghGBhEXZ6ygZsHYFENR4MYCxFUXd/wHNDnOleJyx
|
| zFznw6KI/XoQBrlWFKVtUNNDefHK4cyArLMl7psqAgA16wJpDGTxB39ZZzg3lXp90FI6DVf
|
| et+RCSWNnuORc69BvyBnw==
|
| Reply-To: | djgpp-workers AT delorie DOT com
|
Below is a patch that shall provide the implementation of the [l]lround[f|l]
family of functions for libm.a. Some test cases have been added. I have
compiled and tested the code with gcc473 and with djcross-gcc481.
As usual suggestions, objections and comments are welcome.
Regards,
Juan M. Guerrero
2013-10-19 Juan Manuel Guerrero <juan DOT guerrero AT gmx DOT de>
* djgpp/src/docs/kb/wc204.txi: Info about [l]lround[f|l] family of functions added.
2013-10-06 Juan Manuel Guerrero <juan DOT guerrero AT gmx DOT de>
* djgpp/src/libm/math/lroundl.c: Implementation of lroundl.
* djgpp/src/libm/math/llroundl.c: Implementation of llroundl.
* djgpp/src/libm/math/makefile: [l]lroundl added to goal list.
* djgpp/tests/cygnus/t-llroundl.c: Check for llroundl.
* djgpp/tests/cygnus/t-lroundl.c: Check for lroundl.
* djgpp/tests/cygnus/makefile: [l]lroundl added to goal list.
2013-10-03 Juan Manuel Guerrero <juan DOT guerrero AT gmx DOT de>
* djgpp/src/libm/math/lround.c: Implementation of lround.
* djgpp/src/libm/math/llround.c: Implementation of llround.
* djgpp/src/libm/math/makefile: [l]lround added to goal list.
* djgpp/tests/cygnus/t-llround.c: Check for llround.
* djgpp/tests/cygnus/t-lround.c: Check for lround.
* djgpp/tests/cygnus/makefile: [l]lround added to goal list.
2013-10-02 Juan Manuel Guerrero <juan DOT guerrero AT gmx DOT de>
* djgpp/include/math.h: Prototype of [l]lround[f|l] added.
* djgpp/include/libm/math.h: Prototype of [l]lround[f|l] added.
* djgpp/src/libm/math/lroundf.c: Implementation of lroundf.
* djgpp/src/libm/math/llroundf.c: Implementation of llroundf.
* djgpp/src/libm/math/makefile: [l]lroundf added to goal list.
* djgpp/tests/cygnus/t-lroundf.c: Check for lroundf.
* djgpp/tests/cygnus/makefile: lroundf added to goal list.
diff -aprNU5 djgpp.orig/include/libm/math.h djgpp/include/libm/math.h
--- djgpp.orig/include/libm/math.h 2013-10-14 22:37:58 +0100
+++ djgpp/include/libm/math.h 2013-10-14 22:40:02 +0100
@@ -162,10 +162,16 @@ extern long long int llrintf __P((float)
extern long long int llrint __P((double));
extern long long int llrintl __P((long double));
extern long int lrintf __P((float));
extern long int lrint __P((double));
extern long int lrintl __P((long double));
+extern long long int llroundf __P((float));
+extern long long int llround __P((double));
+extern long long int llroundl __P((long double));
+extern long int lroundf __P((float));
+extern long int lround __P((double));
+extern long int lroundl __P((long double));
extern float roundf __P((float));
extern double round __P((double));
extern long double roundl __P((long double));
extern float truncf __P((float));
extern double trunc __P((double));
diff -aprNU5 djgpp.orig/include/math.h djgpp/include/math.h
--- djgpp.orig/include/math.h 2013-10-14 22:37:58 +0100
+++ djgpp/include/math.h 2013-10-14 22:40:02 +0100
@@ -195,10 +195,14 @@ extern long int lrint(double);
extern long int lrintl(long double);
extern long long int llrint(double);
extern long long int llrintl(long double);
extern double round(double);
extern long double roundl(long double);
+extern long int lround(double);
+extern long int lroundl(long double);
+extern long long int llround(double);
+extern long long int llroundl(long double);
extern float erff(float);
extern float erfcf(float);
extern float hypotf(float, float);
extern float lgammaf(float);
extern float acoshf(float);
@@ -217,10 +221,12 @@ extern float expm1f(float);
extern float log1pf(float);
extern float truncf(float);
extern long int lrintf(float);
extern long long int llrintf(float);
extern float roundf(float);
+extern long int lroundf(float);
+extern long long int llroundf(float);
/* End libm.a. */
#endif /* (__STDC_VERSION__ >= 199901L) || !__STRICT_ANSI__ */
diff -aprNU5 djgpp.orig/src/docs/kb/wc204.txi djgpp/src/docs/kb/wc204.txi
--- djgpp.orig/src/docs/kb/wc204.txi 2013-10-14 22:37:58 +0100
+++ djgpp/src/docs/kb/wc204.txi 2013-10-19 12:39:24 +0100
@@ -1327,5 +1327,16 @@ in @file{libc.a} and one accurate in @fi
@findex roundf AT r{ added}
@findex round AT r{ added}
@findex roundl AT r{ added}
The @acronym{C99} functions @code{roundf}, @code{round} and @code{roundl} were added
to comply with the @acronym{C99} standard.
+
+@cindex @acronym{C99} compliance, @code{math.h}
+@findex lroundf AT r{ added}
+@findex lround AT r{ added}
+@findex lroundl AT r{ added}
+@findex llroundf AT r{ added}
+@findex llround AT r{ added}
+@findex llroundl AT r{ added}
+The @acronym{C99} functions @code{lroundf}, @code{lround}, @code{lroundl}, @code{llroundf},
+@code{llround} and @code{llroundl} were added to comply with the @acronym{C99} standard.
+These functions are available in two versions.
diff -aprNU5 djgpp.orig/src/libm/math/llround.c djgpp/src/libm/math/llround.c
--- djgpp.orig/src/libm/math/llround.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/llround.c 2013-10-17 00:14:46 +0100
@@ -0,0 +1,111 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+/*
+<<llround>>, <<llroundf>>, <<llroundl>>--round to integer, to nearest
+INDEX
+ llround
+INDEX
+ llroundf
+INDEX
+ llroundl
+
+ANSI_SYNOPSIS
+ #include <math.h>
+ long long int llround(double <[x]>);
+ long long int llroundf(float <[x]>);
+ long long int llroundl(long double <[x]>);
+
+DESCRIPTION
+ The <<llround>>, <<llroundl>> and <<llroundl>> functions round their
+ argument to the nearest integer value, rounding halfway cases away from
+ zero, regardless of the current rounding direction. If the rounded
+ value is outside the range of the return type, the numeric result is
+ unspecified (depending upon the floating-point implementation, not the
+ library). A range error may occur if the magnitude of x is too large.
+
+RETURNS
+ <[x]> rounded to an integral value as an integer.
+
+SEEALSO
+ See the <<round>> functions for the return being the same floating-point
+ type as the argument. See also <<llrint>>.
+
+PORTABILITY
+ANSI C, POSIX
+
+*/
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#if defined (__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 8))
+# define __gnuc_extension__ __extension__
+#else
+# define __gnuc_extension__
+#endif
+
+#define DOUBLE_BIAS (0x3FFU)
+#define BIN_DIGITS_IN_FRACTION (52) /* Amount of binary digits in fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAH (20) /* Amount of binary digits in msw of the fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAL (BIN_DIGITS_IN_FRACTION - BIN_DIGITS_IN_MANTISSAH)
+#define ALL_DIGITS_ARE_SIGNIFICANT(exp) ((exp) > (BIN_DIGITS_IN_FRACTION - 1))
+#define NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(exp) ((exp) < BIN_DIGITS_IN_MANTISSAH)
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+
+#define ROUND_MANTISSAH_TO_INTEGER(num, unbiased_exponent) ((long long int)((((uint32_t)(num).dt.mantissah | 0x00100000UL) + (0x0080000UL >> (unbiased_exponent))) >> (BIN_DIGITS_IN_MANTISSAH - (unbiased_exponent))))
+
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long long int)(((((uint64_t)(num).dt.mantissah | 0x00100000UL) << BIN_DIGITS_IN_MANTISSAL) | (uint64_t)(num).dt.mantissal) << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION)))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint32_t rounded_mantissal = (uint32_t)(num).dt.mantissal + (0x80000000UL >> ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH)); \
+ if (result = ((long long int)(num).dt.mantissah | 0x00100000UL), rounded_mantissal < (num).dt.mantissal) result++; \
+ if ((unbiased_exponent) > BIN_DIGITS_IN_MANTISSAH) \
+ result = result << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH) | rounded_mantissal >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent)); \
+ (long long int)result; \
+ }) \
+)
+
+
+#ifdef __STDC__
+long long int
+llround(double x)
+#else
+long long int
+llround(x)
+double x;
+#endif
+{
+ _double_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long long int result;
+
+
+ ieee_value.d = x;
+ unbiased_exponent = ieee_value.dt.exponent - DOUBLE_BIAS;
+ sign = ieee_value.dt.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long long int)x; /* It is left implementation defined what happens. */
+ else if (NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(unbiased_exponent))
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else
+ result = ROUND_MANTISSAH_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+ else
+ {
+ if (ALL_DIGITS_ARE_SIGNIFICANT(unbiased_exponent))
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* >= 2^52 is already an exact integer iff long long int is 64 bit. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/llroundf.c djgpp/src/libm/math/llroundf.c
--- djgpp.orig/src/libm/math/llroundf.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/llroundf.c 2013-10-16 22:34:56 +0100
@@ -0,0 +1,47 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#define FLOAT_BIAS (0x7FU)
+#define BIN_DIGITS_IN_FRACTION (23) /* Amount of binary digits in fraction part of mantissa. */
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long long int)((uint32_t)(num).ft.mantissa | 0x00800000UL) << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long long int)((((uint32_t)(num).ft.mantissa | 0x00800000UL) + (0x00400000UL >> (unbiased_exponent))) >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent))))
+
+
+#ifdef __STDC__
+long long int
+llroundf(float x)
+#else
+long long int
+llroundf(x)
+float x;
+#endif
+{
+ _float_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long long int result;
+
+
+ ieee_value.f = x;
+ unbiased_exponent = ieee_value.ft.exponent - FLOAT_BIAS;
+ sign = ieee_value.ft.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long long int)x; /* It is left implementation defined what happens. */
+ else
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else if (unbiased_exponent > (BIN_DIGITS_IN_FRACTION - 1))
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* >= 2^23 is already an exact integer. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/llroundl.c djgpp/src/libm/math/llroundl.c
--- djgpp.orig/src/libm/math/llroundl.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/llroundl.c 2013-10-19 11:58:36 +0100
@@ -0,0 +1,90 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#if defined (__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 8))
+# define __gnuc_extension__ __extension__
+#else
+# define __gnuc_extension__
+#endif
+
+#define LONG_DOUBLE_BIAS (0x3FFFU)
+#define BIN_DIGITS_IN_FRACTION (63) /* Amount of binary digits in fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAH (31) /* Amount of binary digits in msw of the fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAL (BIN_DIGITS_IN_FRACTION - BIN_DIGITS_IN_MANTISSAH)
+#define ALL_DIGITS_ARE_SIGNIFICANT(exp) ((exp) > (BIN_DIGITS_IN_FRACTION - 1))
+#define NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(exp) ((exp) < BIN_DIGITS_IN_MANTISSAH)
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+
+#define CARRY_INTO_INTEGER_PART(x) ((x) & 0x100000000ULL)
+#define ROUND_MANTISSAH_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint64_t rounded_mantissah = (uint64_t)(num).ldt.mantissah + (0x40000000ULL >> (unbiased_exponent)); \
+ if (CARRY_INTO_INTEGER_PART(rounded_mantissah)) \
+ { \
+ rounded_mantissah >>= 1; \
+ unbiased_exponent++; \
+ } \
+ (long long int)(rounded_mantissah >> (BIN_DIGITS_IN_MANTISSAH - (unbiased_exponent))); \
+ }) \
+)
+
+#define SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) ((uint64_t)(num).ldt.mantissah << BIN_DIGITS_IN_MANTISSAL)
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long long int)((SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) | (uint64_t)(num).ldt.mantissal) << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION)))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint32_t rounded_mantissal = (uint32_t)(num).ldt.mantissal + (0x80000000UL >> ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH)); \
+ if (result = (long long int)(num).ldt.mantissah, rounded_mantissal < (num).ldt.mantissal) \
+ result++; \
+ if ((unbiased_exponent) > BIN_DIGITS_IN_MANTISSAH) \
+ result = result << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH) | rounded_mantissal >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent)); \
+ (long long int)result; \
+ }) \
+)
+
+
+#ifdef __STDC__
+long long int
+llroundl(long double x)
+#else
+long long int
+llroundl(x)
+long double x;
+#endif
+{
+ _longdouble_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long long int result;
+
+
+ ieee_value.ld = x;
+ unbiased_exponent = ieee_value.ldt.exponent - LONG_DOUBLE_BIAS;
+ sign = ieee_value.ldt.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long long int)x; /* It is left implementation defined what happens. */
+ else if (NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(unbiased_exponent))
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else
+ result = ROUND_MANTISSAH_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+ else
+ {
+ if (ALL_DIGITS_ARE_SIGNIFICANT(unbiased_exponent))
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* >= 2^63 is already an exact integer. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/lround.c djgpp/src/libm/math/lround.c
--- djgpp.orig/src/libm/math/lround.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/lround.c 2013-10-16 22:34:56 +0100
@@ -0,0 +1,112 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+/*
+<<lround>>, <<lroundf>>, <<lroundl>>--round to integer, to nearest
+INDEX
+ lround
+INDEX
+ lroundf
+INDEX
+ lroundl
+
+ANSI_SYNOPSIS
+ #include <math.h>
+ long int lround(double <[x]>);
+ long int lroundf(float <[x]>);
+ long int lroundl(long double <[x]>);
+
+DESCRIPTION
+ The <<lround>>, <<lroundl>> and <<lroundl>> functions round their
+ argument to the nearest integer value, rounding halfway cases away from
+ zero, regardless of the current rounding direction. If the rounded
+ value is outside the range of the return type, the numeric result is
+ unspecified (depending upon the floating-point implementation, not the
+ library). A range error may occur if the magnitude of x is too large.
+
+RETURNS
+ <[x]> rounded to an integral value as an integer.
+
+SEEALSO
+ See the <<round>> functions for the return being the same floating-point
+ type as the argument. See also <<lrint>>.
+
+PORTABILITY
+ANSI C, POSIX
+
+*/
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#if defined (__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 8))
+# define __gnuc_extension__ __extension__
+#else
+# define __gnuc_extension__
+#endif
+
+#define DOUBLE_BIAS (0x3FFU)
+#define BIN_DIGITS_IN_FRACTION (52) /* Amount of binary digits in fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAH (20) /* Amount of binary digits in msw of the fraction part of mantissa. */
+#define ALL_DIGITS_ARE_SIGNIFICANT(exp) ((exp) > (BIN_DIGITS_IN_FRACTION - 1))
+#define NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(exp) ((exp) < BIN_DIGITS_IN_MANTISSAH)
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+
+#define ROUND_MANTISSAH_TO_INTEGER(num, unbiased_exponent) ((long int)((((uint32_t)(num).dt.mantissah | 0x00100000UL) + (0x0080000UL >> (unbiased_exponent))) >> (BIN_DIGITS_IN_MANTISSAH - (unbiased_exponent))))
+
+#define SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) (((uint64_t)(num).dt.mantissah | 0x00100000ULL) << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH))
+#define SHIFT_LEFT_MANTISSAL(num, unbiased_exponent) ((uint64_t)(num).dt.mantissal << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION))
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long int)(SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) | SHIFT_LEFT_MANTISSAL(num, unbiased_exponent)))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint32_t rounded_mantissal = (uint32_t)(num).dt.mantissal + (0x80000000UL >> ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH)); \
+ if (result = ((long int)(num).dt.mantissah | 0x00100000UL), rounded_mantissal < (num).dt.mantissal) result++; \
+ if ((unbiased_exponent) > BIN_DIGITS_IN_MANTISSAH) \
+ result = result << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH) | rounded_mantissal >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent)); \
+ (long int)result; \
+ }) \
+)
+
+
+#ifdef __STDC__
+long int
+lround(double x)
+#else
+long int
+lround(x)
+double x;
+#endif
+{
+ _double_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long int result;
+
+
+ ieee_value.d = x;
+ unbiased_exponent = ieee_value.dt.exponent - DOUBLE_BIAS;
+ sign = ieee_value.dt.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long int)x; /* It is left implementation defined what happens. */
+ else if (NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(unbiased_exponent))
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else
+ result = ROUND_MANTISSAH_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+ else
+ {
+ if (ALL_DIGITS_ARE_SIGNIFICANT(unbiased_exponent)) /* >= 2^52 is already an exact integer iff long int is 64 bit. */
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* But this is not the case with djgpp. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/lroundf.c djgpp/src/libm/math/lroundf.c
--- djgpp.orig/src/libm/math/lroundf.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/lroundf.c 2013-10-16 22:34:54 +0100
@@ -0,0 +1,47 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#define FLOAT_BIAS (0x7FU)
+#define BIN_DIGITS_IN_FRACTION (23) /* Amount of binary digits in fraction part of mantissa. */
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long int)((uint32_t)(num).ft.mantissa | 0x00800000UL) << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long int)((((uint32_t)(num).ft.mantissa | 0x00800000UL) + (0x00400000UL >> (unbiased_exponent))) >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent))))
+
+
+#ifdef __STDC__
+long int
+lroundf(float x)
+#else
+long int
+lroundf(x)
+float x;
+#endif
+{
+ _float_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long int result;
+
+
+ ieee_value.f = x;
+ unbiased_exponent = ieee_value.ft.exponent - FLOAT_BIAS;
+ sign = ieee_value.ft.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long int)x; /* It is left implementation defined what happens. */
+ else
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else if (unbiased_exponent > (BIN_DIGITS_IN_FRACTION - 1))
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* >= 2^23 is already an exact integer. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/lroundl.c djgpp/src/libm/math/lroundl.c
--- djgpp.orig/src/libm/math/lroundl.c 1970-01-01 01:00:00 +0100
+++ djgpp/src/libm/math/lroundl.c 2013-10-19 10:19:16 +0100
@@ -0,0 +1,94 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+#include <stdint.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+#if defined (__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 8))
+# define __gnuc_extension__ __extension__
+#else
+# define __gnuc_extension__
+#endif
+
+#define LONG_DOUBLE_BIAS (0x3FFFU)
+#define BIN_DIGITS_IN_FRACTION (63) /* Amount of binary digits in fraction part of mantissa. */
+#define BIN_DIGITS_IN_MANTISSAH (31) /* Amount of binary digits in msw of the fraction part of mantissa. */
+#define ALL_DIGITS_ARE_SIGNIFICANT(exp) ((exp) > (BIN_DIGITS_IN_FRACTION - 1))
+#define NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(exp) ((exp) < BIN_DIGITS_IN_MANTISSAH)
+#define MAGNITUDE_IS_TOO_LARGE(exp) ((exp) > (int)(sizeof(long int) * 8) - 2)
+#define MAGNITUDE_IS_LESS_THAN_ONE(exp) ((exp) < 0)
+#define MAGNITUDE_IS_LESS_THAN_ONE_HALF(exp) ((exp) < -1)
+
+#define CARRY_INTO_INTEGER_PART(x) ((x) & 0x100000000ULL)
+#define ROUND_MANTISSAH_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint64_t rounded_mantissah = (uint64_t)(num).ldt.mantissah + (0x40000000ULL >> (unbiased_exponent)); \
+ if (CARRY_INTO_INTEGER_PART(rounded_mantissah)) \
+ { \
+ rounded_mantissah >>= 1; \
+ (unbiased_exponent)++; \
+ } \
+ (long int)(rounded_mantissah >> (BIN_DIGITS_IN_MANTISSAH - (unbiased_exponent))); \
+ }) \
+)
+
+#define SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) ((uint64_t)(num).ldt.mantissah << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH))
+#define SHIFT_LEFT_MANTISSAL(num, unbiased_exponent) ((uint64_t)(num).ldt.mantissal << ((unbiased_exponent) - BIN_DIGITS_IN_FRACTION))
+#define CONVERT_MANTISSA_TO_INTEGER(num, unbiased_exponent) ((long int)(SHIFT_LEFT_MANTISSAH(num, unbiased_exponent) | SHIFT_LEFT_MANTISSAL(num, unbiased_exponent)))
+#define ROUND_MANTISSA_TO_INTEGER(num, unbiased_exponent) \
+(__gnuc_extension__ \
+ ({ \
+ uint32_t rounded_mantissal = (uint32_t)(num).ldt.mantissal + (0x80000000UL >> ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH)); \
+ if (result = (long int)(num).ldt.mantissah, rounded_mantissal < (num).ldt.mantissal) \
+ { \
+ uint64_t rounded_mantissah = (uint64_t)result + 1; \
+ if (CARRY_INTO_INTEGER_PART(rounded_mantissah)) \
+ result = (long int)(rounded_mantissah >> 1); \
+ } \
+ if ((unbiased_exponent) > BIN_DIGITS_IN_MANTISSAH) \
+ result = result << ((unbiased_exponent) - BIN_DIGITS_IN_MANTISSAH) | rounded_mantissal >> (BIN_DIGITS_IN_FRACTION - (unbiased_exponent)); \
+ (long int)result; \
+ }) \
+)
+
+
+#ifdef __STDC__
+long int
+lroundl(long double x)
+#else
+long int
+lroundl(x)
+long double x;
+#endif
+{
+ _longdouble_union_t ieee_value;
+ int unbiased_exponent, sign;
+ long int result;
+
+
+ ieee_value.ld = x;
+ unbiased_exponent = ieee_value.ldt.exponent - LONG_DOUBLE_BIAS;
+ sign = ieee_value.ldt.sign;
+
+ if (MAGNITUDE_IS_TOO_LARGE(unbiased_exponent)) /* The number is too large. */
+ return (long int)x; /* It is left implementation defined what happens. */
+ else if (NO_SIGNIFICANT_DIGITS_IN_MANTISSAL(unbiased_exponent))
+ {
+ if (MAGNITUDE_IS_LESS_THAN_ONE(unbiased_exponent))
+ result = MAGNITUDE_IS_LESS_THAN_ONE_HALF(unbiased_exponent) ? 0 : 1;
+ else
+ result = ROUND_MANTISSAH_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+ else
+ {
+ if (ALL_DIGITS_ARE_SIGNIFICANT(unbiased_exponent)) /* >= 2^63 is already an exact integer iff long int is 64 bit. */
+ result = CONVERT_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent); /* But this is not the case with djgpp. */
+ else
+ result = ROUND_MANTISSA_TO_INTEGER(ieee_value, unbiased_exponent);
+
+ return sign ? -result : result;
+ }
+}
diff -aprNU5 djgpp.orig/src/libm/math/makefile djgpp/src/libm/math/makefile
--- djgpp.orig/src/libm/math/makefile 2013-10-14 22:37:58 +0100
+++ djgpp/src/libm/math/makefile 2013-10-17 20:37:26 +0100
@@ -174,10 +174,16 @@ SRC += llrint.c
SRC += llrintf.c
SRC += llrintl.c
SRC += round.c
SRC += roundf.c
SRC += roundl.c
+SRC += lround.c
+SRC += lroundf.c
+SRC += lroundl.c
+SRC += llround.c
+SRC += llroundf.c
+SRC += llroundl.c
chobj = w_acos.def w_acosh.def w_asin.def s_asinh.def \
s_atan.def w_atan2.def w_atanh.def w_j0.def \
s_copysign.def w_cosh.def s_erf.def w_exp.def \
s_fabs.def s_floor.def w_fmod.def s_frexp.def \
@@ -185,11 +191,11 @@ chobj = w_acos.def w_acosh.def w_asin.de
w_log10.def s_log1p.def s_matherr.def s_modf.def \
w_pow.def w_remainder.def s_sin.def w_sinh.def \
s_cbrt.def w_sqrt.def s_tan.def s_tanh.def \
s_infinity.def s_isnan.def s_scalbn.def s_nextafter.def \
s_nan.def s_ilogb.def s_expm1.def trunc.def \
- lrint.def llrint.def round.def
+ lrint.def llrint.def round.def lround.def llround.def
CFLAGS = -D_USE_LIBM_MATH_H
EXTRA_FILES = $(TOP)/../../info/libm.info
# chew emits non-fatal warnings, so we redirect them to the void
CHEW = ./chew.exe -f ./doc.str -e /dev/null
diff -aprNU5 djgpp.orig/tests/cygnus/makefile djgpp/tests/cygnus/makefile
--- djgpp.orig/tests/cygnus/makefile 2013-10-14 22:37:58 +0100
+++ djgpp/tests/cygnus/makefile 2013-10-14 22:40:02 +0100
@@ -93,11 +93,11 @@ VEC_OFILES = $(GEN_VEC_FILES:.c=.o)
$(OFILES): CFLAGS = $(DEFS) -fno-builtin -O2 -g -Wall
all: check
check: mtest.exe t-trunc.exe t-truncf.exe t-truncl.exe t-lrint.exe t-lrintf.exe t-lrintl.exe t-llrint.exe t-llrintl.exe \
- t-roundf.exe t-round.exe t-roundl.exe
+ t-roundf.exe t-round.exe t-roundl.exe t-lround.exe t-lroundf.exe t-lroundl.exe t-llround.exe t-llroundf.exe t-llroundl.exe
./mtest.exe > mtest.results
./t-trunc.exe > ttest.results
./t-truncf.exe >> ttest.results
./t-truncl.exe >> ttest.results
./t-lrintf.exe >> ttest.results
@@ -106,10 +106,16 @@ check: mtest.exe t-trunc.exe t-truncf.ex
./t-llrint.exe >> ttest.results
./t-llrintl.exe >> ttest.results
./t-roundf.exe >> ttest.results
./t-round.exe >> ttest.results
./t-roundl.exe >> ttest.results
+ ./t-lroundf.exe >> ttest.results
+ ./t-lround.exe >> ttest.results
+ ./t-lroundl.exe >> ttest.results
+ ./t-llroundf.exe >> ttest.results
+ ./t-llround.exe >> ttest.results
+ ./t-llroundl.exe >> ttest.results
# Pattern rules to generate test vectors. (The funky vec.c=%.c replacement
# is meant to create a pattern rule where actually a normal rule will
# do, since only pattern rules can tell Make that several targets are
# generated all at once. Without this, Make will invoke the vector-
@@ -168,13 +174,31 @@ t-round.exe: t-round.o
$(CC) -o $@ $(LDFLAGS) t-round.o $(LIBS)
t-roundl.exe: t-roundl.o
$(CC) -o $@ $(LDFLAGS) t-roundl.o $(LIBS)
+t-lroundf.exe: t-lroundf.o
+ $(CC) -o $@ $(LDFLAGS) t-lroundf.o $(LIBS)
+
+t-lround.exe: t-lround.o
+ $(CC) -o $@ $(LDFLAGS) t-lround.o $(LIBS)
+
+t-lroundl.exe: t-lroundl.o
+ $(CC) -o $@ $(LDFLAGS) t-lroundl.o $(LIBS)
+
+t-llroundf.exe: t-llroundf.o
+ $(CC) -o $@ $(LDFLAGS) t-llroundf.o $(LIBS)
+
+t-llround.exe: t-llround.o
+ $(CC) -o $@ $(LDFLAGS) t-llround.o $(LIBS)
+
+t-llroundl.exe: t-llroundl.o
+ $(CC) -o $@ $(LDFLAGS) t-llroundl.o $(LIBS)
+
$(OFILES) $(VEC_OFILES) : test.h
clean mostlyclean:
-cd tgen; $(MAKE) $@
- cd $(HERE); $(RM) $(OFILES) $(VEC_OFILES) *~ *.exe mtest.results ttest.results t-trunc*.o t-lrint*.o t-llrint*.o t-round*.o
+ cd $(HERE); $(RM) $(OFILES) $(VEC_OFILES) *~ *.exe mtest.results ttest.results t-trunc*.o t-lrint*.o t-llrint*.o t-round*.o t-lround*.o t-llround*.o
.SECONDARY: $(GEN_PROGS) $(GEN_VEC_FILES)
.PHONY: all check clean mostlyclean
diff -aprNU5 djgpp.orig/tests/cygnus/t-llround.c djgpp/tests/cygnus/t-llround.c
--- djgpp.orig/tests/cygnus/t-llround.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-llround.c 2013-10-16 22:20:50 +0100
@@ -0,0 +1,198 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _double_union_t value; /* test value */
+ const long long int should_be; /* result */
+} entry_t;
+
+static const entry_t tests_double[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.dt = {0x0U, 0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.dt = {0x0U, 0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.dt = {0x1U, 0x0U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.dt = {0x1U, 0x0U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.dt = {0x1U, 0x0U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.dt = {0x1U, 0x0U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.dt = {0xFFFFFFFFU, 0x7FFFFU, 0x7FEU, 0}}, -9.223372036854775808E18}, /* Big number. */
+ {{.dt = {0xFFFFFFFFU, 0x7FFFFU, 0x7FEU, 1}}, -9.223372036854775808E18}, /* Big -number. */
+
+ /* Infs. */
+ {{.dt = {0x0U, 0x0U, 0x7FFU, 0}}, -9.223372036854775808E18}, /* Inf */
+ {{.dt = {0x0U, 0x0U, 0x7FFU, 1}}, -9.223372036854775808E18}, /* -Inf */
+
+ /* NaNs. */
+ {{.dt = {0x1U, 0x0U, 0x7FFU, 0}}, -9.223372036854775808E18}, /* SNaN */
+ {{.dt = {0x1U, 0x0U, 0x7FFU, 1}}, -9.223372036854775808E18}, /* -SNaN */
+ {{.dt = {0x0U, 0xFFFFFU, 0x7FFU, 1}}, -9.223372036854775808E18}, /* QNaN */
+ {{.dt = {0x0U, 0xFFFFFU, 0x7FFU, 0}}, -9.223372036854775808E18}, /* -QNaN */
+
+
+ /* Number. */
+ {{.dt = {0x54442D18U, 0x921FBU, 0x3FFU + 0x001U, 0}}, +3}, /* PI */
+ {{.dt = {0x54442D18U, 0x921FBU, 0x3FFU + 0x001U, 1}}, -3}, /* -PI */
+
+ {{.dt = {0x00000000U, 0xE0000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.875000 */
+ {{.dt = {0x00000000U, 0xE0000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.875000 */
+ {{.dt = {0x00000000U, 0xA0000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.625000 */
+ {{.dt = {0x00000000U, 0xA0000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.625000 */
+ {{.dt = {0x18DEF417U, 0x80002U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.500002 */
+ {{.dt = {0x18DEF417U, 0x80002U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.500002 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.500000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.500000 */
+ {{.dt = {0xE7210BE9U, 0x7F9FDU, 0x3FFU + 0x0U, 0}}, +1}, /* 1.499998 */
+ {{.dt = {0xE7210BE9U, 0x7F9FDU, 0x3FFU + 0x0U, 1}}, -1}, /* -1.499998 */
+ {{.dt = {0x00000000U, 0x60000U, 0x3FFU + 0x0U, 0}}, +1}, /* 1.375000 */
+ {{.dt = {0x00000000U, 0x60000U, 0x3FFU + 0x0U, 1}}, -1}, /* -1.375000 */
+ {{.dt = {0x00000000U, 0x20000U, 0x3FFU + 0x0U, 0}}, +1}, /* 1.125000 */
+ {{.dt = {0x00000000U, 0x20000U, 0x3FFU + 0x0U, 1}}, -1}, /* -1.125000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x016U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x016U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x017U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x017U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x018U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x018U, 1}}, -16777216}, /* -16777216.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 0}}, +2147483647LL}, /* 2147483647.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x020U, 0}}, +4294967296}, /* 4294967296.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x020U, 1}}, -4294967296}, /* -4294967296.000000 */
+ {{.dt = {0xACF13400U, 0x02468U, 0x3FFU + 0x033U, 0}}, 2271815812028928}, /* 2271815812028928.000000 */
+ {{.dt = {0xACF13400U, 0x02468U, 0x3FFU + 0x033U, 1}}, -2271815812028928}, /* -2271815812028928.000000 */
+
+ {{.dt = {0x56789AB0U, 0x01234U, 0x3FFU + 0x034U, 0}}, 4523615625714352}, /* 4523615625714352.000000 */
+ {{.dt = {0x56789AB0U, 0x01234U, 0x3FFU + 0x034U, 1}}, -4523615625714352}, /* -4523615625714352.000000 */
+ {{.dt = {0xA9876543U, 0xFEDCBU, 0x3FFU + 0x034U, 0}}, 8987183256397123}, /* 8987183256397123.000000 */
+ {{.dt = {0xA9876543U, 0xFEDCBU, 0x3FFU + 0x034U, 1}}, -8987183256397123}, /* -8987183256397123.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x035U, 0}}, 9007199254740992}, /* 9007199254740992.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x035U, 1}}, -9007199254740992}, /* -9007199254740992.000000 */
+ {{.dt = {0x6789ABCEU, 0x12345U, 0x3FFU + 0x035U, 0}}, 9647711201744796}, /* 9647711201744796.000000 */
+ {{.dt = {0x6789ABCEU, 0x12345U, 0x3FFU + 0x035U, 1}}, -9647711201744796}, /* -9647711201744796.000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x03FU, 0}}, 9223372036854775808ULL}, /* 9223372036854775808.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x03FU, 1}}, -9.223372036854775808E+18}, /* -9223372036854775808.0000000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x03EU, 0}}, 4611686018427387904}, /* 4611686018427387904.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x03EU, 1}}, -4611686018427387904}, /* -4611686018427387904.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x034U, 0}}, 4503599627370496}, /* 4503599627370496.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x034U, 1}}, -4503599627370496}, /* -4503599627370496.0000000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x033U, 0}}, 2251799813685248}, /* 2251799813685248.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x033U, 1}}, -2251799813685248}, /* -2251799813685248.0000000000 */
+ {{.dt = {0x00000001U, 0x00000U, 0x3FFU + 0x033U, 0}}, 2251799813685249}, /* 2251799813685248.5000000000 */
+ {{.dt = {0x00000001U, 0x00000U, 0x3FFU + 0x033U, 1}}, -2251799813685249}, /* -2251799813685248.5000000000 */
+ {{.dt = {0x00000002U, 0x00000U, 0x3FFU + 0x033U, 0}}, 2251799813685249}, /* 2251799813685249.7500000000 */
+ {{.dt = {0x00000002U, 0x00000U, 0x3FFU + 0x033U, 1}}, -2251799813685249}, /* -2251799813685249.7500000000 */
+
+ {{.dt = {0x00000000U, 0xC0000U, 0x3FFU + 0x000U, 0}}, 2}, /* 1.7500000000 */
+ {{.dt = {0x00000000U, 0xC0000U, 0x3FFU + 0x000U, 1}}, -2}, /* -1.7500000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x000U, 0}}, 2}, /* 1.5000000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x000U, 1}}, -2}, /* -1.5000000000 */
+ {{.dt = {0x00000000U, 0x40000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.2500000000 */
+ {{.dt = {0x00000000U, 0x40000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.2500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.7500000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.7500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.2500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.2500000000 */
+ {{.dt = {0xF8FAB5E6U, 0x948B0U, 0x3FFU + 0xFFFFFFF9U, 0}}, 0}, /* 0.0123456789 */
+ {{.dt = {0xF8FAB5E6U, 0x948B0U, 0x3FFU + 0xFFFFFFF9U, 1}}, -0}, /* -0.0123456789 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFC01U, 0}}, 0}, /* 0.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFC01U, 1}}, -0}, /* -0.0000000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048576}, /* 1048576.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048576}, /* -1048576.0000000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048576}, /* 1048576.2500000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048576}, /* -1048576.2500000000 */
+ {{.dt = {0x80000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.5000000000 */
+ {{.dt = {0x80000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.5000000000 */
+ {{.dt = {0xC0000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.7500000000 */
+ {{.dt = {0xC0000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.7500000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097152}, /* 2097152.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097152}, /* -2097152.0000000000 */
+ {{.dt = {0x20000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097152}, /* 2097152.2500000000 */
+ {{.dt = {0x20000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097152}, /* -2097152.2500000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097153}, /* 2097152.5000000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097153}, /* -2097152.5000000000 */
+ {{.dt = {0x60000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097153}, /* 2097152.7500000000 */
+ {{.dt = {0x60000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097153}, /* -2097152.7500000000 */
+
+ /* Number greater than 2**63 exceeds long int size and will be casted in an implementation defined manner. */
+ {{.dt = {0x00000001U, 0x00000U, 0x3FFU + 0x03FU, 0}}, -9.223372036854775808E+18}, /* 9223372036854777856.0000000000 */
+ {{.dt = {0x00000001U, 0x00000U, 0x3FFU + 0x03FU, 1}}, -9.223372036854775808E+18}, /* -9223372036854777856.0000000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.dt = {0x1FF885F3U, 0xFF793U, 0x3FFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.dt = {0x1FF885F3U, 0xFF793U, 0x3FFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.dt = {0xD2F1A9FCU, 0x0624DU, 0x3FFU + 0xFFFFFFF6U, 0}}, 0}, /* 0.0010000000 */
+ {{.dt = {0xD2F1A9FCU, 0x0624DU, 0x3FFU + 0xFFFFFFF6U, 1}}, -0}, /* -0.0010000000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.dt = {0x8FB86F48U, 0xF9ADCU, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.9876540000 */
+ {{.dt = {0x8FB86F48U, 0xF9ADCU, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.9876540000 */
+ {{.dt = {0x20365653U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000600 */
+ {{.dt = {0x20365653U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000600 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+
+ /* Number greather than 1 and less than 2**63 will be rounded accordingly. */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000000000 */
+ {{.dt = {0x18DEF417U, 0x00002U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000020000 */
+ {{.dt = {0x18DEF417U, 0x00002U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000020000 */
+ {{.dt = {0x00000000U, 0x00003U, 0x3FFU + 0x010U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.dt = {0x00000000U, 0x00003U, 0x3FFU + 0x010U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.dt = {0x00000000U, 0x00008U, 0x3FFU + 0x010U, 0}}, 65537}, /* 65536.5000000000 */
+ {{.dt = {0x00000000U, 0x00008U, 0x3FFU + 0x010U, 1}}, -65537}, /* -65536.5000000000 */
+ {{.dt = {0xA0000000U, 0x00009U, 0x3FFU + 0x010U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.dt = {0xA0000000U, 0x00009U, 0x3FFU + 0x010U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.dt = {0xE0000000U, 0xFFFFFU, 0x3FFU + 0x016U, }}, 8388608}, /* 8388607.5000000000 */
+ {{.dt = {0xE0000000U, 0xFFFFFU, 0x3FFU + 0x016U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.dt = {0xA0000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.dt = {0xA0000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.6250000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 0}}, 2147483647}, /* 2147483647.0000000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 1}}, -2147483647}, /* -2147483647.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x034U, 0}}, 4503599627370496}, /* 4503599627370496.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x034U, 1}}, -4503599627370496}, /* -4503599627370496.0000000000 */
+ {{.dt = {0xFFFFFFFEU, 0xFFFFFU, 0x3FFU + 0x03EU, 0}}, 9223372036854773760}, /* 9223372036854773760.0000000000 */
+ {{.dt = {0xFFFFFFFEU, 0xFFFFFU, 0x3FFU + 0x03EU, 1}}, -9223372036854773760}, /* -9223372036854773760.0000000000 */
+};
+
+static const size_t n_tests_double = sizeof(tests_double) / sizeof(tests_double[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_double; i++)
+ {
+ long long int result = llround(tests_double[i].value.d);
+
+ if (tests_double[i].should_be == result)
+ counter++;
+ else
+ printf("llround test failed: value to round = %.6g result = %lld should be = %lld\n", tests_double[i].value.d, result, tests_double[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_double) ? "llround test failed." : "llround test succeded.");
+
+ return 0;
+}
diff -aprNU5 djgpp.orig/tests/cygnus/t-llroundf.c djgpp/tests/cygnus/t-llroundf.c
--- djgpp.orig/tests/cygnus/t-llroundf.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-llroundf.c 2013-10-14 22:40:02 +0100
@@ -0,0 +1,137 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _float_union_t value; /* test value */
+ const long long int should_be; /* result */
+} entry_t;
+
+static const entry_t tests_float[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.ft = {0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.ft = {0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.ft = {0x1U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.ft = {0x1U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.ft = {0x1U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.ft = {0x1U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.ft = {0xFFFFU, 0xFEU, 0}}, -9223372036854775808ULL}, /* Big number. */
+ {{.ft = {0xFFFFU, 0xFEU, 1}}, -9223372036854775808ULL}, /* Big -number. */
+
+ /* Infs. */
+ {{.ft = {0x0U, 0xFFU, 0}}, -9223372036854775808ULL}, /* Inf */
+ {{.ft = {0x0U, 0xFFU, 1}}, -9223372036854775808ULL}, /* -Inf */
+
+ /* NaNs. */
+ {{.ft = {0x1U, 0xFFU, 0}}, -9223372036854775808ULL}, /* SNaN */
+ {{.ft = {0x1U, 0xFFU, 1}}, -9223372036854775808ULL}, /* -SNaN */
+ {{.ft = {0x7FFFFFU, 0xFFU, 0}}, -9223372036854775808ULL}, /* QNaN */
+ {{.ft = {0x7FFFFFU, 0xFFU, 1}}, -9223372036854775808ULL}, /* -QNaN */
+
+ /* Numbers. */
+ {{.ft = {0x490FDBU, 0x80U, 0}}, +3}, /* PI */
+ {{.ft = {0x490FDBU, 0x80U, 1}}, -3}, /* -PI */
+
+ {{.ft = {0x700000U, 0x7FU, 0}}, +2}, /* 1.875000 */
+ {{.ft = {0x700000U, 0x7FU, 1}}, -2}, /* -1.875000 */
+ {{.ft = {0x500000U, 0x7FU, 0}}, +2}, /* 1.625000 */
+ {{.ft = {0x500000U, 0x7FU, 1}}, -2}, /* -1.625000 */
+ {{.ft = {0x40000FU, 0x7FU, 0}}, +2}, /* 1.500002 */
+ {{.ft = {0x40000FU, 0x7FU, 1}}, -2}, /* -1.500002 */
+ {{.ft = {0x400000U, 0x7FU, 0}}, +2}, /* 1.500000 */
+ {{.ft = {0x400000U, 0x7FU, 1}}, -2}, /* -1.500000 */
+ {{.ft = {0x3FFFF0U, 0x7FU, 0}}, +1}, /* 1.499998 */
+ {{.ft = {0x3FFFF0U, 0x7FU, 1}}, -1}, /* -1.499998 */
+ {{.ft = {0x300000U, 0x7FU, 0}}, +1}, /* 1.375000 */
+ {{.ft = {0x300000U, 0x7FU, 1}}, -1}, /* -1.375000 */
+ {{.ft = {0x100000U, 0x7FU, 0}}, +1}, /* 1.125000 */
+ {{.ft = {0x100000U, 0x7FU, 1}}, -1}, /* -1.125000 */
+
+ {{.ft = {0x000000U, 0x7FU + 0x16U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x16U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x17U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x17U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x18U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x18U, 1}}, -16777216}, /* -16777216.000000 */
+
+ {{.ft = {0x000000U, 0x7FU + 0x1EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 0}}, +2147483648LL}, /* 2147483648.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x20U, 0}}, 4294967296}, /* 4294967296.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x20U, 1}}, -4294967296}, /* -4294967296.000000 */
+
+ /* Different mantissa patterns. */
+ {{.ft = {0x7FFFFFU, 0x96U, 0}}, +16777215}, /* 16777215.000000 */
+ {{.ft = {0x7FF000U, 0x95U, 0}}, +8386560}, /* 8386560.000000 */
+ {{.ft = {0x1555FFU, 0x8DU, 0}}, +19115}, /* 19115.000000 */
+ {{.ft = {0x7FF000U, 0x96U, 1}}, -16773120}, /* -16773120.000000 */
+ {{.ft = {0x7FFFFEU, 0x95U, 1}}, -8388607}, /* -8388607.000000 */
+ {{.ft = {0x1555FFU, 0x8DU, 1}}, -19115}, /* -19115.000000 */
+
+ /* Number greater than 2**63 exceeds long long int size and will be casted in an implementation defined manner. */
+ {{.ft = {0x000001U, 0x7FU + 0x3FU, 0}}, +9223372036854775808ULL}, /* 9223372036854775808.000000 */
+ {{.ft = {0x000001U, 0x7FU + 0x3FU, 1}}, -9223372036854775808ULL}, /* -9223372036854775808.000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.ft = {0x7FBC99U, 0x7FU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.ft = {0x7FBC99U, 0x7FU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.ft = {0x03126FU, 0x7FU + 0xFFFFFFF6U, 0}}, 0}, /* 0.001000 */
+ {{.ft = {0x03126FU, 0x7FU + 0xFFFFFFF6U, 1}}, -0}, /* -0.001000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.ft = {0x7CD6EAU, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.987654 */
+ {{.ft = {0x7CD6EAU, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.987654 */
+ {{.ft = {0x000001U, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.50000006 */
+ {{.ft = {0x000001U, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.50000006 */
+ {{.ft = {0x000000U, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.500000 */
+ {{.ft = {0x000000U, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.500000 */
+
+ /* Number greather than 1 and less than 2**23 will be rounded accordingly. */
+ {{.ft = {0x000000U, 0x7FU + 0x00U, 0}}, 1}, /* 1.0000000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x00U, 1}}, -1}, /* 1.0000000000 */
+ {{.ft = {0x00000FU, 0x7FU + 0x00U, 0}}, 1}, /* 1.000002 */
+ {{.ft = {0x00000FU, 0x7FU + 0x00U, 1}}, -1}, /* 1.000002 */
+ {{.ft = {0x000018U, 0x7FU + 0x10U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.ft = {0x000018U, 0x7FU + 0x10U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.ft = {0x000040U, 0x7FU + 0x10U, 0}}, 65537}, /* 65536.5000000 */
+ {{.ft = {0x000040U, 0x7FU + 0x10U, 1}}, -65537}, /* -65536.5000000 */
+ {{.ft = {0x00004DU, 0x7FU + 0x10U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.ft = {0x00004DU, 0x7FU + 0x10U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.ft = {0x7FFFFFU, 0x7FU + 0x16U, 0}}, 8388608}, /* 8388607.5000000000 */
+ {{.ft = {0x7FFFFFU, 0x7FU + 0x16U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.ft = {0x000005U, 0x7FU + 0x14U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.ft = {0x000005U, 0x7FU + 0x14U, 1}}, -1048577}, /* -1048576.6250000000 */
+};
+
+static const size_t n_tests_float = sizeof(tests_float) / sizeof(tests_float[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_float; i++)
+ {
+ long long int result = llroundf(tests_float[i].value.f);
+
+ if (tests_float[i].should_be == result)
+ counter++;
+ else
+ printf("llroundf test failed: value to round = %.12f result = %lld should be = %lld\n", tests_float[i].value.f, result, tests_float[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_float) ? "llroundf test failed." : "llroundf test succeded.");
+
+ return 0;
+}
diff -aprNU5 djgpp.orig/tests/cygnus/t-llroundl.c djgpp/tests/cygnus/t-llroundl.c
--- djgpp.orig/tests/cygnus/t-llroundl.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-llroundl.c 2013-10-19 12:18:28 +0100
@@ -0,0 +1,147 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _longdouble_union_t value; /* test value */
+ const long long int should_be; /* result */
+} entry_t;
+
+
+static const entry_t tests_long_double[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.ldt = {0x0U, 0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.ldt = {0x0U, 0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.ldt = {0x1U, 0x0U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.ldt = {0x1U, 0x0U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.ldt = {0x0U, 0x80000000U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.ldt = {0x0U, 0x80000000U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.ldt = {0xFFFFFFFFU, 0xFFFFFFFFU, 0x7FFEU, 0}}, -9.223372036854775808E+18}, /* Big number. */
+ {{.ldt = {0xFFFFFFFFU, 0xFFFFFFFFU, 0x7FFEU, 1}}, -9.223372036854775808E+18}, /* Big -number. */
+
+ /* Infs. */
+ {{.ldt = {0x0U, 0x80000000U, 0x7FFFU, 0}}, -9.223372036854775808E+18}, /* Inf */
+ {{.ldt = {0x0U, 0x80000000U, 0x7FFFU, 1}}, -9.223372036854775808E+18}, /* -Inf */
+
+ /* NaNs. */
+ {{.ldt = {0x1U, 0x80000000U, 0x7FFFU, 0}}, -9.223372036854775808E+18}, /* SNaN */
+ {{.ldt = {0x1U, 0x80000000U, 0x7FFFU, 1}}, -9.223372036854775808E+18}, /* -SNaN */
+ {{.ldt = {0x0U, 0xFFFFFFFFU, 0x7FFFU, 0}}, -9.223372036854775808E+18}, /* QNaN */
+ {{.ldt = {0x0U, 0xFFFFFFFFU, 0x7FFFU, 1}}, -9.223372036854775808E+18}, /* -QNaN */
+
+
+ /* Number. */
+ {{.ldt = {0x2168C000U, 0xC90FDAA2U, 0x3FFFU + 0x0001U, 0}}, +3}, /* PI */
+ {{.ldt = {0x2168C000U, 0xC90FDAA2U, 0x3FFFU + 0x0001U, 1}}, -3}, /* -PI */
+
+ {{.ldt = {0x00000000U, 0xF0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.875000 */
+ {{.ldt = {0x00000000U, 0xF0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.875000 */
+ {{.ldt = {0x00000000U, 0xD0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.625000 */
+ {{.ldt = {0x00000000U, 0xD0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.625000 */
+ {{.ldt = {0xF7A0B800U, 0xC00010C6U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.500002 */
+ {{.ldt = {0xF7A0B800U, 0xC00010C6U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.500002 */
+ {{.ldt = {0x00000000U, 0xC0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.500000 */
+ {{.ldt = {0x00000000U, 0xC0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.500000 */
+ {{.ldt = {0x085F4800U, 0xBFFFEF39U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.499998 */
+ {{.ldt = {0x085F4800U, 0xBFFFEF39U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.499998 */
+ {{.ldt = {0x00000000U, 0xB0000000U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.375000 */
+ {{.ldt = {0x00000000U, 0xB0000000U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.375000 */
+ {{.ldt = {0x00000000U, 0x90000000U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.125000 */
+ {{.ldt = {0x00000000U, 0x90000000U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.125000 */
+
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0016U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0016U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0017U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0017U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0018U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0018U, 1}}, -16777216}, /* -16777216.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 0}}, +2147483647LL}, /* 2147483647.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0020U, 0}}, +4294967296}, /* 4294967296.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0020U, 1}}, -4294967296}, /* -4294967296.000000 */
+
+ /* Number greater than 2**63 exceeds long long int size and will be casted in an implementation defined manner. */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x003FU, 0}}, -9.223372036854775808E18}, /* 9223372036854775808.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x003FU, 1}}, -9.223372036854775808E18}, /* -9223372036854775808.000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.ldt = {0xC42F9800U, 0xFFBC98FFU, 0x3FFFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.ldt = {0xC42F9800U, 0xFFBC98FFU, 0x3FFFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.ldt = {0x8D4FE000U, 0x83126E97U, 0x3FFFU + 0xFFFFFFF6U, 0}}, 0}, /* 0.0010000000 */
+ {{.ldt = {0x8D4FE000U, 0x83126E97U, 0x3FFFU + 0xFFFFFFF6U, 1}}, -0}, /* -0.0010000000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.ldt = {0xC37A4000U, 0xFCD6E47DU, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.9876540000 */
+ {{.ldt = {0xC37A4000U, 0xFCD6E47DU, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.9876540000 */
+ {{.ldt = {0xB2B29800U, 0x80000101U, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000600 */
+ {{.ldt = {0xB2B29800U, 0x80000101U, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000600 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+
+ /* Number greather than 1 and less than 2**31 will be rounded accordingly. */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0000U, 0}}, 1}, /* 1.0000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0000U, 1}}, -1}, /* -1.0000000000 */
+ {{.ldt = {0xF7A0B800U, 0x800010C6U, 0x3FFFU + 0x0000U, 0}}, 1}, /* 1.0000020000 */
+ {{.ldt = {0xF7A0B800U, 0x800010C6U, 0x3FFFU + 0x0000U, 1}}, -1}, /* -1.0000020000 */
+ {{.ldt = {0x00000000U, 0x80001800U, 0x3FFFU + 0x0010U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.ldt = {0x00000000U, 0x80001800U, 0x3FFFU + 0x0010U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.ldt = {0x00000000U, 0x80004000U, 0x3FFFU + 0x0010U, 0}}, 65537}, /* 65536.5000000000 */
+ {{.ldt = {0x00000000U, 0x80004000U, 0x3FFFU + 0x0010U, 1}}, -65537}, /* -65536.5000000000 */
+ {{.ldt = {0x00000000U, 0x80004D00U, 0x3FFFU + 0x0010U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.ldt = {0x00000000U, 0x80004D00U, 0x3FFFU + 0x0010U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.ldt = {0x00000000U, 0xFFFFFF00U, 0x3FFFU + 0x0016U, 0}}, 8388608}, /* 8388607.5000000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFF00U, 0x3FFFU + 0x0016U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.ldt = {0x00000000U, 0x80000500U, 0x3FFFU + 0x0014U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.ldt = {0x00000000U, 0x80000500U, 0x3FFFU + 0x0014U, 1}}, -1048577}, /* -1048576.6250000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 0}}, 2147483647}, /* 2147483647.0000000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 1}}, -2147483647}, /* -2147483647.0000000000 */
+
+ /* Number greather than 2**31 and less than 2**63 will be rounded accordingly. */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 0}}, 2147483648LL}, /* 2147483648.0000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 1}}, -2147483648LL}, /* -2147483648.0000000000 */
+ {{.ldt = {0xFFFFFFFEU, 0xFFFFFFFFU, 0x3FFFU + 0x003EU, 0}}, 9223372036854775807}, /* 9223372036854775807.0000000000 */
+ {{.ldt = {0xFFFFFFFEU, 0xFFFFFFFFU, 0x3FFFU + 0x003EU, 1}}, -9223372036854775807}, /* -9223372036854775807.0000000000 */
+ {{.ldt = {0x00000067U, 0x80000000U, 0x3FFFU + 0x0037U, 0}}, 36028797018963968}, /* 36028797018963968.4023437500 */
+ {{.ldt = {0x00000067U, 0x80000000U, 0x3FFFU + 0x0037U, 1}}, -36028797018963968}, /* -36028797018963968.4023437500 */
+ {{.ldt = {0x00000080U, 0x80000000U, 0x3FFFU + 0x0037U, 0}}, 36028797018963969}, /* 36028797018963968.5000000000000000 */
+ {{.ldt = {0x00000080U, 0x80000000U, 0x3FFFU + 0x0037U, 1}}, -36028797018963969}, /* -36028797018963968.5000000000000000 */
+ {{.ldt = {0x00000083U, 0x80000000U, 0x3FFFU + 0x0037U, 0}}, 36028797018963969}, /* 36028797018963968.5117187500000000 */
+ {{.ldt = {0x00000083U, 0x80000000U, 0x3FFFU + 0x0037U, 1}}, -36028797018963969}, /* -36028797018963968.5117187500000000 */
+
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x003FU, 0}}, -9223372036854775808ULL}, /* 9223372036854775808.0000000000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x003FU, 1}}, -9223372036854775808ULL}, /* 9223372036854775808.0000000000000000 */
+};
+
+static const size_t n_tests_long_double = sizeof(tests_long_double) / sizeof(tests_long_double[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_long_double; i++)
+ {
+ long long int result = llroundl(tests_long_double[i].value.ld);
+
+ if (tests_long_double[i].should_be == result)
+ counter++;
+ else
+ printf("llroundl test failed: value to round = %.12Lg result = %lld should be = %lld\n", tests_long_double[i].value.ld, result, tests_long_double[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_long_double) ? "llroundl test failed." : "llroundl test succeded.");
+
+ return 0;
+}
diff -aprNU5 djgpp.orig/tests/cygnus/t-lround.c djgpp/tests/cygnus/t-lround.c
--- djgpp.orig/tests/cygnus/t-lround.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-lround.c 2013-10-14 22:40:02 +0100
@@ -0,0 +1,180 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _double_union_t value; /* test value */
+ const long int should_be; /* result */
+} entry_t;
+
+
+static const entry_t tests_double[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.dt = {0x0U, 0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.dt = {0x0U, 0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.dt = {0x1U, 0x0U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.dt = {0x1U, 0x0U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.dt = {0x1U, 0x0U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.dt = {0x1U, 0x0U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.dt = {0xFFFFFFFFU, 0x7FFFFU, 0x7FEU, 0}}, -2147483648LL}, /* Big number. */
+ {{.dt = {0xFFFFFFFFU, 0x7FFFFU, 0x7FEU, 1}}, -2147483648LL}, /* Big -number. */
+
+ /* Infs. */
+ {{.dt = {0x0U, 0x0U, 0x7FFU, 0}}, -2147483648LL}, /* Inf */
+ {{.dt = {0x0U, 0x0U, 0x7FFU, 1}}, -2147483648LL}, /* -Inf */
+
+ /* NaNs. */
+ {{.dt = {0x1U, 0x0U, 0x7FFU, 0}}, -2147483648LL}, /* SNaN */
+ {{.dt = {0x1U, 0x0U, 0x7FFU, 1}}, -2147483648LL}, /* -SNaN */
+ {{.dt = {0x0U, 0xFFFFFU, 0x7FFU, 1}}, -2147483648LL}, /* QNaN */
+ {{.dt = {0x0U, 0xFFFFFU, 0x7FFU, 0}}, -2147483648LL}, /* -QNaN */
+
+
+ /* Number. */
+ {{.dt = {0x54442D18U, 0x921FBU, 0x3FFU + 0x001U, 0}}, +3}, /* PI */
+ {{.dt = {0x54442D18U, 0x921FBU, 0x3FFU + 0x001U, 1}}, -3}, /* -PI */
+
+ {{.dt = {0x00000000U, 0xE0000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.875000 */
+ {{.dt = {0x00000000U, 0xE0000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.875000 */
+ {{.dt = {0x00000000U, 0xA0000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.625000 */
+ {{.dt = {0x00000000U, 0xA0000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.625000 */
+ {{.dt = {0x18DEF417U, 0x80002U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.500002 */
+ {{.dt = {0x18DEF417U, 0x80002U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.500002 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x0U, 0}}, +2}, /* 1.500000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x0U, 1}}, -2}, /* -1.500000 */
+ {{.dt = {0xE7210BE9U, 0x7F9FDU, 0x3FFU + 0x0U, 0}}, +1}, /* 1.499998 */
+ {{.dt = {0xE7210BE9U, 0x7F9FDU, 0x3FFU + 0x0U, 1}}, -1}, /* -1.499998 */
+ {{.dt = {0x00000000U, 0x60000U, 0x3FFU + 0x0U, 0}}, +1}, /* 1.375000 */
+ {{.dt = {0x00000000U, 0x60000U, 0x3FFU + 0x0U, 1}}, -1}, /* -1.375000 */
+ {{.dt = {0x00000000U, 0x20000U, 0x3FFU + 0x0U, 0}}, +1}, /* 1.125000 */
+ {{.dt = {0x00000000U, 0x20000U, 0x3FFU + 0x0U, 1}}, -1}, /* -1.125000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x016U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x016U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x017U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x017U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x018U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x018U, 1}}, -16777216}, /* -16777216.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 0}}, +2147483647LL}, /* 2147483647.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x020U, 0}}, -2147483648LL}, /* 4294967296.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x020U, 1}}, -2147483648LL}, /* -4294967296.000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 0}}, -2147483648LL}, /* 2147483648.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.dt = {0xF9ADD3C1U, 0x00001U, 0x3FFU + 0x013U, 0}}, 524289}, /* 524288.9876543210 */
+ {{.dt = {0xF9ADD3C1U, 0x00001U, 0x3FFU + 0x013U, 1}}, -524289}, /* -524288.9876543210 */
+ {{.dt = {0x22C27E88U, 0x00001U, 0x3FFU + 0x013U, 0}}, 524289}, /* 524288.5678901235 */
+ {{.dt = {0x22C27E88U, 0x00001U, 0x3FFU + 0x013U, 1}}, -524289}, /* -524288.5678901235 */
+ {{.dt = {0xE9E03FDAU, 0x00000U, 0x3FFU + 0x013U, 0}}, 524288}, /* 524288.4567890123 */
+ {{.dt = {0xE9E03FDAU, 0x00000U, 0x3FFU + 0x013U, 1}}, -524288}, /* -524288.4567890123 */
+ {{.dt = {0x06522C3EU, 0x00000U, 0x3FFU + 0x013U, 0}}, 524288}, /* 524288.0123456789 */
+ {{.dt = {0x06522C3EU, 0x00000U, 0x3FFU + 0x013U, 1}}, -524288}, /*- 524288.0123456789 */
+
+ {{.dt = {0x00000000U, 0xC0000U, 0x3FFU + 0x000U, 0}}, 2}, /* 1.7500000000 */
+ {{.dt = {0x00000000U, 0xC0000U, 0x3FFU + 0x000U, 1}}, -2}, /* -1.7500000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x000U, 0}}, 2}, /* 1.5000000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0x000U, 1}}, -2}, /* -1.5000000000 */
+ {{.dt = {0x00000000U, 0x40000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.2500000000 */
+ {{.dt = {0x00000000U, 0x40000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.2500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.7500000000 */
+ {{.dt = {0x00000000U, 0x80000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.7500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.2500000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.2500000000 */
+ {{.dt = {0xF8FAB5E6U, 0x948B0U, 0x3FFU + 0xFFFFFFF9U, 0}}, 0}, /* 0.0123456789 */
+ {{.dt = {0xF8FAB5E6U, 0x948B0U, 0x3FFU + 0xFFFFFFF9U, 1}}, -0}, /* -0.0123456789 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFC01U, 0}}, 0}, /* 0.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFC01U, 1}}, -0}, /* -0.0000000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048576}, /* 1048576.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048576}, /* -1048576.0000000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048576}, /* 1048576.2500000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048576}, /* -1048576.2500000000 */
+ {{.dt = {0x80000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.5000000000 */
+ {{.dt = {0x80000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.5000000000 */
+ {{.dt = {0xC0000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.7500000000 */
+ {{.dt = {0xC0000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.7500000000 */
+
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097152}, /* 2097152.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097152}, /* -2097152.0000000000 */
+ {{.dt = {0x20000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097152}, /* 2097152.2500000000 */
+ {{.dt = {0x20000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097152}, /* -2097152.2500000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097153}, /* 2097152.5000000000 */
+ {{.dt = {0x40000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097153}, /* -2097152.5000000000 */
+ {{.dt = {0x60000000U, 0x00000U, 0x3FFU + 0x015U, 0}}, 2097153}, /* 2097152.7500000000 */
+ {{.dt = {0x60000000U, 0x00000U, 0x3FFU + 0x015U, 1}}, -2097153}, /* -2097152.7500000000 */
+
+ /* Number greater than 2**31 exceeds long int size and will be casted in an implementation defined manner. */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 0}}, -2147483648LL}, /* 2147483648.000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x01FU, 0}}, -2147483648LL}, /* -2147483648.000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.dt = {0x1FF885F3U, 0xFF793U, 0x3FFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.dt = {0x1FF885F3U, 0xFF793U, 0x3FFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.dt = {0xD2F1A9FCU, 0x0624DU, 0x3FFU + 0xFFFFFFF6U, 0}}, 0}, /* 0.0010000000 */
+ {{.dt = {0xD2F1A9FCU, 0x0624DU, 0x3FFU + 0xFFFFFFF6U, 1}}, -0}, /* -0.0010000000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.dt = {0x8FB86F48U, 0xF9ADCU, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.9876540000 */
+ {{.dt = {0x8FB86F48U, 0xF9ADCU, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.9876540000 */
+ {{.dt = {0x20365653U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000600 */
+ {{.dt = {0x20365653U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000600 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 0}}, 1}, /* -0.5000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+
+ /* Number greather than 1 and less than 2**31 will be rounded accordingly. */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000000000 */
+ {{.dt = {0x00000000U, 0x00000U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000000000 */
+ {{.dt = {0x18DEF417U, 0x00002U, 0x3FFU + 0x000U, 0}}, 1}, /* 1.0000020000 */
+ {{.dt = {0x18DEF417U, 0x00002U, 0x3FFU + 0x000U, 1}}, -1}, /* -1.0000020000 */
+ {{.dt = {0x00000000U, 0x00003U, 0x3FFU + 0x010U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.dt = {0x00000000U, 0x00003U, 0x3FFU + 0x010U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.dt = {0x00000000U, 0x00008U, 0x3FFU + 0x010U, 0}}, 65537}, /* 65536.5000000000 */
+ {{.dt = {0x00000000U, 0x00008U, 0x3FFU + 0x010U, 1}}, -65537}, /* -65536.5000000000 */
+ {{.dt = {0xA0000000U, 0x00009U, 0x3FFU + 0x010U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.dt = {0xA0000000U, 0x00009U, 0x3FFU + 0x010U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.dt = {0xE0000000U, 0xFFFFFU, 0x3FFU + 0x016U, }}, 8388608}, /* 8388607.5000000000 */
+ {{.dt = {0xE0000000U, 0xFFFFFU, 0x3FFU + 0x016U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.dt = {0xA0000000U, 0x00000U, 0x3FFU + 0x014U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.dt = {0xA0000000U, 0x00000U, 0x3FFU + 0x014U, 1}}, -1048577}, /* -1048576.6250000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 0}}, 2147483647}, /* 2147483647.0000000000 */
+ {{.dt = {0xFFC00000U, 0xFFFFFU, 0x3FFU + 0x01EU, 1}}, -2147483647}, /* -2147483647.0000000000 */
+};
+
+static const size_t n_tests_double = sizeof(tests_double) / sizeof(tests_double[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_double; i++)
+ {
+ long int result = lround(tests_double[i].value.d);
+
+ if (tests_double[i].should_be == result)
+ counter++;
+ else
+ printf("lround test failed: value to round = %.6g result = %ld should be = %ld\n", tests_double[i].value.d, result, tests_double[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_double) ? "lround test failed." : "lround test succeded.");
+
+ return 0;
+}
diff -aprNU5 djgpp.orig/tests/cygnus/t-lroundf.c djgpp/tests/cygnus/t-lroundf.c
--- djgpp.orig/tests/cygnus/t-lroundf.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-lroundf.c 2013-10-14 22:40:02 +0100
@@ -0,0 +1,137 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _float_union_t value; /* test value */
+ const long int should_be; /* result */
+} entry_t;
+
+static const entry_t tests_float[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.ft = {0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.ft = {0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.ft = {0x1U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.ft = {0x1U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.ft = {0x1U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.ft = {0x1U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.ft = {0xFFFFU, 0xFEU, 0}}, -2147483648LL}, /* Big number. */
+ {{.ft = {0xFFFFU, 0xFEU, 1}}, -2147483648LL}, /* Big -number. */
+
+ /* Infs. */
+ {{.ft = {0x0U, 0xFFU, 0}}, -2147483648LL}, /* Inf */
+ {{.ft = {0x0U, 0xFFU, 1}}, -2147483648LL}, /* -Inf */
+
+ /* NaNs. */
+ {{.ft = {0x1U, 0xFFU, 0}}, -2147483648LL}, /* SNaN */
+ {{.ft = {0x1U, 0xFFU, 1}}, -2147483648LL}, /* -SNaN */
+ {{.ft = {0x7FFFFFU, 0xFFU, 0}}, -2147483648LL}, /* QNaN */
+ {{.ft = {0x7FFFFFU, 0xFFU, 1}}, -2147483648LL}, /* -QNaN */
+
+ /* Numbers. */
+ {{.ft = {0x490FDBU, 0x80U, 0}}, +3}, /* PI */
+ {{.ft = {0x490FDBU, 0x80U, 1}}, -3}, /* -PI */
+
+ {{.ft = {0x700000U, 0x7FU, 0}}, +2}, /* 1.875000 */
+ {{.ft = {0x700000U, 0x7FU, 1}}, -2}, /* -1.875000 */
+ {{.ft = {0x500000U, 0x7FU, 0}}, +2}, /* 1.625000 */
+ {{.ft = {0x500000U, 0x7FU, 1}}, -2}, /* -1.625000 */
+ {{.ft = {0x40000FU, 0x7FU, 0}}, +2}, /* 1.500002 */
+ {{.ft = {0x40000FU, 0x7FU, 1}}, -2}, /* -1.500002 */
+ {{.ft = {0x400000U, 0x7FU, 0}}, +2}, /* 1.500000 */
+ {{.ft = {0x400000U, 0x7FU, 1}}, -2}, /* -1.500000 */
+ {{.ft = {0x3FFFF0U, 0x7FU, 0}}, +1}, /* 1.499998 */
+ {{.ft = {0x3FFFF0U, 0x7FU, 1}}, -1}, /* -1.499998 */
+ {{.ft = {0x300000U, 0x7FU, 0}}, +1}, /* 1.375000 */
+ {{.ft = {0x300000U, 0x7FU, 1}}, -1}, /* -1.375000 */
+ {{.ft = {0x100000U, 0x7FU, 0}}, +1}, /* 1.125000 */
+ {{.ft = {0x100000U, 0x7FU, 1}}, -1}, /* -1.125000 */
+
+ {{.ft = {0x000000U, 0x7FU + 0x16U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x16U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x17U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x17U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x18U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x18U, 1}}, -16777216}, /* -16777216.000000 */
+
+ {{.ft = {0x000000U, 0x7FU + 0x1EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 0}}, +2147483648LL}, /* 2147483648.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x20U, 0}}, -2147483648LL}, /* 4294967296.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x20U, 1}}, -2147483648LL}, /* -4294967296.000000 */
+
+ /* Different mantissa patterns. */
+ {{.ft = {0x7FFFFFU, 0x96U, 0}}, +16777215}, /* 16777215.000000 */
+ {{.ft = {0x7FF000U, 0x95U, 0}}, +8386560}, /* 8386560.000000 */
+ {{.ft = {0x1555FFU, 0x8DU, 0}}, +19115}, /* 19115.000000 */
+ {{.ft = {0x7FF000U, 0x96U, 1}}, -16773120}, /* -16773120.000000 */
+ {{.ft = {0x7FFFFEU, 0x95U, 1}}, -8388607}, /* -8388607.000000 */
+ {{.ft = {0x1555FFU, 0x8DU, 1}}, -19115}, /* -19115.000000 */
+
+ /* Number greater than 2**31 exceeds long int size and will be casted in an implementation defined manner. */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 0}}, +2147483648LL}, /* 2147483648.000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x1FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.ft = {0x7FBC99U, 0x7FU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.ft = {0x7FBC99U, 0x7FU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.ft = {0x03126FU, 0x7FU + 0xFFFFFFF6U, 0}}, 0}, /* 0.001000 */
+ {{.ft = {0x03126FU, 0x7FU + 0xFFFFFFF6U, 1}}, -0}, /* -0.001000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.ft = {0x7CD6EAU, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.987654 */
+ {{.ft = {0x7CD6EAU, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.987654 */
+ {{.ft = {0x000001U, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.50000006 */
+ {{.ft = {0x000001U, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.50000006 */
+ {{.ft = {0x000000U, 0x7FU + 0xFFFFFFFFU, 0}}, 1}, /* 0.500000 */
+ {{.ft = {0x000000U, 0x7FU + 0xFFFFFFFFU, 1}}, -1}, /* -0.500000 */
+
+ /* Number greather than 1 and less than 2**23 will be rounded accordingly. */
+ {{.ft = {0x000000U, 0x7FU + 0x00U, 0}}, 1}, /* 1.0000000000 */
+ {{.ft = {0x000000U, 0x7FU + 0x00U, 1}}, -1}, /* 1.0000000000 */
+ {{.ft = {0x00000FU, 0x7FU + 0x00U, 0}}, 1}, /* 1.000002 */
+ {{.ft = {0x00000FU, 0x7FU + 0x00U, 1}}, -1}, /* 1.000002 */
+ {{.ft = {0x000018U, 0x7FU + 0x10U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.ft = {0x000018U, 0x7FU + 0x10U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.ft = {0x000040U, 0x7FU + 0x10U, 0}}, 65537}, /* 65536.5000000 */
+ {{.ft = {0x000040U, 0x7FU + 0x10U, 1}}, -65537}, /* -65536.5000000 */
+ {{.ft = {0x00004DU, 0x7FU + 0x10U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.ft = {0x00004DU, 0x7FU + 0x10U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.ft = {0x7FFFFFU, 0x7FU + 0x16U, 0}}, 8388608}, /* 8388607.5000000000 */
+ {{.ft = {0x7FFFFFU, 0x7FU + 0x16U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.ft = {0x000005U, 0x7FU + 0x14U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.ft = {0x000005U, 0x7FU + 0x14U, 1}}, -1048577}, /* -1048576.6250000000 */
+};
+
+static const size_t n_tests_float = sizeof(tests_float) / sizeof(tests_float[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_float; i++)
+ {
+ long int result = lroundf(tests_float[i].value.f);
+
+ if (tests_float[i].should_be == result)
+ counter++;
+ else
+ printf("lroundf test failed: value to round = %.12f result = %ld should be = %ld\n", tests_float[i].value.f, result, tests_float[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_float) ? "lroundf test failed." : "lroundf test succeded.");
+
+ return 0;
+}
diff -aprNU5 djgpp.orig/tests/cygnus/t-lroundl.c djgpp/tests/cygnus/t-lroundl.c
--- djgpp.orig/tests/cygnus/t-lroundl.c 1970-01-01 01:00:00 +0100
+++ djgpp/tests/cygnus/t-lroundl.c 2013-10-14 22:40:02 +0100
@@ -0,0 +1,131 @@
+/* Copyright (C) 2013 DJ Delorie, see COPYING.DJ for details */
+
+
+#include <stdio.h>
+#include <math.h>
+#include <libc/ieee.h>
+
+typedef struct {
+ const _longdouble_union_t value; /* test value */
+ const long int should_be; /* result */
+} entry_t;
+
+
+static const entry_t tests_long_double[] =
+{
+ /* test value */
+ /* value should be */
+
+ /* Zeros. */
+ {{.ldt = {0x0U, 0x0U, 0x0U, 0}}, 0}, /* 0.0 */
+ {{.ldt = {0x0U, 0x0U, 0x0U, 1}}, 0}, /* -0.0 */
+
+ /* Subnormals aka denormals. */
+ {{.ldt = {0x1U, 0x0U, 0x0U, 0}}, 0}, /* Very small number. */
+ {{.ldt = {0x1U, 0x0U, 0x0U, 1}}, 0}, /* Very small -number. */
+
+ /* Normals. */
+ {{.ldt = {0x0U, 0x80000000U, 0x1U, 0}}, 0}, /* Small number. */
+ {{.ldt = {0x0U, 0x80000000U, 0x1U, 1}}, 0}, /* Small -number. */
+ {{.ldt = {0xFFFFFFFFU, 0xFFFFFFFFU, 0x7FFEU, 0}}, -2147483648LL}, /* Big number. */
+ {{.ldt = {0xFFFFFFFFU, 0xFFFFFFFFU, 0x7FFEU, 1}}, -2147483648LL}, /* Big -number. */
+
+ /* Infs. */
+ {{.ldt = {0x0U, 0x80000000U, 0x7FFFU, 0}}, -2147483648LL}, /* Inf */
+ {{.ldt = {0x0U, 0x80000000U, 0x7FFFU, 1}}, -2147483648LL}, /* -Inf */
+
+ /* NaNs. */
+ {{.ldt = {0x1U, 0x80000000U, 0x7FFFU, 0}}, -2147483648LL}, /* SNaN */
+ {{.ldt = {0x1U, 0x80000000U, 0x7FFFU, 1}}, -2147483648LL}, /* -SNaN */
+ {{.ldt = {0x0U, 0xFFFFFFFFU, 0x7FFFU, 0}}, -2147483648LL}, /* QNaN */
+ {{.ldt = {0x0U, 0xFFFFFFFFU, 0x7FFFU, 1}}, -2147483648LL}, /* -QNaN */
+
+ /* Number. */
+ {{.ldt = {0x2168C000U, 0xC90FDAA2U, 0x3FFFU + 0x0001U, 0}}, +3}, /* PI */
+ {{.ldt = {0x2168C000U, 0xC90FDAA2U, 0x3FFFU + 0x0001U, 1}}, -3}, /* -PI */
+
+ {{.ldt = {0x00000000U, 0xF0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.875000 */
+ {{.ldt = {0x00000000U, 0xF0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.875000 */
+ {{.ldt = {0x00000000U, 0xD0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.625000 */
+ {{.ldt = {0x00000000U, 0xD0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.625000 */
+ {{.ldt = {0xF7A0B800U, 0xC00010C6U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.500002 */
+ {{.ldt = {0xF7A0B800U, 0xC00010C6U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.500002 */
+ {{.ldt = {0x00000000U, 0xC0000000U, 0x3FFFU + 0x0U, 0}}, +2}, /* 1.500000 */
+ {{.ldt = {0x00000000U, 0xC0000000U, 0x3FFFU + 0x0U, 1}}, -2}, /* -1.500000 */
+ {{.ldt = {0x085F4800U, 0xBFFFEF39U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.499998 */
+ {{.ldt = {0x085F4800U, 0xBFFFEF39U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.499998 */
+ {{.ldt = {0x00000000U, 0xB0000000U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.375000 */
+ {{.ldt = {0x00000000U, 0xB0000000U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.375000 */
+ {{.ldt = {0x00000000U, 0x90000000U, 0x3FFFU + 0x0U, 0}}, +1}, /* 1.125000 */
+ {{.ldt = {0x00000000U, 0x90000000U, 0x3FFFU + 0x0U, 1}}, -1}, /* -1.125000 */
+
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0016U, 0}}, +4194304}, /* 4194304.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0016U, 1}}, -4194304}, /* -4194304.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0017U, 0}}, +8388608}, /* 8388608.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0017U, 1}}, -8388608}, /* -8388608.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0018U, 0}}, +16777216}, /* 16777216.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0018U, 1}}, -16777216}, /* -16777216.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001EU, 0}}, +1073741824}, /* 1073741824.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001EU, 1}}, -1073741824}, /* -1073741824.000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 0}}, +2147483647LL}, /* 2147483647.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 1}}, -2147483648LL}, /* -2147483648.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0020U, 0}}, -2147483648LL}, /* 4294967296.000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0020U, 1}}, -2147483648LL}, /* -4294967296.000000 */
+
+ /* Number greater than 2**31 exceeds long int size and will be casted in an implementation defined manner. */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 0}}, -2147483648LL}, /* 2147483648.0000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x001FU, 1}}, -2147483648LL}, /* -2147483648.0000000000 */
+
+ /* Number less than 0.5 will be truncated to 0. */
+ {{.ldt = {0xC42F9800U, 0xFFBC98FFU, 0x3FFFU + 0xFFFFFFFEU, 0}}, 0}, /* 0.4994857609 */
+ {{.ldt = {0xC42F9800U, 0xFFBC98FFU, 0x3FFFU + 0xFFFFFFFEU, 1}}, -0}, /* -0.4994857609 */
+ {{.ldt = {0x8D4FE000U, 0x83126E97U, 0x3FFFU + 0xFFFFFFF6U, 0}}, 0}, /* 0.0010000000 */
+ {{.ldt = {0x8D4FE000U, 0x83126E97U, 0x3FFFU + 0xFFFFFFF6U, 1}}, -0}, /* -0.0010000000 */
+
+ /* Number greater than 0.5 and less than 1 will be rounded to 1. */
+ {{.ldt = {0xC37A4000U, 0xFCD6E47DU, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.9876540000 */
+ {{.ldt = {0xC37A4000U, 0xFCD6E47DU, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.9876540000 */
+ {{.ldt = {0xB2B29800U, 0x80000101U, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000600 */
+ {{.ldt = {0xB2B29800U, 0x80000101U, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000600 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0xFFFFFFFFU, 0}}, 1}, /* 0.5000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0xFFFFFFFFU, 1}}, -1}, /* -0.5000000000 */
+
+ /* Number greather than 1 and less than 2**31 will be rounded accordingly. */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0000U, 0}}, 1}, /* 1.0000000000 */
+ {{.ldt = {0x00000000U, 0x80000000U, 0x3FFFU + 0x0000U, 1}}, -1}, /* -1.0000000000 */
+ {{.ldt = {0xF7A0B800U, 0x800010C6U, 0x3FFFU + 0x0000U, 0}}, 1}, /* 1.0000020000 */
+ {{.ldt = {0xF7A0B800U, 0x800010C6U, 0x3FFFU + 0x0000U, 1}}, -1}, /* -1.0000020000 */
+ {{.ldt = {0x00000000U, 0x80001800U, 0x3FFFU + 0x0010U, 0}}, 65536}, /* 65536.1875000000 */
+ {{.ldt = {0x00000000U, 0x80001800U, 0x3FFFU + 0x0010U, 1}}, -65536}, /* -65536.1875000000 */
+ {{.ldt = {0x00000000U, 0x80004000U, 0x3FFFU + 0x0010U, 0}}, 65537}, /* 65536.5000000000 */
+ {{.ldt = {0x00000000U, 0x80004000U, 0x3FFFU + 0x0010U, 1}}, -65537}, /* -65536.5000000000 */
+ {{.ldt = {0x00000000U, 0x80004D00U, 0x3FFFU + 0x0010U, 0}}, 65537}, /* 65536.6015625000 */
+ {{.ldt = {0x00000000U, 0x80004D00U, 0x3FFFU + 0x0010U, 1}}, -65537}, /* -65536.6015625000 */
+ {{.ldt = {0x00000000U, 0xFFFFFF00U, 0x3FFFU + 0x0016U, 0}}, 8388608}, /* 8388607.5000000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFF00U, 0x3FFFU + 0x0016U, 1}}, -8388608}, /* -8388607.5000000000 */
+ {{.ldt = {0x00000000U, 0x80000500U, 0x3FFFU + 0x0014U, 0}}, 1048577}, /* 1048576.6250000000 */
+ {{.ldt = {0x00000000U, 0x80000500U, 0x3FFFU + 0x0014U, 1}}, -1048577}, /* -1048576.6250000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 0}}, 2147483647}, /* 2147483647.0000000000 */
+ {{.ldt = {0x00000000U, 0xFFFFFFFEU, 0x3FFFU + 0x001EU, 1}}, -2147483647}, /* -2147483647.0000000000 */
+};
+
+static const size_t n_tests_long_double = sizeof(tests_long_double) / sizeof(tests_long_double[0]);
+
+
+int main(void)
+{
+ unsigned int i, counter;
+
+ for (counter = i = 0; i < n_tests_long_double; i++)
+ {
+ long int result = lroundl(tests_long_double[i].value.ld);
+
+ if (tests_long_double[i].should_be == result)
+ counter++;
+ else
+ printf("lroundl test failed: value to round = %.12Lf result = %ld should be = %ld\n", tests_long_double[i].value.ld, result, tests_long_double[i].should_be);
+ }
+ printf("%s\n", (counter < n_tests_long_double) ? "lroundl test failed." : "lroundl test succeded.");
+
+ return 0;
+}
- Raw text -