def METHOD_NAME(center):
    # do the axis sum with a dot
    # this gets called a LOT so worth optimizing
    radii_sq = np.dot((points - center) ** 2, ones)
    # residuals are difference between mean
    # use our sum mean vs .mean() as it is slightly faster
    return radii_sq - (radii_sq.sum() / len(radii_sq))