From c25bd5bc78f841f04adf6de6384653331e5f6ced Mon Sep 17 00:00:00 2001
From: Stefano D'Angelo <stefano.dangelo@orastron.com>
Date: Wed, 9 Aug 2023 13:11:09 +0200
Subject: [PATCH] bw_math wip + settling on c++11 + some tests

---
 examples/common/android/android.mk |   2 +-
 examples/common/vst3/vst3.mk       |   2 +-
 include/bw_math.h                  | 171 +++++----
 test/Makefile                      |  16 +
 test/bw_math.c                     | 536 +++++++++++++++++++++++++++++
 5 files changed, 655 insertions(+), 72 deletions(-)
 create mode 100644 test/Makefile
 create mode 100644 test/bw_math.c

diff --git a/examples/common/android/android.mk b/examples/common/android/android.mk
index 1b5077e..465dcde 100644
--- a/examples/common/android/android.mk
+++ b/examples/common/android/android.mk
@@ -39,7 +39,7 @@ CXXFLAGS := \
 	-Wall \
 	-Wextra \
 	-Wpedantic \
-	-std=c++20
+	-std=c++11
 LDFLAGS := \
 	-shared \
 	-static-libstdc++ \
diff --git a/examples/common/vst3/vst3.mk b/examples/common/vst3/vst3.mk
index c1d79a3..567947c 100644
--- a/examples/common/vst3/vst3.mk
+++ b/examples/common/vst3/vst3.mk
@@ -17,7 +17,7 @@ CXXFLAGS := \
 	-Wall \
 	-Wextra \
 	-Wpedantic \
-	-std=c++20
+	-std=c++11
 LDFLAGS := \
 	-shared
 
diff --git a/include/bw_math.h b/include/bw_math.h
index 83a1883..f89401b 100644
--- a/include/bw_math.h
+++ b/include/bw_math.h
@@ -20,17 +20,20 @@
 
 /*!
  *  module_type {{{ utility }}}
- *  version {{{ 0.6.0 }}}
+ *  version {{{ 1.0.0 }}}
  *  requires {{{ bw_common }}}
  *  description {{{
  *    A collection of mathematical routines that strive to be better suited to
  *    DSP than, e.g., those supplied by your C standard library.
  * 
  *    Such a goal is hopefully accomplished by:
+ *
  *    * being as branchless as reasonable/convenient;
- *    * not handling uninteresting corner cases, such as invalid,
- *      infinity, and NaN input values;
- *    * returning approximated results (indicated in this documentation).
+ *    * not handling uninteresting corner cases, such as out-of-range,
+ *      NaN, and sometimes infinity input values (out-of-range and NaN inputs
+ *      are always considered invalid and lead to undefined behavior);
+ *    * returning approximated results (indicated in this documentation);
+ *    * making no distinction between `0.f` and `-0.f`.
  * 
  *    In practice they should guarantee fast and consistent performance, but
  *    always do your own benchmarking.
@@ -41,6 +44,12 @@
  *  }}}
  *  changelog {{{
  *    <ul>
+ *      <li>Version <strong>1.0.0</strong>:
+ *        <ul>
+ *          <li>WRITEME: document changes!!!</li>
+ *          <li>Avoided using reserved identifiers.</li>
+ *        </ul>
+ *      </li>
  *      <li>Version <strong>0.6.0</strong>:
  *        <ul>
  *          <li>Added debugging code.</li>
@@ -82,8 +91,8 @@
  *  }}}
  */
 
-#ifndef _BW_MATH_H
-#define _BW_MATH_H
+#ifndef BW_MATH_H
+#define BW_MATH_H
 
 #ifdef __cplusplus
 extern "C" {
@@ -154,22 +163,18 @@ static inline float bw_absf(float x);
 /*! <<<```
  *    Returns the absolute value of `x`.
  *
- *    #### bw_min0xf()
+ *    #### bw_min0f()
  *  ```>>> */
-static inline float bw_min0xf(float x);
+static inline float bw_min0f(float x);
 /*! <<<```
  *    Returns the minimum of `0.f` and `x`.
  *
- *    `x` must be finite.
- *
- *    #### bw_max0xf()
+ *    #### bw_max0f()
  *  ```>>> */
-static inline float bw_max0xf(float x);
+static inline float bw_max0f(float x);
 /*! <<<```
  *    Returns the maximum of `0.f` and `x`.
  *
- *    `x` must be finite.
- *
  *    #### bw_minf()
  *  ```>>> */
 static inline float bw_minf(float a, float b);
@@ -193,7 +198,7 @@ static inline float bw_clipf(float x, float m, float M);
  *    Returns `x` unless it is smaller than `m`, in which case it returns `m`,
  *    or bigger than `M`, in which case it returns `M`.
  *
- *    `x`, `m`, and `M` must be finite.
+ *    `M` must be greater than or equal to `m`.
  *
  *    #### bw_truncf()
  *  ```>>> */
@@ -202,6 +207,8 @@ static inline float bw_truncf(float x);
  *    Returns `x` with its fractional part set to zero (i.e., rounded towards
  *    zero).
  *
+ *    `x` must be finite.
+ *
  *    #### bw_roundf()
  *  ```>>> */
 static inline float bw_roundf(float x);
@@ -211,6 +218,8 @@ static inline float bw_roundf(float x);
  *    Halfway cases are rounded away from zero. E.g., `bw_roundf(0.5f)` gives
  *    `1.f` and `bw_roundf(-0.5f)` gives `-1.f`.
  *
+ *    `x` must be finite.
+ *
  *    #### bw_floorf()
  *  ```>>> */
 static inline float bw_floorf(float x);
@@ -218,6 +227,8 @@ static inline float bw_floorf(float x);
  *    Returns the biggest integer less or equal than `x` (i.e., `x` is rounded
  *    down).
  *
+ *    `x` must be finite.
+ *
  *    #### bw_ceilf()
  *  ```>>> */
 static inline float bw_ceilf(float x);
@@ -225,6 +236,8 @@ static inline float bw_ceilf(float x);
  *    Returns the smallest integer greater or equal than `x` (i.e., `x` is
  *    rounded up).
  *
+ *    `x` must be finite.
+ *
  *    #### bw_intfracf()
  *  ```>>> */
 static inline void bw_intfracf(float x, float *i, float *f);
@@ -427,12 +440,6 @@ static inline float bw_acoshf_3(float x);
 /* WARNING: This part of the file is not part of the public API. Its content may
  * change at any time in future versions. Please, do not use it directly. */
 
-typedef union {
-	float f;
-	int32_t i;
-	uint32_t u;
-} _bw_floatint;
-
 // I hope the target architecture and compiler will use conditional ops here
 
 static inline int32_t bw_signfilli32(int32_t x) {
@@ -467,87 +474,105 @@ static inline uint32_t bw_clipu32(uint32_t x, uint32_t m, uint32_t M) {
 // floating point numbers
 
 static inline float bw_copysignf(float x, float y) {
-	_bw_floatint v = {.f = x};
-	_bw_floatint s = {.f = y};
-	v.i = (v.i & 0x7fffffff) | (s.i & 0x80000000);
+	BW_ASSERT(!bw_is_nan(x));
+	BW_ASSERT(!bw_is_nan(y));
+	union { float f; uint32_t u; } v, s;
+	v.f = x;
+	s.f = y;
+	v.u = (v.u & 0x7fffffffu) | (s.u & 0x80000000u);
+	BW_ASSERT(!bw_is_nan(v.f));
 	return v.f;
 }
 
 static inline float bw_signf(float x) {
+	BW_ASSERT(!bw_is_nan(x));
 	static const float y[4] = { 0.f, 1.f, 0.f, -1.f };
-	_bw_floatint v = {.f = x};
-	return y[bw_minu32(v.u & 0x7fffffff, 1) | ((v.u >> 30) & 0x2)];
+	union { float f; uint32_t u; } v;
+	v.f = x;
+	const float r = y[bw_minu32(v.u & 0x7fffffffu, 1) | ((v.u >> 30) & 0x2)];
+	BW_ASSERT(!bw_is_nan(r));
+	return r;
 }
 
 static inline float bw_absf(float x) {
-	_bw_floatint v = {.f = x};
-	v.i = v.i & 0x7fffffff;
+	BW_ASSERT(!bw_is_nan(x));
+	union { float f; uint32_t u; } v;
+	v.f = x;
+	v.u = v.u & 0x7fffffffu;
+	BW_ASSERT(!bw_is_nan(v.f));
 	return v.f;
 }
 
-static inline float bw_min0xf(float x) {
-	BW_ASSERT(bw_is_finite(x));
-	const float y = 0.5f * (x - bw_absf(x));
-	BW_ASSERT(bw_is_finite(y));
-	return y;
+static inline float bw_min0f(float x) {
+	BW_ASSERT(!bw_is_nan(x));
+	union { float f; int32_t i; } v;
+	v.f = x;
+	v.i = bw_mini32(0, v.i);
+	BW_ASSERT(!bw_is_nan(v.f));
+	return v.f;
 }
 
-static inline float bw_max0xf(float x) {
-	BW_ASSERT(bw_is_finite(x));
-	const float y = 0.5f * (x + bw_absf(x));
-	BW_ASSERT(bw_is_finite(y));
-	return y;
+static inline float bw_max0f(float x) {
+	BW_ASSERT(!bw_is_nan(x));
+	union { float f; int32_t i; } v;
+	v.f = x;
+	v.i = bw_maxi32(0, v.i);
+	BW_ASSERT(!bw_is_nan(v.f));
+	return v.f;
 }
 
 static inline float bw_minf(float a, float b) {
-	BW_ASSERT(bw_is_finite(a));
-	BW_ASSERT(bw_is_finite(b));
-	const float y = a + bw_min0xf(b - a);
-	BW_ASSERT(bw_is_finite(y));
+	BW_ASSERT(!bw_is_nan(a));
+	BW_ASSERT(!bw_is_nan(b));
+	const float y = a < b ? a : b;
+	BW_ASSERT(!bw_is_nan(y));
 	return y;
 }
 
 static inline float bw_maxf(float a, float b) {
-	BW_ASSERT(bw_is_finite(a));
-	BW_ASSERT(bw_is_finite(b));
-	const float y = a + bw_max0xf(b - a);
-	BW_ASSERT(bw_is_finite(y));
+	BW_ASSERT(!bw_is_nan(a));
+	BW_ASSERT(!bw_is_nan(b));
+	const float y = a > b ? a : b;
+	BW_ASSERT(!bw_is_nan(y));
 	return y;
 }
 
 static inline float bw_clipf(float x, float m, float M) {
-	BW_ASSERT(bw_is_finite(x));
-	BW_ASSERT(bw_is_finite(m));
-	BW_ASSERT(bw_is_finite(M));
+	BW_ASSERT(!bw_is_nan(x));
+	BW_ASSERT(!bw_is_nan(m));
+	BW_ASSERT(!bw_is_nan(M));
+	BW_ASSERT(M >= m);
 	const float y = bw_minf(bw_maxf(x, m), M);
-	BW_ASSERT(bw_is_finite(y));
+	BW_ASSERT(!bw_is_nan(y));
 	return y;
 }
 
 static inline float bw_truncf(float x) {
 	BW_ASSERT(bw_is_finite(x));
-	_bw_floatint v = {.f = x};
-	int32_t ex = (v.i & 0x7f800000) >> 23;
+	union { float f; uint32_t u; } v;
+	v.f = x;
+	const int32_t ex = (v.u & 0x7f800000u) >> 23;
 	int32_t m = (~0u) << bw_clipi32(150 - ex, 0, 23);
 	m &= bw_signfilli32(126 - ex) | 0x80000000;
-	v.i &= m;
+	v.u &= m;
 	BW_ASSERT(bw_is_finite(v.f));
 	return v.f;
 }
 
 static inline float bw_roundf(float x) {
 	BW_ASSERT(bw_is_finite(x));
-	_bw_floatint v = {.f = x};
-	int32_t ex = (v.i & 0x7f800000) >> 23;
-	int32_t sh = bw_clipi32(150 - ex, 0, 23);
+	union { float f; uint32_t u; } v, s;
+	v.f = x;
+	const int32_t ex = (v.u & 0x7f800000u) >> 23;
+	const int32_t sh = bw_clipi32(150 - ex, 0, 24);
 	int32_t mt = (~0u) << sh;
 	mt &= bw_signfilli32(126 - ex) | 0x80000000;
 	int32_t mr = (1 << sh) >> 1;
 	mr &= bw_signfilli32(125 - ex);
-	_bw_floatint s = {.f = bw_copysignf(1.f, x)};
-	int32_t ms = bw_signfilli32((v.i & mr) << (32 - sh));
-	v.i &= mt;
-	s.i &= ms;
+	s.f = bw_copysignf(1.f, x);
+	int32_t ms = bw_signfilli32(((v.u | 0x00800000u) & mr) << (32 - sh));
+	v.u &= mt;
+	s.u &= ms;
 	const float y = v.f + s.f;
 	BW_ASSERT(bw_is_finite(y));
 	return y;
@@ -555,9 +580,10 @@ static inline float bw_roundf(float x) {
 
 static inline float bw_floorf(float x) {
 	BW_ASSERT(bw_is_finite(x));
-	_bw_floatint t = {.f = bw_truncf(x)}; // first bit set when t < 0
-	_bw_floatint y = {.f = x - t.f}; // first bit set when t > x
-	_bw_floatint s = {.f = 1.f};
+	union { float f; int32_t i; } t, y, s;
+	t.f = bw_truncf(x); // first bit set when t < 0
+	y.f = x - t.f; // first bit set when t > x
+	s.f = 1.f;
 	s.i &= bw_signfilli32(t.i & y.i);
 	const float r = t.f - s.f;
 	BW_ASSERT(bw_is_finite(r));
@@ -566,9 +592,10 @@ static inline float bw_floorf(float x) {
 
 static inline float bw_ceilf(float x) {
 	BW_ASSERT(bw_is_finite(x));
-	_bw_floatint t = {.f = bw_truncf(x)}; // first bit set when t < 0
-	_bw_floatint y = {.f = x - t.f}; // first bit set when t > x
-	_bw_floatint s = {.f = 1.f};
+	union { float f; int32_t i; } t, y, s;
+	t.f = bw_truncf(x); // first bit set when t < 0
+	y.f = t.f - x; // first bit set when t < x
+	s.f = 1.f;
 	s.i &= bw_signfilli32(~t.i & y.i);
 	const float r = t.f + s.f;
 	BW_ASSERT(bw_is_finite(r));
@@ -584,7 +611,8 @@ static inline void bw_intfracf(float x, float *i, float *f) {
 }
 
 static inline float bw_rcpf_2(float x) {
-	_bw_floatint v = {.f = x};
+	union { float f; int32_t i; } v;
+	v.f = x;
 	v.i = 0x7ef0e840 - v.i;
 	v.f = v.f + v.f - x * v.f * v.f;
 	v.f = v.f + v.f - x * v.f * v.f;
@@ -645,7 +673,8 @@ static inline float bw_tanf_3(float x) {
 static inline float bw_log2f_3(float x) {
 	BW_ASSERT(bw_is_finite(x));
 	BW_ASSERT(x >= 1.175494350822287e-38f);
-	_bw_floatint v = {.f = x};
+	union { float f; int32_t i; } v;
+	v.f = x;
 	int e = v.i >> 23;
 	v.i = (v.i & 0x007fffff) | 0x3f800000;
 	const float y = (float)e - 129.213475204444817f + v.f * (3.148297929334117f + v.f * (-1.098865286222744f + v.f * 0.1640425613334452f));
@@ -674,7 +703,8 @@ static inline float bw_pow2f_3(float x) {
 	BW_ASSERT(x <= 127.999f);
 	if (x < -126.f)
 		return 0.f;
-	_bw_floatint v = {.f = x};
+	union { float f; int32_t i; } v;
+	v.f = x;
 	int xi = (int)x;
 	int l = xi - ((v.i >> 31) & 1);
 	float f = x - (float)l;
@@ -724,7 +754,8 @@ static inline float bw_sqrtf_2(float x) {
 		BW_ASSERT(bw_is_finite(y));
 		return y;
 	}
-	_bw_floatint v = {.f = x};
+	union { float f; uint32_t u; } v;
+	v.f = x;
 	v.u = (((v.u - 0x3f82a127) >> 1) + 0x3f7d8fc7) & 0x7fffffff;
 	float r = bw_rcpf_2(x);
 	v.f = v.f + v.f * (0.5f - 0.5f * r * v.f * v.f);
diff --git a/test/Makefile b/test/Makefile
new file mode 100644
index 0000000..0916997
--- /dev/null
+++ b/test/Makefile
@@ -0,0 +1,16 @@
+CC := gcc
+CFLAGS := -I../include
+
+all: build/bw_math
+	./build/bw_math
+
+build/bw_math: bw_math.c ../include/bw_math.h | build
+	${CC} ${CFLAGS} $^ -o $@
+
+build:
+	mkdir -p $@
+
+clean:
+	rm -fr build/
+
+.PHONY: all
diff --git a/test/bw_math.c b/test/bw_math.c
new file mode 100644
index 0000000..fc60231
--- /dev/null
+++ b/test/bw_math.c
@@ -0,0 +1,536 @@
+#include <stdlib.h>
+#include <stdio.h>
+
+#include <bw_math.h>
+
+int n_ok = 0;
+int n_ko = 0;
+
+#define TEST_F(expr, expected) \
+{ \
+	union { float f; uint32_t u; } v_expr, v_expected; \
+	v_expr.f = expr; \
+	v_expected.f = expected; \
+	if (expr == v_expected.f) { \
+		printf("[OK] %s = %g\n", #expr, v_expected.f); \
+		n_ok++; \
+	} else { \
+		printf("[KO] %s = %g [0x%x] (expected %g [0x%x]) - line %d\n", #expr, v_expr.f, v_expr.u, v_expected.f, v_expected.u, __LINE__); \
+		n_ko++; \
+	} \
+}
+
+#define TEST_BW_INTFRAC(val, expected_i, expected_f) \
+{ \
+	union { float f; uint32_t u; } v_res_i, v_res_f, v_expected_i, v_expected_f; \
+	v_expected_i.f = expected_i; \
+	v_expected_f.f = expected_f; \
+	bw_intfracf(val, &v_res_i.f, &v_res_f.f); \
+	if (v_res_i.f == v_expected_i.f && v_res_f.f == v_expected_f.f) { \
+		printf("[OK] bw_intfracf(%g) = %g, %g\n", val, expected_i, expected_f); \
+		n_ok++; \
+	} else { \
+		printf("[KO] bw_intfracf(%g) = %g [0x%x], %g [0x%x] (expected %g [0x%x], %g [0x%x]) - line %d\n", val, v_res_i.f, v_res_i.u, v_res_f.f, v_res_f.u, v_expected_i.f, v_expected_i.u, v_expected_f.f, v_expected_f.u, __LINE__); \
+		n_ko++; \
+	} \
+}
+
+int main() {
+	printf("\nbw_math unit tests\n");
+	printf("------------------\n\n");
+
+	TEST_F(bw_copysignf(INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_copysignf(INFINITY, -INFINITY), -INFINITY);
+	TEST_F(bw_copysignf(INFINITY, 1e3f), INFINITY);
+	TEST_F(bw_copysignf(INFINITY, -1e3f), -INFINITY);
+	TEST_F(bw_copysignf(INFINITY, 1.f), INFINITY);
+	TEST_F(bw_copysignf(INFINITY, -1.f), -INFINITY);
+	TEST_F(bw_copysignf(INFINITY, 1e-3f), INFINITY);
+	TEST_F(bw_copysignf(INFINITY, -1e-3f), -INFINITY);
+	TEST_F(bw_copysignf(INFINITY, 0.f), INFINITY);
+	TEST_F(bw_copysignf(INFINITY, -0.f), -INFINITY);
+	TEST_F(bw_copysignf(1e3f, INFINITY), 1e3f);
+	TEST_F(bw_copysignf(1e3f, -INFINITY), -1e3f);
+	TEST_F(bw_copysignf(1e3f, 1e3f), 1e3f);
+	TEST_F(bw_copysignf(1e3f, -1e3f), -1e3f);
+	TEST_F(bw_copysignf(1e3f, 1.f), 1e3f);
+	TEST_F(bw_copysignf(1e3f, -1.f), -1e3f);
+	TEST_F(bw_copysignf(1e3f, 1e-3f), 1e3f);
+	TEST_F(bw_copysignf(1e3f, -1e-3f), -1e3f);
+	TEST_F(bw_copysignf(1e3f, 0.f), 1e3f);
+	TEST_F(bw_copysignf(1e3f, -0.f), -1e3f);
+	TEST_F(bw_copysignf(1.f, INFINITY), 1.f);
+	TEST_F(bw_copysignf(1.f, -INFINITY), -1.f);
+	TEST_F(bw_copysignf(1.f, 1e3f), 1.f);
+	TEST_F(bw_copysignf(1.f, -1e3f), -1.f);
+	TEST_F(bw_copysignf(1.f, 1.f), 1.f);
+	TEST_F(bw_copysignf(1.f, -1.f), -1.f);
+	TEST_F(bw_copysignf(1.f, 1e-3f), 1.f);
+	TEST_F(bw_copysignf(1.f, -1e-3f), -1.f);
+	TEST_F(bw_copysignf(1.f, 0.f), 1.f);
+	TEST_F(bw_copysignf(1.f, -0.f), -1.f);
+	TEST_F(bw_copysignf(1e-3f, INFINITY), 1e-3f);
+	TEST_F(bw_copysignf(1e-3f, -INFINITY), -1e-3f);
+	TEST_F(bw_copysignf(1e-3f, 1e3f), 1e-3f);
+	TEST_F(bw_copysignf(1e-3f, -1e3f), -1e-3f);
+	TEST_F(bw_copysignf(1e-3f, 1.f), 1e-3f);
+	TEST_F(bw_copysignf(1e-3f, -1.f), -1e-3f);
+	TEST_F(bw_copysignf(1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_copysignf(1e-3f, -1e-3f), -1e-3f);
+	TEST_F(bw_copysignf(1e-3f, 0.f), 1e-3f);
+	TEST_F(bw_copysignf(1e-3f, -0.f), -1e-3f);
+	TEST_F(bw_copysignf(0.f, INFINITY), 0.f);
+	TEST_F(bw_copysignf(0.f, -INFINITY), -0.f);
+	TEST_F(bw_copysignf(0.f, 1e3f), 0.f);
+	TEST_F(bw_copysignf(0.f, -1e3f), -0.f);
+	TEST_F(bw_copysignf(0.f, 1.f), 0.f);
+	TEST_F(bw_copysignf(0.f, -1.f), -0.f);
+	TEST_F(bw_copysignf(0.f, 1e-3f), 0.f);
+	TEST_F(bw_copysignf(0.f, -1e-3f), -0.f);
+	TEST_F(bw_copysignf(0.f, 0.f), 0.f);
+	TEST_F(bw_copysignf(0.f, -0.f), -0.f);
+	TEST_F(bw_copysignf(-INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, -INFINITY), -INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, 1e3f), INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, -1e3f), -INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, 1.f), INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, -1.f), -INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, 1e-3f), INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, -1e-3f), -INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, 0.f), INFINITY);
+	TEST_F(bw_copysignf(-INFINITY, -0.f), -INFINITY);
+	TEST_F(bw_copysignf(-1e3f, INFINITY), 1e3f);
+	TEST_F(bw_copysignf(-1e3f, -INFINITY), -1e3f);
+	TEST_F(bw_copysignf(-1e3f, 1e3f), 1e3f);
+	TEST_F(bw_copysignf(-1e3f, -1e3f), -1e3f);
+	TEST_F(bw_copysignf(-1e3f, 1.f), 1e3f);
+	TEST_F(bw_copysignf(-1e3f, -1.f), -1e3f);
+	TEST_F(bw_copysignf(-1e3f, 1e-3f), 1e3f);
+	TEST_F(bw_copysignf(-1e3f, -1e-3f), -1e3f);
+	TEST_F(bw_copysignf(-1e3f, 0.f), 1e3f);
+	TEST_F(bw_copysignf(-1e3f, -0.f), -1e3f);
+	TEST_F(bw_copysignf(-1.f, INFINITY), 1.f);
+	TEST_F(bw_copysignf(-1.f, -INFINITY), -1.f);
+	TEST_F(bw_copysignf(-1.f, 1e3f), 1.f);
+	TEST_F(bw_copysignf(-1.f, -1e3f), -1.f);
+	TEST_F(bw_copysignf(-1.f, 1.f), 1.f);
+	TEST_F(bw_copysignf(-1.f, -1.f), -1.f);
+	TEST_F(bw_copysignf(-1.f, 1e-3f), 1.f);
+	TEST_F(bw_copysignf(-1.f, -1e-3f), -1.f);
+	TEST_F(bw_copysignf(-1.f, 0.f), 1.f);
+	TEST_F(bw_copysignf(-1.f, -0.f), -1.f);
+	TEST_F(bw_copysignf(-1e-3f, INFINITY), 1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, -INFINITY), -1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, 1e3f), 1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, -1e3f), -1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, 1.f), 1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, -1.f), -1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, -1e-3f), -1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, 0.f), 1e-3f);
+	TEST_F(bw_copysignf(-1e-3f, -0.f), -1e-3f);
+	TEST_F(bw_copysignf(-0.f, INFINITY), 0.f);
+	TEST_F(bw_copysignf(-0.f, -INFINITY), -0.f);
+	TEST_F(bw_copysignf(-0.f, 1e3f), 0.f);
+	TEST_F(bw_copysignf(-0.f, -1e3f), -0.f);
+	TEST_F(bw_copysignf(-0.f, 1.f), 0.f);
+	TEST_F(bw_copysignf(-0.f, -1.f), -0.f);
+	TEST_F(bw_copysignf(-0.f, 1e-3f), 0.f);
+	TEST_F(bw_copysignf(-0.f, -1e-3f), -0.f);
+	TEST_F(bw_copysignf(-0.f, 0.f), 0.f);
+	TEST_F(bw_copysignf(-0.f, -0.f), -0.f);
+	
+	TEST_F(bw_signf(INFINITY), 1.f);
+	TEST_F(bw_signf(-INFINITY), -1.f);
+	TEST_F(bw_signf(1e3f), 1.f);
+	TEST_F(bw_signf(-1e3f), -1.f);
+	TEST_F(bw_signf(1.f), 1.f);
+	TEST_F(bw_signf(-1.f), -1.f);
+	TEST_F(bw_signf(1e-3f), 1.f);
+	TEST_F(bw_signf(-1e-3f), -1.f);
+	TEST_F(bw_signf(0.f), 0.f);
+	TEST_F(bw_signf(-0.f), 0.f);
+	
+	TEST_F(bw_absf(INFINITY), INFINITY);
+	TEST_F(bw_absf(-INFINITY), INFINITY);
+	TEST_F(bw_absf(1e3f), 1e3f);
+	TEST_F(bw_absf(-1e3f), 1e3f);
+	TEST_F(bw_absf(1.f), 1.f);
+	TEST_F(bw_absf(-1.f), 1.f);
+	TEST_F(bw_absf(1e-3f), 1e-3f);
+	TEST_F(bw_absf(-1e-3f), 1e-3f);
+	TEST_F(bw_absf(0.f), 0.f);
+	TEST_F(bw_absf(-0.f), 0.f);
+	
+	TEST_F(bw_min0f(INFINITY), 0.f);
+	TEST_F(bw_min0f(-INFINITY), -INFINITY);
+	TEST_F(bw_min0f(1e3f), 0.f);
+	TEST_F(bw_min0f(-1e3f), -1e3f);
+	TEST_F(bw_min0f(1.f), 0.f);
+	TEST_F(bw_min0f(-1.f), -1.f);
+	TEST_F(bw_min0f(1e-3f), 0.f);
+	TEST_F(bw_min0f(-1e-3f), -1e-3f);
+	TEST_F(bw_min0f(0.f), 0.f);
+	TEST_F(bw_min0f(-0.f), -0.f);
+	
+	TEST_F(bw_max0f(INFINITY), INFINITY);
+	TEST_F(bw_max0f(-INFINITY), 0.f);
+	TEST_F(bw_max0f(1e3f), 1e3f);
+	TEST_F(bw_max0f(-1e3f), 0.f);
+	TEST_F(bw_max0f(1.f), 1.f);
+	TEST_F(bw_max0f(-1.f), 0.f);
+	TEST_F(bw_max0f(1e-3f), 1e-3f);
+	TEST_F(bw_max0f(-1e-3f), 0.f);
+	TEST_F(bw_max0f(0.f), 0.f);
+	TEST_F(bw_max0f(-0.f), 0.f);
+	
+	TEST_F(bw_minf(INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_minf(INFINITY, 1e3f), 1e3f);
+	TEST_F(bw_minf(INFINITY, 1.f), 1.f);
+	TEST_F(bw_minf(INFINITY, 1e-3f), 1e-3f);
+	TEST_F(bw_minf(INFINITY, 0.f), 0.f);
+	TEST_F(bw_minf(INFINITY, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(INFINITY, -1.f), -1.f);
+	TEST_F(bw_minf(INFINITY, -1e3f), -1e3f);
+	TEST_F(bw_minf(INFINITY, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(1e3f, INFINITY), 1e3f);
+	TEST_F(bw_minf(1e3f, 1e3f), 1e3f);
+	TEST_F(bw_minf(1e3f, 1.f), 1.f);
+	TEST_F(bw_minf(1e3f, 1e-3f), 1e-3f);
+	TEST_F(bw_minf(1e3f, 0.f), 0.f);
+	TEST_F(bw_minf(1e3f, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(1e3f, -1.f), -1.f);
+	TEST_F(bw_minf(1e3f, -1e3f), -1e3f);
+	TEST_F(bw_minf(1e3f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(1.f, INFINITY), 1.f);
+	TEST_F(bw_minf(1.f, 1e3f), 1.f);
+	TEST_F(bw_minf(1.f, 1.f), 1.f);
+	TEST_F(bw_minf(1.f, 1e-3f), 1e-3f);
+	TEST_F(bw_minf(1.f, 0.f), 0.f);
+	TEST_F(bw_minf(1.f, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(1.f, -1.f), -1.f);
+	TEST_F(bw_minf(1.f, -1e3f), -1e3f);
+	TEST_F(bw_minf(1.f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(1e-3f, INFINITY), 1e-3f);
+	TEST_F(bw_minf(1e-3f, 1e3f), 1e-3f);
+	TEST_F(bw_minf(1e-3f, 1.f), 1e-3f);
+	TEST_F(bw_minf(1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_minf(1e-3f, 0.f), 0.f);
+	TEST_F(bw_minf(1e-3f, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(1e-3f, -1.f), -1.f);
+	TEST_F(bw_minf(1e-3f, -1e3f), -1e3f);
+	TEST_F(bw_minf(1e-3f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(0.f, INFINITY), 0.f);
+	TEST_F(bw_minf(0.f, 1e3f), 0.f);
+	TEST_F(bw_minf(0.f, 1.f), 0.f);
+	TEST_F(bw_minf(0.f, 1e-3f), 0.f);
+	TEST_F(bw_minf(0.f, 0.f), 0.f);
+	TEST_F(bw_minf(0.f, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(0.f, -1.f), -1.f);
+	TEST_F(bw_minf(0.f, -1e3f), -1e3f);
+	TEST_F(bw_minf(0.f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(-1e-3f, INFINITY), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, 1e3f), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, 1.f), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, 1e-3f), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, 0.f), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, -1e-3f), -1e-3f);
+	TEST_F(bw_minf(-1e-3f, -1.f), -1.f);
+	TEST_F(bw_minf(-1e-3f, -1e3f), -1e3f);
+	TEST_F(bw_minf(-1e-3f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(-1.f, INFINITY), -1.f);
+	TEST_F(bw_minf(-1.f, 1e3f), -1.f);
+	TEST_F(bw_minf(-1.f, 1.f), -1.f);
+	TEST_F(bw_minf(-1.f, 1e-3f), -1.f);
+	TEST_F(bw_minf(-1.f, 0.f), -1.f);
+	TEST_F(bw_minf(-1.f, -1e-3f), -1.f);
+	TEST_F(bw_minf(-1.f, -1.f), -1.f);
+	TEST_F(bw_minf(-1.f, -1e3f), -1e3f);
+	TEST_F(bw_minf(-1.f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(-1e3f, INFINITY), -1e3f);
+	TEST_F(bw_minf(-1e3f, 1e3f), -1e3f);
+	TEST_F(bw_minf(-1e3f, 1.f), -1e3f);
+	TEST_F(bw_minf(-1e3f, 1e-3f), -1e3f);
+	TEST_F(bw_minf(-1e3f, 0.f), -1e3f);
+	TEST_F(bw_minf(-1e3f, -1e-3f), -1e3f);
+	TEST_F(bw_minf(-1e3f, -1.f), -1e3f);
+	TEST_F(bw_minf(-1e3f, -1e3f), -1e3f);
+	TEST_F(bw_minf(-1e3f, -INFINITY), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, INFINITY), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, 1e3f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, 1.f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, 1e-3f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, 0.f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, -1e-3f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, -1.f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, -1e3f), -INFINITY);
+	TEST_F(bw_minf(-INFINITY, -INFINITY), -INFINITY);
+	
+	TEST_F(bw_maxf(INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_maxf(INFINITY, 1e3f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, 1.f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, 1e-3f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, 0.f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, -1e-3f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, -1.f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, -1e3f), INFINITY);
+	TEST_F(bw_maxf(INFINITY, -INFINITY), INFINITY);
+	TEST_F(bw_maxf(1e3f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(1e3f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(1e3f, 1.f), 1e3f);
+	TEST_F(bw_maxf(1e3f, 1e-3f), 1e3f);
+	TEST_F(bw_maxf(1e3f, 0.f), 1e3f);
+	TEST_F(bw_maxf(1e3f, -1e-3f), 1e3f);
+	TEST_F(bw_maxf(1e3f, -1.f), 1e3f);
+	TEST_F(bw_maxf(1e3f, -1e3f), 1e3f);
+	TEST_F(bw_maxf(1e3f, -INFINITY), 1e3f);
+	TEST_F(bw_maxf(1.f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(1.f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(1.f, 1.f), 1.f);
+	TEST_F(bw_maxf(1.f, 1e-3f), 1.f);
+	TEST_F(bw_maxf(1.f, 0.f), 1.f);
+	TEST_F(bw_maxf(1.f, -1e-3f), 1.f);
+	TEST_F(bw_maxf(1.f, -1.f), 1.f);
+	TEST_F(bw_maxf(1.f, -1e3f), 1.f);
+	TEST_F(bw_maxf(1.f, -INFINITY), 1.f);
+	TEST_F(bw_maxf(1e-3f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(1e-3f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(1e-3f, 1.f), 1.f);
+	TEST_F(bw_maxf(1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(1e-3f, 0.f), 1e-3f);
+	TEST_F(bw_maxf(1e-3f, -1e-3f), 1e-3f);
+	TEST_F(bw_maxf(1e-3f, -1.f), 1e-3f);
+	TEST_F(bw_maxf(1e-3f, -1e3f), 1e-3f);
+	TEST_F(bw_maxf(1e-3f, -INFINITY), 1e-3f);
+	TEST_F(bw_maxf(0.f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(0.f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(0.f, 1.f), 1.f);
+	TEST_F(bw_maxf(0.f, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(0.f, 0.f), 0.f);
+	TEST_F(bw_maxf(0.f, -1e-3f), 0.f);
+	TEST_F(bw_maxf(0.f, -1.f), 0.f);
+	TEST_F(bw_maxf(0.f, -1e3f), 0.f);
+	TEST_F(bw_maxf(0.f, -INFINITY), 0.f);
+	TEST_F(bw_maxf(-1e-3f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(-1e-3f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(-1e-3f, 1.f), 1.f);
+	TEST_F(bw_maxf(-1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(-1e-3f, 0.f), 0.f);
+	TEST_F(bw_maxf(-1e-3f, -1e-3f), -1e-3f);
+	TEST_F(bw_maxf(-1e-3f, -1.f), -1e-3f);
+	TEST_F(bw_maxf(-1e-3f, -1e3f), -1e-3f);
+	TEST_F(bw_maxf(-1e-3f, -INFINITY), -1e-3f);
+	TEST_F(bw_maxf(-1.f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(-1.f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(-1.f, 1.f), 1.f);
+	TEST_F(bw_maxf(-1.f, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(-1.f, 0.f), 0.f);
+	TEST_F(bw_maxf(-1.f, -1e-3f), -1e-3f);
+	TEST_F(bw_maxf(-1.f, -1.f), -1.f);
+	TEST_F(bw_maxf(-1.f, -1e3f), -1.f);
+	TEST_F(bw_maxf(-1.f, -INFINITY), -1.f);
+	TEST_F(bw_maxf(-1e3f, INFINITY), INFINITY);
+	TEST_F(bw_maxf(-1e3f, 1e3f), 1e3f);
+	TEST_F(bw_maxf(-1e3f, 1.f), 1.f);
+	TEST_F(bw_maxf(-1e3f, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(-1e3f, 0.f), 0.f);
+	TEST_F(bw_maxf(-1e3f, -1e-3f), -1e-3f);
+	TEST_F(bw_maxf(-1e3f, -1.f), -1.f);
+	TEST_F(bw_maxf(-1e3f, -1e3f), -1e3f);
+	TEST_F(bw_maxf(-1e3f, -INFINITY), -1e3f);
+	TEST_F(bw_maxf(-INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_maxf(-INFINITY, 1e3f), 1e3f);
+	TEST_F(bw_maxf(-INFINITY, 1.f), 1.f);
+	TEST_F(bw_maxf(-INFINITY, 1e-3f), 1e-3f);
+	TEST_F(bw_maxf(-INFINITY, 0.f), 0.f);
+	TEST_F(bw_maxf(-INFINITY, -1e-3f), -1e-3f);
+	TEST_F(bw_maxf(-INFINITY, -1.f), -1.f);
+	TEST_F(bw_maxf(-INFINITY, -1e3f), -1e3f);
+	TEST_F(bw_maxf(-INFINITY, -INFINITY), -INFINITY);
+	
+	TEST_F(bw_clipf(INFINITY, -INFINITY, INFINITY), INFINITY);
+	TEST_F(bw_clipf(INFINITY, -1e3f, 1e3f), 1e3f);
+	TEST_F(bw_clipf(INFINITY, -1.f, 1.f), 1.f);
+	TEST_F(bw_clipf(INFINITY, -1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_clipf(INFINITY, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(1e3f, -INFINITY, INFINITY), 1e3f);
+	TEST_F(bw_clipf(1e3f, -1e3f, 1e3f), 1e3f);
+	TEST_F(bw_clipf(1e3f, -1.f, 1.f), 1.f);
+	TEST_F(bw_clipf(1e3f, -1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_clipf(1e3f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(1.f, -INFINITY, INFINITY), 1.f);
+	TEST_F(bw_clipf(1.f, -1e3f, 1e3f), 1.f);
+	TEST_F(bw_clipf(1.f, -1.f, 1.f), 1.f);
+	TEST_F(bw_clipf(1.f, -1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_clipf(1.f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(1e-3f, -INFINITY, INFINITY), 1e-3f);
+	TEST_F(bw_clipf(1e-3f, -1e3f, 1e3f), 1e-3f);
+	TEST_F(bw_clipf(1e-3f, -1.f, 1.f), 1e-3f);
+	TEST_F(bw_clipf(1e-3f, -1e-3f, 1e-3f), 1e-3f);
+	TEST_F(bw_clipf(1e-3f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(0.f, -INFINITY, INFINITY), 0.f);
+	TEST_F(bw_clipf(0.f, -1e3f, 1e3f), 0.f);
+	TEST_F(bw_clipf(0.f, -1.f, 1.f), 0.f);
+	TEST_F(bw_clipf(0.f, -1e-3f, 1e-3f), 0.f);
+	TEST_F(bw_clipf(0.f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(-1e-3f, -INFINITY, INFINITY), -1e-3f);
+	TEST_F(bw_clipf(-1e-3f, -1e3f, 1e3f), -1e-3f);
+	TEST_F(bw_clipf(-1e-3f, -1.f, 1.f), -1e-3f);
+	TEST_F(bw_clipf(-1e-3f, -1e-3f, 1e-3f), -1e-3f);
+	TEST_F(bw_clipf(-1e-3f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(-1.f, -INFINITY, INFINITY), -1.f);
+	TEST_F(bw_clipf(-1.f, -1e3f, 1e3f), -1.f);
+	TEST_F(bw_clipf(-1.f, -1.f, 1.f), -1.f);
+	TEST_F(bw_clipf(-1.f, -1e-3f, 1e-3f), -1e-3f);
+	TEST_F(bw_clipf(-1.f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(-1e3f, -INFINITY, INFINITY), -1e3f);
+	TEST_F(bw_clipf(-1e3f, -1e3f, 1e3f), -1e3f);
+	TEST_F(bw_clipf(-1e3f, -1.f, 1.f), -1.f);
+	TEST_F(bw_clipf(-1e3f, -1e-3f, 1e-3f), -1e-3f);
+	TEST_F(bw_clipf(-1e3f, 0.f, 0.f), 0.f);
+	TEST_F(bw_clipf(-INFINITY, -INFINITY, INFINITY), -INFINITY);
+	TEST_F(bw_clipf(-INFINITY, -1e3f, 1e3f), -1e3f);
+	TEST_F(bw_clipf(-INFINITY, -1.f, 1.f), -1.f);
+	TEST_F(bw_clipf(-INFINITY, -1e-3f, 1e-3f), -1e-3f);
+	TEST_F(bw_clipf(-INFINITY, 0.f, 0.f), 0.f);
+	
+	TEST_F(bw_truncf(1.234e38f), 1.234e38f);
+	TEST_F(bw_truncf(1001.f), 1001.f);
+	TEST_F(bw_truncf(1000.9f), 1000.f);
+	TEST_F(bw_truncf(1000.5f), 1000.f);
+	TEST_F(bw_truncf(1000.1f), 1000.f);
+	TEST_F(bw_truncf(1000.0f), 1000.f);
+	TEST_F(bw_truncf(999.9f), 999.f);
+	TEST_F(bw_truncf(999.5f), 999.f);
+	TEST_F(bw_truncf(999.1f), 999.f);
+	TEST_F(bw_truncf(999.f), 999.f);
+	TEST_F(bw_truncf(1.5f), 1.f);
+	TEST_F(bw_truncf(1.f), 1.f);
+	TEST_F(bw_truncf(0.9f), 0.f);
+	TEST_F(bw_truncf(0.5f), 0.f);
+	TEST_F(bw_truncf(0.1f), 0.f);
+	TEST_F(bw_truncf(0.0f), 0.f);
+	TEST_F(bw_truncf(-0.0f), 0.f);
+	TEST_F(bw_truncf(-0.1f), 0.f);
+	TEST_F(bw_truncf(-0.5f), 0.f);
+	TEST_F(bw_truncf(-0.9f), 0.f);
+	TEST_F(bw_truncf(-1.f), -1.f);
+	TEST_F(bw_truncf(-1.5f), -1.f);
+	TEST_F(bw_truncf(-999.f), -999.f);
+	TEST_F(bw_truncf(-999.1f), -999.f);
+	TEST_F(bw_truncf(-999.5f), -999.f);
+	TEST_F(bw_truncf(-999.9f), -999.f);
+	TEST_F(bw_truncf(-1000.0f), -1000.f);
+	TEST_F(bw_truncf(-1000.1f), -1000.f);
+	TEST_F(bw_truncf(-1000.5f), -1000.f);
+	TEST_F(bw_truncf(-1000.9f), -1000.f);
+	TEST_F(bw_truncf(-1001.f), -1001.f);
+	TEST_F(bw_truncf(-1.234e38f), -1.234e38f);
+	
+	TEST_F(bw_roundf(1.234e38f), 1.234e38f);
+	TEST_F(bw_roundf(1001.f), 1001.f);
+	TEST_F(bw_roundf(1000.9f), 1001.f);
+	TEST_F(bw_roundf(1000.5f), 1001.f);
+	TEST_F(bw_roundf(1000.1f), 1000.f);
+	TEST_F(bw_roundf(1000.0f), 1000.f);
+	TEST_F(bw_roundf(999.9f), 1000.f);
+	TEST_F(bw_roundf(999.5f), 1000.f);
+	TEST_F(bw_roundf(999.1f), 999.f);
+	TEST_F(bw_roundf(999.f), 999.f);
+	TEST_F(bw_roundf(1.5f), 2.f);
+	TEST_F(bw_roundf(1.f), 1.f);
+	TEST_F(bw_roundf(0.9f), 1.f);
+	TEST_F(bw_roundf(0.5f), 1.f);
+	TEST_F(bw_roundf(0.1f), 0.f);
+	TEST_F(bw_roundf(0.0f), 0.f);
+	TEST_F(bw_roundf(-0.0f), 0.f);
+	TEST_F(bw_roundf(-0.1f), 0.f);
+	TEST_F(bw_roundf(-0.5f), -1.f);
+	TEST_F(bw_roundf(-0.9f), -1.f);
+	TEST_F(bw_roundf(-1.f), -1.f);
+	TEST_F(bw_roundf(-1.5f), -2.f);
+	TEST_F(bw_roundf(-999.f), -999.f);
+	TEST_F(bw_roundf(-999.1f), -999.f);
+	TEST_F(bw_roundf(-999.5f), -1000.f);
+	TEST_F(bw_roundf(-999.9f), -1000.f);
+	TEST_F(bw_roundf(-1000.0f), -1000.f);
+	TEST_F(bw_roundf(-1000.1f), -1000.f);
+	TEST_F(bw_roundf(-1000.5f), -1001.f);
+	TEST_F(bw_roundf(-1000.9f), -1001.f);
+	TEST_F(bw_roundf(-1001.f), -1001.f);
+	TEST_F(bw_roundf(-1.234e38f), -1.234e38f);
+	
+	TEST_F(bw_floorf(1.234e38f), 1.234e38f);
+	TEST_F(bw_floorf(1001.f), 1001.f);
+	TEST_F(bw_floorf(1000.9f), 1000.f);
+	TEST_F(bw_floorf(1000.5f), 1000.f);
+	TEST_F(bw_floorf(1000.1f), 1000.f);
+	TEST_F(bw_floorf(1000.0f), 1000.f);
+	TEST_F(bw_floorf(999.9f), 999.f);
+	TEST_F(bw_floorf(999.5f), 999.f);
+	TEST_F(bw_floorf(999.1f), 999.f);
+	TEST_F(bw_floorf(999.f), 999.f);
+	TEST_F(bw_floorf(1.5f), 1.f);
+	TEST_F(bw_floorf(1.f), 1.f);
+	TEST_F(bw_floorf(0.9f), 0.f);
+	TEST_F(bw_floorf(0.5f), 0.f);
+	TEST_F(bw_floorf(0.1f), 0.f);
+	TEST_F(bw_floorf(0.0f), 0.f);
+	TEST_F(bw_floorf(-0.0f), 0.f);
+	TEST_F(bw_floorf(-0.1f), -1.f);
+	TEST_F(bw_floorf(-0.5f), -1.f);
+	TEST_F(bw_floorf(-0.9f), -1.f);
+	TEST_F(bw_floorf(-1.f), -1.f);
+	TEST_F(bw_floorf(-1.5f), -2.f);
+	TEST_F(bw_floorf(-999.f), -999.f);
+	TEST_F(bw_floorf(-999.1f), -1000.f);
+	TEST_F(bw_floorf(-999.5f), -1000.f);
+	TEST_F(bw_floorf(-999.9f), -1000.f);
+	TEST_F(bw_floorf(-1000.0f), -1000.f);
+	TEST_F(bw_floorf(-1000.1f), -1001.f);
+	TEST_F(bw_floorf(-1000.5f), -1001.f);
+	TEST_F(bw_floorf(-1000.9f), -1001.f);
+	TEST_F(bw_floorf(-1001.f), -1001.f);
+	TEST_F(bw_floorf(-1.234e38f), -1.234e38f);
+	
+	TEST_F(bw_ceilf(1.234e38f), 1.234e38f);
+	TEST_F(bw_ceilf(1001.f), 1001.f);
+	TEST_F(bw_ceilf(1000.9f), 1001.f);
+	TEST_F(bw_ceilf(1000.5f), 1001.f);
+	TEST_F(bw_ceilf(1000.1f), 1001.f);
+	TEST_F(bw_ceilf(1000.0f), 1000.f);
+	TEST_F(bw_ceilf(999.9f), 1000.f);
+	TEST_F(bw_ceilf(999.5f), 1000.f);
+	TEST_F(bw_ceilf(999.1f), 1000.f);
+	TEST_F(bw_ceilf(999.f), 999.f);
+	TEST_F(bw_ceilf(1.5f), 2.f);
+	TEST_F(bw_ceilf(1.f), 1.f);
+	TEST_F(bw_ceilf(0.9f), 1.f);
+	TEST_F(bw_ceilf(0.5f), 1.f);
+	TEST_F(bw_ceilf(0.1f), 1.f);
+	TEST_F(bw_ceilf(0.0f), 0.f);
+	TEST_F(bw_ceilf(-0.0f), 0.f);
+	TEST_F(bw_ceilf(-0.1f), 0.f);
+	TEST_F(bw_ceilf(-0.5f), 0.f);
+	TEST_F(bw_ceilf(-0.9f), 0.f);
+	TEST_F(bw_ceilf(-1.f), -1.f);
+	TEST_F(bw_ceilf(-1.5f), -1.f);
+	TEST_F(bw_ceilf(-999.f), -999.f);
+	TEST_F(bw_ceilf(-999.1f), -999.f);
+	TEST_F(bw_ceilf(-999.5f), -999.f);
+	TEST_F(bw_ceilf(-999.9f), -999.f);
+	TEST_F(bw_ceilf(-1000.0f), -1000.f);
+	TEST_F(bw_ceilf(-1000.1f), -1000.f);
+	TEST_F(bw_ceilf(-1000.5f), -1000.f);
+	TEST_F(bw_ceilf(-1000.9f), -1000.f);
+	TEST_F(bw_ceilf(-1001.f), -1001.f);
+	TEST_F(bw_ceilf(-1.234e38f), -1.234e38f);
+	
+	TEST_BW_INTFRAC(1.f, 1.f, 0.f);
+	TEST_BW_INTFRAC(1.123f, 1.f, (1.123f - 1.f));
+	TEST_BW_INTFRAC(1.999f, 1.f, (1.999f - 1.f));
+
+	printf("\nsuceeded: %d, failed: %d\n\n", n_ok, n_ko);
+	
+	return n_ko ? EXIT_FAILURE : EXIT_SUCCESS;
+}