TEST3

pytim/interface.py

@@ -190,7 +190,7 @@ def _define_cluster_group(self):
             labels = np.array(labels)
 
             # counts is not necessarily ordered by size of cluster.
-            sorting = np.argsort(counts)[::-1]
+            sorting = np.argsort(counts,stable=True)[::-1]
             # labels for atoms in each cluster starting from the largest
             unique_labels = np.sort(np.unique(labels[labels > -1]))
             # by default, all elements of the cluster_group are in

pytim/itim.py

@@ -380,7 +380,7 @@ def _assign_layers(self):
         _y = utilities.get_y(self.cluster_group, self.normal)
         _z = utilities.get_z(self.cluster_group, self.normal)
 
-        sort = np.argsort(_z + _radius * np.sign(_z))
+        sort = np.argsort(_z + _radius * np.sign(_z),stable=True)
         # NOTE: np.argsort returns the sorted *indices*
 
         if self.multiproc and ('win' not in self.system.lower()):

pytim/observables/contactangle.py

@@ -496,7 +496,7 @@ def _fit_ellipse_fitzgibbon(x, y):
         C = np.zeros((6, 6), dtype=float)
         C[0, 2], C[2, 0], C[1, 1] = 2, 2, -1
         eigval, eigvec = np.linalg.eig(np.linalg.solve(S, C))
-        sort = np.argsort(eigval)
+        sort = np.argsort(eigval,stable=True)
         eigval, eigvec = eigval[sort], eigvec[sort]
         lam = np.nonzero(np.logical_and(
             eigval > 0, np.isfinite(eigval)))[0][-1]
@@ -515,7 +515,7 @@ def _fit_ellipse_flusser(x, y):
         C = np.array(((0, 0, 2), (0, -1, 0), (2, 0, 0)), dtype=float)
         M = np.linalg.solve(C, M)
         eigval, eigvec = np.linalg.eig(M)
-        sort = np.argsort(eigval)[::-1]
+        sort = np.argsort(eigval,stable=True)[::-1]
         eigval, eigvec = eigval[sort], eigvec[sort]
         con = 4 * eigvec[0] * eigvec[2] - eigvec[1]**2 > 0
         ak = eigvec[:, np.nonzero(con)[0]]