def METHOD_NAME(
    init_distance_threshold, threshold_ratio, scales_to_test
):
    coords, reference, evaluation, _ = get_dummy_gamma_set()
    all_scales = np.concatenate([[1], scales_to_test])
    distance_thresholds_to_test = init_distance_threshold * all_scales
    dose_thresholds_to_test = distance_thresholds_to_test * threshold_ratio
    gamma_results = []
    for dose, distance in zip(dose_thresholds_to_test, distance_thresholds_to_test):
        gamma_results.append(
            pymedphys.gamma(
                coords,
                reference,
                coords,
                evaluation,
                dose,
                distance,
                lower_percent_dose_cutoff=0,
            )
        )
    for i, scale in enumerate(scales_to_test):
        print(scale)
        abs_diff = np.abs(gamma_results[i + 1] - gamma_results[0] / scale)
        print(np.max(abs_diff))
        assert np.all(abs_diff <= 0.1)