blob: 6b2dfb39be457f3109ce3480ef17d681367e476a (
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
31
32
33
34
35
36
37
38
|
/*
* ====================================================
* x87 FP implementation contributed to Newlib by
* Dave Korn, November 2007. This file is placed in the
* public domain. Permission to use, copy, modify, and
* distribute this software is freely granted.
* ====================================================
*/
#ifdef __GNUC__
#if !defined(_SOFT_FLOAT)
#include <math.h>
/*
* Fast math version of llrintf(x)
* Return x rounded to integral value according to the prevailing
* rounding mode.
* Method:
* Using inline x87 asms.
* Exception:
* Governed by x87 FPCR.
*/
long long int _f_llrintf (float x)
{
long long int _result;
asm ("fistpll %0" : "=m" (_result) : "t" (x) : "st");
return _result;
}
/* For now, we only have the fast math version. */
long long int llrintf (float x) {
return _f_llrintf(x);
}
#endif /* !_SOFT_FLOAT */
#endif /* __GNUC__ */
|