Changes to logic for determining order of output array in elementwise functions, clip, matmul (#1678)

When multiple input arrays are both C- and F-contiguous and need to be cast, the output defaults to C-contiguous instead of F-contiguous

Author: ndgrigorian

dpctl/tensor/_clip.py

@@ -791,17 +791,17 @@ def clip(x, /, min=None, max=None, out=None, order="K"):
 
         if order == "K":
             if (
-                x.flags.f_contiguous
-                and a_min.flags.f_contiguous
-                and a_max.flags.f_contiguous
-            ):
-                order = "F"
-            elif (
                 x.flags.c_contiguous
                 and a_min.flags.c_contiguous
                 and a_max.flags.c_contiguous
             ):
                 order = "C"
+            elif (
+                x.flags.f_contiguous
+                and a_min.flags.f_contiguous
+                and a_max.flags.f_contiguous
+            ):
+                order = "F"
         if order == "K":
             buf1 = _empty_like_orderK(a_min, buf1_dt)
         else:

dpctl/tensor/_elementwise_common.py

@@ -858,10 +858,10 @@ def __call__(self, o1, o2, /, *, out=None, order="K"):
             return out
 
         if order == "K":
-            if src1.flags.f_contiguous and src2.flags.f_contiguous:
-                order = "F"
-            elif src1.flags.c_contiguous and src2.flags.c_contiguous:
+            if src1.flags.c_contiguous and src2.flags.c_contiguous:
                 order = "C"
+            elif src1.flags.f_contiguous and src2.flags.f_contiguous:
+                order = "F"
         if order == "K":
             buf1 = _empty_like_orderK(src1, buf1_dt)
         else:

dpctl/tensor/_linear_algebra_functions.py

@@ -942,10 +942,10 @@ def matmul(x1, x2, out=None, dtype=None, order="K"):
         return out
 
     if order == "K":
-        if x1.flags.f_contiguous and x2.flags.f_contiguous:
-            order = "F"
-        elif x1.flags.c_contiguous and x2.flags.c_contiguous:
+        if x1.flags.c_contiguous and x2.flags.c_contiguous:
             order = "C"
+        elif x1.flags.f_contiguous and x2.flags.f_contiguous:
+            order = "F"
     if order == "K":
         buf1 = _empty_like_orderK(x1, buf1_dt)
     else: