def METHOD_NAME(model):
    t = make_fake_toas_uniform(56000, 59000, 10, model, obs="@")
    t.table["error"] = 1 * u.us
    delta_f = (1 / (t.last_MJD - t.first_MJD)).to(u.Hz)
    m_2 = deepcopy(model)
    m_2.F0.quantity += 2 * delta_f
    m_2.TRACK.value = "-2"
    r = Residuals(t, m_2)
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) > 1
    r = Residuals(t, m_2, track_mode="nearest")
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) <= 1
    r = Residuals(t, m_2, track_mode="use_pulse_numbers")
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) > 1
    m_2.TRACK.value = "0"
    r = Residuals(t, m_2)
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) <= 1
    r = Residuals(t, m_2, track_mode="nearest")
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) <= 1
    r = Residuals(t, m_2, track_mode="use_pulse_numbers")
    assert np.amax(r.phase_resids) - np.amin(r.phase_resids) > 1