You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
first_in =findfirst((X->begin#= In[14]:11 =#set(X) ⊆ P
end), vsol.Xk) =888783all((X->begin#= In[14]:14 =#set(X) ⊆ P
end), vsol.Xk[first_in:end]) =truetspan(vsol[first_in]) = [0.0887895, 0.0887897]
max values of x
ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.04850572743425879ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:end])) =0.0492794261847055ρ([1.0, 0.0], Flowpipe(vsol.Xk[end:end])) =0.0492794261847055
min values of x
-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.048000002027475806-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:end]))) =0.048000002027475806-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[end:end]))) =0.04854423385309506
max values of v
ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.0007994567734929507ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:end])) =0.0007994567734929507ρ([0.0, 1.0], Flowpipe(vsol.Xk[end:end])) =0.0005874379912341613
min values of v
-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.0006327526764030105-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:end]))) =0.0003450783028747313-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[end:end]))) =0.00034507830287473130.0003450783028747313
no pv, no jitter, delta = 1e-8, GLGM06
prob =embrake_no_pv(ζ=0.0, Tsample=1e-4);
# GLGM06, order 1, static#@time sol = solve(prob, alg=GLGM06(δ=1e-8, max_order=1, static=true, dim=4, ngens=4), max_jumps=2);@time sol =solve(prob, alg=GLGM06(δ=1e-8, max_order=1, static=true, dim=4, ngens=4), max_jumps=1000);
final_diams(sol)
0.414480 seconds (271.44 k allocations:1.691 GiB, 1.54% gc time)
final diameter of current I =17.746604024816392
final diameter of position x =0.0009518262957221718safety_property(sol)
first_in =findfirst((X->begin#= In[56]:11 =#set(X) ⊆ P
end), vsol.Xk) =9015845all((X->begin#= In[56]:14 =#set(X) ⊆ P
end), vsol.Xk[first_in:end]) =truetspan(vsol[first_in]) = [0.0900594, 0.0900595]
max values of x
ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.04868287562433759ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:end])) =0.04938757641401575ρ([1.0, 0.0], Flowpipe(vsol.Xk[end:end])) =0.04938757641401575
min values of x
-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.04800000018332063-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:end]))) =0.04800000018332063-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[end:end]))) =0.048435750118293576
max values of v
ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.0007950689857855407ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:end])) =0.0007950689857855407ρ([0.0, 1.0], Flowpipe(vsol.Xk[end:end])) =0.0006232115341202552
min values of v
-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.0005699644062819035-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:end]))) =0.00030943639703057016-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[end:end]))) =0.000309436397030570160.00030943639703057016
no pv, no jitter, delta = 1e-9, GLGM06
no pv, with jitter, delta = 1e-8, ASB07
prob =embrake_no_pv(ζ=[-1e-8, 1e-7], Tsample=1e-4);
# GLGM06, order 1, static@time sol =solve(prob, alg=GLGM06(δ=1e-8, max_order=1, static=true, dim=4, ngens=4), max_jumps=1000);
final_diams(sol)
2.109413 seconds (795.14 k allocations:1.716 GiB, 0.67% gc time)
final diameter of current I =17.746604024816392
final diameter of position x =0.0009518262957221718safety_property(sol, ngens=4)
first_in =findfirst((X->begin#= In[56]:11 =#set(X) ⊆ P
end), vsol.Xk) =9015845all((X->begin#= In[56]:14 =#set(X) ⊆ P
end), vsol.Xk[first_in:end]) =truetspan(vsol[first_in]) = [0.0900594, 0.0900595]
max values of x
ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.04868287562433759ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:end])) =0.04938757641401575ρ([1.0, 0.0], Flowpipe(vsol.Xk[end:end])) =0.04938757641401575
min values of x
-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.04800000018332063-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:end]))) =0.04800000018332063-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[end:end]))) =0.048435750118293576
max values of v
ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.0007950689857855407ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:end])) =0.0007950689857855407ρ([0.0, 1.0], Flowpipe(vsol.Xk[end:end])) =0.0006232115341202552
min values of v
-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.0005699644062819035-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:end]))) =0.00030943639703057016-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[end:end]))) =0.000309436397030570160.00030943639703057016
pv1, no jitter, delta = 1e-8, ASB07
prob =embrake_pv_1(ζ=0.0, Tsample=1e-4);
# order 3@time sol =solve(prob, alg=ASB07(δ=1e-8, max_order=3, dim=4, ngens=12, static=true), max_jumps=1_000);
final_diams(sol)
final diameter of current I =2.9414439401694423
final diameter of position x =0.0001229578152434907safety_property(sol, ngens=12)
first_in =findfirst((X->begin#= In[56]:11 =#set(X) ⊆ P
end), vsol.Xk) =8654084all((X->begin#= In[56]:14 =#set(X) ⊆ P
end), vsol.Xk[first_in:end]) =truetspan(vsol[first_in]) = [0.0865408, 0.0865409]
max values of x
ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.048161937383798405ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:end])) =0.04897071950798819ρ([1.0, 0.0], Flowpipe(vsol.Xk[end:end])) =0.04897071950798819
min values of x
-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.048000000058297514-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:end]))) =0.048000000058297514-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[end:end]))) =0.0488477616927447
max values of v
ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.0008165062857196601ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:end])) =0.0008165062857196601ρ([0.0, 1.0], Flowpipe(vsol.Xk[end:end])) =0.0004932708432558988
min values of v
-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.0007460799328546838-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:end]))) =0.0004412635986992693-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[end:end]))) =0.00044126359869926930.0004412635986992693
pv1, with jitter, delta = 1e-8, ASB07
prob =embrake_pv_1(ζ=[-1e-8, 1e-7], Tsample=1e-4);
@time sol =solve(prob, alg=ASB07(δ=1e-8, max_order=1, dim=4, ngens=4, static=true), max_jumps=1);
@time sol =solve(prob, alg=ASB07(δ=1e-8, max_order=1, dim=4, ngens=4, static=true), max_jumps=1000);
final_diams(sol)
0.722996 seconds (1.23 M allocations:71.672 MiB)
16.972784 seconds (110.52 M allocations:10.243 GiB, 43.87% gc time)
final diameter of current I =154.207482098115
final diameter of position x =0.008210036200221565safety_property(sol, ngens=4, ε=0.005)
first_in =findfirst((X->begin#= In[56]:11 =#set(X) ⊆ P
end), vsol.Xk) =7251941all((X->begin#= In[56]:14 =#set(X) ⊆ P
end), vsol.Xk[first_in:end]) =falsetspan(vsol[first_in]) = [0.0724397, 0.0724398]
max values of x
ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.04825960654834148ρ([1.0, 0.0], Flowpipe(vsol.Xk[first_in:end])) =0.05301736140115831ρ([1.0, 0.0], Flowpipe(vsol.Xk[end:end])) =0.05301736140115831
min values of x
-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.04500000058208545-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[first_in:end]))) =0.04480732520093674-(ρ([-1.0, 0.0], Flowpipe(vsol.Xk[end:end]))) =0.04480732520093674
max values of v
ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:first_in])) =0.0018794473599369806ρ([0.0, 1.0], Flowpipe(vsol.Xk[first_in:end])) =0.0018794473599369806ρ([0.0, 1.0], Flowpipe(vsol.Xk[end:end])) =0.0018293238731492435
min values of v
-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:first_in]))) =0.0007834008334339722-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[first_in:end]))) =-0.0008971962975791285-(ρ([0.0, -1.0], Flowpipe(vsol.Xk[end:end]))) =-0.0008971962975791285-0.0008971962975791285