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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
#ifndef FP_EMULATION_H_INCLUDED
#define FP_EMULATION_H_INCLUDED
#include <math.h>
int is_snan(float x)
{
union32f_ud fb;
fb.f = x;
return __builtin_isnan(x) && (fb.u & 0x00400000) == 0;
}
int is_qnan(float x)
{
return __builtin_isnan(x) && !is_snan(x);
}
#define INTEL_SSE_MATH_OP(op, a, b) \
({ \
union32f_ud tmp; \
float result; \
if (is_qnan(a)) \
result = a; \
else if (is_snan(a)) \
{ \
tmp.f = a; tmp.u |= 0x400000; result = tmp.f; \
} else if (is_snan(b)) \
{ \
tmp.f = b; tmp.u |= 0x400000; result = tmp.f; \
} else if (is_qnan(b)) \
result = b; \
else \
result = a op b; \
result; \
})
float
intel_sse_math_add(float a, float b)
{
return INTEL_SSE_MATH_OP(+, a, b);
}
float
intel_sse_math_sub(float a, float b)
{
return INTEL_SSE_MATH_OP(-, a, b);
}
float
intel_sse_math_mul(float a, float b)
{
return INTEL_SSE_MATH_OP(*, a, b);
}
float
intel_sse_math_div(float a, float b)
{
return INTEL_SSE_MATH_OP(/, a, b);
}
/* +-3 == +-inf,
+-2 == +-Nan, use -2 only,
+-1 == +- normal number,
0 == undefined/not intialized */
int
state_handler(int src0, int src1, __int128 a, __int128 b, char op)
{
/* Nan */
if (src0 == -2 || src1 == -2)
return -2;
if (abs(src0) == 3 || abs(src1) == 3)
{
/* +INF + +INF = +INF, -INF + -INF = -INF */
if (src0 == src1)
return src0;
/* Positive result */
if ((op == '*' && ((src0 < 0 && b < 0) || (src1 < 0 && a < 0))) ||
(op == '*' && ((src0 > 0 && b > 0) || (src1 > 0 && a > 0))))
return 3;
/* -INF * (positive normal) = -INF */
if (op == '*' && ((src0 < 0 || src1 < 0) && (a >= 0 && b >= 0)))
return -3;
/* INF * 0 = NaN */
if (a == 0 || b == 0)
return -2;
/* -INF + +INF = NaN, else +INF or -INF */
return (op == '+' && src0 != src1) ? -2 : 3 * ((src0 / 3) + (src1 / 3));
}
/* Normal number case */
__int128 result = op == '+' ? a + b : a * b;
return result < 0 ? -1 : 1;
}
int
state_handler_float (float src)
{
if (isnan (src))
return -2;
else if (isinf (src))
return src > 0 ? 3 : -3;
else
return src >= 0 ? 1 : -1;
}
#endif
|
