return np.sqrt(np.sum((self.coord[:, np.newaxis] - self.coord) ** 2, axis=2))