summaryrefslogtreecommitdiff
path: root/src/rtapi/rtapi_math_i386.h
blob: 74b5825a909740980aa1881aae02c184d9c8570f (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
39
40
41
42
43
44
45
46
47
/* GPL, etc -- based loosely on mathinline.h */
#ifndef RTAPI_MATH_I386_H
#define RTAPI_MATH_I386_H

#define MATHOP(f,i) extern __inline double f(double __x) { \
    register double __result; \
    __asm __volatile__ (i : "=t" (__result) : "0" (__x)); \
    return __result; \
}

MATHOP(sin, "fsin")
MATHOP(cos, "fcos")
MATHOP(fabs, "fabs")
MATHOP(sqrt, "fsqrt")

extern __inline double tan (double __x) {
    register long double __value;
    register long double __value2 __attribute__ ((__unused__));
    __asm __volatile__ ("fptan" : "=t" (__value2), "=u" (__value) : "0" (__x));
    return __value;
}

extern __inline double atan2 (double __y, double __x) {
    register long double __value;
    __asm __volatile__ ("fpatan" : "=t" (__value) : "0" (__x), "u" (__y) : "st(1)");
    return __value;
}

extern __inline double atan (double __y) {
    return atan2(__y, 1.);
}

extern __inline double asin (double __x) {
    return atan2(__x, sqrt (1.0 - __x * __x));
}

extern __inline double acos (double __x) {
    return atan2(sqrt(1.0 - __x * __x), __x);
}

extern __inline double fmax(double __y, double __x) {
    return __y > __x ? __y : __x;
}
extern __inline double fmin(double __y, double __x) {
    return __y < __x ? __y : __x;
}
#endif /* RTAPI_MATH_I386_H */