blob: 91e993d7e447d861740693a1d3732e303291879e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
/*
* Written by Michael Shalayeff. Public Domain
*/
#include <sys/types.h>
#include <machine/ieeefp.h>
#include "math.h"
float
roundf(float x)
{
u_int64_t ofpsr, fpsr;
__asm__ __volatile__("fstds %%fr0,0(%1)" : "=m" (ofpsr) : "r" (&ofpsr));
fpsr = (ofpsr & ~((u_int64_t)FP_RM << (9 + 32))) |
((u_int64_t)FP_RN << (9 + 32));
__asm__ __volatile__("fldds 0(%0), %%fr0" :: "r" (&fpsr));
__asm__ __volatile__("frnd,sgl %0,%0" : "+f" (x));
__asm__ __volatile__("fldds 0(%0), %%fr0" :: "r" (&ofpsr));
return (x);
}
|
