StructuralIdentifiability.jl icon indicating copy to clipboard operation
StructuralIdentifiability.jl copied to clipboard

ModelingToolkit extension completely freezing Julia

Open isaacsas opened this issue 4 months ago • 8 comments

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.

isaacsas avatar Oct 03 '24 21:10 isaacsas