qemu-e2k/tests/tcg/i386/test-i386-sse-exceptions.c

814 lines
28 KiB
C
Raw Normal View History

target/i386: fix IEEE SSE floating-point exception raising The SSE instruction implementations all fail to raise the expected IEEE floating-point exceptions because they do nothing to convert the exception state from the softfloat machinery into the exception flags in MXCSR. Fix this by adding such conversions. Unlike for x87, emulated SSE floating-point operations might be optimized using hardware floating point on the host, and so a different approach is taken that is compatible with such optimizations. The required invariant is that all exceptions set in env->sse_status (other than "denormal operand", for which the SSE semantics are different from those in the softfloat code) are ones that are set in the MXCSR; the emulated MXCSR is updated lazily when code reads MXCSR, while when code sets MXCSR, the exceptions in env->sse_status are set accordingly. A few instructions do not raise all the exceptions that would be raised by the softfloat code, and those instructions are made to save and restore the softfloat exception state accordingly. Nothing is done about "denormal operand"; setting that (only for the case when input denormals are *not* flushed to zero, the opposite of the logic in the softfloat code for such an exception) will require custom code for relevant instructions, or else architecture-specific conditionals in the softfloat code for when to set such an exception together with custom code for various SSE conversion and rounding instructions that do not set that exception. Nothing is done about trapping exceptions (for which there is minimal and largely broken support in QEMU's emulation in the x87 case and no support at all in the SSE case). Signed-off-by: Joseph Myers <joseph@codesourcery.com> Message-Id: <alpine.DEB.2.21.2006252358000.3832@digraph.polyomino.org.uk> Signed-off-by: Paolo Bonzini <pbonzini@redhat.com>
2020-06-26 01:58:31 +02:00
/* Test SSE exceptions. */
#include <float.h>
#include <stdint.h>
#include <stdio.h>
volatile float f_res;
volatile double d_res;
volatile float f_snan = __builtin_nansf("");
volatile float f_half = 0.5f;
volatile float f_third = 1.0f / 3.0f;
volatile float f_nan = __builtin_nanl("");
volatile float f_inf = __builtin_inff();
volatile float f_ninf = -__builtin_inff();
volatile float f_one = 1.0f;
volatile float f_two = 2.0f;
volatile float f_zero = 0.0f;
volatile float f_nzero = -0.0f;
volatile float f_min = FLT_MIN;
volatile float f_true_min = 0x1p-149f;
volatile float f_max = FLT_MAX;
volatile float f_nmax = -FLT_MAX;
volatile double d_snan = __builtin_nans("");
volatile double d_half = 0.5;
volatile double d_third = 1.0 / 3.0;
volatile double d_nan = __builtin_nan("");
volatile double d_inf = __builtin_inf();
volatile double d_ninf = -__builtin_inf();
volatile double d_one = 1.0;
volatile double d_two = 2.0;
volatile double d_zero = 0.0;
volatile double d_nzero = -0.0;
volatile double d_min = DBL_MIN;
volatile double d_true_min = 0x1p-1074;
volatile double d_max = DBL_MAX;
volatile double d_nmax = -DBL_MAX;
volatile int32_t i32_max = INT32_MAX;
#define IE (1 << 0)
#define ZE (1 << 2)
#define OE (1 << 3)
#define UE (1 << 4)
#define PE (1 << 5)
#define EXC (IE | ZE | OE | UE | PE)
uint32_t mxcsr_default = 0x1f80;
uint32_t mxcsr_ftz = 0x9f80;
int main(void)
{
uint32_t mxcsr;
int32_t i32_res;
int ret = 0;
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = f_snan;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: widen float snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = d_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: narrow float underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = d_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: narrow float overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: narrow float inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = d_snan;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: narrow float snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundss $4, %0, %0" : "=x" (f_res) : "0" (f_min));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: roundss min\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundss $12, %0, %0" : "=x" (f_res) : "0" (f_min));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: roundss no-inexact min\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundss $4, %0, %0" : "=x" (f_res) : "0" (f_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: roundss snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundss $12, %0, %0" : "=x" (f_res) : "0" (f_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: roundss no-inexact snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundsd $4, %0, %0" : "=x" (d_res) : "0" (d_min));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: roundsd min\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundsd $12, %0, %0" : "=x" (d_res) : "0" (d_min));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: roundsd no-inexact min\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundsd $4, %0, %0" : "=x" (d_res) : "0" (d_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: roundsd snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("roundsd $12, %0, %0" : "=x" (d_res) : "0" (d_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: roundsd no-inexact snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("comiss %1, %0" : : "x" (f_nan), "x" (f_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: comiss nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("ucomiss %1, %0" : : "x" (f_nan), "x" (f_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: ucomiss nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("ucomiss %1, %0" : : "x" (f_snan), "x" (f_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: ucomiss snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("comisd %1, %0" : : "x" (d_nan), "x" (d_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: comisd nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("ucomisd %1, %0" : : "x" (d_nan), "x" (d_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: ucomisd nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("ucomisd %1, %0" : : "x" (d_snan), "x" (d_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: ucomisd snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max + f_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: float add overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max + f_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: float add inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_inf + f_ninf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float add inf -inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_snan + f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float add snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
f_res = f_true_min + f_true_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float add FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max + d_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: double add overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max + d_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: double add inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_inf + d_ninf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double add inf -inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_snan + d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double add snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
d_res = d_true_min + d_true_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double add FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max - f_nmax;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: float sub overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max - f_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: float sub inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_inf - f_inf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float sub inf inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_snan - f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float sub snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
f_res = f_min - f_true_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float sub FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max - d_nmax;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: double sub overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max - d_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: double sub inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_inf - d_inf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double sub inf inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_snan - d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double sub snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
d_res = d_min - d_true_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double sub FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max * f_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: float mul overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_third * f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: float mul inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_min * f_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float mul underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_inf * f_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float mul inf 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_snan * f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float mul snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
f_res = f_min * f_half;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float mul FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max * d_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: double mul overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_third * d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: double mul inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_min * d_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double mul underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_inf * d_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double mul inf 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_snan * d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double mul snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
d_res = d_min * d_half;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double mul FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_max / f_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: float div overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_one / f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: float div inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_min / f_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float div underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_one / f_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != ZE) {
printf("FAIL: float div 1 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_inf / f_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: float div inf 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_nan / f_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: float div nan 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_zero / f_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float div 0 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_inf / f_inf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float div inf inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
f_res = f_snan / f_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: float div snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
f_res = f_min / f_two;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: float div FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_max / d_min;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (OE | PE)) {
printf("FAIL: double div overflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_one / d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: double div inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_min / d_max;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double div underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_one / d_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != ZE) {
printf("FAIL: double div 1 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_inf / d_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: double div inf 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_nan / d_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: double div nan 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_zero / d_zero;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double div 0 0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_inf / d_inf;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double div inf inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
d_res = d_snan / d_third;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: double div snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_ftz));
d_res = d_min / d_two;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != (UE | PE)) {
printf("FAIL: double div FTZ underflow\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) : "0" (f_max));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: sqrtss inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) : "0" (f_nmax));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtss -max\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) : "0" (f_ninf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtss -inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) : "0" (f_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtss snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) : "0" (f_nzero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: sqrtss -0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtss %0, %0" : "=x" (f_res) :
"0" (-__builtin_nanf("")));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: sqrtss -nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) : "0" (d_max));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: sqrtsd inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) : "0" (d_nmax));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtsd -max\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) : "0" (d_ninf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtsd -inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) : "0" (d_snan));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: sqrtsd snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) : "0" (d_nzero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: sqrtsd -0\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("sqrtsd %0, %0" : "=x" (d_res) :
"0" (-__builtin_nan("")));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: sqrtsd -nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("maxss %1, %0" : : "x" (f_nan), "x" (f_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: maxss nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("minss %1, %0" : : "x" (f_nan), "x" (f_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: minss nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("maxsd %1, %0" : : "x" (d_nan), "x" (d_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: maxsd nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("minsd %1, %0" : : "x" (d_nan), "x" (d_zero));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: minsd nan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtsi2ss %1, %0" : "=x" (f_res) : "m" (i32_max));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: cvtsi2ss inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtsi2sd %1, %0" : "=x" (d_res) : "m" (i32_max));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: cvtsi2sd exact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtss2si %1, %0" : "=r" (i32_res) : "x" (1.5f));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: cvtss2si inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtss2si %1, %0" : "=r" (i32_res) : "x" (0x1p31f));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvtss2si 0x1p31\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtss2si %1, %0" : "=r" (i32_res) : "x" (f_inf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvtss2si inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtsd2si %1, %0" : "=r" (i32_res) : "x" (1.5));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: cvtsd2si inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtsd2si %1, %0" : "=r" (i32_res) : "x" (0x1p31));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvtsd2si 0x1p31\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvtsd2si %1, %0" : "=r" (i32_res) : "x" (d_inf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvtsd2si inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttss2si %1, %0" : "=r" (i32_res) : "x" (1.5f));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: cvttss2si inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttss2si %1, %0" : "=r" (i32_res) : "x" (0x1p31f));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvttss2si 0x1p31\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttss2si %1, %0" : "=r" (i32_res) : "x" (f_inf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvttss2si inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttsd2si %1, %0" : "=r" (i32_res) : "x" (1.5));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != PE) {
printf("FAIL: cvttsd2si inexact\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttsd2si %1, %0" : "=r" (i32_res) : "x" (0x1p31));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvttsd2si 0x1p31\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("cvttsd2si %1, %0" : "=r" (i32_res) : "x" (d_inf));
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != IE) {
printf("FAIL: cvttsd2si inf\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("rcpss %0, %0" : "=x" (f_res) : "0" (f_snan));
f_res += f_one;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: rcpss snan\n");
ret = 1;
}
__asm__ volatile ("ldmxcsr %0" : : "m" (mxcsr_default));
__asm__ volatile ("rsqrtss %0, %0" : "=x" (f_res) : "0" (f_snan));
f_res += f_one;
__asm__ volatile ("stmxcsr %0" : "=m" (mxcsr));
if ((mxcsr & EXC) != 0) {
printf("FAIL: rsqrtss snan\n");
ret = 1;
}
return ret;
}