Try different min/max implementation

Author: juntyr

quaddtype/numpy_quaddtype/src/ops.hpp

@@ -514,53 +514,33 @@ quad_mod(const Sleef_quad *a, const Sleef_quad *b)
 static inline Sleef_quad
 quad_minimum(const Sleef_quad *in1, const Sleef_quad *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);
+    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in1, *in1) ? *in1 : *in2)
+        : Sleef_icmpltq1(*in1, *in2) ? *in1
+                                    : *in2;
 }
 
 static inline Sleef_quad
 quad_maximum(const Sleef_quad *in1, const Sleef_quad *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);
+    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in1, *in1) ? *in1 : *in2)
+        : Sleef_icmpgtq1(*in1, *in2) ? *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);
+    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in2, *in2) ? *in1 : *in2)
+           : Sleef_icmpleq1(*in1, *in2) ? *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);
+    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in2, *in2) ? *in1 : *in2)
+        : Sleef_icmpgeq1(*in1, *in2) ? *in1
+                                    : *in2;
 }
 
 static inline Sleef_quad