def METHOD_NAME(toas):
    temp_par = """
            F0    10 1 0.0001
            JUMP -fe L-wide 0.01 1 0.001
            JUMP -fe S-wide 0.02 1 0.001
            JUMP mjd 55000 56000 0.03 1 0.001
            JUMP freq 1400 1440 0.04 1 0.001
            JUMP tel ao 0.07 1 0.001
               """
    model = get_model(StringIO(temp_par))
    assert len(model.jumps) == 5
    assert len(model.jump_phase(toas, 0.0)) == toas.ntoas
    jump_phase0 = model.F0.quantity * (
        model.JUMP1.quantity + model.JUMP4.quantity + model.JUMP5.quantity
    )
    assert (model.jump_phase(toas, 0.0)[0] - jump_phase0) < 1e-16 * jump_phase0.unit