From dc4d89e2b91012334855bbbf184df32aa5d6d982 Mon Sep 17 00:00:00 2001 From: David Monniaux Date: Sat, 13 Apr 2019 22:11:06 +0200 Subject: some more examples --- test/monniaux/math/exceptions.c | 43 +++++++++++++++++++++++++++++++++++++---- 1 file changed, 39 insertions(+), 4 deletions(-) (limited to 'test') diff --git a/test/monniaux/math/exceptions.c b/test/monniaux/math/exceptions.c index d4a6a26b..05423e62 100644 --- a/test/monniaux/math/exceptions.c +++ b/test/monniaux/math/exceptions.c @@ -2,8 +2,9 @@ #include #include -#ifdef __K1C__ -#include +#pragma STDC FENV_ACCESS ON + +#ifdef __GNUC__ 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_k1_get(K1_SFR_CS); @@ -17,8 +18,6 @@ int feclearexcept(int excepts) { } #endif -#pragma STDC FENV_ACCESS ON - double add(double x, double y) { return x+y; } @@ -27,15 +26,51 @@ 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); } -- cgit