def METHOD_NAME():
    up = np.array([0, 0, 1])
    query_point = [FITTER.x0 + (0.2 * R_SPHERE), FITTER.y0, FITTER.z0]
    dist = FITTER.distance_to_shell_along_vector_from_point(up, query_point, np.arange(1, 50))
    ground_truth = np.sqrt((R_SPHERE ** 2) - ((0.2 * R_SPHERE + FITTER.x0 - FITTER.x0) ** 2))
    np.testing.assert_almost_equal(ground_truth, dist, decimal=0)
    print(ground_truth - dist)
    query_point = [FITTER.x0 + (0.5 * R_SPHERE), FITTER.y0, FITTER.z0]
    dist = FITTER.distance_to_shell_along_vector_from_point(up, query_point, np.arange(1, 50))
    ground_truth = np.sqrt((R_SPHERE ** 2) - ((0.5 * R_SPHERE + FITTER.x0 - FITTER.x0) ** 2))
    np.testing.assert_almost_equal(ground_truth, dist, decimal=0)
    query_point = [FITTER.x0 + (0.75 * R_SPHERE), FITTER.y0, FITTER.z0]
    dist = FITTER.distance_to_shell_along_vector_from_point(up, query_point, np.arange(1, 50))
    ground_truth = np.sqrt((R_SPHERE ** 2) - ((0.75 * R_SPHERE + FITTER.x0 - FITTER.x0) ** 2))
    np.testing.assert_almost_equal(ground_truth, dist, decimal=0)
    query_point = [FITTER.x0 + (0.9 * R_SPHERE), FITTER.y0, FITTER.z0]
    dist = FITTER.distance_to_shell_along_vector_from_point(up, query_point, np.arange(1, 50))
    ground_truth = np.sqrt((R_SPHERE ** 2) - ((0.9 * R_SPHERE + FITTER.x0 - FITTER.x0) ** 2))
    np.testing.assert_almost_equal(ground_truth, dist, decimal=0)