Explicitly handle min/max for -0.0 and +0.0

Author: juntyr

quaddtype/numpy_quaddtype/src/ops.hpp

@@ -514,26 +514,52 @@ quad_mod(const Sleef_quad *a, const Sleef_quad *b)
 static inline Sleef_quad
 quad_minimum(const Sleef_quad *in1, const Sleef_quad *in2)
 {
-    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in1, *in1) ? *in1 : *in2)
-           : Sleef_fminq1(*in1, *in2);
+    if (Sleef_iunordq1(*in1, *in2)) {
+        return Sleef_iunordq1(*in1, *in1) ? *in1 : *in2;
+    }
+    // minimum(-0.0, +0.0) = -0.0
+    if (Sleef_icmpeqq1(*in1, QUAD_ZERO) && Sleef_icmpeqq1(*in2, QUAD_ZERO)) {
+        return Sleef_icmpleq1(Sleef_copysignq1(QUAD_ONE, *in1), Sleef_copysignq1(QUAD_ONE, *in2)) ? *in1 : *in2;
+    }
+    return Sleef_fminq1(*in1, *in2);
 }
 
 static inline Sleef_quad
 quad_maximum(const Sleef_quad *in1, const Sleef_quad *in2)
 {
-    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in1, *in1) ? *in1 : *in2)
-           : Sleef_fmaxq1(*in1, *in2);
+    if (Sleef_iunordq1(*in1, *in2)) {
+        return Sleef_iunordq1(*in1, *in1) ? *in1 : *in2;
+    }
+    // maximum(-0.0, +0.0) = +0.0
+    if (Sleef_icmpeqq1(*in1, QUAD_ZERO) && Sleef_icmpeqq1(*in2, QUAD_ZERO)) {
+        return Sleef_icmpgeq1(Sleef_copysignq1(QUAD_ONE, *in1), Sleef_copysignq1(QUAD_ONE, *in2)) ? *in1 : *in2;
+    }
+    return Sleef_fmaxq1(*in1, *in2);
 }
 
 static inline Sleef_quad
 quad_fmin(const Sleef_quad *in1, const Sleef_quad *in2)
 {
+    if (Sleef_iunordq1(*in1, *in2)) {
+        return Sleef_iunordq1(*in2, *in2) ? *in1 : *in2;
+    }
+    // fmin(-0.0, +0.0) = -0.0
+    if (Sleef_icmpeqq1(*in1, QUAD_ZERO) && Sleef_icmpeqq1(*in2, QUAD_ZERO)) {
+        return Sleef_icmpleq1(Sleef_copysignq1(QUAD_ONE, *in1), Sleef_copysignq1(QUAD_ONE, *in2)) ? *in1 : *in2;
+    }
     return Sleef_fminq1(*in1, *in2);
 }
 
 static inline Sleef_quad
 quad_fmax(const Sleef_quad *in1, const Sleef_quad *in2)
 {
+    if (Sleef_iunordq1(*in1, *in2)) {
+        return Sleef_iunordq1(*in2, *in2) ? *in1 : *in2;
+    }
+    // maximum(-0.0, +0.0) = +0.0
+    if (Sleef_icmpeqq1(*in1, QUAD_ZERO) && Sleef_icmpeqq1(*in2, QUAD_ZERO)) {
+        return Sleef_icmpgeq1(Sleef_copysignq1(QUAD_ONE, *in1), Sleef_copysignq1(QUAD_ONE, *in2)) ? *in1 : *in2;
+    }
     return Sleef_fmaxq1(*in1, *in2);
 }