Try to fix min/max functions

Author: juntyr

quaddtype/numpy_quaddtype/src/ops.hpp

@@ -515,32 +515,26 @@ 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_icmpleq1(*in1, *in2) ? *in1
-                                        : *in2;
+           : 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_icmpgeq1(*in1, *in2) ? *in1
-                                        : *in2;
+           : Sleef_fmaxq1(*in1, *in2);
 }
 
 static inline Sleef_quad
 quad_fmin(const Sleef_quad *in1, const Sleef_quad *in2)
 {
-    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in2, *in2) ? *in1 : *in2)
-           : Sleef_icmpleq1(*in1, *in2) ? *in1
-                                        : *in2;
+    return Sleef_fminq1(*in1, *in2);
 }
 
 static inline Sleef_quad
 quad_fmax(const Sleef_quad *in1, const Sleef_quad *in2)
 {
-    return Sleef_iunordq1(*in1, *in2)   ? (Sleef_iunordq1(*in2, *in2) ? *in1 : *in2)
-           : Sleef_icmpgeq1(*in1, *in2) ? *in1
-                                        : *in2;
+    return Sleef_fmaxq1(*in1, *in2);
 }
 
 static inline Sleef_quad

quaddtype/tests/test_quaddtype.py

@@ -525,7 +525,7 @@ def test_hyperbolic_functions(op, val):
 
     # For finite non-zero results
     # Use relative tolerance for exponential functions due to their rapid growth
-    rtol = 1e-14 if abs(float_result) < 1e100 else 1e-10
+    rtol = 1e-13 if abs(float_result) < 1e100 else 1e-10
     np.testing.assert_allclose(float(quad_result), float_result, rtol=rtol, atol=1e-15,
                                err_msg=f"Value mismatch for {op}({val})")