diff options
Diffstat (limited to 'test/monniaux/math')
-rw-r--r-- | test/monniaux/math/exceptions.c | 41 | ||||
-rw-r--r-- | test/monniaux/math/rounding.c | 36 |
2 files changed, 77 insertions, 0 deletions
diff --git a/test/monniaux/math/exceptions.c b/test/monniaux/math/exceptions.c new file mode 100644 index 00000000..d4a6a26b --- /dev/null +++ b/test/monniaux/math/exceptions.c @@ -0,0 +1,41 @@ +#include <stdio.h> +#include <fenv.h> +#include <float.h> + +#ifdef __K1C__ +#include <mppa_bare_runtime/k1c/registers.h> +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); + 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_k1_wfxl(K1_SFR_CS, mask); + return 0; +} +#endif + +#pragma STDC FENV_ACCESS ON + +double add(double x, double y) { + return x+y; +} + +double mul(double x, double y) { + return x*y; +} + +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)); +} diff --git a/test/monniaux/math/rounding.c b/test/monniaux/math/rounding.c new file mode 100644 index 00000000..c2ce85e3 --- /dev/null +++ b/test/monniaux/math/rounding.c @@ -0,0 +1,36 @@ +#include <stdio.h> +#include <fenv.h> + +#ifdef __K1C__ +#include <mppa_bare_runtime/k1c/registers.h> +int fesetround(int rounding_mode) { + if (rounding_mode < 0 || rounding_mode > 3) return 1; + unsigned long long cs = __builtin_k1_get(K1_SFR_CS); + cs = (cs & ~(3 << 16)) | (rounding_mode << 16); + __builtin_k1_set(K1_SFR_CS, cs); + return 0; +} + +int fegetround(void) { + unsigned long long cs = __builtin_k1_get(K1_SFR_CS); + return (cs >> 16) & 3; +} +#endif + +#pragma STDC FENV_ACCESS ON + +double add(double x) { + return x+0.1; +} + +int main() { + if (fesetround(FE_DOWNWARD)) { + printf("fesetround(FE_DOWNWARD) unsuccessful\n"); + } + double low = add(3.0); + if (fesetround(FE_UPWARD)) { + printf("fesetround(FE_UPWARD) unsuccessful\n"); + } + double high = add(3.0); + printf("%d %a %d\n", low==high, high-low, fegetround()); +} |