Files
mstorsjo-llvm-mingw/test/crt-test.c
T
Martin Storsjö 61ddd976fe Revert "test: crt-test: Skip a test for cos(-NAN) == -NAN when optimizing"
This reverts commit 1d8c107d49.

After 49e14021be, when using F(-NAN)
instead of -F(NAN) to hide the negation, this test passes correctly.
I hadn't realized that GCC and Clang optimized cos(-x) = cos(x), when
I made 1d8c107d49 (and at the time, the
F() macro didn't work properly with -NAN on all toolchains).
2026-04-24 15:39:58 +03:00

2390 lines
91 KiB
C

/*
* Copyright (c) 2018 Martin Storsjo
*
* This file is part of llvm-mingw.
*
* Permission to use, copy, modify, and/or distribute this software for any
* purpose with or without fee is hereby granted, provided that the above
* copyright notice and this permission notice appear in all copies.
*
* THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
* WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
* MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
* ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
* WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
* ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
* OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
*/
#if defined(__linux__) && !defined(_GNU_SOURCE)
#define _GNU_SOURCE
#endif
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <stdlib.h>
#include <fenv.h>
#include <inttypes.h>
#include <stdarg.h>
#include <errno.h>
#include <limits.h>
#include <float.h>
#include <wchar.h>
#ifdef _WIN32
#include <windows.h>
#endif
#if defined(__MINGW32__) && defined(__MSVCRT_VERSION__) && __MSVCRT_VERSION__ < 0xE00
#define MSVCRT_DLL
#endif
#if defined(MSVCRT_DLL) && defined(__USE_MINGW_ANSI_STDIO) && __USE_MINGW_ANSI_STDIO == 0
#define MSVCRT_DLL_NOANSI
#endif
#ifndef _WIN32
extern char **environ;
#endif
#ifdef _WIN32
static void invalid_parameter(const wchar_t *expression, const wchar_t *function, const wchar_t *file, unsigned int line, uintptr_t pReserved) {
}
#endif
int tests = 0, fails = 0;
const char *context = "";
#define TEST(x) do { \
tests++; \
if (!(x)) { \
fails++; \
printf("%s:%d: %s\"%s\" failed\n", __FILE__, __LINE__, context, #x); \
} \
} while (0)
#define TEST_STR(x, expect) do { \
tests++; \
if (strcmp((x), (expect))) { \
fails++; \
printf("%s:%d: %sexpected \"%s\", got \"%s\"\n", __FILE__, __LINE__, context, (expect), (x)); \
} \
} while (0)
#define TEST_XSTR(strcmp, type, x, expect, fmt) do { \
const type *str = (x); \
tests++; \
if (!str || strcmp(str, (expect))) { \
fails++; \
printf("%s:%d: %sexpected \"" fmt "\", got \"" fmt "\"\n", __FILE__, __LINE__, context, (expect), str); \
} \
} while (0)
#define TEST_WSTR(x, expect) TEST_XSTR(wcscmp, wchar_t, x, expect, "%ls")
#define TEST_FLT(x, expect) do { \
tests++; \
if ((x) != (expect)) { \
fails++; \
printf("%s:%d: %s%s failed, expected %f, got %f\n", __FILE__, __LINE__, context, #x, (double)(expect), (double)(x)); \
} \
} while (0)
#define TEST_FLT_EXPR(x, expr) do { \
tests++; \
if (!(expr)) { \
fails++; \
printf("%s:%d: %s%s failed, got %f\n", __FILE__, __LINE__, context, #expr, (double)(x)); \
} \
} while (0)
#define TEST_FLT_NAN_ANY(x) do { \
tests++; \
if (!isnan(x)) { \
fails++; \
printf("%s:%d: %s%s failed, expected any NAN, got %f\n", __FILE__, __LINE__, context, #x, (double)(x)); \
} \
} while (0)
#if defined(__MINGW32__) || defined(TEST_NAN_PRESERVATION)
// Only check for NAN sign presevation on mingw (where the tests currently
// pass on all platforms and we want to keep checking that) or if
// TEST_NAN_PRESERVATION is defined. This aspect does fail for some
// functions on some libcs - but ignore those failures as it's a rather
// overly pedantic test.
#define NAN_MATCHING_SIGNS(a, b) (!!signbit(a) == !!signbit(b))
#else
#define NAN_MATCHING_SIGNS(a, b) (1)
#endif
#define TEST_FLT_NAN(x, expect) do { \
tests++; \
long double val = (x); \
long double expval = (expect); \
if (!isnan(val) || !NAN_MATCHING_SIGNS(val, expval)) { \
fails++; \
printf("%s:%d: %s%s failed, expected %f (sign %d), got %f (sign %d)\n", __FILE__, __LINE__, context, #x, (double)expval, signbit(expval), (double)val, signbit(val)); \
} \
} while (0)
#define TEST_FLT_ACCURACY(x, expect, accuracy) do { \
long double val = (x); \
long double diff = fabsl(val - (expect)); \
tests++; \
if (diff <= (accuracy)) { \
/* All ok, not NAN */ \
} else { \
fails++; \
printf("%s:%d: %s%s failed, expected %f, got %f (diff %f > %f)\n", __FILE__, __LINE__, context, #x, (double)(expect), (double)val, (double)diff, (double)(accuracy)); \
} \
} while (0)
#define TEST_FLT_SIGN(x, expect) do { \
tests++; \
long double val = (x); \
long double expval = (expect); \
if (val != expval || !!signbit(val) != !!signbit(expval)) { \
fails++; \
printf("%s:%d: %s%s failed, expected %f, got %f\n", __FILE__, __LINE__, context, #x, (double)expval, (double)val); \
} \
} while (0)
#define TEST_INT(x, expect) do { \
tests++; \
if ((x) != (expect)) { \
fails++; \
printf("%s:%d: %s%s failed, expected %lld, got %lld\n", __FILE__, __LINE__, context, #x, (long long)(expect), (long long)(x)); \
} \
} while (0)
#define TEST_PTR(x, expect) do { \
tests++; \
if ((x) != (expect)) { \
fails++; \
printf("%s:%d: %s%s failed, expected %p, got %p\n", __FILE__, __LINE__, context, #x, (expect), (x)); \
} \
} while (0)
char char_wrap_impl(char c) {
return c;
}
double double_wrap_impl(double f) {
return f;
}
long long_wrap_impl(long l) {
return l;
}
unsigned long ulong_wrap_impl(unsigned long u) {
return u;
}
long long longlong_wrap_impl(long long l) {
return l;
}
unsigned long long ulonglong_wrap_impl(unsigned long long u) {
return u;
}
const char *str_wrap_impl(const char *str) {
return str;
}
char (*char_wrap)(char c) = char_wrap_impl;
double (*double_wrap)(double f) = double_wrap_impl;
long (*long_wrap)(long l) = long_wrap_impl;
unsigned long (*ulong_wrap)(unsigned long u) = ulong_wrap_impl;
long long (*longlong_wrap)(long long l) = longlong_wrap_impl;
unsigned long long (*ulonglong_wrap)(unsigned long long l) = ulonglong_wrap_impl;
const char *(*str_wrap)(const char *str) = str_wrap_impl;
#define C(x) char_wrap(x)
#define F(x) double_wrap(x)
#define L(x) long_wrap(x)
#define UL(x) ulong_wrap(x)
#define LL(x) longlong_wrap(x)
#define ULL(x) ulonglong_wrap(x ## ULL)
#define S(x) str_wrap(x)
int vsscanf_wrap(const char *str, const char *fmt, ...) {
va_list ap;
int ret;
va_start(ap, fmt);
ret = vsscanf(str, fmt, ap);
va_end(ap);
return ret;
}
int vswprintf_wrap(wchar_t *buf, size_t n, const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = vswprintf(buf, n, fmt, ap);
va_end(ap);
return ret;
}
#ifdef _WIN32
#ifdef __cplusplus
int vswprintf_overload_wrap(wchar_t *buf, const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = vswprintf(buf, fmt, ap);
va_end(ap);
return ret;
}
#endif
int _vscprintf_wrap(const char *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = _vscprintf(fmt, ap);
va_end(ap);
return ret;
}
int _vswprintf_wrap(wchar_t *buf, const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = _vswprintf(buf, fmt, ap);
va_end(ap);
return ret;
}
int _vsnwprintf_wrap(wchar_t *buf, size_t n, const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = _vsnwprintf(buf, n, fmt, ap);
va_end(ap);
return ret;
}
#ifdef __MINGW32__
int vsnwprintf_wrap(wchar_t *buf, size_t n, const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = vsnwprintf(buf, n, fmt, ap);
va_end(ap);
return ret;
}
#endif
int _vscwprintf_wrap(const wchar_t *fmt, ...) {
va_list ap;
va_start(ap, fmt);
int ret = _vscwprintf(fmt, ap);
va_end(ap);
return ret;
}
#endif
double int_to_double(uint64_t i) {
union {
uint64_t i;
double d;
} u;
u.i = i;
return u.d;
}
void fill_wchar(wchar_t *buf, wchar_t c, size_t n) {
for (size_t i = 0; i < n; i++)
buf[i] = c;
}
#define ARRAYLEN(x) (sizeof(x)/sizeof(*x))
void test_strings() {
char buf[200];
int i;
uint64_t myconst = 0xbaadf00dcafe;
void *retptr;
memset(buf, '#', sizeof(buf));
retptr = memcpy(buf, S("foo"), L(4));
TEST(retptr == buf);
TEST_STR(buf, "foo");
TEST_INT(buf[5], '#');
#if defined(__GLIBC__) || defined(__MINGW32__)
memset(buf, '#', sizeof(buf));
retptr = mempcpy(buf, S("foo"), L(4));
TEST(retptr == buf + 4);
TEST_STR(buf, "foo");
TEST_INT(buf[5], '#');
#endif
memset(buf, '#', sizeof(buf));
memcpy(buf, "foobar", 7);
retptr = memmove(buf + 2, S(buf), L(3));
TEST(retptr == buf + 2);
TEST_STR(buf, "fofoor");
TEST_INT(buf[8], '#');
memset(buf, '#', sizeof(buf));
memcpy(buf, "foobar", 7);
retptr = memmove(buf, S(buf + 2), L(3));
TEST(retptr == buf);
TEST_STR(buf, "obabar");
TEST_INT(buf[8], '#');
retptr = memset(buf, C('#'), sizeof(buf) + L(0));
TEST(retptr == buf);
TEST_INT(buf[0], '#');
TEST_INT(buf[sizeof(buf)-1], '#');
memset(buf, '#', sizeof(buf));
retptr = strcpy(buf, S("foo"));
TEST(retptr == buf);
TEST_STR(buf, "foo");
TEST_INT(buf[5], '#');
memset(buf, '#', sizeof(buf));
retptr = strncpy(buf, S("foobar"), L(3));
TEST(retptr == buf);
TEST_INT(buf[3], '#');
buf[3] = '\0';
TEST_STR(buf, "foo");
memset(buf, '#', sizeof(buf));
retptr = strncpy(buf, S("foobar"), sizeof(buf) + L(0));
TEST(retptr == buf);
TEST_STR(buf, "foobar");
TEST_INT(buf[sizeof(buf)-1], '\0');
memset(buf, '#', sizeof(buf));
strcpy(buf, "foo");
retptr = strcat(buf, S("bar"));
TEST(retptr == buf);
TEST_STR(buf, "foobar");
memset(buf, '#', sizeof(buf));
strcpy(buf, "foo");
retptr = strncat(buf, S("bar"), L(5));
TEST(retptr == buf);
TEST_STR(buf, "foobar");
memset(buf, '#', sizeof(buf));
strcpy(buf, "foo");
retptr = strncat(buf, S("bar"), L(2));
TEST(retptr == buf);
TEST_STR(buf, "fooba");
memset(buf, '#', sizeof(buf));
snprintf(buf, sizeof(buf), "%f", 3.141592654);
#ifndef MSVCRT_DLL_NOANSI
TEST_STR(buf, "3.141593");
#endif
TEST_INT(buf[sizeof(buf)-1], '#');
snprintf(buf, sizeof(buf), "%e", 42.0);
#ifndef MSVCRT_DLL_NOANSI
TEST_STR(buf, "4.200000e+01");
#endif
snprintf(buf, sizeof(buf), "%a", 42.0);
// Different implementations of printf differ in formatting of %a
// (with differing number of trailing zeros). Additionally, with
// __USE_MINGW_ANSI_STDIO, this produces 0xa.8p+2 instead of 0x1.5p+5,
// but it at least gets parsed back to the right value.
TEST_FLT(strtod(buf, NULL), 42.0);
snprintf(buf, sizeof(buf), "%g", 42.0);
TEST_STR(buf, "42");
snprintf(buf, sizeof(buf), "%g", 0.000061035156250);
#ifndef MSVCRT_DLL_NOANSI
TEST_STR(buf, "6.10352e-05");
const char formats[] = "feagFEAG";
for (int i = 0; formats[i]; i++) {
char formatbuf[3] = { '%', formats[i], '\0' };
snprintf(buf, sizeof(buf), formatbuf, INFINITY);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "INF");
else
TEST_STR(buf, "inf");
snprintf(buf, sizeof(buf), formatbuf, -INFINITY);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "-INF");
else
TEST_STR(buf, "-inf");
#ifndef _MSC_VER
// MSVC prints these differently by default, printing e.g. "-nan(ind)"
snprintf(buf, sizeof(buf), formatbuf, NAN);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "NAN");
else
TEST_STR(buf, "nan");
#endif
}
#endif /* !MSVCRT_DLL_NOANSI */
#if !defined(__MINGW32__) || (!defined(__i386__) && !defined(__x86_64__)) || (defined(__USE_MINGW_ANSI_STDIO) && __USE_MINGW_ANSI_STDIO)
// On mingw, long double formatting on x86 only works if __USE_MINGW_ANSI_STDIO is defined.
long double print_val_ld = 42.0L;
snprintf(buf, sizeof(buf), "%Lf", print_val_ld);
TEST_STR(buf, "42.000000");
snprintf(buf, sizeof(buf), "%Le", print_val_ld);
#ifndef MSVCRT_DLL_NOANSI
TEST_STR(buf, "4.200000e+01");
#endif
snprintf(buf, sizeof(buf), "%La", print_val_ld);
TEST_FLT(strtod(buf, NULL), 42.0);
snprintf(buf, sizeof(buf), "%Lg", print_val_ld);
TEST_STR(buf, "42");
#ifndef MSVCRT_DLL_NOANSI
for (int i = 0; formats[i]; i++) {
long double inf = INFINITY;
long double nan = NAN;
char formatbuf[4] = { '%', 'L', formats[i], '\0' };
snprintf(buf, sizeof(buf), formatbuf, inf);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "INF");
else
TEST_STR(buf, "inf");
snprintf(buf, sizeof(buf), formatbuf, -inf);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "-INF");
else
TEST_STR(buf, "-inf");
#ifndef _MSC_VER
// MSVC prints these differently by default, printing e.g. "-nan(ind)"
snprintf(buf, sizeof(buf), formatbuf, nan);
if (formats[i] >= 'A' && formats[i] <= 'Z')
TEST_STR(buf, "NAN");
else
TEST_STR(buf, "nan");
#endif
}
#endif /* !MSVCRT_DLL_NOANSI */
#endif /* long double tests */
snprintf(buf, sizeof(buf), "%" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64" %" PRIx64, myconst + 0, myconst + 1, myconst + 2, myconst + 3, myconst + 4, myconst + 5, myconst + 6, myconst + 7, myconst + 8, myconst + 9);
TEST_STR(buf, "baadf00dcafe baadf00dcaff baadf00dcb00 baadf00dcb01 baadf00dcb02 baadf00dcb03 baadf00dcb04 baadf00dcb05 baadf00dcb06 baadf00dcb07");
char fmt[10] = { '%', '+', '0', '5', 'd', '\0' };
snprintf(buf, sizeof(buf), fmt, 42);
TEST_STR(buf, "+0042");
fmt[4] = 'u';
snprintf(buf, sizeof(buf), fmt, 42);
TEST_STR(buf, "00042");
fmt[4] = 'o';
snprintf(buf, sizeof(buf), fmt, 42);
TEST_STR(buf, "00052");
fmt[4] = 'x';
snprintf(buf, sizeof(buf), fmt, 42);
TEST_STR(buf, "0002a");
#ifdef _WIN32
TEST_INT(_scprintf("%d", 12345), 5);
TEST_INT(_vscprintf_wrap("%d", 12345), 5);
#endif
wchar_t wbuf[200], wbuf2[5];
int ret;
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = swprintf(wbuf, ARRAYLEN(wbuf), L"%d", 42);
TEST_INT(ret, 2);
TEST_WSTR(wbuf, L"42");
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = vswprintf_wrap(wbuf, ARRAYLEN(wbuf), L"%d", 42);
TEST_INT(ret, 2);
TEST_WSTR(wbuf, L"42");
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = swprintf(wbuf2, ARRAYLEN(wbuf2), L"%d", 12345);
TEST_INT(ret, -1);
#ifndef __GLIBC__
// On overflow, glibc doesn't seem to null terminate the array.
// https://sourceware.org/bugzilla/show_bug.cgi?id=27857, fixed in
// glibc 2.37.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '\0');
#endif
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
#ifdef _WIN32
#ifdef __cplusplus
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = swprintf(wbuf, L"%d", 42);
TEST_WSTR(wbuf, L"42");
TEST_INT(ret, 2);
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = vswprintf_overload_wrap(wbuf, L"%d", 42);
TEST_WSTR(wbuf, L"42");
TEST_INT(ret, 2);
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
#endif
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = _swprintf(wbuf, L"%d", 42);
TEST_WSTR(wbuf, L"42");
TEST_INT(ret, 2);
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
fill_wchar(wbuf, '#', ARRAYLEN(wbuf));
ret = _vswprintf_wrap(wbuf, L"%d", 42);
TEST_WSTR(wbuf, L"42");
TEST_INT(ret, 2);
TEST_INT(wbuf[ARRAYLEN(wbuf)-1], '#');
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = _snwprintf(wbuf2, ARRAYLEN(wbuf2), L"%d", 12345);
TEST_INT(ret, 5);
// On overflow, _snwprintf doesn't null terminate.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '5');
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = _snwprintf(wbuf2, ARRAYLEN(wbuf2), L"%d", 123456);
TEST_INT(ret, -1);
// On overflow, _snwprintf doesn't null terminate.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '5');
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = _vsnwprintf_wrap(wbuf2, ARRAYLEN(wbuf2), L"%d", 12345);
TEST_INT(ret, 5);
// On overflow, _vsnwprintf doesn't null terminate.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '5');
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
#ifdef __MINGW32__
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = snwprintf(wbuf2, ARRAYLEN(wbuf2), L"%d", 123456);
TEST_INT(ret, 6);
// On overflow, snwprintf null terminates.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '\0');
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
TEST_INT(snwprintf(NULL, 0, L"%d", 123456), 6);
fill_wchar(wbuf2, '#', ARRAYLEN(wbuf2));
ret = vsnwprintf_wrap(wbuf2, ARRAYLEN(wbuf2), L"%d", 123456);
TEST_INT(ret, 6);
// On overflow, vsnwprintf null terminates.
TEST_INT(wbuf2[ARRAYLEN(wbuf2)-1], '\0');
wbuf2[ARRAYLEN(wbuf2)-1] = '\0';
TEST_WSTR(wbuf2, L"1234");
TEST_INT(vsnwprintf_wrap(NULL, 0, L"%d", 123456), 6);
#endif
TEST_INT(_scwprintf(L"%d", 123456), 6);
TEST_INT(_vscwprintf_wrap(L"%d", 123456), 6);
#endif
uint64_t val0, val1, val2, val3, val4, val5, val6, val7, val8, val9;
if (sscanf("baadf00dcafe baadf00dcaff baadf00dcb00 baadf00dcb01 baadf00dcb02 baadf00dcb03 baadf00dcb04 baadf00dcb05 baadf00dcb06 baadf00dcb07", "%" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64, &val0, &val1, &val2, &val3, &val4, &val5, &val6, &val7, &val8, &val9) != 10) {
fails++;
printf("sscanf failed\n");
} else {
int64_t diff = 0;
diff += llabs((int64_t)(val0 - 0 - myconst));
diff += llabs((int64_t)(val1 - 1 - myconst));
diff += llabs((int64_t)(val2 - 2 - myconst));
diff += llabs((int64_t)(val3 - 3 - myconst));
diff += llabs((int64_t)(val4 - 4 - myconst));
diff += llabs((int64_t)(val5 - 5 - myconst));
diff += llabs((int64_t)(val6 - 6 - myconst));
diff += llabs((int64_t)(val7 - 7 - myconst));
diff += llabs((int64_t)(val8 - 8 - myconst));
diff += llabs((int64_t)(val9 - 9 - myconst));
if (diff != 0) {
fails++;
printf("sscanf output failed\n");
}
}
tests++;
val0 = val1 = val2 = val3 = val4 = val5 = val6 = val7 = val8 = val9 = 0xff;
if (vsscanf_wrap("baadf00dcafe baadf00dcaff baadf00dcb00 baadf00dcb01 baadf00dcb02 baadf00dcb03 baadf00dcb04 baadf00dcb05 baadf00dcb06 baadf00dcb07", "%" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64" %" SCNx64, &val0, &val1, &val2, &val3, &val4, &val5, &val6, &val7, &val8, &val9) != 10) {
fails++;
printf("vsscanf failed\n");
} else {
int64_t diff = 0;
diff += llabs((int64_t)(val0 - 0 - myconst));
diff += llabs((int64_t)(val1 - 1 - myconst));
diff += llabs((int64_t)(val2 - 2 - myconst));
diff += llabs((int64_t)(val3 - 3 - myconst));
diff += llabs((int64_t)(val4 - 4 - myconst));
diff += llabs((int64_t)(val5 - 5 - myconst));
diff += llabs((int64_t)(val6 - 6 - myconst));
diff += llabs((int64_t)(val7 - 7 - myconst));
diff += llabs((int64_t)(val8 - 8 - myconst));
diff += llabs((int64_t)(val9 - 9 - myconst));
if (diff != 0) {
fails++;
printf("vsscanf output failed\n");
}
}
tests++;
float val_f;
double val_d;
if (sscanf("0.8 0.8", "%f %lf", &val_f, &val_d) != 2) {
fails++;
printf("sscanf for floats failed\n");
} else {
double diff = 0;
diff += fabs(val_f - 0.8);
diff += fabs(val_d - 0.8);
// Testing for the success case, to make NaNs trigger the failure case
if (diff < 0.0001) {
// All ok
} else {
fails++;
printf("sscanf output for floats failed, %f %f - diff %f\n", val_f, val_d, diff);
}
}
tests++;
#if !defined(__MINGW32__) || (!defined(__i386__) && !defined(__x86_64__)) || (defined(__USE_MINGW_ANSI_STDIO) && __USE_MINGW_ANSI_STDIO)
long double val_ld;
if (sscanf("0.8", "%Lf", &val_ld) != 1) {
fails++;
printf("sscanf for long double failed\n");
} else {
// Testing for the success case, to make NaNs trigger the failure case
if (fabsl(val_ld - 0.8) < 0.0001) {
// All ok
} else {
fails++;
printf("sscanf output for long double failed, %f vs %f\n", (double) val_ld, 0.8);
}
}
tests++;
#endif
}
void test_tokenize() {
#define TEST_STRTOK(strtok, type, prefix, strcmp, fmt) do { \
type buf[] = prefix ## " foo bar baz "; \
TEST_XSTR(strcmp, type, strtok(buf, prefix ## " "), prefix ## "foo", fmt); \
TEST_XSTR(strcmp, type, strtok(NULL, prefix ## " "), prefix ## "bar", fmt); \
TEST_XSTR(strcmp, type, strtok(NULL, prefix ## " "), prefix ## "baz", fmt); \
TEST(strtok(NULL, prefix ## " ") == NULL); \
} while (0)
#define TEST_STRTOK_R(strtok_r, type, prefix, strcmp, fmt) do { \
type buf[] = prefix ## " foo bar baz "; \
type *ptr; \
TEST_XSTR(strcmp, type, strtok_r(buf, prefix ## " ", &ptr), prefix ## "foo", fmt); \
TEST_XSTR(strcmp, type, strtok_r(NULL, prefix ## " ", &ptr), prefix ## "bar", fmt); \
TEST_XSTR(strcmp, type, strtok_r(NULL, prefix ## " ", &ptr), prefix ## "baz", fmt); \
TEST(strtok_r(NULL, prefix ## " ", &ptr) == NULL); \
} while (0)
TEST_STRTOK(strtok, char, , strcmp, "%s");
#ifndef _MSC_VER
TEST_STRTOK_R(strtok_r, char, , strcmp, "%s");
#endif
TEST_STRTOK_R(wcstok, wchar_t, L, wcscmp, "%ls");
#ifdef _WIN32
TEST_STRTOK(_wcstok, wchar_t, L, wcscmp, "%ls");
TEST_STRTOK_R(wcstok_s, wchar_t, L, wcscmp, "%ls");
#ifdef __cplusplus
TEST_STRTOK(wcstok, wchar_t, L, wcscmp, "%ls");
#endif
#endif
}
void test_parse_numbers() {
#ifdef _WIN32
_set_invalid_parameter_handler(invalid_parameter);
#endif
errno = 0;
TEST_INT(strtol("foo", NULL, 100), 0);
TEST_INT(errno, EINVAL);
#define TEST_STRTOD(strtod, type, prefix) do { \
const type *curstr = prefix ## "1.e5z"; \
type *end; \
errno = 0; \
TEST_FLT(strtod(curstr, &end), 1.e5); \
TEST_INT(errno, 0); \
TEST_INT((int)(end - curstr), 4); \
curstr = prefix ## "0x125p-1z"; \
errno = 0; \
TEST_FLT(strtod(curstr, &end), 146.5); \
TEST_INT(errno, 0); \
TEST_INT((int)(end - curstr), 8); \
} while (0)
TEST_STRTOD(strtod, char, );
TEST_STRTOD(strtof, char, );
TEST_STRTOD(strtold, char, );
#ifndef MSVCRT_DLL_NOANSI
TEST_STRTOD(wcstod, wchar_t, L);
#endif
TEST_STRTOD(wcstof, wchar_t, L);
TEST_STRTOD(wcstold, wchar_t, L);
#define TEST_STRTOD_VALUE(strtod, prefix, val, expect) \
errno = 0; \
TEST_FLT_SIGN(strtod(prefix ## #val, NULL), expect); \
TEST_INT(errno, 0)
#define TEST_STRTOD_RANGE_EXPECT(strtod, prefix, val, expect) \
errno = 0; \
TEST_FLT_SIGN(strtod(prefix ## #val, NULL), expect); \
TEST_INT(errno, ERANGE)
#define TEST_STRTOD_STRTOF_RANGE_EXPECT(strtod, strtof, prefix, val, expect) \
TEST_STRTOD_VALUE(strtod, prefix, val, val); \
TEST_STRTOD_RANGE_EXPECT(strtof, prefix, val, expect)
#define TEST_STRTOF_32B_RANGE(strtod, strtof, prefix) \
TEST_STRTOD_STRTOF_RANGE_EXPECT(strtod, strtof, prefix, 1.e40, HUGE_VALF); \
TEST_STRTOD_STRTOF_RANGE_EXPECT(strtod, strtof, prefix, -1.e40, -HUGE_VALF); \
TEST_STRTOD_STRTOF_RANGE_EXPECT(strtod, strtof, prefix, 1.e-60, 0.0); \
TEST_STRTOD_STRTOF_RANGE_EXPECT(strtod, strtof, prefix, -1.e-60, -0.0)
TEST_STRTOF_32B_RANGE(strtod, strtof, );
TEST_STRTOF_32B_RANGE(wcstod, wcstof, L);
#define TEST_STRTOD_64B_RANGE(strtod, prefix) \
TEST_STRTOD_RANGE_EXPECT(strtod, prefix, 1.e310, HUGE_VAL); \
TEST_STRTOD_RANGE_EXPECT(strtod, prefix, -1.e310, -HUGE_VAL); \
TEST_STRTOD_RANGE_EXPECT(strtod, prefix, 1.e-400, 0.0); \
TEST_STRTOD_RANGE_EXPECT(strtod, prefix, -1.e-400, -0.0)
TEST_STRTOD_64B_RANGE(strtod, );
#ifndef MSVCRT_DLL_NOANSI
// This test fails on msvcrt.dll when not using -D__USE_MINGW_ANSI_STDIO=1
// (which sets __USE_MINGW_STRTOX), hence the ifdef.
//
// When using the mingw provided wcstod on x86, this test succeeds (and
// thus is included). But when running such a binary emulated on ARM, this
// particular test fails, as the range checks in __mingw_wcstod rely on
// 80 bit long doubles really having more precision than 64 bit doubles.
TEST_STRTOD_64B_RANGE(wcstod, L);
#endif
#if !defined(_MSC_VER) && (__SIZEOF_LONG_DOUBLE__ > __SIZEOF_DOUBLE__)
#define TEST_STRTOLD_80B_RANGE(strtold, prefix) \
TEST_STRTOD_VALUE(strtold, prefix, 1.e310, 1e310L); \
TEST_STRTOD_VALUE(strtold, prefix, -1.e310, -1e310L); \
TEST_STRTOD_VALUE(strtold, prefix, 1.e-310, 1e-310L); \
TEST_STRTOD_VALUE(strtold, prefix, -1.e-310, -1e-310L)
TEST_STRTOLD_80B_RANGE(strtold, );
TEST_STRTOLD_80B_RANGE(wcstold, L);
#else
TEST_STRTOD_64B_RANGE(strtold, );
TEST_STRTOD_64B_RANGE(wcstold, L);
#endif
}
void test_environment() {
int env_ok = 0;
char buf1[] = "CRT_TEST_VAR=1";
putenv(buf1);
for (char **ptr = environ; *ptr; ptr++)
if (!strcmp(*ptr, "CRT_TEST_VAR=1"))
env_ok = 1;
if (!env_ok) {
fails++;
printf("Variable set by putenv not found found in environ\n");
}
tests++;
env_ok = 0;
char buf2[] = "CRT_TEST_VAR=2";
putenv(buf2);
for (char **ptr = environ; *ptr; ptr++)
if (!strcmp(*ptr, "CRT_TEST_VAR=2"))
env_ok = 1;
if (!env_ok) {
fails++;
printf("Variable updated by putenv not found found in environ\n");
}
tests++;
}
void test_math_rounding() {
int i;
#define TEST_FLOOR(floor) \
TEST_FLT(floor(F(3.9)), 3.0); \
TEST_FLT(floor(F(-3.3)), -4.0); \
TEST_FLT(floor(F(-3.9)), -4.0); \
TEST_FLT(floor(F(17179869184.0)), 17179869184.0); \
TEST_FLT(floor(F(1329227995784915872903807060280344576.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT(floor(F(INFINITY)), INFINITY); \
TEST_FLT(floor(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(floor(F(NAN)), F(NAN)); \
TEST_FLT_NAN(floor(F(-NAN)), F(-NAN))
TEST_FLOOR(floor);
TEST_FLOOR(floorf);
TEST_FLOOR(floorl);
#define TEST_CEIL(ceil) \
TEST_FLT(ceil(F(3.9)), 4.0); \
TEST_FLT(ceil(F(-3.3)), -3.0); \
TEST_FLT(ceil(F(-3.9)), -3.0); \
TEST_FLT(ceil(F(17179869184.0)), 17179869184.0); \
TEST_FLT(ceil(F(1329227995784915872903807060280344576.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT(ceil(F(INFINITY)), INFINITY); \
TEST_FLT(ceil(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(ceil(F(NAN)), F(NAN)); \
TEST_FLT_NAN(ceil(F(-NAN)), F(-NAN))
TEST_CEIL(ceil);
TEST_CEIL(ceilf);
TEST_CEIL(ceill);
#define TEST_TRUNC(trunc) \
TEST_FLT(trunc(F(3.9)), 3.0); \
TEST_FLT(trunc(F(-3.3)), -3.0); \
TEST_FLT(trunc(F(-3.9)), -3.0); \
TEST_FLT(trunc(F(17179869184.0)), 17179869184.0); \
TEST_FLT(trunc(F(1329227995784915872903807060280344576.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT(trunc(F(INFINITY)), INFINITY); \
TEST_FLT(trunc(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(trunc(F(NAN)), F(NAN)); \
TEST_FLT_NAN(trunc(F(-NAN)), F(-NAN))
TEST_TRUNC(trunc);
TEST_TRUNC(truncf);
TEST_TRUNC(truncl);
for (i = 0; i < 2; i++) {
if (i == 0) {
// Use the default env in the first round here
context = "FE_DFL_ENV ";
} else {
fesetround(FE_TONEAREST); // Only set it on the second round
context = "FE_TONEAREST ";
}
#define TEST_LRINT_NEAREST(lrint) \
TEST_INT(lrint(F(3.3)), 3); \
TEST_INT(lrint(F(3.6)), 4); \
TEST_INT(lrint(F(3.5)), 4); \
TEST_INT(lrint(F(4.5)), 4); \
TEST_INT(lrint(F(1073741824.0)), 1073741824); \
TEST_INT(lrint(F(-3.3)), -3); \
TEST_INT(lrint(F(-3.6)), -4); \
TEST_INT(lrint(F(-3.5)), -4); \
TEST_INT(lrint(F(-4.5)), -4)
#define TEST_LLRINT_NEAREST(llrint) \
TEST_INT(llrint(F(17179869184.0)), 17179869184); \
TEST_INT(llrint(F(1152921504606846976.0)), 1152921504606846976)
TEST_LRINT_NEAREST(llrint);
TEST_LRINT_NEAREST(llrintf);
TEST_LRINT_NEAREST(llrintl);
TEST_LRINT_NEAREST(lrint);
TEST_LRINT_NEAREST(lrintf);
TEST_LRINT_NEAREST(lrintl);
TEST_LLRINT_NEAREST(llrint);
TEST_LLRINT_NEAREST(llrintf);
TEST_LLRINT_NEAREST(llrintl);
#define TEST_RINT_NEAREST(rint) \
TEST_FLT(rint(F(3.3)), 3.0); \
TEST_FLT(rint(F(3.6)), 4.0); \
TEST_FLT(rint(F(3.5)), 4.0); \
TEST_FLT(rint(F(4.5)), 4.0); \
TEST_FLT(rint(F(17179869184.0)), 17179869184.0); \
TEST_FLT(rint(F(1329227995784915872903807060280344576.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT_NAN(rint(F(NAN)), F(NAN)); \
TEST_FLT(rint(F(-3.3)), -3.0); \
TEST_FLT(rint(F(-3.6)), -4.0); \
TEST_FLT(rint(F(-3.5)), -4.0); \
TEST_FLT(rint(F(-4.5)), -4.0); \
TEST_FLT_NAN(rint(F(-NAN)), F(-NAN))
TEST_RINT_NEAREST(rint);
TEST_RINT_NEAREST(rintf);
TEST_RINT_NEAREST(rintl);
TEST_RINT_NEAREST(nearbyint);
TEST_RINT_NEAREST(nearbyintf);
TEST_RINT_NEAREST(nearbyintl);
}
fesetround(FE_TOWARDZERO);
context = "FE_TOWARDZERO ";
#define TEST_LRINT_TOWARDZERO(lrint) \
TEST_INT(lrint(F(3.3)), 3); \
TEST_INT(lrint(F(3.6)), 3); \
TEST_INT(lrint(F(-3.3)), -3); \
TEST_INT(lrint(F(-3.6)), -3)
TEST_LRINT_TOWARDZERO(llrint);
TEST_LRINT_TOWARDZERO(llrintf);
TEST_LRINT_TOWARDZERO(llrintl);
TEST_LRINT_TOWARDZERO(lrint);
TEST_LRINT_TOWARDZERO(lrintf);
TEST_LRINT_TOWARDZERO(lrintl);
#define TEST_RINT_TOWARDZERO(rint) \
TEST_FLT(rint(F(3.3)), 3.0); \
TEST_FLT(rint(F(3.6)), 3.0); \
TEST_FLT(rint(F(-3.3)), -3.0); \
TEST_FLT(rint(F(-3.6)), -3.0)
TEST_RINT_TOWARDZERO(rint);
TEST_RINT_TOWARDZERO(rintf);
TEST_RINT_TOWARDZERO(rintl);
TEST_RINT_TOWARDZERO(nearbyint);
TEST_RINT_TOWARDZERO(nearbyintf);
TEST_RINT_TOWARDZERO(nearbyintl);
fesetround(FE_DOWNWARD);
context = "FE_DOWNWARD ";
#define TEST_LRINT_DOWNWARD(lrint) \
TEST_INT(lrint(F(3.3)), 3); \
TEST_INT(lrint(F(3.6)), 3); \
TEST_INT(lrint(F(-3.3)), -4); \
TEST_INT(lrint(F(-3.6)), -4)
TEST_LRINT_DOWNWARD(llrint);
TEST_LRINT_DOWNWARD(llrintf);
TEST_LRINT_DOWNWARD(llrintl);
TEST_LRINT_DOWNWARD(lrint);
TEST_LRINT_DOWNWARD(lrintf);
TEST_LRINT_DOWNWARD(lrintl);
#define TEST_RINT_DOWNWARD(rint) \
TEST_FLT(rint(F(3.3)), 3.0); \
TEST_FLT(rint(F(3.6)), 3.0); \
TEST_FLT(rint(F(-3.3)), -4.0); \
TEST_FLT(rint(F(-3.6)), -4.0)
TEST_RINT_DOWNWARD(rint);
TEST_RINT_DOWNWARD(rintf);
TEST_RINT_DOWNWARD(rintl);
TEST_RINT_DOWNWARD(nearbyint);
TEST_RINT_DOWNWARD(nearbyintf);
TEST_RINT_DOWNWARD(nearbyintl);
fesetround(FE_UPWARD);
context = "FE_UPWARD ";
#define TEST_LRINT_UPWARD(lrint) \
TEST_INT(lrint(F(3.3)), 4); \
TEST_INT(lrint(F(3.6)), 4); \
TEST_INT(lrint(F(-3.3)), -3); \
TEST_INT(lrint(F(-3.6)), -3)
TEST_LRINT_UPWARD(llrint);
TEST_LRINT_UPWARD(llrintf);
TEST_LRINT_UPWARD(llrintl);
TEST_LRINT_UPWARD(lrint);
TEST_LRINT_UPWARD(lrintf);
TEST_LRINT_UPWARD(lrintl);
#define TEST_RINT_UPWARD(rint) \
TEST_FLT(rint(F(3.3)), 4.0); \
TEST_FLT(rint(F(3.6)), 4.0); \
TEST_FLT(rint(F(-3.3)), -3.0); \
TEST_FLT(rint(F(-3.6)), -3.0)
TEST_RINT_UPWARD(rint);
TEST_RINT_UPWARD(rintf);
TEST_RINT_UPWARD(rintl);
TEST_RINT_UPWARD(nearbyint);
TEST_RINT_UPWARD(nearbyintf);
TEST_RINT_UPWARD(nearbyintl);
context = "";
fesetround(FE_TONEAREST);
#define TEST_LROUND(lround) \
TEST_INT(lround(F(3.3)), 3); \
TEST_INT(lround(F(3.6)), 4); \
TEST_INT(lround(F(3.5)), 4); \
TEST_INT(lround(F(4.5)), 5); \
TEST_INT(lround(F(1073741824.0)), 1073741824); \
TEST_INT(lround(F(-3.3)), -3); \
TEST_INT(lround(F(-3.6)), -4); \
TEST_INT(lround(F(-3.5)), -4); \
TEST_INT(lround(F(-4.5)), -5)
#define TEST_LLROUND(llrint) \
TEST_INT(llround(F(17179869184.0)), 17179869184); \
TEST_INT(llround(F(1152921504606846976.0)), 1152921504606846976)
TEST_LROUND(llround);
TEST_LROUND(llroundf);
TEST_LROUND(llroundl);
TEST_LROUND(lround);
TEST_LROUND(lroundf);
TEST_LROUND(lroundl);
TEST_LLROUND(llround);
TEST_LLROUND(llroundf);
TEST_LLROUND(llroundl);
#define TEST_ROUND(round) \
TEST_FLT(round(F(3.3)), 3.0); \
TEST_FLT(round(F(3.6)), 4.0); \
TEST_FLT(round(F(3.5)), 4.0); \
TEST_FLT(round(F(4.5)), 5.0); \
TEST_FLT(round(F(17179869184.0)), 17179869184.0); \
TEST_FLT(round(F(1329227995784915872903807060280344576.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT(round(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(round(F(NAN)), F(NAN)); \
TEST_FLT(round(F(-3.3)), -3.0); \
TEST_FLT(round(F(-3.6)), -4.0); \
TEST_FLT(round(F(-3.5)), -4.0); \
TEST_FLT(round(F(-4.5)), -5.0); \
TEST_FLT(round(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(round(F(-NAN)), F(-NAN))
TEST_ROUND(round);
TEST_ROUND(roundf);
TEST_ROUND(roundl);
}
void test_math_roots() {
#define TEST_SQRT(sqrt) \
TEST_FLT(sqrt(F(9)), 3.0); \
TEST_FLT(sqrt(F(0.25)), 0.5); \
TEST_FLT(sqrt(F(INFINITY)), INFINITY); \
TEST_FLT_NAN_ANY(sqrt(F(-1.0))); \
TEST_FLT_NAN_ANY(sqrt(F(-INFINITY))); \
TEST_FLT_NAN(sqrt(F(NAN)), F(NAN)); \
TEST_FLT_NAN(sqrt(F(-NAN)), F(-NAN))
TEST_SQRT(sqrt);
TEST_SQRT(sqrtf);
TEST_SQRT(sqrtl);
#define TEST_CBRT(cbrt) \
TEST_FLT_ACCURACY(cbrt(F(27)), 3.0, 0.001); \
TEST_FLT_ACCURACY(cbrt(F(-27)), -3.0, 0.001); \
TEST_FLT_ACCURACY(cbrt(F(0.125)), 0.5, 0.001); \
TEST_FLT_ACCURACY(cbrt(F(-0.125)), -0.5, 0.001); \
TEST_FLT(cbrt(F(INFINITY)), INFINITY); \
TEST_FLT(cbrt(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(cbrt(F(NAN)), F(NAN)); \
TEST_FLT_NAN(cbrt(F(-NAN)), F(-NAN))
TEST_CBRT(cbrt);
TEST_CBRT(cbrtf);
TEST_CBRT(cbrtl);
#define TEST_HYPOT(hypot) \
TEST_FLT_ACCURACY(hypot(F(1.0), F(1.0)), 1.414214, 0.001); \
TEST_FLT_ACCURACY(hypot(F(-1.0), F(1.0)), 1.414214, 0.001); \
TEST_FLT_ACCURACY(hypot(F(1.0), F(-1.0)), 1.414214, 0.001); \
TEST_FLT_ACCURACY(hypot(F(-1.0), F(-1.0)), 1.414214, 0.001); \
TEST_FLT(hypot(F(INFINITY), F(0.0)), INFINITY); \
TEST_FLT(hypot(F(-INFINITY), F(0.0)), INFINITY); \
TEST_FLT(hypot(F(0.0), F(INFINITY)), INFINITY); \
TEST_FLT(hypot(F(0.0), F(-INFINITY)), INFINITY); \
TEST_FLT_NAN_ANY(hypot(F(NAN), F(0.0))); \
TEST_FLT_NAN_ANY(hypot(F(0.0), F(NAN)))
TEST_HYPOT(hypot);
TEST_HYPOT(hypotf);
TEST_HYPOT(hypotl);
}
void test_math_mod() {
double retd;
float retf;
long double retl;
#define TEST_MODF(modf, retd) \
TEST_FLT_ACCURACY(modf(F(2.1), &retd), 0.1, 0.001); \
TEST_FLT(retd, 2); \
TEST_FLT_ACCURACY(modf(F(-2.1), &retd), -0.1, 0.001); \
TEST_FLT(retd, -2); \
TEST_FLT(modf(F(17179869184.0), &retd), 0); \
TEST_FLT(retd, 17179869184.0); \
TEST_FLT(modf(F(1329227995784915872903807060280344576.0), &retd), 0); \
TEST_FLT(retd, 1329227995784915872903807060280344576.0); \
TEST_FLT(modf(F(INFINITY), &retd), 0); \
TEST_FLT(retd, INFINITY); \
TEST_FLT(modf(F(-INFINITY), &retd), 0); \
TEST_FLT(retd, -INFINITY); \
TEST_FLT_NAN(modf(F(NAN), &retd), F(NAN)); \
TEST_FLT_NAN(retd, F(NAN)); \
TEST_FLT_NAN(modf(F(-NAN), &retd), F(-NAN)); \
TEST_FLT_NAN(retd, F(-NAN))
TEST_MODF(modf, retd);
TEST_MODF(modff, retf);
TEST_MODF(modfl, retl);
#define TEST_FMOD(fmod) \
TEST_FLT_ACCURACY(fmod(F(3.9), F(4.0)), 3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(7.9), F(4.0)), 3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(-3.9), F(4.0)), -3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(3.9), F(-4.0)), 3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(7.9), F(-4.0)), 3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(-3.9), F(-4.0)), -3.9, 0.001); \
TEST_FLT(fmod(F(17179869184.0), F(17180917760.0)), 17179869184.0); \
TEST_FLT(fmod(F(17179869184.0), F(1.0)), 0.0); \
TEST_FLT(fmod(F(1329227995784915872903807060280344576.0), F(1330526069999549579810939684362649600.0)), 1329227995784915872903807060280344576.0); \
TEST_FLT(fmod(F(1329227995784915872903807060280344576.0), F(1.0)), 0.0); \
TEST_FLT_NAN_ANY(fmod(F(INFINITY), F(4.0))); \
TEST_FLT_NAN_ANY(fmod(F(-INFINITY), F(4.0))); \
TEST_FLT_NAN(fmod(F(0), F(NAN)), F(NAN)); \
TEST_FLT_NAN(fmod(F(0), F(-NAN)), F(-NAN)); \
TEST_FLT_NAN(fmod(F(NAN), F(1)), F(NAN)); \
TEST_FLT_NAN(fmod(F(-NAN), F(1)), F(-NAN)); \
TEST_FLT_NAN_ANY(fmod(F(3.9), F(0))); \
TEST_FLT_ACCURACY(fmod(F(3.9), F(INFINITY)), 3.9, 0.001); \
TEST_FLT_ACCURACY(fmod(F(3.9), F(-INFINITY)), 3.9, 0.001)
TEST_FMOD(fmod);
TEST_FMOD(fmodf);
TEST_FMOD(fmodl);
#define TEST_REMAINDER(remainder) \
TEST_FLT_ACCURACY(remainder(F(1.9), F(4.0)), 1.9, 0.001); \
TEST_FLT(remainder(F(2.0), F(4.0)), 2.0); \
TEST_FLT(remainder(F(6.0), F(4.0)), -2.0); \
TEST_FLT(remainder(F(-6.0), F(4.0)), 2.0); \
TEST_FLT_ACCURACY(remainder(F(3.9), F(4.0)), -0.1, 0.001); \
TEST_FLT_ACCURACY(remainder(F(-2.0), F(4.0)), -2.0, 0.001); \
TEST_FLT_ACCURACY(remainder(F(-3.9), F(4.0)), 0.1, 0.001); \
TEST_FLT_ACCURACY(remainder(F(-4.1), F(4.0)), -0.1, 0.001); \
TEST_FLT_ACCURACY(remainder(F(3.9), F(-4.0)), -0.1, 0.001); \
TEST_FLT_ACCURACY(remainder(F(-3.9), F(-4.0)), 0.1, 0.001); \
TEST_FLT(remainder(F(17179869184.0), F(17180917760.0)), -1048576.0); \
TEST_FLT(remainder(F(17179869184.0), F(1.0)), 0.0); \
TEST_FLT(remainder(F(1329227995784915872903807060280344576.0), F(1330526069999549579810939684362649600.0)), -1298074214633706907132624082305024.0); \
TEST_FLT(remainder(F(1329227995784915872903807060280344576.0), F(1.0)), 0.0); \
TEST_FLT_NAN_ANY(remainder(F(INFINITY), F(4.0))); \
TEST_FLT_NAN_ANY(remainder(F(-INFINITY), F(4.0))); \
TEST_FLT_NAN(remainder(F(0), F(NAN)), F(NAN)); \
TEST_FLT_NAN(remainder(F(0), F(-NAN)), F(-NAN)); \
TEST_FLT_NAN(remainder(F(NAN), F(1)), F(NAN)); \
TEST_FLT_NAN(remainder(F(-NAN), F(1)), F(-NAN)); \
TEST_FLT_NAN_ANY(remainder(F(1.9), F(0)))
TEST_REMAINDER(remainder);
TEST_REMAINDER(remainderf);
TEST_REMAINDER(remainderl);
int quo = 42;
#define TEST_REMQUO(remquo) \
TEST_FLT_ACCURACY(remquo(F(1.9), F(4.0), &quo), 1.9, 0.001); \
TEST_INT(quo, 0); \
TEST_FLT(remquo(F(2.0), F(4.0), &quo), 2.0); \
TEST_INT(quo, 0); \
TEST_FLT(remquo(F(6.0), F(4.0), &quo), -2.0); \
TEST_INT(quo, 2); \
TEST_FLT(remquo(F(-6.0), F(4.0), &quo), 2.0); \
TEST_INT(quo, -2); \
TEST_FLT(remquo(F(17179869184.0), F(17180917760.0), &quo), -1048576.0); \
TEST_INT(quo, 1); \
TEST_FLT(remquo(F(1329227995784915872903807060280344576.0), F(1330526069999549579810939684362649600.0), &quo), -1298074214633706907132624082305024.0); \
TEST_INT(quo, 1); \
TEST_FLT_ACCURACY(remquo(F(3.9), F(4.0), &quo), -0.1, 0.001); \
TEST_INT(quo, 1); \
TEST_FLT_ACCURACY(remquo(F(-2.0), F(4.0), &quo), -2.0, 0.001); \
TEST_INT(quo, 0); \
TEST_FLT_ACCURACY(remquo(F(-3.9), F(4.0), &quo), 0.1, 0.001); \
TEST_INT(quo, -1); \
TEST_FLT_ACCURACY(remquo(F(-4.1), F(4.0), &quo), -0.1, 0.001); \
TEST_INT(quo, -1); \
TEST_FLT_ACCURACY(remquo(F(3.9), F(-4.0), &quo), -0.1, 0.001); \
TEST_INT(quo, -1); \
TEST_FLT_ACCURACY(remquo(F(-3.9), F(-4.0), &quo), 0.1, 0.001); \
TEST_INT(quo, 1); \
TEST_FLT_NAN_ANY(remquo(F(INFINITY), F(4.0), &quo)); \
TEST_FLT_NAN_ANY(remquo(F(-INFINITY), F(4.0), &quo)); \
TEST_FLT_NAN(remquo(F(0), F(NAN), &quo), F(NAN)); \
TEST_FLT_NAN(remquo(F(0), F(-NAN), &quo), F(-NAN)); \
TEST_FLT_NAN(remquo(F(NAN), F(0), &quo), F(NAN)); \
TEST_FLT_NAN(remquo(F(-NAN), F(0), &quo), F(-NAN)); \
TEST_FLT_NAN_ANY(remquo(F(1.9), F(0), &quo));
TEST_REMQUO(remquo);
TEST_REMQUO(remquof);
TEST_REMQUO(remquol);
}
void test_math_log_exp() {
#define TEST_LOG(log, HUGE_VAL) \
TEST_FLT_ACCURACY(log(F(1.0)), 0.0, 0.001); \
TEST_FLT_ACCURACY(log(F(2.7182818)), 1.0, 0.001); \
TEST_FLT_ACCURACY(log(F(7.3890561)), 2.0, 0.001); \
TEST_FLT_ACCURACY(log(F(0.3678794)), -1.0, 0.001); \
TEST_FLT(log(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(log(F(NAN)), F(NAN)); \
TEST_FLT_NAN(log(F(-NAN)), F(-NAN)); \
TEST_FLT(log(F(0)), -HUGE_VAL); \
TEST_FLT_NAN_ANY(log(F(-1.0))); \
TEST_FLT_NAN_ANY(log(F(-INFINITY)))
TEST_LOG(log, HUGE_VAL);
TEST_LOG(logf, HUGE_VALF);
TEST_LOG(logl, HUGE_VALL);
#define TEST_LOG2(log2, HUGE_VAL) \
TEST_FLT_ACCURACY(log2(F(1.0)), 0.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(8.0)), 3.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(1024.0)), 10.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(1048576.0)), 20.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(4294967296.0)), 32.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(0.5)), -1.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(0.125)), -3.0, 0.001); \
TEST_FLT_ACCURACY(log2(F(9.7656e-04)), -10, 0.001); \
TEST_FLT_ACCURACY(log2(F(9.5367e-07)), -20, 0.001); \
TEST_FLT_ACCURACY(log2(F(3.5527e-15)), -48, 0.001); \
TEST_FLT_ACCURACY(log2(F(7.8886e-31)), -100, 0.001); \
TEST_FLT_ACCURACY(log2(F(7.3468e-40)), -130, 0.001); \
TEST_FLT_ACCURACY(log2(F(1.225000)), 0.292782, 0.001); /* This, with log2f, crashes the mingw-w64 softfloat implementation */ \
TEST_FLT(log2(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(log2(F(NAN)), F(NAN)); \
TEST_FLT_NAN(log2(F(-NAN)), F(-NAN)); \
TEST_FLT(log2(F(0)), -HUGE_VAL); \
TEST_FLT_NAN_ANY(log2(F(-1.0))); \
TEST_FLT_NAN_ANY(log2(F(-INFINITY)))
TEST_LOG2(log2, HUGE_VAL);
TEST_LOG2(log2f, HUGE_VALF);
TEST_LOG2(log2l, HUGE_VALL);
TEST_FLT_ACCURACY(log2(F(9.8813e-324)), -1073, 0.001);
TEST_FLT_ACCURACY(log2f(F(7.1746e-43)), -140, 0.001);
TEST_FLT_ACCURACY(log2l(F(7.1746e-43)), -140, 0.001);
#define TEST_LOG10(log10, HUGE_VAL) \
TEST_FLT_ACCURACY(log10(F(1.0)), 0.0, 0.001); \
TEST_FLT_ACCURACY(log10(F(10.0)), 1.0, 0.001); \
TEST_FLT_ACCURACY(log10(F(100.0)), 2.0, 0.001); \
TEST_FLT_ACCURACY(log10(F(0.1)), -1.0, 0.001); \
TEST_FLT(log10(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(log10(F(NAN)), F(NAN)); \
TEST_FLT_NAN(log10(F(-NAN)), F(-NAN)); \
TEST_FLT(log10(F(0)), -HUGE_VAL); \
TEST_FLT_NAN_ANY(log10(F(-1.0))); \
TEST_FLT_NAN_ANY(log10(F(-INFINITY)))
TEST_LOG10(log10, HUGE_VAL);
TEST_LOG10(log10f, HUGE_VALF);
TEST_LOG10(log10l, HUGE_VALL);
#define TEST_LOG1P(log1p, HUGE_VAL) \
TEST_FLT_ACCURACY(log1p(F(0.0)), 0.0, 0.001); \
TEST_FLT_ACCURACY(log1p(F(1.718282)), 1.0, 0.001); \
TEST_FLT_ACCURACY(log1p(F(-0.632120)), -1.0, 0.001); \
TEST_FLT(log1p(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(log1p(F(NAN)), F(NAN)); \
TEST_FLT_NAN(log1p(F(-NAN)), F(-NAN)); \
TEST_FLT(log1p(F(-1.0)), -HUGE_VAL); \
TEST_FLT_NAN_ANY(log1p(F(-2.0))); \
TEST_FLT_NAN_ANY(log1p(F(-INFINITY)))
TEST_LOG1P(log1p, HUGE_VAL);
TEST_LOG1P(log1pf, HUGE_VALF);
TEST_LOG1P(log1pl, HUGE_VALL);
#define TEST_EXP(exp) \
TEST_FLT_ACCURACY(exp(F(0.0)), 1.0, 0.001); \
TEST_FLT_ACCURACY(exp(F(1.0)), 2.7182818, 0.001); \
TEST_FLT_ACCURACY(exp(F(2.0)), 7.3890561, 0.001); \
TEST_FLT_ACCURACY(exp(F(-1.0)), 0.3678794, 0.001); \
TEST_FLT(exp(F(INFINITY)), INFINITY); \
TEST_FLT(exp(F(-INFINITY)), 0); \
TEST_FLT_NAN(exp(F(NAN)), F(NAN)); \
TEST_FLT_NAN(exp(F(-NAN)), F(-NAN))
TEST_EXP(exp);
TEST_EXP(expf);
TEST_EXP(expl);
#define TEST_EXP2(exp2) \
TEST_FLT_ACCURACY(exp2(F(0.0)), 1.0, 0.001); \
TEST_FLT_ACCURACY(exp2(F(3.0)), 8.0, 0.001); \
TEST_FLT_ACCURACY(exp2(F(10.0)), 1024.0, 0.001); \
TEST_FLT_ACCURACY(exp2(F(20.0)), 1048576.0, 0.001); \
TEST_FLT_ACCURACY(exp2(F(32.0)), 4294967296.0, 0.001); \
TEST_FLT_ACCURACY(exp2(F(-2.0)), 0.25, 0.001); \
TEST_FLT(exp2(F(INFINITY)), INFINITY); \
TEST_FLT(exp2(F(-INFINITY)), 0); \
TEST_FLT_NAN(exp2(F(NAN)), F(NAN)); \
TEST_FLT_NAN(exp2(F(-NAN)), F(-NAN))
TEST_EXP2(exp2);
TEST_EXP2(exp2f);
TEST_EXP2(exp2l);
#define TEST_EXPM1(expm1) \
TEST_FLT_ACCURACY(expm1(F(0.0)), 0.0, 0.001); \
TEST_FLT_ACCURACY(expm1(F(1.0)), 1.718282, 0.001); \
TEST_FLT_ACCURACY(expm1(F(-1.0)), -0.632120, 0.001); \
TEST_FLT(expm1(F(INFINITY)), INFINITY); \
TEST_FLT(expm1(F(-INFINITY)), -1.0); \
TEST_FLT_NAN(expm1(F(NAN)), F(NAN)); \
TEST_FLT_NAN(expm1(F(-NAN)), F(-NAN))
TEST_EXPM1(expm1);
TEST_EXPM1(expm1f);
TEST_EXPM1(expm1l);
}
void test_math_inspect_manipulate() {
#define TEST_LDEXP(ldexp) \
TEST_FLT_ACCURACY(ldexp(F(0.0), 1), 0.0, 0.001); \
TEST_FLT_ACCURACY(ldexp(F(2.0), 2), 8.0, 0.001); \
TEST_FLT_ACCURACY(ldexp(F(2.0), -2), 0.5, 0.001); \
TEST_FLT(ldexp(F(INFINITY), -42), INFINITY); \
TEST_FLT(ldexp(F(-INFINITY), 42), -INFINITY); \
TEST_FLT_NAN(ldexp(F(NAN), 42), F(NAN)); \
TEST_FLT_NAN(ldexp(F(-NAN), 42), F(-NAN))
TEST_LDEXP(ldexp);
TEST_LDEXP(ldexpf);
TEST_LDEXP(ldexpl);
#define TEST_SCALBN(scalbn) \
TEST_FLT_ACCURACY(scalbn(F(0.0), 1), 0.0, 0.001); \
TEST_FLT_ACCURACY(scalbn(F(2.0), 2), 8.0, 0.001); \
TEST_FLT_ACCURACY(scalbn(F(2.0), -2), 0.5, 0.001); \
TEST_FLT(scalbn(F(INFINITY), -42), INFINITY); \
TEST_FLT(scalbn(F(-INFINITY), 42), -INFINITY); \
TEST_FLT_NAN(scalbn(F(NAN), 42), F(NAN)); \
TEST_FLT_NAN(scalbn(F(-NAN), 42), F(-NAN))
TEST_SCALBN(scalbn);
TEST_SCALBN(scalbnf);
TEST_SCALBN(scalbnf);
TEST_SCALBN(scalbln);
TEST_SCALBN(scalblnf);
TEST_SCALBN(scalblnf);
int iret;
#define TEST_FREXP(frexp) \
TEST_FLT_NAN(frexp(F(NAN), &iret), F(NAN)); \
TEST_FLT_NAN(frexp(F(-NAN), &iret), F(-NAN)); \
iret = 0; \
TEST_FLT(frexp(F(0x1.4p+42), &iret), 0.625); \
TEST_INT(iret, 43)
#define TEST_FREXP_INF(frexp) \
TEST_FLT(frexp(F(INFINITY), &iret), INFINITY); \
TEST_FLT(frexp(F(-INFINITY), &iret), -INFINITY)
TEST_FREXP(frexp);
TEST_FREXP(frexpf);
TEST_FREXP(frexpl);
#if !defined(MSVCRT_DLL) || !(defined(__arm__) || defined(__aarch64__) || defined(__arm64ec__))
// On msvcrt.dll on arm, frexp*(INFINITY) returns NAN, not INFINITY.
TEST_FREXP_INF(frexp);
TEST_FREXP_INF(frexpf);
TEST_FREXP_INF(frexpl);
#endif
#define TEST_ILOGB(ilogb) \
TEST_INT(ilogb(F(1.0)), 0); \
TEST_INT(ilogb(F(0.25)), -2); \
TEST_INT(ilogb(F(-0.25)), -2); \
TEST_INT(ilogb(F(0.0)), FP_ILOGB0); \
TEST_INT(ilogb(F(INFINITY)), INT_MAX); \
TEST_INT(ilogb(F(-INFINITY)), INT_MAX); \
TEST_INT(ilogb(F(NAN)), FP_ILOGBNAN); \
TEST_INT(ilogb(F(-NAN)), FP_ILOGBNAN)
TEST_ILOGB(ilogb);
TEST_ILOGB(ilogbf);
TEST_ILOGB(ilogbl);
TEST_INT(ilogb(3.49514e-308), -1022);
TEST_INT(ilogb(1.74757e-308), -1023);
TEST_INT(ilogb(9.8813e-324), -1073);
TEST_INT(ilogbf(F(3.69292e-38)), -125);
TEST_INT(ilogbf(F(4.61616e-39)), -128);
TEST_INT(ilogbf(F(1.4013e-45)), -149);
TEST_INT(ilogbl(3.49514e-308), -1022);
TEST_INT(ilogbl(1.74757e-308), -1023);
TEST_INT(ilogbl(9.8813e-324), -1073);
#define TEST_LOGB(logb) \
TEST_FLT(logb(F(1.0)), 0.0); \
TEST_FLT(logb(F(0.25)), -2.0); \
TEST_FLT(logb(F(-0.25)), -2.0); \
TEST_FLT(logb(F(0.0)), -INFINITY); \
TEST_FLT(logb(F(INFINITY)), INFINITY); \
TEST_FLT(logb(F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(logb(F(NAN)), F(NAN)); \
TEST_FLT_NAN(logb(F(-NAN)), F(-NAN))
TEST_LOGB(logb);
TEST_LOGB(logbf);
TEST_LOGB(logbl);
TEST_FLT(logb(3.49514e-308), -1022.0);
TEST_FLT(logb(1.74757e-308), -1023.0);
TEST_FLT(logb(9.8813e-324), -1073.0);
TEST_FLT(logbf(F(3.69292e-38)), -125.0);
TEST_FLT(logbf(F(4.61616e-39)), -128.0);
TEST_FLT(logbf(F(1.4013e-45)), -149.0);
TEST_FLT(logbl(3.49514e-308), -1022.0);
TEST_FLT(logbl(1.74757e-308), -1023.0);
TEST_FLT(logbl(9.8813e-324), -1073.0);
}
void test_math_pow() {
#define TEST_POW(pow) \
TEST_FLT(pow(F(2.0), F(0.0)), 1.0); \
TEST_FLT(pow(F(2.0), F(0.0)), 1.0); \
TEST_FLT(pow(F(10.0), F(0.0)), 1.0); \
TEST_FLT(pow(F(10.0), F(1.0)), 10.0); \
TEST_FLT_ACCURACY(pow(F(10.0), F(0.5)), 3.162278, 0.01); \
TEST_FLT_NAN_ANY(pow(F(-1.0), F(1.5))); \
TEST_FLT_SIGN(pow(F(0.0), F(3.0)), 0.0); \
TEST_FLT_SIGN(pow(F(-0.0), F(3.0)), -0.0); \
TEST_FLT_SIGN(pow(F(0.0), F(4.2)), 0.0); \
TEST_FLT_SIGN(pow(F(-0.0), F(4.2)), 0.0); \
TEST_FLT_SIGN(pow(F(INFINITY), F(-0.5)), 0.0); \
TEST_FLT(pow(F(INFINITY), F(0.5)), INFINITY); \
TEST_FLT_SIGN(pow(F(-INFINITY), F(-3)), -0.0); \
TEST_FLT_SIGN(pow(F(-INFINITY), F(-0.5)), 0.0); \
TEST_FLT(pow(F(-INFINITY), F(3.0)), -INFINITY); \
TEST_FLT(pow(F(-INFINITY), F(2.5)), INFINITY); \
TEST_FLT(pow(F(2.0), F(INFINITY)), INFINITY); \
TEST_FLT(pow(F(1.0), F(INFINITY)), 1.0); \
TEST_FLT_SIGN(pow(F(0.5), F(INFINITY)), 0.0); \
TEST_FLT_SIGN(pow(F(2.0), F(-INFINITY)), 0.0); \
TEST_FLT(pow(F(1.0), F(-INFINITY)), 1.0); \
TEST_FLT(pow(F(0.5), F(-INFINITY)), INFINITY); \
TEST_FLT(pow(F(-2.0), F(INFINITY)), INFINITY); \
TEST_FLT(pow(F(-1.0), F(INFINITY)), 1.0); \
TEST_FLT_SIGN(pow(F(-0.5), F(INFINITY)), 0.0); \
TEST_FLT_SIGN(pow(F(-2.0), F(-INFINITY)), 0.0); \
TEST_FLT(pow(F(-1.0), F(-INFINITY)), 1.0); \
TEST_FLT(pow(F(-0.5), F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(pow(F(NAN), F(2.0)), F(NAN)); \
TEST_FLT_NAN(pow(F(-NAN), F(2.0)), F(-NAN)); \
TEST_FLT_NAN(pow(F(2.0), F(NAN)), F(NAN)); \
TEST_FLT_NAN(pow(F(2.0), F(-NAN)), F(-NAN)); \
TEST_FLT(pow(F(1.0), F(NAN)), 1.0); \
TEST_FLT(pow(F(1.0), F(INFINITY)), 1.0); \
TEST_FLT(pow(F(1.0), F(-INFINITY)), 1.0); \
TEST_FLT(pow(F(NAN), F(0.0)), 1.0); \
TEST_FLT(pow(F(INFINITY), F(0.0)), 1.0); \
TEST_FLT(pow(F(-INFINITY), F(0.0)), 1.0)
TEST_POW(pow);
TEST_POW(powf);
TEST_POW(powl);
}
void test_math_trig() {
#define TEST_COS(cos) \
TEST_FLT_ACCURACY(cos(F(0.0)), 1.0, 0.01); \
TEST_FLT_ACCURACY(cos(F(3.141592654)/2), 0.0, 0.01); \
TEST_FLT_ACCURACY(cos(F(3.141592654)), -1.0, 0.01); \
TEST_FLT_ACCURACY(cos(3*F(3.141592654)/2), 0.0, 0.01); \
TEST_FLT_ACCURACY(cos(2*F(3.141592654)), 1.0, 0.01); \
TEST_FLT_NAN_ANY(cos(F(INFINITY))); \
TEST_FLT_NAN_ANY(cos(F(-INFINITY))); \
TEST_FLT_NAN(cos(F(NAN)), F(NAN)); \
TEST_FLT_NAN(cos(F(-NAN)), F(-NAN))
TEST_COS(cos);
TEST_COS(cosf);
TEST_COS(cosl);
#define TEST_SIN(sin) \
TEST_FLT_ACCURACY(sin(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(sin(F(3.141592654)/2), 1.0, 0.01); \
TEST_FLT_ACCURACY(sin(F(3.141592654)), 0.0, 0.01); \
TEST_FLT_ACCURACY(sin(3*F(3.141592654)/2), -1.0, 0.01); \
TEST_FLT_ACCURACY(sin(2*F(3.141592654)), 0.0, 0.01); \
TEST_FLT_NAN_ANY(sin(F(INFINITY))); \
TEST_FLT_NAN_ANY(sin(F(-INFINITY))); \
TEST_FLT_NAN(sin(F(NAN)), F(NAN)); \
TEST_FLT_NAN(sin(F(-NAN)), F(-NAN))
TEST_SIN(sin);
TEST_SIN(sinf);
TEST_SIN(sinl);
#define TEST_TAN(tan) \
TEST_FLT_ACCURACY(tan(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(tan(F(1.0)), 1.557408, 0.01); \
TEST_FLT_ACCURACY(tan(F(3.141592654)/4), 1.0, 0.01); \
TEST_FLT_ACCURACY(tan(3*F(3.141592654)/4), -1.0, 0.01); \
TEST_FLT_ACCURACY(tan(5*F(3.141592654)/4), 1.0, 0.01); \
TEST_FLT_ACCURACY(tan(7*F(3.141592654)/4), -1.0, 0.01); \
TEST_FLT_NAN_ANY(tan(F(INFINITY))); \
TEST_FLT_NAN_ANY(tan(F(-INFINITY))); \
TEST_FLT_NAN(tan(F(NAN)), F(NAN)); \
TEST_FLT_NAN(tan(F(-NAN)), F(-NAN))
TEST_TAN(tan);
TEST_TAN(tanf);
TEST_TAN(tanl);
#define TEST_ACOS(acos) \
TEST_FLT_ACCURACY(acos(F(1.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(acos(F(0.0)), 3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(acos(F(-1.0)), 3.141592654, 0.01); \
TEST_FLT_NAN_ANY(acos(F(1.1))); \
TEST_FLT_NAN_ANY(acos(F(-1.1))); \
TEST_FLT_NAN_ANY(acos(F(INFINITY))); \
TEST_FLT_NAN_ANY(acos(F(-INFINITY))); \
TEST_FLT_NAN(acos(F(NAN)), F(NAN)); \
TEST_FLT_NAN(acos(F(-NAN)), F(-NAN))
TEST_ACOS(acos);
TEST_ACOS(acosf);
TEST_ACOS(acosl);
#define TEST_ASIN(asin) \
TEST_FLT_ACCURACY(asin(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(asin(F(1.0)), 3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(asin(F(-1.0)), -3.141592654/2, 0.01); \
TEST_FLT_NAN_ANY(asin(F(1.1))); \
TEST_FLT_NAN_ANY(asin(F(-1.1))); \
TEST_FLT_NAN_ANY(asin(F(INFINITY))); \
TEST_FLT_NAN_ANY(asin(F(-INFINITY))); \
TEST_FLT_NAN(asin(F(NAN)), F(NAN)); \
TEST_FLT_NAN(asin(F(-NAN)), F(-NAN))
TEST_ASIN(asin);
TEST_ASIN(asinf);
TEST_ASIN(asinl);
#define TEST_ATAN(atan) \
TEST_FLT_ACCURACY(atan(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(atan(F(1.0)), 3.141592654/4, 0.01); \
TEST_FLT_ACCURACY(atan(F(-1.0)), -3.141592654/4, 0.01); \
TEST_FLT_ACCURACY(atan(F(INFINITY)), 3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(atan(F(-INFINITY)), -3.141592654/2, 0.01); \
TEST_FLT_NAN(atan(F(NAN)), F(NAN)); \
TEST_FLT_NAN(atan(F(-NAN)), F(-NAN))
TEST_ATAN(atan);
TEST_ATAN(atanf);
TEST_ATAN(atanl);
#define TEST_ATAN2(atan2) \
TEST_FLT_ACCURACY(atan2(F(0.0), F(-1.0)), 3.141592654, 0.01); \
TEST_FLT_SIGN(atan2(F(0.0), F(1.0)), 0.0); \
TEST_FLT_SIGN(atan2(F(-0.0), F(1.0)), -0.0); \
TEST_FLT_ACCURACY(atan2(F(-1.0), F(0.0)), -3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(atan2(F(1.0), F(0.0)), 3.141592654/2, 0.01); \
TEST_FLT_SIGN(atan2(F(0.0), F(0.0)), 0.0); \
TEST_FLT_ACCURACY(atan2(F(0.0), F(-0.0)), 3.141592654, 0.01); \
TEST_FLT_SIGN(atan2(F(-0.0), F(0.0)), -0.0); \
TEST_FLT_ACCURACY(atan2(F(-0.0), F(-0.0)), -3.141592654, 0.01); \
TEST_FLT_ACCURACY(atan2(F(1.0), F(-INFINITY)), 3.141592654, 0.01); \
TEST_FLT_ACCURACY(atan2(F(-1.0), F(-INFINITY)), -3.141592654, 0.01); \
TEST_FLT_ACCURACY(atan2(F(1.0), F(INFINITY)), 0.0, 0.01); \
TEST_FLT_ACCURACY(atan2(F(INFINITY), F(1.0)), 3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(atan2(F(-INFINITY), F(1.0)), -3.141592654/2, 0.01); \
TEST_FLT_ACCURACY(atan2(F(INFINITY), F(-INFINITY)), 3*3.141592654/4, 0.01); \
TEST_FLT_ACCURACY(atan2(F(-INFINITY), F(-INFINITY)), -3*3.141592654/4, 0.01); \
TEST_FLT_ACCURACY(atan2(F(INFINITY), F(INFINITY)), 3.141592654/4, 0.01); \
TEST_FLT_ACCURACY(atan2(F(-INFINITY), F(INFINITY)), -3.141592654/4, 0.01); \
TEST_FLT_NAN_ANY(atan2(F(NAN), F(1.0))); \
TEST_FLT_NAN_ANY(atan2(F(1.0), F(NAN)))
TEST_ATAN2(atan2);
TEST_ATAN2(atan2f);
TEST_ATAN2(atan2l);
#if defined(__linux__) || defined(__MINGW32__)
double outSin = 42.0, outCos = 42.0;
float outSinf = 42.0, outCosf = 42.0;
long double outSinl = 42.0, outCosl = 42.0;
#define TEST_SINCOS(sincos, outSin, outCos) \
sincos(F(0.0), &outSin, &outCos); \
TEST_FLT_ACCURACY(outSin, 0.0, 0.01); \
TEST_FLT_ACCURACY(outCos, 1.0, 0.01)
TEST_SINCOS(sincos, outSin, outCos);
TEST_SINCOS(sincosf, outSinf, outCosf);
TEST_SINCOS(sincosl, outSinl, outCosl);
#endif
}
void test_math_trig_hyp() {
#define TEST_ACOSH(acosh) \
TEST_FLT_ACCURACY(acosh(F(1.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(acosh(F(2.0)), 1.316958, 0.01); \
TEST_FLT_NAN_ANY(acosh(F(0.0))); \
TEST_FLT_NAN_ANY(acosh(F(-4.0))); \
TEST_FLT_NAN_ANY(acosh(F(-INFINITY))); \
TEST_FLT(acosh(F(INFINITY)), INFINITY); \
TEST_FLT_NAN(acosh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(acosh(F(-NAN)), F(-NAN))
TEST_ACOSH(acosh);
TEST_ACOSH(acoshf);
TEST_ACOSH(acoshl);
#define TEST_ASINH(asinh) \
TEST_FLT_ACCURACY(asinh(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(asinh(F(1.0)), 0.881374, 0.01); \
TEST_FLT_ACCURACY(asinh(F(2.0)), 1.443636, 0.01); \
TEST_FLT_ACCURACY(asinh(F(-1.0)), -0.881374, 0.01); \
TEST_FLT_ACCURACY(asinh(F(-2.0)), -1.443636, 0.01); \
TEST_FLT(asinh(F(INFINITY)), INFINITY); \
TEST_FLT(asinh(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(asinh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(asinh(F(-NAN)), F(-NAN))
TEST_ASINH(asinh);
TEST_ASINH(asinhf);
TEST_ASINH(asinhl);
#define TEST_ATANH(atanh) \
TEST_FLT_ACCURACY(atanh(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(atanh(F(0.5)), 0.549307, 0.01); \
TEST_FLT_ACCURACY(atanh(F(-0.5)), -0.549307, 0.01); \
TEST_FLT(atanh(F(1.0)), INFINITY); \
TEST_FLT(atanh(F(-1.0)), -INFINITY); \
TEST_FLT_NAN_ANY(atanh(F(2.0))); \
TEST_FLT_NAN_ANY(atanh(F(-2.0))); \
TEST_FLT_NAN(atanh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(atanh(F(-NAN)), F(-NAN))
TEST_ATANH(atanh);
TEST_ATANH(atanhf);
TEST_ATANH(atanhl);
#define TEST_COSH(cosh) \
TEST_FLT_ACCURACY(cosh(F(0.0)), 1.0, 0.01); \
TEST_FLT_ACCURACY(cosh(F(1.316958)), 2.0, 0.01); \
TEST_FLT_ACCURACY(cosh(F(-1.316958)), 2.0, 0.01); \
TEST_FLT(cosh(F(INFINITY)), INFINITY); \
TEST_FLT(cosh(F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(cosh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(cosh(F(-NAN)), F(-NAN))
TEST_COSH(cosh);
TEST_COSH(coshf);
TEST_COSH(coshl);
#define TEST_SINH(sinh) \
TEST_FLT_ACCURACY(sinh(F(0.0)), 0.0, 0.01); \
TEST_FLT_ACCURACY(sinh(F(0.881374)), 1.0, 0.01); \
TEST_FLT_ACCURACY(sinh(F(1.443636)), 2.0, 0.01); \
TEST_FLT_ACCURACY(sinh(F(-0.881374)), -1.0, 0.01); \
TEST_FLT_ACCURACY(sinh(F(-1.443636)), -2.0, 0.01); \
TEST_FLT(sinh(F(INFINITY)), INFINITY); \
TEST_FLT(sinh(F(-INFINITY)), -INFINITY); \
TEST_FLT_NAN(sinh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(sinh(F(-NAN)), F(-NAN))
TEST_SINH(sinh);
TEST_SINH(sinhf);
TEST_SINH(sinhl);
#define TEST_TANH(tanh) \
TEST_FLT(tanh(F(0.0)), 0.0); \
TEST_FLT_ACCURACY(tanh(F(0.549307)), 0.5, 0.01); \
TEST_FLT_ACCURACY(tanh(F(-0.549307)), -0.5, 0.01); \
TEST_FLT(tanh(F(INFINITY)), 1.0); \
TEST_FLT(tanh(F(-INFINITY)), -1.0); \
TEST_FLT_NAN(tanh(F(NAN)), F(NAN)); \
TEST_FLT_NAN(tanh(F(-NAN)), F(-NAN))
TEST_TANH(tanh);
TEST_TANH(tanhf);
TEST_TANH(tanhl);
}
void test_math_misc() {
// The plain "NAN" constant in MSVC is negative, while it is positive
// in other environments.
double pNAN = int_to_double(0x7ff8000000000000ULL);
double nNAN = int_to_double(0xfff8000000000000ULL);
#define TEST_FMA(fma) \
TEST_FLT(fma(F(2), F(3), F(4)), 10); \
TEST_FLT_NAN(fma(F(NAN), F(3), F(4)), F(NAN)); \
TEST_FLT_NAN(fma(F(2), F(NAN), F(4)), F(NAN)); \
TEST_FLT_NAN(fma(F(2), F(3), F(NAN)), F(NAN))
TEST_FMA(fma);
TEST_FMA(fmaf);
TEST_FMA(fmal);
#define TEST_FABS(fabs, double) \
TEST_FLT_SIGN(fabs((double)F(0.0)), 0.0); \
TEST_FLT_SIGN(fabs((double)F(-0.0)), 0.0); \
TEST_FLT(fabs((double)F(3.125)), 3.125); \
TEST_FLT(fabs((double)F(-3.125)), 3.125); \
TEST_FLT(fabs((double)F(INFINITY)), INFINITY); \
TEST_FLT(fabs((double)F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(fabs((double)F(pNAN)), F(pNAN)); \
TEST_FLT_NAN(fabs((double)F(nNAN)), F(pNAN))
TEST_FABS(fabs, double);
TEST_FABS(fabsf, float);
TEST_FABS(fabsl, long double);
#define TEST_ERF(erf) \
TEST_FLT(erf(F(0.0)), 0.0); \
TEST_FLT_ACCURACY(erf(F(1.0)), 0.842701, 0.001); \
TEST_FLT_ACCURACY(erf(F(-1.0)), -0.842701, 0.001); \
TEST_FLT_ACCURACY(erf(F(2.0)), 0.995322, 0.001); \
TEST_FLT_ACCURACY(erf(F(-2.0)), -0.995322, 0.001); \
TEST_FLT(erf(F(INFINITY)), 1.0); \
TEST_FLT(erf(F(-INFINITY)), -1.0); \
TEST_FLT_NAN(erf(F(NAN)), F(NAN)); \
TEST_FLT_NAN(erf(F(-NAN)), F(-NAN))
TEST_ERF(erf);
TEST_ERF(erff);
TEST_ERF(erfl);
#define TEST_ERFC(erfc) \
TEST_FLT(erfc(F(0.0)), 1.0); \
TEST_FLT_ACCURACY(erfc(F(1.0)), 0.157299, 0.001); \
TEST_FLT_ACCURACY(erfc(F(-1.0)), 1.842701, 0.001); \
TEST_FLT_ACCURACY(erfc(F(2.0)), 0.004678, 0.001); \
TEST_FLT_ACCURACY(erfc(F(-2.0)), 1.995322, 0.001); \
TEST_FLT(erfc(F(INFINITY)), 0.0); \
TEST_FLT(erfc(F(-INFINITY)), 2.0); \
TEST_FLT_NAN(erfc(F(NAN)), F(NAN)); \
TEST_FLT_NAN(erfc(F(-NAN)), F(-NAN))
TEST_ERFC(erfc);
TEST_ERFC(erfcf);
TEST_ERFC(erfcl);
#define TEST_TGAMMA(tgamma, HUGE_VAL) \
TEST_FLT(tgamma(F(0.0)), HUGE_VAL); \
TEST_FLT(tgamma(F(-0.0)), -HUGE_VAL); \
TEST_FLT_ACCURACY(tgamma(F(0.5)), 1.772454, 0.001); \
TEST_FLT(tgamma(F(1.0)), 1.0); \
TEST_FLT_ACCURACY(tgamma(F(1.5)), 0.886227, 0.001); \
TEST_FLT(tgamma(F(2.0)), 1.0); \
TEST_FLT_ACCURACY(tgamma(F(3.3)), 2.683437, 0.001); \
TEST_FLT(tgamma(F(5.0)), 24.0); \
TEST_FLT_ACCURACY(tgamma(F(-0.5)), -3.544908, 0.001); \
TEST_FLT_NAN_ANY(tgamma(F(-1.0))); \
TEST_FLT_ACCURACY(tgamma(F(-1.5)), 2.363272, 0.001); \
TEST_FLT(tgamma(F(INFINITY)), INFINITY); \
TEST_FLT_NAN_ANY(tgamma(F(-INFINITY))); \
TEST_FLT_NAN(tgamma(F(NAN)), F(NAN)); \
TEST_FLT_NAN(tgamma(F(-NAN)), F(-NAN))
TEST_TGAMMA(tgamma, HUGE_VAL);
TEST_TGAMMA(tgammaf, HUGE_VALF);
TEST_TGAMMA(tgammal, HUGE_VALL);
#ifdef _MSC_VER
#define TEST_SIGNGAM(x) do { } while (0)
#else
#define TEST_SIGNGAM(x) do { x; signgam = 42; } while (0)
signgam = 42;
#endif
#define TEST_LGAMMA(lgamma, HUGE_VAL) \
TEST_FLT(lgamma(F(0.0)), INFINITY); \
TEST_FLT(lgamma(F(-0.0)), INFINITY); \
TEST_FLT_ACCURACY(lgamma(F(0.5)), 0.572365, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT(lgamma(F(1.0)), 0.0); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT_ACCURACY(lgamma(F(1.5)), -0.120782, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT(lgamma(F(2.0)), 0.0); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT_ACCURACY(lgamma(F(3.3)), 0.987099, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT_ACCURACY(lgamma(F(5.0)), 3.178054, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT_ACCURACY(lgamma(F(-0.5)), 1.265512, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, -1)); \
TEST_FLT(lgamma(F(-1.0)), HUGE_VAL); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT_ACCURACY(lgamma(F(-1.5)), 0.860047, 0.001); \
TEST_SIGNGAM(TEST_INT(signgam, 1)); \
TEST_FLT(lgamma(F(INFINITY)), INFINITY); \
TEST_FLT(lgamma(F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(lgamma(F(NAN)), F(NAN)); \
TEST_FLT_NAN(lgamma(F(-NAN)), F(-NAN))
TEST_LGAMMA(lgamma, HUGE_VAL);
TEST_LGAMMA(lgammaf, HUGE_VALF);
TEST_LGAMMA(lgammal, HUGE_VALL);
TEST_FLT_NAN_ANY(nan("foo"));
TEST_FLT_NAN_ANY(nanf("foo"));
TEST_FLT_NAN_ANY(nanl("foo"));
#define TEST_NEXTAFTER(nextafter, DBL_MAX, DBL_EPSILON) \
TEST_FLT(nextafter(F(1.0), F(2.0)), 1.0 + DBL_EPSILON); \
TEST_FLT(nextafter(F(INFINITY), F(INFINITY)), INFINITY); \
TEST_FLT(nextafter(F(INFINITY), F(-INFINITY)), DBL_MAX); \
TEST_FLT(nextafter(F(-INFINITY), F(-INFINITY)), -INFINITY); \
TEST_FLT(nextafter(F(-INFINITY), F(INFINITY)), -DBL_MAX); \
TEST_FLT_NAN(nextafter(F(NAN), F(0.0)), F(NAN)); \
TEST_FLT_NAN(nextafter(F(0.0), F(NAN)), F(NAN)); \
TEST_FLT_NAN(nextafter(F(-NAN), F(0.0)), F(-NAN)); \
TEST_FLT_NAN(nextafter(F(0.0), F(-NAN)), F(-NAN))
TEST_NEXTAFTER(nextafter, DBL_MAX, DBL_EPSILON);
TEST_NEXTAFTER(nextafterf, FLT_MAX, FLT_EPSILON);
TEST_NEXTAFTER(nextafterl, LDBL_MAX, LDBL_EPSILON);
TEST_NEXTAFTER(nexttoward, DBL_MAX, DBL_EPSILON);
TEST_NEXTAFTER(nexttowardf, FLT_MAX, FLT_EPSILON);
TEST_NEXTAFTER(nexttowardl, LDBL_MAX, LDBL_EPSILON);
#define TEST_FDIM(fdim) \
TEST_FLT(fdim(F(2.0), F(1.0)), 1.0); \
TEST_FLT(fdim(F(1.0), F(2.0)), 0.0); \
TEST_FLT(fdim(F(INFINITY), F(1.0)), INFINITY); \
TEST_FLT(fdim(F(1.0), F(-INFINITY)), INFINITY); \
TEST_FLT(fdim(F(-1.0), F(INFINITY)), 0.0); \
TEST_FLT(fdim(F(-INFINITY), F(1.0)), 0.0); \
TEST_FLT(fdim(F(-INFINITY), F(INFINITY)), 0.0); \
TEST_FLT(fdim(F(INFINITY), F(-INFINITY)), INFINITY); \
TEST_FLT_NAN(fdim(F(NAN), F(0.0)), F(NAN)); \
TEST_FLT_NAN(fdim(F(0.0), F(NAN)), F(NAN)); \
TEST_FLT_NAN(fdim(F(-NAN), F(0.0)), F(-NAN)); \
TEST_FLT_NAN(fdim(F(0.0), F(-NAN)), F(-NAN))
TEST_FDIM(fdim);
TEST_FDIM(fdimf);
TEST_FDIM(fdiml);
#define TEST_FMAX(fmax) \
TEST_FLT(fmax(F(1.0), F(0.0)), 1.0); \
TEST_FLT(fmax(F(0.0), F(1.0)), 1.0); \
TEST_FLT(fmax(F(INFINITY), F(1.0)), INFINITY); \
TEST_FLT(fmax(F(-INFINITY), F(1.0)), 1.0); \
TEST_FLT(fmax(F(1.0), F(INFINITY)), INFINITY); \
TEST_FLT(fmax(F(1.0), F(-INFINITY)), 1.0); \
TEST_FLT(fmax(F(1.0), F(NAN)), 1.0); \
TEST_FLT(fmax(F(NAN), F(1.0)), 1.0); \
TEST_FLT_NAN_ANY(fmax(F(NAN), F(-NAN)))
TEST_FMAX(fmax);
TEST_FMAX(fmaxf);
TEST_FMAX(fmaxl);
#define TEST_FMIN(fmin) \
TEST_FLT(fmin(F(1.0), F(0.0)), 0.0); \
TEST_FLT(fmin(F(0.0), F(1.0)), 0.0); \
TEST_FLT(fmin(F(0.0), F(-1.0)), -1.0); \
TEST_FLT(fmin(F(-1.0), F(0.0)), -1.0); \
TEST_FLT(fmin(F(INFINITY), F(1.0)), 1.0); \
TEST_FLT(fmin(F(-INFINITY), F(1.0)), -INFINITY); \
TEST_FLT(fmin(F(1.0), F(INFINITY)), 1.0); \
TEST_FLT(fmin(F(1.0), F(-INFINITY)), -INFINITY); \
TEST_FLT(fmin(F(1.0), F(NAN)), 1.0); \
TEST_FLT(fmin(F(NAN), F(1.0)), 1.0); \
TEST_FLT_NAN_ANY(fmin(F(NAN), F(-NAN)))
TEST_FMIN(fmin);
TEST_FMIN(fminf);
TEST_FMIN(fminl);
TEST_INT(isgreater(F(0.0), F(0.0)), 0);
TEST_INT(isgreater(F(1.0), F(0.0)), 1);
TEST_INT(isgreater(F(0.0), F(1.0)), 0);
TEST_INT(isgreater(F(INFINITY), F(0.0)), 1);
TEST_INT(isgreater(F(-INFINITY), F(0.0)), 0);
TEST_INT(isgreater(F(0.0), F(INFINITY)), 0);
TEST_INT(isgreater(F(0.0), F(-INFINITY)), 1);
TEST_INT(isgreater(F(0.0), F(NAN)), 0);
TEST_INT(isgreater(F(NAN), F(0.0)), 0);
TEST_INT(isgreater(F(NAN), F(NAN)), 0);
TEST_INT(isgreaterequal(F(0.0), F(0.0)), 1);
TEST_INT(isgreaterequal(F(1.0), F(0.0)), 1);
TEST_INT(isgreaterequal(F(0.0), F(1.0)), 0);
TEST_INT(isgreaterequal(F(INFINITY), F(0.0)), 1);
TEST_INT(isgreaterequal(F(-INFINITY), F(0.0)), 0);
TEST_INT(isgreaterequal(F(0.0), F(INFINITY)), 0);
TEST_INT(isgreaterequal(F(0.0), F(-INFINITY)), 1);
TEST_INT(isgreaterequal(F(0.0), F(NAN)), 0);
TEST_INT(isgreaterequal(F(NAN), F(0.0)), 0);
TEST_INT(isgreaterequal(F(NAN), F(NAN)), 0);
TEST_INT(isless(F(0.0), F(0.0)), 0);
TEST_INT(isless(F(1.0), F(0.0)), 0);
TEST_INT(isless(F(0.0), F(1.0)), 1);
TEST_INT(isless(F(INFINITY), F(0.0)), 0);
TEST_INT(isless(F(-INFINITY), F(0.0)), 1);
TEST_INT(isless(F(0.0), F(INFINITY)), 1);
TEST_INT(isless(F(0.0), F(-INFINITY)), 0);
TEST_INT(isless(F(0.0), F(NAN)), 0);
TEST_INT(isless(F(NAN), F(0.0)), 0);
TEST_INT(isless(F(NAN), F(NAN)), 0);
TEST_INT(islessequal(F(0.0), F(0.0)), 1);
TEST_INT(islessequal(F(1.0), F(0.0)), 0);
TEST_INT(islessequal(F(0.0), F(1.0)), 1);
TEST_INT(islessequal(F(INFINITY), F(0.0)), 0);
TEST_INT(islessequal(F(-INFINITY), F(0.0)), 1);
TEST_INT(islessequal(F(0.0), F(INFINITY)), 1);
TEST_INT(islessequal(F(0.0), F(-INFINITY)), 0);
TEST_INT(islessequal(F(0.0), F(NAN)), 0);
TEST_INT(islessequal(F(NAN), F(0.0)), 0);
TEST_INT(islessequal(F(NAN), F(NAN)), 0);
TEST_INT(islessgreater(F(0.0), F(0.0)), 0);
TEST_INT(islessgreater(F(1.0), F(0.0)), 1);
TEST_INT(islessgreater(F(0.0), F(1.0)), 1);
TEST_INT(islessgreater(F(INFINITY), F(0.0)), 1);
TEST_INT(islessgreater(F(-INFINITY), F(0.0)), 1);
TEST_INT(islessgreater(F(0.0), F(INFINITY)), 1);
TEST_INT(islessgreater(F(0.0), F(-INFINITY)), 1);
TEST_INT(islessgreater(F(0.0), F(NAN)), 0);
TEST_INT(islessgreater(F(NAN), F(0.0)), 0);
TEST_INT(islessgreater(F(NAN), F(NAN)), 0);
TEST_INT(isunordered(F(0.0), F(0.0)), 0);
TEST_INT(isunordered(F(1.0), F(0.0)), 0);
TEST_INT(isunordered(F(0.0), F(1.0)), 0);
TEST_INT(isunordered(F(INFINITY), F(0.0)), 0);
TEST_INT(isunordered(F(-INFINITY), F(0.0)), 0);
TEST_INT(isunordered(F(0.0), F(INFINITY)), 0);
TEST_INT(isunordered(F(0.0), F(-INFINITY)), 0);
TEST_INT(isunordered(F(0.0), F(NAN)), 1);
TEST_INT(isunordered(F(NAN), F(0.0)), 1);
TEST_INT(isunordered(F(NAN), F(NAN)), 1);
#define TEST_COPYSIGN(copysign) \
TEST_FLT_ACCURACY(copysign(F(3.125), F(1)), 3.125, 0.0001); \
TEST_FLT_ACCURACY(copysign(F(3.125), F(-1)), -3.125, 0.0001); \
TEST_FLT_ACCURACY(copysign(F(-3.125), F(-1)), -3.125, 0.0001); \
TEST_FLT_ACCURACY(copysign(F(-3.125), F(1)), 3.125, 0.0001); \
TEST_FLT_ACCURACY(copysign(F(3.125), F(nNAN)), -3.125, 0.0001); \
TEST_FLT(copysign(F(INFINITY), F(1)), INFINITY); \
TEST_FLT(copysign(F(INFINITY), F(-1)), -INFINITY); \
TEST_FLT(copysign(F(-INFINITY), F(-1)), -INFINITY); \
TEST_FLT(copysign(F(-INFINITY), F(1)), INFINITY); \
TEST_FLT_NAN(copysign(F(pNAN), F(-1)), F(nNAN)); \
TEST_FLT_NAN(copysign(F(nNAN), F(pNAN)), F(pNAN))
TEST_COPYSIGN(copysign);
TEST_COPYSIGN(copysignf);
TEST_COPYSIGN(copysignl);
#ifdef _WIN32
TEST_COPYSIGN(_copysign);
TEST_COPYSIGN(_copysignf);
TEST_COPYSIGN(_copysignl);
TEST_FLT_ACCURACY(_chgsignl(F(3.125)), -3.125, 0.0001);
TEST_FLT_ACCURACY(_chgsignl(F(-3.125)), 3.125, 0.0001);
TEST_FLT(_chgsignl(F(INFINITY)), -INFINITY);
TEST_FLT(_chgsignl(F(-INFINITY)), INFINITY);
TEST_FLT_NAN(_chgsignl(F(NAN)), F(-NAN));
TEST_FLT_NAN(_chgsignl(F(-NAN)), F(NAN));
#endif
}
void test_compiler_helpers() {
TEST_INT(L(7) / L(7), 1); // __rt_sdiv
TEST_INT(L(-7) / L(7), -1); // __rt_sdiv
TEST_INT(L(-7) / L(-7), 1); // __rt_sdiv
TEST_INT(L(7) / L(-7), -1); // __rt_sdiv
TEST_INT(L(1073741824) / L(3), 357913941); // __rt_sdiv
TEST_INT(L(0) / L(3), 0); // __rt_sdiv
TEST_INT(L(0) / L(-3), 0); // __rt_sdiv
TEST_INT(L(1024) / L(357913941), 0); // __rt_sdiv
TEST_INT(L(1073741824) / L(357913941), 3); // __rt_sdiv
TEST_INT(L(2147483647) / L(1), 2147483647); // __rt_sdiv
TEST_INT(L(2147483647) / L(-1), -2147483647); // __rt_sdiv
TEST_INT(L(-2147483648) / L(1), (long) -2147483648LL); // __rt_sdiv
TEST_INT(UL(7) / L(7), 1); // __rt_udiv
TEST_INT(UL(4294967289) / L(7), 613566755); // __rt_udiv
TEST_INT(UL(4294967289) / L(1), 4294967289UL); // __rt_udiv
TEST_INT(UL(1073741824) / L(3), 357913941); // __rt_udiv
TEST_INT(UL(0) / L(3), 0); // __rt_udiv
TEST_INT(UL(1024) / L(357913941), 0); // __rt_udiv
TEST_INT(UL(1073741824) / L(357913941), 3); // __rt_udiv
TEST_INT(UL(2147483647) / L(1), 2147483647); // __rt_udiv
TEST_INT(LL(7) / 7, 1); // __rt_sdiv64
TEST_INT(LL(-7) / 7, -1); // __rt_sdiv64
TEST_INT(LL(-7) / -7, 1); // __rt_sdiv64
TEST_INT(LL(7) / -7, -1); // __rt_sdiv64
TEST_INT(LL(1073741824) / 3, 357913941); // __rt_sdiv64
TEST_INT(LL(0) / 3, 0); // __rt_sdiv64
TEST_INT(LL(0) / -3, 0); // __rt_sdiv64
TEST_INT(LL(1024) / 357913941, 0); // __rt_sdiv64
TEST_INT(LL(1073741824) / 357913941, 3); // __rt_sdiv64
TEST_INT(LL(2147483647) / LL(1), 2147483647); // __rt_sdiv64
TEST_INT(LL(2147483647) / LL(-1), -2147483647); // __rt_sdiv64
TEST_INT(LL(-2147483648LL) / LL(1), -2147483648LL); // __rt_sdiv64
TEST_INT(LL(0) / LL(2305843009213693952), 0); // __rt_sdiv64
TEST_INT(LL(0) / LL(2305843009213693953), 0); // __rt_sdiv64
TEST_INT(LL(0) / LL(2147483648), 0); // __rt_sdiv64
TEST_INT(LL(0) / LL(4294967296), 0); // __rt_sdiv64
TEST_INT(LL(4294967296) / LL(4294967296), 1); // __rt_sdiv64
TEST_INT(LL(4294967295) / LL(4294967296), 0); // __rt_sdiv64
TEST_INT(LL(883547321287490176) / LL(128), 6902713447558517LL); // __rt_sdiv64
TEST_INT(ULL(7) / 7, 1); // __rt_udiv64
TEST_INT(ULL(4294967289) / LL(7), 613566755); // __rt_udiv64
TEST_INT(ULL(4294967289) / LL(1), 4294967289ULL); // __rt_udiv64
TEST_INT(ULL(1073741824) / LL(3), 357913941); // __rt_udiv64
TEST_INT(ULL(0) / LL(3), 0); // __rt_udiv64
TEST_INT(ULL(1024) / LL(357913941), 0); // __rt_udiv64
TEST_INT(ULL(1073741824) / LL(357913941), 3); // __rt_udiv64
TEST_INT(ULL(2147483647) / LL(1), 2147483647); // __rt_udiv64
TEST_INT(ULL(18446744073709551615) / LL(1), 18446744073709551615ULL); // __rt_udiv64
TEST_INT(ULL(0) / ULL(2305843009213693952), 0); // __rt_udiv64
TEST_INT(ULL(0) / ULL(2305843009213693953), 0); // __rt_udiv64
TEST_INT(ULL(0) / ULL(2147483648), 0); // __rt_udiv64
TEST_INT(ULL(0) / ULL(4294967296), 0); // __rt_udiv64
TEST_INT(ULL(4294967296) / ULL(4294967296), 1); // __rt_udiv64
TEST_INT(ULL(4294967297) / ULL(8589934593), 0); // __rt_udiv64
TEST_INT(ULL(883547321287490176) / ULL(128), 6902713447558517ULL); // __rt_udiv64
TEST_INT(L(7) % L(7), 0); // __rt_sdiv
TEST_INT(L(-7) % L(7), 0); // __rt_sdiv
TEST_INT(L(-7) % L(-7), 0); // __rt_sdiv
TEST_INT(L(7) % L(-7), 0); // __rt_sdiv
TEST_INT(L(1073741824) % L(3), 1); // __rt_sdiv
TEST_INT(L(0) % L(3), 0); // __rt_sdiv
TEST_INT(L(0) % L(-3), 0); // __rt_sdiv
TEST_INT(L(1024) % L(357913941), 1024); // __rt_sdiv
TEST_INT(L(1073741824) % L(357913941), 1); // __rt_sdiv
TEST_INT(L(2147483647) % L(1), 0); // __rt_sdiv
TEST_INT(L(2147483647) % L(-1), 0); // __rt_sdiv
TEST_INT(L(-2147483648) % L(1), 0); // __rt_sdiv
TEST_INT(UL(7) % L(7), 0); // __rt_udiv
TEST_INT(UL(4294967289) % L(7), 4); // __rt_udiv
TEST_INT(UL(4294967289) % L(1), 0); // __rt_udiv
TEST_INT(UL(1073741824) % L(3), 1); // __rt_udiv
TEST_INT(UL(0) % L(3), 0); // __rt_udiv
TEST_INT(UL(1024) % L(357913941), 1024); // __rt_udiv
TEST_INT(UL(1073741824) % L(357913941), 1); // __rt_udiv
TEST_INT(UL(2147483647) % L(1), 0); // __rt_udiv
TEST_INT(LL(7) % 7, 0); // __rt_sdiv64
TEST_INT(LL(-7) % 7, 0); // __rt_sdiv64
TEST_INT(LL(-7) % -7, 0); // __rt_sdiv64
TEST_INT(LL(7) % -7, 0); // __rt_sdiv64
TEST_INT(LL(1073741824) % 3, 1); // __rt_sdiv64
TEST_INT(LL(0) % 3, 0); // __rt_sdiv64
TEST_INT(LL(0) % -3, 0); // __rt_sdiv64
TEST_INT(LL(1024) % 357913941, 1024); // __rt_sdiv64
TEST_INT(LL(1073741824) % 357913941, 1); // __rt_sdiv64
TEST_INT(LL(2147483647) % LL(1), 0); // __rt_sdiv64
TEST_INT(LL(2147483647) % LL(-1), 0); // __rt_sdiv64
TEST_INT(LL(-2147483648) % LL(1), 0); // __rt_sdiv64
TEST_INT(LL(0) % LL(2305843009213693952), 0); // __rt_sdiv64
TEST_INT(LL(0) % LL(2305843009213693953), 0); // __rt_sdiv64
TEST_INT(LL(0) % LL(2147483648), 0); // __rt_sdiv64
TEST_INT(LL(0) % LL(4294967296), 0); // __rt_sdiv64
TEST_INT(LL(4294967296) % LL(4294967296), 0); // __rt_sdiv64
TEST_INT(LL(4294967295) % LL(4294967296), 4294967295LL); // __rt_sdiv64
TEST_INT(ULL(7) % 7, 0); // __rt_udiv64
TEST_INT(ULL(4294967289) % LL(7), 4); // __rt_udiv64
TEST_INT(ULL(4294967289) % LL(1), 0); // __rt_udiv64
TEST_INT(ULL(1073741824) % LL(3), 1); // __rt_udiv64
TEST_INT(ULL(0) % LL(3), 0); // __rt_udiv64
TEST_INT(ULL(1024) % LL(357913941), 1024); // __rt_udiv64
TEST_INT(ULL(1073741824) % LL(357913941), 1); // __rt_udiv64
TEST_INT(ULL(2147483647) % LL(1), 0); // __rt_udiv64
TEST_INT(ULL(18446744073709551615) % LL(1), 0); // __rt_udiv64
TEST_INT(ULL(0) % ULL(2305843009213693952), 0); // __rt_udiv64
TEST_INT(ULL(0) % ULL(2305843009213693953), 0); // __rt_udiv64
TEST_INT(ULL(0) % ULL(2147483648), 0); // __rt_udiv64
TEST_INT(ULL(0) % ULL(4294967296), 0); // __rt_udiv64
TEST_INT(ULL(4294967296) % ULL(4294967296), 0); // __rt_udiv64
TEST_INT(ULL(4294967297) % ULL(8589934593), 4294967297ULL); // __rt_udiv64
TEST_INT((unsigned long long)F(4.2), 4);
TEST_INT((signed long long)F(4.2), 4);
TEST_INT((unsigned long long)F(123456789012345678.0), 123456789012345680ULL);
TEST_INT((signed long long)F(123456789012345678.0), 123456789012345680ULL);
TEST_INT((signed long long)F(-123456789012345.0), -123456789012345LL);
TEST_INT((unsigned long long)(float)F(4.2), 4);
TEST_INT((signed long long)(float)F(4.2), 4);
TEST_INT((unsigned long long)(float)F(274877906944), 274877906944ULL);
TEST_INT((signed long long)(float)F(274877906944), 274877906944ULL);
TEST_INT((signed long long)(float)F(-274877906944), -274877906944LL);
TEST_INT((unsigned long long)(long double)F(4.2), 4);
TEST_INT((signed long long)(long double)F(4.2), 4);
TEST_INT((unsigned long long)(long double)F(274877906944), 274877906944ULL);
TEST_INT((signed long long)(long double)F(274877906944), 274877906944ULL);
TEST_INT((signed long long)(long double)F(-274877906944), -274877906944LL);
TEST_FLT((double)LL(4), 4.0);
TEST_FLT((float)LL(4), 4.0);
TEST_FLT((long double)LL(4), 4.0);
TEST_FLT((double)LL(123456789012345), 123456789012345.0);
TEST_FLT((double)LL(-123456789012345), -123456789012345.0);
TEST_FLT((float)LL(274877906944), 274877906944.0);
TEST_FLT((float)LL(-274877906944), -274877906944.0);
TEST_FLT((long double)LL(274877906944), 274877906944.0);
TEST_FLT((long double)LL(-274877906944), -274877906944.0);
TEST_FLT((double)ULL(4), 4.0);
TEST_FLT((float)ULL(4), 4.0);
TEST_FLT((long double)ULL(4), 4.0);
TEST_FLT((double)ULL(17293822569102704640), 17293822569102704640.0);
TEST_FLT((float)ULL(17293822569102704640), 17293822569102704640.0);
TEST_FLT((long double)ULL(17293822569102704640), 17293822569102704640.0);
#ifdef __SIZEOF_INT128__
TEST_INT((__uint128_t)F(4.2), 4);
TEST_INT((__int128_t)F(4.2), 4);
TEST_INT((__uint128_t)F(123456789012345678.0), 123456789012345680ULL);
TEST_INT((__int128_t)F(123456789012345678.0), 123456789012345680ULL);
TEST_INT((__int128_t)F(-123456789012345.0), -123456789012345LL);
TEST_INT((__uint128_t)(float)F(4.2), 4);
TEST_INT((__int128_t)(float)F(4.2), 4);
TEST_INT((__uint128_t)(float)F(274877906944), 274877906944ULL);
TEST_INT((__int128_t)(float)F(274877906944), 274877906944ULL);
TEST_INT((__int128_t)(float)F(-274877906944), -274877906944LL);
TEST_INT((__uint128_t)(long double)F(4.2), 4);
TEST_INT((__int128_t)(long double)F(4.2), 4);
TEST_INT((__uint128_t)(long double)F(274877906944), 274877906944ULL);
TEST_INT((__int128_t)(long double)F(274877906944), 274877906944ULL);
TEST_INT((__int128_t)(long double)F(-274877906944), -274877906944LL);
TEST_FLT((double)(__int128_t)LL(4), 4.0);
TEST_FLT((float)(__int128_t)LL(4), 4.0);
TEST_FLT((long double)(__int128_t)LL(4), 4.0);
TEST_FLT((double)(__int128_t)LL(123456789012345), 123456789012345.0);
TEST_FLT((double)(__int128_t)LL(-123456789012345), -123456789012345.0);
TEST_FLT((float)(__int128_t)LL(274877906944), 274877906944.0);
TEST_FLT((float)(__int128_t)LL(-274877906944), -274877906944.0);
TEST_FLT((long double)(__int128_t)LL(274877906944), 274877906944.0);
TEST_FLT((long double)(__int128_t)LL(-274877906944), -274877906944.0);
TEST_FLT((double)(__uint128_t)ULL(4), 4.0);
TEST_FLT((float)(__uint128_t)ULL(4), 4.0);
TEST_FLT((long double)(__uint128_t)ULL(4), 4.0);
TEST_FLT((double)(__uint128_t)ULL(17293822569102704640), 17293822569102704640.0);
TEST_FLT((float)(__uint128_t)ULL(17293822569102704640), 17293822569102704640.0);
TEST_FLT((long double)(__uint128_t)ULL(17293822569102704640), 17293822569102704640.0);
#endif
}
void test_win32_intrinsics() {
#ifdef _WIN32
long value = 0;
__int64 ret;
__int64 value64 = 0;
void *ptr = NULL;
void *ptr1 = &value;
void *ptr2 = &value64;
void *ret_ptr;
#define TEST_FUNC(expr, var, expected, expected_ret) do { \
ret = expr; \
TEST_INT(var, expected); \
TEST_INT(ret, expected_ret); \
var = expected; \
} while (0)
#define TEST_FUNC_PTR(expr, var, expected, expected_ret) do { \
ret_ptr = expr; \
TEST_PTR(var, expected); \
TEST_PTR(ret_ptr, expected_ret); \
var = expected; \
} while (0)
TEST_FUNC(InterlockedBitTestAndSet(&value, 0), value, 1, 0);
TEST_FUNC(InterlockedBitTestAndSet(&value, 2), value, 5, 0);
TEST_FUNC(InterlockedBitTestAndSet(&value, 2), value, 5, 1);
TEST_FUNC(InterlockedBitTestAndReset(&value, 2), value, 1, 1);
TEST_FUNC(InterlockedBitTestAndReset(&value, 2), value, 1, 0);
TEST_FUNC(InterlockedBitTestAndReset(&value, 0), value, 0, 1);
#ifdef _WIN64
TEST_FUNC(InterlockedBitTestAndSet64(&value64, 0), value64, 1, 0);
TEST_FUNC(InterlockedBitTestAndSet64(&value64, 2), value64, 5, 0);
TEST_FUNC(InterlockedBitTestAndSet64(&value64, 2), value64, 5, 1);
TEST_FUNC(InterlockedBitTestAndSet64(&value64, 40), value64, 0x10000000005, 0);
TEST_FUNC(InterlockedBitTestAndReset64(&value64, 40), value64, 5, 1);
TEST_FUNC(InterlockedBitTestAndReset64(&value64, 2), value64, 1, 1);
TEST_FUNC(InterlockedBitTestAndReset64(&value64, 2), value64, 1, 0);
TEST_FUNC(InterlockedBitTestAndReset64(&value64, 0), value64, 0, 1);
#endif
TEST_FUNC(InterlockedIncrement(&value), value, 1, 1);
TEST_FUNC(InterlockedDecrement(&value), value, 0, 0);
TEST_FUNC(InterlockedAdd(&value, 7), value, 7, 7);
TEST_FUNC(InterlockedAdd(&value, -2), value, 5, 5);
TEST_FUNC(InterlockedAdd64(&value64, 7), value64, 7, 7);
TEST_FUNC(InterlockedAdd64(&value64, 0x10000000000), value64, 0x10000000007, 0x10000000007);
TEST_FUNC(InterlockedIncrement64(&value64), value64, 0x10000000008, 0x10000000008);
TEST_FUNC(InterlockedDecrement64(&value64), value64, 0x10000000007, 0x10000000007);
TEST_FUNC(InterlockedAdd64(&value64, -0x10000000002), value64, 5, 5);
// Exchange functions return the previous value
TEST_FUNC(InterlockedExchangeAdd(&value, 1), value, 6, 5);
TEST_FUNC(InterlockedExchange(&value, 2), value, 2, 6);
TEST_FUNC(InterlockedCompareExchange(&value, 7, 1), value, 2, 2);
TEST_FUNC(InterlockedCompareExchange(&value, 5, 2), value, 5, 2);
TEST_FUNC_PTR(InterlockedExchangePointer(&ptr, ptr1), ptr, ptr1, (void*)NULL);
TEST_FUNC_PTR(InterlockedExchangePointer(&ptr, ptr2), ptr, ptr2, ptr1);
TEST_FUNC_PTR(InterlockedCompareExchangePointer(&ptr, NULL, ptr1), ptr, ptr2, ptr2);
TEST_FUNC_PTR(InterlockedCompareExchangePointer(&ptr, NULL, ptr2), ptr, (void*)NULL, ptr2);
TEST_FUNC(InterlockedExchangeAdd64(&value64, 0x10000000000), value64, 0x10000000005, 5);
TEST_FUNC(InterlockedExchange64(&value64, 0x10000000000), value64, 0x10000000000, 0x10000000005);
TEST_FUNC(InterlockedCompareExchange64(&value64, 7, 1), value64, 0x10000000000, 0x10000000000);
TEST_FUNC(InterlockedCompareExchange64(&value64, 0x20000000005, 0x10000000000), value64, 0x20000000005, 0x10000000000);
// Logical operations returns the previous value
TEST_FUNC(InterlockedOr(&value, 2), value, 7, 5);
TEST_FUNC(InterlockedOr(&value, 2), value, 7, 7);
TEST_FUNC(InterlockedAnd(&value, 2), value, 2, 7);
TEST_FUNC(InterlockedAnd(&value, 2), value, 2, 2);
TEST_FUNC(InterlockedXor(&value, 2), value, 0, 2);
TEST_FUNC(InterlockedXor(&value, 2), value, 2, 0);
TEST_FUNC(InterlockedXor(&value, 2), value, 0, 2);
TEST_FUNC(InterlockedOr64(&value64, 2), value64, 0x20000000007, 0x20000000005);
TEST_FUNC(InterlockedOr64(&value64, 0x10000000000), value64, 0x30000000007, 0x20000000007);
TEST_FUNC(InterlockedAnd64(&value64, 0x20000000000), value64, 0x20000000000, 0x30000000007);
TEST_FUNC(InterlockedAnd64(&value64, 0x20000000000), value64, 0x20000000000, 0x20000000000);
TEST_FUNC(InterlockedXor64(&value64, 0x20000000000), value64, 0, 0x20000000000);
TEST_FUNC(InterlockedXor64(&value64, 0x20000000000), value64, 0x20000000000, 0);
TEST_FUNC(InterlockedXor64(&value64, 0x20000000000), value64, 0, 0x20000000000);
unsigned long idx = 42;
// If no bit is set, idx is set to an undefined value.
TEST_INT(BitScanForward(&idx, UL(0)), 0);
TEST_FUNC(BitScanForward(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(BitScanForward(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(BitScanForward(&idx, UL(0x80000001)), idx, 0, 1);
TEST_INT(BitScanReverse(&idx, UL(0)), 0);
TEST_FUNC(BitScanReverse(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(BitScanReverse(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(BitScanReverse(&idx, UL(0x80000001)), idx, 31, 1);
#if !defined(_M_ARM) && !defined(__arm__) && !defined(__i386__)
// These seem to be unavailable on 32 bit arm even in MSVC. They're also missing
// on i386 mingw.
TEST_INT(BitScanForward64(&idx, UL(0)), 0);
TEST_FUNC(BitScanForward64(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(BitScanForward64(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(BitScanForward64(&idx, UL(0x80000001)), idx, 0, 1);
TEST_FUNC(BitScanForward64(&idx, ULL(0x8000000000000000)), idx, 63, 1);
TEST_INT(BitScanReverse64(&idx, UL(0)), 0);
TEST_FUNC(BitScanReverse64(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(BitScanReverse64(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(BitScanReverse64(&idx, UL(0x80000001)), idx, 31, 1);
TEST_FUNC(BitScanReverse64(&idx, ULL(0x8000000000000000)), idx, 63, 1);
#endif
// Test intrinsics versions. Not all combinations are available.
TEST_FUNC(_interlockedbittestandset(&value, 0), value, 1, 0);
TEST_FUNC(_interlockedbittestandset(&value, 2), value, 5, 0);
TEST_FUNC(_interlockedbittestandset(&value, 2), value, 5, 1);
TEST_FUNC(_interlockedbittestandreset(&value, 2), value, 1, 1);
TEST_FUNC(_interlockedbittestandreset(&value, 2), value, 1, 0);
TEST_FUNC(_interlockedbittestandreset(&value, 0), value, 0, 1);
TEST_FUNC(_InterlockedIncrement(&value), value, 1, 1);
TEST_FUNC(_InterlockedDecrement(&value), value, 0, 0);
TEST_FUNC(_bittestandset(&value, 0), value, 1, 0);
TEST_FUNC(_bittestandset(&value, 2), value, 5, 0);
TEST_FUNC(_bittestandset(&value, 2), value, 5, 1);
TEST_INT(_bittest(&value, 0), 1);
TEST_INT(_bittest(&value, 1), 0);
TEST_FUNC(_bittestandreset(&value, 2), value, 1, 1);
TEST_FUNC(_bittestandreset(&value, 2), value, 1, 0);
TEST_FUNC(_bittestandreset(&value, 0), value, 0, 1);
TEST_FUNC(_bittestandcomplement(&value, 2), value, 4, 0);
TEST_FUNC(_bittestandcomplement(&value, 2), value, 0, 1);
#ifdef _WIN64
TEST_FUNC(_interlockedbittestandset64(&value64, 0), value64, 1, 0);
TEST_FUNC(_interlockedbittestandset64(&value64, 2), value64, 5, 0);
TEST_FUNC(_interlockedbittestandset64(&value64, 2), value64, 5, 1);
TEST_FUNC(_interlockedbittestandset64(&value64, 40), value64, 0x10000000005, 0);
TEST_FUNC(_interlockedbittestandset64(&value64, 41), value64, 0x30000000005, 0);
TEST_FUNC(_interlockedbittestandreset64(&value64, 40), value64, 0x20000000005, 1);
TEST_FUNC(_interlockedbittestandreset64(&value64, 2), value64, 0x20000000001, 1);
TEST_FUNC(_interlockedbittestandreset64(&value64, 2), value64, 0x20000000001, 0);
TEST_FUNC(_interlockedbittestandreset64(&value64, 0), value64, 0x20000000000, 1);
TEST_FUNC(_InterlockedIncrement64(&value64), value64, 0x20000000001, 0x20000000001);
TEST_FUNC(_InterlockedDecrement64(&value64), value64, 0x20000000000, 0x20000000000);
TEST_FUNC(_interlockedbittestandreset64(&value64, 41), value64, 0, 1);
TEST_FUNC(_bittestandset64(&value64, 0), value64, 1, 0);
TEST_FUNC(_bittestandset64(&value64, 2), value64, 5, 0);
TEST_FUNC(_bittestandset64(&value64, 2), value64, 5, 1);
TEST_FUNC(_bittestandset64(&value64, 40), value64, 0x10000000005, 0);
TEST_FUNC(_bittestandset64(&value64, 41), value64, 0x30000000005, 0);
TEST_INT(_bittest64(&value64, 0), 1);
TEST_INT(_bittest64(&value64, 1), 0);
TEST_INT(_bittest64(&value64, 41), 1);
TEST_INT(_bittest64(&value64, 42), 0);
TEST_FUNC(_bittestandreset64(&value64, 40), value64, 0x20000000005, 1);
TEST_FUNC(_bittestandreset64(&value64, 2), value64, 0x20000000001, 1);
TEST_FUNC(_bittestandreset64(&value64, 2), value64, 0x20000000001, 0);
TEST_FUNC(_bittestandreset64(&value64, 0), value64, 0x20000000000, 1);
TEST_FUNC(_bittestandcomplement64(&value64, 2), value64, 0x20000000004, 0);
TEST_FUNC(_bittestandcomplement64(&value64, 2), value64, 0x20000000000, 1);
TEST_FUNC(_bittestandcomplement64(&value64, 40), value64, 0x30000000000, 0);
TEST_FUNC(_bittestandcomplement64(&value64, 40), value64, 0x20000000000, 1);
TEST_FUNC(_bittestandreset64(&value64, 41), value64, 0, 1);
#endif
TEST_FUNC(_InterlockedExchangeAdd(&value, 1), value, 1, 0);
TEST_FUNC(_InterlockedExchange(&value, 2), value, 2, 1);
TEST_FUNC(_InterlockedCompareExchange(&value, 7, 1), value, 2, 2);
TEST_FUNC(_InterlockedCompareExchange(&value, 0, 2), value, 0, 2);
TEST_FUNC_PTR(_InterlockedExchangePointer(&ptr, ptr1), ptr, ptr1, (void*)NULL);
TEST_FUNC_PTR(_InterlockedExchangePointer(&ptr, ptr2), ptr, ptr2, ptr1);
TEST_FUNC_PTR(_InterlockedCompareExchangePointer(&ptr, NULL, ptr1), ptr, ptr2, ptr2);
TEST_FUNC_PTR(_InterlockedCompareExchangePointer(&ptr, NULL, ptr2), ptr, (void*)NULL, ptr2);
#ifdef _WIN64
TEST_FUNC(_InterlockedExchangeAdd64(&value64, 0x20000000000), value64, 0x20000000000, 0);
TEST_FUNC(_InterlockedExchange64(&value64, 0x10000000000), value64, 0x10000000000, 0x20000000000);
TEST_FUNC(_InterlockedCompareExchange64(&value64, 7, 1), value64, 0x10000000000, 0x10000000000);
TEST_FUNC(_InterlockedCompareExchange64(&value64, 0x20000000000, 0x10000000000), value64, 0x20000000000, 0x10000000000);
#endif
TEST_FUNC(_InterlockedOr(&value, 2), value, 2, 0);
TEST_FUNC(_InterlockedOr(&value, 5), value, 7, 2);
TEST_FUNC(_InterlockedAnd(&value, 2), value, 2, 7);
TEST_FUNC(_InterlockedAnd(&value, 2), value, 2, 2);
TEST_FUNC(_InterlockedXor(&value, 2), value, 0, 2);
TEST_FUNC(_InterlockedXor(&value, 2), value, 2, 0);
TEST_FUNC(_InterlockedXor(&value, 2), value, 0, 2);
#ifdef _WIN64
TEST_FUNC(_InterlockedOr64(&value64, 0x10000000000), value64, 0x30000000000, 0x20000000000);
TEST_FUNC(_InterlockedAnd64(&value64, 0x10000000000), value64, 0x10000000000, 0x30000000000);
TEST_FUNC(_InterlockedXor64(&value64, 0x10000000000), value64, 0, 0x10000000000);
TEST_FUNC(_InterlockedXor64(&value64, 0x10000000000), value64, 0x10000000000, 0);
TEST_FUNC(_InterlockedXor64(&value64, 0x10000000000), value64, 0, 0x10000000000);
#endif
TEST_INT(_BitScanForward(&idx, UL(0)), 0);
TEST_FUNC(_BitScanForward(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(_BitScanForward(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(_BitScanForward(&idx, UL(0x80000001)), idx, 0, 1);
TEST_INT(_BitScanReverse(&idx, UL(0)), 0);
TEST_FUNC(_BitScanReverse(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(_BitScanReverse(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(_BitScanReverse(&idx, UL(0x80000001)), idx, 31, 1);
#ifdef _WIN64
TEST_INT(_BitScanForward64(&idx, UL(0)), 0);
TEST_FUNC(_BitScanForward64(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(_BitScanForward64(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(_BitScanForward64(&idx, UL(0x80000001)), idx, 0, 1);
TEST_FUNC(_BitScanForward64(&idx, ULL(0x8000000000000000)), idx, 63, 1);
TEST_INT(_BitScanReverse64(&idx, UL(0)), 0);
TEST_FUNC(_BitScanReverse64(&idx, UL(1)), idx, 0, 1);
TEST_FUNC(_BitScanReverse64(&idx, UL(0x80000000)), idx, 31, 1);
TEST_FUNC(_BitScanReverse64(&idx, UL(0x80000001)), idx, 31, 1);
TEST_FUNC(_BitScanReverse64(&idx, ULL(0x8000000000000000)), idx, 63, 1);
#endif
#endif
}
int main(int argc, char *argv[]) {
test_strings();
test_tokenize();
test_parse_numbers();
test_environment();
test_math_rounding();
test_math_roots();
test_math_mod();
test_math_log_exp();
test_math_inspect_manipulate();
test_math_pow();
test_math_trig();
test_math_trig_hyp();
test_math_misc();
test_compiler_helpers();
test_win32_intrinsics();
printf("%d tests, %d failures\n", tests, fails);
return fails > 0;
}