numpy
diff --git a/‎doc/release/upcoming_changes/27661.compatibility.rst
Copy file name to clipboard
+5Lines changed: 5 additions & 0 deletions b/‎doc/release/upcoming_changes/27661.compatibility.rst
Copy file name to clipboard
+5Lines changed: 5 additions & 0 deletions
diff --git a/‎numpy/lib/_function_base_impl.py
Copy file name to clipboardExpand all lines: numpy/lib/_function_base_impl.py
+1-1Lines changed: 1 addition & 1 deletion b/‎numpy/lib/_function_base_impl.py
Copy file name to clipboardExpand all lines: numpy/lib/_function_base_impl.py
+1-1Lines changed: 1 addition & 1 deletion
diff --git a/‎numpy/lib/tests/test_function_base.py
Copy file name to clipboardExpand all lines: numpy/lib/tests/test_function_base.py
+6Lines changed: 6 additions & 0 deletions b/‎numpy/lib/tests/test_function_base.py
Copy file name to clipboardExpand all lines: numpy/lib/tests/test_function_base.py
+6Lines changed: 6 additions & 0 deletions
@@ -0,0 +1,5 @@
+* `numpy.cov` now properly transposes single-row (2d array) design matrices
+  when ``rowvar=False``. Previously, single-row design matrices would
+  return a scalar in this scenario, which is not correct, so this
+  is a behavior change and an array of the appropriate shape will
+  now be returned.
@@ -2736,7 +2736,7 @@ def cov(m, y=None, rowvar=True, bias=False, ddof=None, fweights=None,
             dtype = np.result_type(m, y, np.float64)
 
     X = array(m, ndmin=2, dtype=dtype)
-    if not rowvar and X.shape[0] != 1:
+    if not rowvar and m.ndim != 1:
         X = X.T
     if X.shape[0] == 0:
         return np.array([]).reshape(0, 0)
 
@@ -2509,6 +2509,12 @@ def test_cov_dtype(self, test_type):
         res = cov(cast_x1, dtype=test_type)
         assert test_type == res.dtype
 
+    def test_gh_27658(self):
+        x = np.ones((3, 1))
+        expected = np.cov(x, ddof=0, rowvar=True)
+        actual = np.cov(x.T, ddof=0, rowvar=False)
+        assert_allclose(actual, expected, strict=True)
+
 
 class Test_I0: