More cases to test17, and added -NAN generation to divpd ([DYNAREC] too)

2 files changed, 28 insertions, 4 deletions

diff --git a/src/dynarec/arm64/dynarec_arm64_660f.c b/src/dynarec/arm64/dynarec_arm64_660f.c
index abd12c04..b90c27cb 100755
--- a/src/dynarec/arm64/dynarec_arm64_660f.c
+++ b/src/dynarec/arm64/dynarec_arm64_660f.c
@@ -769,8 +769,21 @@ uintptr_t dynarec64_660F(dynarec_arm_t* dyn, uintptr_t addr, uintptr_t ip, int n
             INST_NAME("DIVPD Gx, Ex");

             nextop = F8;

             GETEX(q0, 0);

-            GETGX(v0);

-            VFDIVQD(v0, v0, q0);

+            GETGX(q1);

+            if(!box64_dynarec_fastnan) {

+                v0 = fpu_get_scratch(dyn);

+                v1 = fpu_get_scratch(dyn);

+                // check if any input value was NAN

+                VFMAXQD(v0, q0, q1);    // propagate NAN

+                FCMEQQD(v0, v0, v0);    // 0 if NAN, 1 if not NAN

+            }

+            VFDIVQD(q1, q1, q0);

+            if(!box64_dynarec_fastnan) {

+                FCMEQQD(v1, q1, q1);    // 0 => out is NAN

+                VBICQ(v1, v0, v1);      // forget it in any input was a NAN already

+                VSHLQ_64(v1, v1, 63);   // only keep the sign bit

+                VORRQ(q1, q1, v1);      // NAN -> -NAN

+            }

             break;

         case 0x5F:

             INST_NAME("MAXPD Gx, Ex");

diff --git a/src/emu/x64run660f.c b/src/emu/x64run660f.c
index 5db9994a..ada264d0 100644
--- a/src/emu/x64run660f.c
+++ b/src/emu/x64run660f.c
@@ -55,6 +55,9 @@ int Run660F(x64emu_t *emu, rex_t rex)
     uint32_t tmp32u;

     uint64_t tmp64u;

     float tmpf;

+    #ifndef NOALIGN

+    int is_nan;

+    #endif

     reg64_t *oped, *opgd;

     sse_regs_t *opex, *opgx, eax1, *opex2;

     mmx87_regs_t *opem, *opgm;

@@ -912,8 +915,16 @@ int Run660F(x64emu_t *emu, rex_t rex)
         nextop = F8;

         GETEX(0);

         GETGX;

-        GX->d[0] /= EX->d[0];

-        GX->d[1] /= EX->d[1];

+        for (int i=0; i<2; ++i) {

+            #ifndef NOALIGN

+            is_nan = isnan(GX->d[i]) || isnan(EX->d[i]);

+            #endif

+            GX->d[i] /= EX->d[i];

+            #ifndef NOALIGN

+            if(!is_nan && isnan(GX->d[i]))

+                GX->d[i] = -NAN;

+            #endif

+        }

         break;

     case 0x5F:                      /* MAXPD Gx, Ex */

         nextop = F8;