Fix outlier handling in qap when len(seeds)==n

graspologic/match/qap.py

@@ -400,6 +400,8 @@ def _quadratic_assignment_faq(
 
     # check outlier cases
     if n == 0 or partial_match.shape[0] == n:
+        # Cannot assume partial_match is sorted.
+        partial_match = np.row_stack(sorted(partial_match, key=lambda x: x[0]))
         score = _calc_score(A, B, S, partial_match[:, 1])
         res = {"col_ind": partial_match[:, 1], "fun": score, "nit": 0}
         return OptimizeResult(res)

tests/test_match.py

@@ -141,6 +141,24 @@ def test_barycenter_SGM(self):
         score = chr12c.score_
         assert 11156 == score
 
+        W1 = np.array(range(n))
+        W2 = [pi[z] for z in W1]
+        chr12c = self.barycenter.fit(A, B, W1, W2)
+        score = chr12c.score_
+        assert 11156 == score
+
+        W1 = np.array(range(n))
+        W2 = [pi[z] for z in W1]
+        # Shuffle seed pairs.
+        pairs = [i for i in zip(W1, W2)]
+        random.shuffle(pairs)
+        W1 = list(list(zip(*pairs))[0])
+        W2 = list(list(zip(*pairs))[1])
+        chr12c = self.barycenter.fit(A, B, W1, W2)
+        score = chr12c.score_
+        assert 11156 == score
+
+
     def test_rand_SGM(self):
         A, B = self._get_AB()
         chr12c = self.rand.fit(A, B)