StructuralIdentifiability.jl
StructuralIdentifiability.jl copied to clipboard
ModelingToolkit extension completely freezing Julia
using ModelingToolkit
using StructuralIdentifiability
const SI = StructuralIdentifiability
function testSI()
t = ModelingToolkit.t_nounits
D = ModelingToolkit.D_nounits
@variables X(t)
@parameters p d
eqs = [
D(X) ~ p - d*X
]
@mtkbuild osys = ODESystem(eqs, t)
measured_quantities = [X]
funcs_to_check = [osys.p]
SI.assess_identifiability(osys; measured_quantities, funcs_to_check)
end
causes Julia to completely hang (as reported on Slack by @TorkelE).
Note that changing to funcs_to_check = [p]
fixes the problem.
Environment:
(StructuralIdentifiability) pkg> st
Project StructuralIdentifiability v0.5.9
Status `~/.julia/dev/StructuralIdentifiability/Project.toml`
⌅ [c3fe647b] AbstractAlgebra v0.41.11
[861a8166] Combinatorics v1.0.2
[864edb3b] DataStructures v0.18.20
⌅ [0b43b601] Groebner v0.7.5
[c8e1da08] IterTools v1.10.0
[1914dd2f] MacroTools v0.5.13
⌅ [2edaba10] Nemo v0.45.7
⌅ [3e851597] ParamPunPam v0.4.1
[aea7be01] PrecompileTools v1.2.1
[27ebfcd6] Primes v0.5.6
[a759f4b9] TimerOutputs v0.5.24
[ade2ca70] Dates
[37e2e46d] LinearAlgebra
[56ddb016] Logging
[9a3f8284] Random
Info Packages marked with ⌅ have new versions available but compatibility constraints restrict them from upgrading. To see why use `status --outdated`
on Julia 1.10.5.