def _sphere_constraint(rd, r0, R_adj):
 return (R_adj - np.sqrt(np.sum(((rd - r0) ** 2))))

def _sphere_constraint ( rd , r0 , R_adj ) : return ( R_adj - np . sqrt ( np . sum ( ( ( rd - r0 ) ** 2 ) ) ) )
sphere fitting constraint .