diff options
Diffstat (limited to 'tests32/test19.c')
| -rw-r--r-- | tests32/test19.c | 158 |
1 files changed, 108 insertions, 50 deletions
diff --git a/tests32/test19.c b/tests32/test19.c index edb1f755..5d4c412f 100644 --- a/tests32/test19.c +++ b/tests32/test19.c @@ -32,7 +32,7 @@ uint64_t _fucompp_(double a, double b) :"=a" (ret):"m"(a), "m"(b):"cc"); return ret; } -uint64_t _fistpw_(double a) +uint64_t _fistpw_(double a, uint8_t rd) { uint16_t ret; uint16_t t1, t2; @@ -40,15 +40,31 @@ uint64_t _fistpw_(double a) "fldl %1\n" "fstcw %2\n" "mov %2, %%ax\n" - "mov $0x0c, %%ah\n" + "mov %4, %%ah\n" "mov %%ax, %3\n" "fldcw %3\n" "fistp %0\n" "fldcw %2\n" - :"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc"); + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); return ret; } -uint64_t _fistpl_(double a) +uint64_t _fisttpw_(double a, uint8_t rd) +{ + uint16_t ret; + uint16_t t1, t2; + asm volatile ( + "fldl %1\n" + "fstcw %2\n" + "mov %2, %%ax\n" + "mov %4, %%ah\n" + "mov %%ax, %3\n" + "fldcw %3\n" + "fisttp %0\n" + "fldcw %2\n" + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); + return ret; +} +uint64_t _fistpl_(double a, uint8_t rd) { uint32_t ret; uint16_t t1, t2; @@ -56,15 +72,31 @@ uint64_t _fistpl_(double a) "fldl %1\n" "fstcw %2\n" "mov %2, %%ax\n" - "mov $0x0c, %%ah\n" + "mov %4, %%ah\n" "mov %%ax, %3\n" "fldcw %3\n" "fistpl %0\n" "fldcw %2\n" - :"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc"); + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); + return ret; +} +uint64_t _fisttpl_(double a, uint8_t rd) +{ + uint32_t ret; + uint16_t t1, t2; + asm volatile ( + "fldl %1\n" + "fstcw %2\n" + "mov %2, %%ax\n" + "mov %4, %%ah\n" + "mov %%ax, %3\n" + "fldcw %3\n" + "fisttpl %0\n" + "fldcw %2\n" + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); return ret; } -uint64_t _fistpq_(double a) +uint64_t _fistpq_(double a, uint8_t rd) { uint64_t ret; uint16_t t1, t2; @@ -72,15 +104,31 @@ uint64_t _fistpq_(double a) "fldl %1\n" "fstcw %2\n" "mov %2, %%ax\n" - "mov $0x0c, %%ah\n" + "mov %4, %%ah\n" "mov %%ax, %3\n" "fldcw %3\n" "fistpq %0\n" "fldcw %2\n" - :"=m" (ret):"m"(a), "m"(t1), "m"(t2):"cc"); + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); + return ret; +} +uint64_t _fisttpq_(double a, uint8_t rd) +{ + uint64_t ret; + uint16_t t1, t2; + asm volatile ( + "fldl %1\n" + "fstcw %2\n" + "mov %2, %%ax\n" + "mov %4, %%ah\n" + "mov %%ax, %3\n" + "fldcw %3\n" + "fisttpq %0\n" + "fldcw %2\n" + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); return ret; } -uint64_t _frndint_(double a) +uint64_t _frndint_(double a, uint8_t rd) { uint64_t ret; uint16_t t1, t2; @@ -88,13 +136,13 @@ uint64_t _frndint_(double a) "fldl %1\n" "fstcw %2\n" "mov %2, %%ax\n" - "mov $0x0c, %%ah\n" + "mov %4, %%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"); + :"=m" (ret):"m"(a), "m"(t1), "m"(t2), "m"(rd):"cc"); return ret; } #endif @@ -105,44 +153,54 @@ int main(int argc, const char** argv) 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)); + uint64_t tests[] = { + 0x0, + 0x4082c00000000000LL, + 0x8000000000000000LL, + 0x40dfffc000000000LL, + 0x40e0002000000000LL, // 0x7fff and 0x8001 as double + 0xc0e0002000000000LL, + 0xc0dfffc000000000LL, // -0x8001 and -0x7fff as double + 0x3ff0000000000000LL, + 0x3fe89d9000000000LL, + 0x7ff0000000000000LL, + 0xfff0000000000000LL, + 0x3ff0002ca0000000LL, + 0xaeff000025000000LL, + 0x3ff0000050000000LL, + 0xc082c00000000000LL, + 0x7ff8000000000001LL, +}; +int n = sizeof(tests)/sizeof(tests[0]); +for(int i=0; i<n; ++i) + for(int j=0; j<n; ++j) { + *(uint64_t*)&a = tests[i]; + *(uint64_t*)&b = tests[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); + } +for(int i=0; i<n; ++i) { + *(uint64_t*)&a = tests[i]; + for(int rd=0; rd<3; ++rd) { + printf("Rounding %d\n", rd); + printf(" FRNDINT 0x%llx => 0x%llx\n", *(uint64_t*)&a, _frndint_(a, rd<<2)); + printf(" FISTP 0x%llx => word: %x\n", *(uint64_t*)&a, _fistpw_(a, rd<<2)); + printf(" FISTP 0x%llx => long: %x\n", *(uint64_t*)&a, _fistpl_(a, rd<<2)); + printf(" FISTP 0x%llx => quad: %llx\n", *(uint64_t*)&a, _fistpq_(a, rd<<2)); + printf(" FISTTP 0x%llx => word: %x\n", *(uint64_t*)&a, _fisttpw_(a, rd<<2)); + printf(" FISTTP 0x%llx => long: %x\n", *(uint64_t*)&a, _fisttpl_(a, rd<<2)); + printf(" FISTTP 0x%llx => quad: %llx\n", *(uint64_t*)&a, _fisttpq_(a, rd<<2)); + } } printf("\nDone\n"); } |