Try to fix fmin-fmax when isnan(in2)

juntyr · juntyr · commit 720d2ba5ce80 · 2025-07-26T14:42:38.000Z
diff --git a/quaddtype/numpy_quaddtype/src/ops.hpp b/quaddtype/numpy_quaddtype/src/ops.hpp
@@ -405,15 +405,15 @@ quad_maximum(const Sleef_quad *in1, const Sleef_quad *in2)
 static inline Sleef_quad
 quad_fmin(const Sleef_quad *in1, const Sleef_quad *in2)
 {
-    return Sleef_iunordq1(*in1, *in1) ? (
+    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)
 {
-    return Sleef_iunordq1(*in1, *in1) ? (
+    return Sleef_iunordq1(*in1, *in2) ? (
         Sleef_iunordq1(*in2, *in2) ? *in1 : *in2
     ) : Sleef_icmpgeq1(*in1, *in2) ? *in1 : *in2;
 }

Original file line number	Diff line number	Diff line change
`@@ -405,15 +405,15 @@ quad_maximum(const Sleef_quad in1, const Sleef_quad in2)`
`405`	`405`	`static inline Sleef_quad`
`406`	`406`	`quad_fmin(const Sleef_quad in1, const Sleef_quad in2)`
`407`	`407`	`{`
`408`		`- return Sleef_iunordq1(in1, in1) ? (`
	`408`	`+ return Sleef_iunordq1(in1, in2) ? (`
`409`	`409`	`Sleef_iunordq1(in2, in2) ? in1 : in2`
`410`	`410`	`) : Sleef_icmpleq1(in1, in2) ? in1 : in2;`
`411`	`411`	`}`
`412`	`412`
`413`	`413`	`static inline Sleef_quad`
`414`	`414`	`quad_fmax(const Sleef_quad in1, const Sleef_quad in2)`
`415`	`415`	`{`
`416`		`- return Sleef_iunordq1(in1, in1) ? (`
	`416`	`+ return Sleef_iunordq1(in1, in2) ? (`
`417`	`417`	`Sleef_iunordq1(in2, in2) ? in1 : in2`
`418`	`418`	`) : Sleef_icmpgeq1(in1, in2) ? in1 : in2;`
`419`	`419`	`}`