def METHOD_NAME():
    # Dummy points in 3d
    xyzs = np.array([[0., 0., 0.],
                     [1., 0., 0.],
                     [0., 1., 0.],
                     [0., 1., 1.],
                     [0., 0., 0.]])
    dist_mat = np.sqrt(np.sum(
        (xyzs[np.newaxis, :, :] - xyzs[:, np.newaxis, :])**2, axis=-1))
    # Calculate geometric center
    center = np.mean(xyzs, axis=0)
    # Distances to the center:
    dists_ref = np.sqrt(np.sum((xyzs - center)**2, axis=-1))
    dists, status = GetCMDists(dist_mat, len(xyzs))
    assert status
    np.testing.assert_allclose(dists, dists_ref)