def METHOD_NAME(model, toas):
    toas.compute_pulse_numbers(model)
    toas["delta_pulse_number"][:10] = 1
    toas["delta_pulse_number"][10:20] = -1
    outtim = StringIO()
    toas.write_TOA_file(outtim)
    outtim.seek(0)
    newtoas = get_TOAs(outtim)
    assert (newtoas[:10]["delta_pulse_number"] == 1).all()
    assert (newtoas[10:20]["delta_pulse_number"] == -1).all()
    assert (newtoas[20:]["delta_pulse_number"] == 0).all()
    assert (toas["delta_pulse_number"] == newtoas["delta_pulse_number"]).all()