numpy
diff --git a/‎numpy/core/src/umath/matmul.c.src
Lines changed: 7 additions & 6 deletions b/‎numpy/core/src/umath/matmul.c.src
Lines changed: 7 additions & 6 deletions
diff --git a/‎numpy/core/tests/test_multiarray.py
Lines changed: 12 additions & 2 deletions b/‎numpy/core/tests/test_multiarray.py
Lines changed: 12 additions & 2 deletions
@@ -282,18 +282,19 @@ BOOL_matmul_inner_noblas(void *_ip1, npy_intp is1_m, npy_intp is1_n,
 
     for (m = 0; m < dm; m++) {
         for (p = 0; p < dp; p++) {
+            npy_bool *ip1tmp = ip1;
+            npy_bool *ip2tmp = ip2;
             *(npy_bool *)op = NPY_FALSE;
             for (n = 0; n < dn; n++) {
-                npy_bool val1 = (*(npy_bool *)ip1);
-                npy_bool val2 = (*(npy_bool *)ip2);
+                npy_bool val1 = (*(npy_bool *)ip1tmp);
+                npy_bool val2 = (*(npy_bool *)ip2tmp);
                 if (val1 != 0 && val2 != 0) {
                     *(npy_bool *)op = NPY_TRUE;
+                    break;
                 }
-                ip2 += is2_n;
-                ip1 += is1_n;
+                ip2tmp += is2_n;
+                ip1tmp += is1_n;
             }
-            ip1 -= ib1_n;
-            ip2 -= ib2_n;
             op  +=  os_p;
             ip2 += is2_p;
         }
 
@@ -6254,10 +6254,20 @@ def __mul__(self, other):
     def test_matmul_bool(self):
         # gh-14439
         a = np.array([[1, 0],[1, 1]], dtype=bool)
-        assert np.max(a.view(np.int8)) == 1
+        assert np.max(a.view(np.uint8)) == 1
         b = np.matmul(a, a)
         # matmul with boolean output should always be 0, 1
-        assert np.max(b.view(np.int8)) == 1
+        assert np.max(b.view(np.uint8)) == 1
+
+        rg = np.random.default_rng(np.random.PCG64(43))
+        d = rg.integers(2, size=4*5, dtype=np.int8)
+        d = d.reshape(4, 5) > 0
+        out1 = np.matmul(d, d.reshape(5, 4))
+        out2 = np.dot(d, d.reshape(5, 4))
+        assert_equal(out1, out2)
+
+        c = np.matmul(np.zeros((2, 0), dtype=bool), np.zeros(0, dtype=bool))
+        assert not np.any(c)
 
 
 if sys.version_info[:2] >= (3, 5):