1 files changed, 84 insertions, 0 deletions
diff --git a/test/monniaux/math/exceptions.c b/test/monniaux/math/exceptions.c
new file mode 100644
index 00000000..84ed54db
--- /dev/null
+++ b/test/monniaux/math/exceptions.c
@@ -0,0 +1,84 @@
+#include <stdio.h>
+#include <fenv.h>
+#include <float.h>
+
+#pragma STDC FENV_ACCESS ON
+
+#if defined(__KVX__) && !defined(__COMPCERT__)
+int fetestexcept(int excepts) {
+  int mask = (K1_SFR_CS_IO_MASK | K1_SFR_CS_DZ_MASK | K1_SFR_CS_OV_MASK | K1_SFR_CS_UN_MASK | K1_SFR_CS_IN_MASK) & excepts;
+  unsigned long long cs = __builtin_kvx_get(K1_SFR_CS);
+  return cs & mask;
+}
+
+int feclearexcept(int excepts) {
+  int mask = (K1_SFR_CS_IO_MASK | K1_SFR_CS_DZ_MASK | K1_SFR_CS_OV_MASK | K1_SFR_CS_UN_MASK | K1_SFR_CS_IN_MASK) & excepts;
+  __builtin_kvx_wfxl(K1_SFR_CS, mask);
+  return 0;
+}
+#endif
+
+double add(double x, double y) {
+  return x+y;
+}
+
+double mul(double x, double y) {
+  return x*y;
+}
+
+float double2float(double x) {
+  return x;
+}
+
+float uint2float(unsigned x) {
+  return x;
+}
+
+double ulong2double(unsigned long x) {
+  return x;
+}
+
+unsigned double2uint(double x) {
+  return x;
+}
+
+int main() {
+  printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
+
+  double v1 = add(3.0, 0.1);
+  printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
+  feclearexcept(FE_INEXACT);
+
+  printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
+  double v2 = mul(DBL_MAX, DBL_MAX);
+  printf("%g %x\n", v2, fetestexcept(FE_ALL_EXCEPT));
+  feclearexcept(FE_ALL_EXCEPT);
+
+  double v3 = mul(DBL_MIN, DBL_MIN);
+  printf("%g %x\n", v3, fetestexcept(FE_ALL_EXCEPT));
+  feclearexcept(FE_ALL_EXCEPT);
+
+  double v4 = double2float(DBL_MAX);
+  printf("%g %x\n", v4, fetestexcept(FE_ALL_EXCEPT));
+  feclearexcept(FE_ALL_EXCEPT);
+
+  float v5 = uint2float(0xC07FDFFFU);
+  printf("%g %x\n", v5, fetestexcept(FE_ALL_EXCEPT)); // BUG 0 should have INEXACT
+  feclearexcept(FE_ALL_EXCEPT);
+
+  double v6 = ulong2double(0x11217763AFF77C7CUL);
+  printf("%g %x\n", v6, fetestexcept(FE_ALL_EXCEPT)); // BUG 0 should have INEXACT
+  feclearexcept(FE_ALL_EXCEPT);
+
+  unsigned v7 = double2uint(-0.25); // softfloat says "0 and inexact" but here we have "0 and overflow" (due to negative input for unsigned?)
+  printf("%u %x\n", v7, fetestexcept(FE_ALL_EXCEPT));
+  feclearexcept(FE_ALL_EXCEPT);
+
+  // +41F.307672C5496EF
+  double d8 = 0x1.307672C5496EFp32;
+  unsigned v8 = double2uint(d8);
+  printf("%g %x %x\n", d8, v8, fetestexcept(FE_ALL_EXCEPT));
+  // BUG reports 307672C5 and inexact, but should report overflow
+  // bug comes from x86 https://gcc.gnu.org/bugzilla/show_bug.cgi?id=89175
+  feclearexcept(FE_ALL_EXCEPT);
+}