1#define _GNU_SOURCE
2#include <fenv.h>
3#include <features.h>
4
5static inline double get_fpscr_f(void)
6{
7 double d;
8 __asm__ __volatile__("mffs %0" : "=d"(d));
9 return d;
10}
11
12static inline long get_fpscr(void)
13{
14 return (union {double f; long i;}) {get_fpscr_f()}.i;
15}
16
17static inline void set_fpscr_f(double fpscr)
18{
19 __asm__ __volatile__("mtfsf 255, %0" : : "d"(fpscr));
20}
21
22static void set_fpscr(long fpscr)
23{
24 set_fpscr_f((union {long i; double f;}) {fpscr}.f);
25}
26
27int feclearexcept(int mask)
28{
29 mask &= FE_ALL_EXCEPT;
30 if (mask & FE_INVALID) mask |= FE_ALL_INVALID;
31 set_fpscr(get_fpscr() & ~mask);
32 return 0;
33}
34
35int feraiseexcept(int mask)
36{
37 mask &= FE_ALL_EXCEPT;
38 if (mask & FE_INVALID) mask |= FE_INVALID_SOFTWARE;
39 set_fpscr(get_fpscr() | mask);
40 return 0;
41}
42
43int fetestexcept(int mask)
44{
45 return get_fpscr() & mask & FE_ALL_EXCEPT;
46}
47
48int fegetround(void)
49{
50 return get_fpscr() & 3;
51}
52
53hidden int __fesetround(int r)
54{
55 set_fpscr(get_fpscr() & ~3L | r);
56 return 0;
57}
58
59int fegetenv(fenv_t *envp)
60{
61 *envp = get_fpscr_f();
62 return 0;
63}
64
65int fesetenv(const fenv_t *envp)
66{
67 set_fpscr_f(envp != FE_DFL_ENV ? *envp : 0);
68 return 0;
69}