blob: de3334a626d881f0b46e10071d7b036cf3faa248 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
#include <fenv.h>
long double
roundl (long double x) {
long double retval;
unsigned short saved_cw, _cw;
__asm__ (
"fnstcw %0;"
: "=m" (saved_cw)
); /* save control word */
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */
__asm__ (
"fldcw %0;"
:
: "m" (_cw)
); /* load the rounding control */
__asm__ (
"frndint;"
: "=t" (retval)
: "0" (x)
); /* do the rounding */
__asm__ (
"fldcw %0;"
:
: "m" (saved_cw)
); /* restore control word */
return retval;
}
|