At least on my system the following test code crashes with a Julia bug (Unreachable reached, signal (4): Illegal instruction).
|
A_itv = IntervalMatrix([-1.1..0.9 -4.1.. -3.9; 3.9..4.1 -1.1..0.9]) |
|
X0 = BallInf(zeros(2), 0.1) |
|
sys_homog = IVP(CLCS(A_itv, Universe(2)), X0) |
|
B = IntervalMatrix([1.0..1.0 0.0..0.0; 0.0..0.0 1.0..1.0]) |
|
U = ConstantInput(BallInf(ones(2), 0.5)) |
|
sys_heterog = IVP(CLCCS(A_itv, B, Universe(2), U), X0) |
|
algorithm = "interval_matrix" |
|
|
|
for cont_sys in [sys_homog, sys_heterog] |
|
# default order |
|
discr_sys = discretize(cont_sys, δ; algorithm=algorithm) |
I tried to debug it and apparently there is a problem when we create the result system:
|
return IVP(CLCDS(ϕ, B, stateset(𝑆.s), Ud), Ω0) |
The system itself can be created without problems if you outsorce this line to a new helper function. But for some reason the compiler seems surprised that it can return that result and then crashes.
At least on my system the following test code crashes with a Julia bug (Unreachable reached, signal (4): Illegal instruction).
Reachability.jl/test/ReachSets/unit_discretization.jl
Lines 151 to 161 in 6fb396d
I tried to debug it and apparently there is a problem when we create the result system:
Reachability.jl/src/ReachSets/discretize.jl
Line 832 in 6fb396d
The system itself can be created without problems if you outsorce this line to a new helper function. But for some reason the compiler seems surprised that it can return that result and then crashes.