1 files changed, 124 insertions, 0 deletions
diff --git a/test/test-fp_Infcases.c b/test/test-fp_Infcases.c
new file mode 100644
index 00000000..88a89ead
--- /dev/null
+++ b/test/test-fp_Infcases.c
@@ -0,0 +1,124 @@
+/* test cases for floating point comparison transformations
+ * compile with -DFLOAT_TYPE=float
+ *          or  -DFLOAT_TYPE=double
+ *          or  -DFLOAT_TYPE="long double"
+ */
+
+#include <assert.h>
+#define _GNU_SOURCE
+#include <math.h> /* for NaNs and infinity values */
+
+int main() {
+
+  volatile FLOAT_TYPE a, b;
+
+#ifdef INFINITY
+  FLOAT_TYPE inf = (FLOAT_TYPE) INFINITY;
+#else
+  FLOAT_TYPE inf = 1.0 / 0.0; /* produces infinity */
+#endif
+  FLOAT_TYPE negZero = 1.0 / -inf;
+  FLOAT_TYPE posZero = 0.0;
+
+  /* plus infinity */
+  a = (1.0 / 0.0);	/* positive infinity */
+  b = (1.0 / 0.0);	/* positive infinity */
+  assert(!(a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert((a >= b));
+  assert(!(a != b));
+  assert((a == b));
+
+  b = -(1.0 / 0.0);	/* negative infinity */
+  assert(!(a < b));
+  assert(!(a <= b));
+  assert((a > b));
+  assert((a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = 1.0 / -(1.0 / 0.0);	/* negative 0 */
+  assert(!(a < b));
+  assert(!(a <= b));
+  assert((a > b));
+  assert((a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = 0.0;	/* positive 0 */
+  assert(!(a < b));
+  assert(!(a <= b));
+  assert((a > b));
+  assert((a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = -42.0;
+  assert(!(a < b));
+  assert(!(a <= b));
+  assert((a > b));
+  assert((a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = 42.0;
+  assert(!(a < b));
+  assert(!(a <= b));
+  assert((a > b));
+  assert((a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  /* negative infinity */
+  a = -(1.0 / 0.0);	
+  b = (1.0 / 0.0);	/* positive infinity */
+  assert((a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert(!(a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = -(1.0 / 0.0);	/* negative infinity */
+  assert(!(a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert((a >= b));
+  assert(!(a != b));
+  assert((a == b));
+
+  b = 1.0 / -(1.0 / 0.0);	/* negative 0 */
+  assert((a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert(!(a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = 0.0;	/* positive 0 */
+  assert((a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert(!(a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = -42.0;
+  assert((a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert(!(a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+  b = 42.0;
+  assert((a < b));
+  assert((a <= b));
+  assert(!(a > b));
+  assert(!(a >= b));
+  assert((a != b));
+  assert(!(a == b));
+
+}
+