def METHOD_NAME(centers_sphere):
    """Make sure all spheres are entirely contained within the domain."""
    d = np.linalg.norm(centers_sphere - [1, 2, 3], axis=1)
    r_max = max(d) + _RADIUS
    assert r_max < 1 or r_max == pytest.approx(1)