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
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
#include <string.h>
#include <stdio.h>
#include <stddef.h>
#include <stdlib.h>
#include <stdint.h>
#include <math.h>
#if defined(__x86_64__)
#error Nope!
#else
uint64_t _fucomip_(double a, double b)
{
uint32_t ret;
asm volatile (
"fldl %2\n"
"fldl %1\n"
"fucomip %%st(1)\n"
"fstp %%st(0)\n"
"pushfl\n"
"popl %%eax\n"
:"=a" (ret):"m"(a), "m"(b):"cc");
return ret;
}
uint64_t _fucompp_(double a, double b)
{
uint32_t ret;
asm volatile (
"fldl %2\n"
"fldl %1\n"
"fucompp\n"
"fstsw %%ax\n"
:"=a" (ret):"m"(a), "m"(b):"cc");
return ret;
}
uint64_t _fistpw_(double a)
{
uint16_t ret;
uint16_t t1, t2;
asm volatile (
"fldl %1\n"
"fstcw %2\n"
"mov %2, %%ax\n"
"mov $0x0c, %%ah\n"
"mov %%ax, %3\n"
"fldcw %3\n"
"fistp %0\n"
"fldcw %2\n"
:"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc");
return ret;
}
uint64_t _fistpl_(double a)
{
uint32_t ret;
uint16_t t1, t2;
asm volatile (
"fldl %1\n"
"fstcw %2\n"
"mov %2, %%ax\n"
"mov $0x0c, %%ah\n"
"mov %%ax, %3\n"
"fldcw %3\n"
"fistpl %0\n"
"fldcw %2\n"
:"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc");
return ret;
}
uint64_t _fistpq_(double a)
{
uint64_t ret;
uint16_t t1, t2;
asm volatile (
"fldl %1\n"
"fstcw %2\n"
"mov %2, %%ax\n"
"mov $0x0c, %%ah\n"
"mov %%ax, %3\n"
"fldcw %3\n"
"fistpq %0\n"
"fldcw %2\n"
:"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc");
return ret;
}
uint64_t _frndint_(double a)
{
uint64_t ret;
uint16_t t1, t2;
asm volatile (
"fldl %1\n"
"fstcw %2\n"
"mov %2, %%ax\n"
"mov $0x0c, %%ah\n"
"mov %%ax, %3\n"
"fldcw %3\n"
"frndint\n"
"fstpl %0\n"
"fldcw %2\n"
:"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc");
return ret;
}
#endif
int main(int argc, const char** argv)
{
double a, b;
char pf, cf, zf;
char pa, ca, za;
uint64_t flags;
uint64_t tests[][2] = {
{0x0, 0x4082c00000000000LL},
{0x8000000000000000LL, 0x4082c00000000000LL},
{0x8000000000000000LL, 0x0},
{0x8000000000000000LL, 0x3ff0000000000000LL},
{0x3ff0000000000000LL, 0x3fe89d9000000000LL},
{0x3ff0000000000000LL, 0x7ff0000000000000LL},
{0xfff0000000000000LL, 0x7ff0000000000000LL},
{0x3ff0002ca0000000LL, 0xaeff000025000000LL},
{0x3ff0000050000000LL, 0xc082c00000000000LL},
{0x0, 0x7ff8000000000001LL},
{0x7ff8000000000001LL, 0x7ff8000000000001LL},
};
int n = sizeof(tests)/sizeof(tests[0]);
for(int i=0; i<n; ++i)
for(int j=0; j<2; ++j) {
*(uint64_t*)&a = tests[i][0+j];
*(uint64_t*)&b = tests[i][1-j];
printf("FUCOMI 0x%llx 0x%llx ", *(uint64_t*)&a, *(uint64_t*)&b);
flags = _fucomip_(a, b);
ca = (flags>>0)&1?'C':'-';
za = (flags>>(0+6))&1?'Z':'-';
pa = (flags>>(0+2))&1?'P':'-';
printf("%c%c%c\n", za, pa, ca);
printf("FUCOMPP 0x%llx 0x%llx ", *(uint64_t*)&a, *(uint64_t*)&b);
flags = _fucompp_(a, b);
ca = (flags>>8)&1?'C':'-';
za = (flags>>(8+6))&1?'Z':'-';
pa = (flags>>(8+2))&1?'P':'-';
printf("%c%c%c\n", za, pa, ca);
printf("FRNDINT 0x%llx => 0x%llx\n", *(uint64_t*)&a, _frndint_(a));
printf("FRNDINT 0x%llx => 0x%llx\n", *(uint64_t*)&b, _frndint_(b));
printf("FISTP 0x%llx => word: %x\n", *(uint64_t*)&a, _fistpw_(a));
printf("FISTP 0x%llx => word: %x\n", *(uint64_t*)&b, _fistpw_(b));
printf("FISTP 0x%llx => long: %x\n", *(uint64_t*)&a, _fistpl_(a));
printf("FISTP 0x%llx => long: %x\n", *(uint64_t*)&b, _fistpl_(b));
printf("FISTP 0x%llx => quad: %llx\n", *(uint64_t*)&a, _fistpq_(a));
printf("FISTP 0x%llx => quad: %llx\n", *(uint64_t*)&b, _fistpq_(b));
}
printf("\nDone\n");
}
|