def dist_diff(self, geom): diff = self.coords[:, None, :]-geom.coords[None, :, :] dist = np.sqrt(np.sum(diff**2, 2)) dist[np.diag_indices(len(self))] = np.inf return dist, diff