Error when simulating ODESystem with difference equations #1304
Description
Trying out the status of #894
Here's a reproducing example generated with ODEProblemExpr
using ModelingToolkit, OrdinaryDiffEq
prob = begin
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:624 =#
f = begin
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:417 =#
var"##f#569" = (ModelingToolkit.ODEFunctionClosure)(function (ˍ₋arg1, ˍ₋arg2, t)
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:284 =#
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:285 =#
let var"model3₊mass1₊vel(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[1]), var"model3₊mass2₊vel(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[2]), var"model3₊mass1₊pos(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[3]), var"model3₊mass2₊pos(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[4]), var"discrete_pid₊x2(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[5]), var"discrete_pid₊x1(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[6]), model3₊mass1₊m = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[1]), model3₊sd₊damper₊c = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[2]), model3₊sd₊spring₊k = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[3]), model3₊mass2₊m = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[4])
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:377 =#
(SymbolicUtils.Code.create_array)(typeof(ˍ₋arg1), nothing, Val{1}(), Val{(6,)}(), (+)((+)((*)(-3.3333333333333335e-5, var"model3₊mass2₊pos(t)"), (*)(3.3333333333333335e-5, (>=)(t, 1)), (*)(0.006666666666666666, var"discrete_pid₊x2(t)"), (*)(-1, (Difference(t; dt=0.01, update=false))(var"discrete_pid₊x1(t)"))), (+)(var"discrete_pid₊x1(t)")), (+)((*)(-0.006666666666666666, var"model3₊mass2₊pos(t)"), (*)(0.006666666666666666, (>=)(t, 1)), (*)(0.3333333333333333, var"discrete_pid₊x2(t)"), (*)(-1, (Difference(t; dt=0.01, update=false))(var"discrete_pid₊x2(t)"))), (/)((+)((+)((*)(-1, (+)(var"model3₊mass1₊pos(t)", (*)(-1//1, var"model3₊mass2₊pos(t)")), model3₊sd₊spring₊k), (*)(-133266.33333333334, var"discrete_pid₊x2(t)"), (*)(100.0, var"discrete_pid₊x1(t)"), (*)(-1333.6683333333333, var"model3₊mass2₊pos(t)")), (+)((*)(1333.6683333333333, (>=)(t, 1)), (*)(-1, model3₊sd₊damper₊c, (+)((*)(-1//1, var"model3₊mass2₊vel(t)"), var"model3₊mass1₊vel(t)")))), model3₊mass1₊m), (/)((+)((*)((+)(var"model3₊mass1₊pos(t)", (*)(-1//1, var"model3₊mass2₊pos(t)")), model3₊sd₊spring₊k), (*)(2, (>=)(t, 10)), (*)(model3₊sd₊damper₊c, (+)((*)(-1//1, var"model3₊mass2₊vel(t)"), var"model3₊mass1₊vel(t)"))), model3₊mass2₊m), var"model3₊mass1₊vel(t)", var"model3₊mass2₊vel(t)")
end
end, function (ˍ₋out, ˍ₋arg1, ˍ₋arg2, t)
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:284 =#
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:285 =#
let var"model3₊mass1₊vel(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[1]), var"model3₊mass2₊vel(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[2]), var"model3₊mass1₊pos(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[3]), var"model3₊mass2₊pos(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[4]), var"discrete_pid₊x2(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[5]), var"discrete_pid₊x1(t)" = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg1[6]), model3₊mass1₊m = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[1]), model3₊sd₊damper₊c = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[2]), model3₊sd₊spring₊k = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[3]), model3₊mass2₊m = #= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:171 =# @inbounds(ˍ₋arg2[4])
#= /home/fredrikb/.julia/packages/Symbolics/wJd48/src/build_function.jl:373 =#
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:331 =# @inbounds begin
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:327 =#
ˍ₋out[1] = (+)((+)((*)(-3.3333333333333335e-5, var"model3₊mass2₊pos(t)"), (*)(3.3333333333333335e-5, (>=)(t, 1)), (*)(0.006666666666666666, var"discrete_pid₊x2(t)"), (*)(-1, (Difference(t; dt=0.01, update=false))(var"discrete_pid₊x1(t)"))), (+)(var"discrete_pid₊x1(t)"))
ˍ₋out[2] = (+)((*)(-0.006666666666666666, var"model3₊mass2₊pos(t)"), (*)(0.006666666666666666, (>=)(t, 1)), (*)(0.3333333333333333, var"discrete_pid₊x2(t)"), (*)(-1, (Difference(t; dt=0.01, update=false))(var"discrete_pid₊x2(t)")))
ˍ₋out[3] = (/)((+)((+)((*)(-1, (+)(var"model3₊mass1₊pos(t)", (*)(-1//1, var"model3₊mass2₊pos(t)")), model3₊sd₊spring₊k), (*)(-133266.33333333334, var"discrete_pid₊x2(t)"), (*)(100.0, var"discrete_pid₊x1(t)"), (*)(-1333.6683333333333, var"model3₊mass2₊pos(t)")), (+)((*)(1333.6683333333333, (>=)(t, 1)), (*)(-1, model3₊sd₊damper₊c, (+)((*)(-1//1, var"model3₊mass2₊vel(t)"), var"model3₊mass1₊vel(t)")))), model3₊mass1₊m)
ˍ₋out[4] = (/)((+)((*)((+)(var"model3₊mass1₊pos(t)", (*)(-1//1, var"model3₊mass2₊pos(t)")), model3₊sd₊spring₊k), (*)(2, (>=)(t, 10)), (*)(model3₊sd₊damper₊c, (+)((*)(-1//1, var"model3₊mass2₊vel(t)"), var"model3₊mass1₊vel(t)"))), model3₊mass2₊m)
ˍ₋out[5] = var"model3₊mass1₊vel(t)"
ˍ₋out[6] = var"model3₊mass2₊vel(t)"
#= /home/fredrikb/.julia/packages/SymbolicUtils/vKBkE/src/code.jl:329 =#
nothing
end
end
end)
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:418 =#
var"##tgrad#570" = nothing
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:419 =#
var"##jac#571" = nothing
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:420 =#
M = [0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0; 1.0 0.0 0.0 0.0 0.0 0.0; 0.0 1.0 0.0 0.0 0.0 0.0; 0.0 0.0 1.0 0.0 0.0 0.0; 0.0 0.0 0.0 1.0 0.0 0.0]
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:421 =#
ODEFunction{true}(var"##f#569", jac = var"##jac#571", tgrad = var"##tgrad#570", mass_matrix = M, jac_prototype = nothing, syms = [Symbol("model3₊mass1₊vel(t)"), Symbol("model3₊mass2₊vel(t)"), Symbol("model3₊mass1₊pos(t)"), Symbol("model3₊mass2₊pos(t)"), Symbol("discrete_pid₊x2(t)"), Symbol("discrete_pid₊x1(t)")], indepsym = :t)
end
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:625 =#
u0 = [0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:626 =#
tspan = (0.0, 20.0)
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:627 =#
p = [1.0, 10.0, 1000.0, 1.0]
#= /home/fredrikb/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:628 =#
ODEProblem(f, u0, tspan, p; )
end
solve(prob, Rosenbrock23())julia> solve(prob, Rosenbrock23())
ERROR: MethodError: Cannot `convert` an object of type SymbolicUtils.Mul{Float64, Int64, Dict{Any, Number}, Nothing} to an object of type Float64
Closest candidates are:
convert(::Type{T}, ::VectorizationBase.AbstractSIMD) where T<:Union{Bool, Float16, Float32, Float64, Int16, Int32, Int64, Int8, UInt16, UInt32, UInt64, UInt8, SIMDTypes.Bit} at /home/fredrikb/.julia/packages/VectorizationBase/bOEgn/src/base_defs.jl:154
convert(::Type{T}, ::Intervals.AnchoredInterval{P, T, L, R} where {L<:Intervals.Bounded, R<:Intervals.Bounded}) where {P, T} at /home/fredrikb/.julia/packages/Intervals/ua9cq/src/anchoredinterval.jl:181
convert(::Type{T}, ::Intervals.Interval{T, L, R} where {L<:Intervals.Bound, R<:Intervals.Bound}) where T at /home/fredrikb/.julia/packages/Intervals/ua9cq/src/interval.jl:253
...
Stacktrace:
[1] setindex!(A::Vector{Float64}, x::SymbolicUtils.Mul{Float64, Int64, Dict{Any, Number}, Nothing}, i1::Int64)
@ Base ./array.jl:839
[2] (::var"#8#10")(ˍ₋out::Vector{Float64}, ˍ₋arg1::Vector{Float64}, ˍ₋arg2::Vector{Float64}, t::Float64)
@ Main ./REPL[2]:19
[3] ODEFunctionClosure
@ ~/.julia/dev/ModelingToolkit/src/systems/diffeqs/abstractodesystem.jl:374 [inlined]
[4] ODEFunction
@ ~/.julia/packages/SciMLBase/h4Gxc/src/scimlfunctions.jl:334 [inlined]
[5] _initialize_dae!(integrator::OrdinaryDiffEq.ODEIntegrator{Rosenbrock23{6, true, DefaultLinSolve, Val{:forward}}, true, Vector{Float64}, Nothing, Float64, Vector{Float64}, Float64, Float64, Float64, Float64, Vector{Vector{Float64}}, ODESolution{Float64, 2, Vector{Vector{Float64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{Float64}}}, ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, Rosenbrock23{6, true, DefaultLinSolve, Val{:forward}}, OrdinaryDiffEq.InterpolationData{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.Rosenbrock23Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, OrdinaryDiffEq.Rosenbrock23Tableau{Float64}, SciMLBase.TimeGradientWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Float64}, Vector{Float64}}, SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, DefaultLinSolve, SparseDiffTools.ForwardColorJacCache{Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, Float64}, Float64, 6}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, Float64}, Float64, 6}}, Vector{Float64}, Vector{Vector{NTuple{6, Float64}}}, UnitRange{Int64}, Nothing}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.TimeGradientWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Float64}, Vector{Float64}}, Float64}, Float64, 1}}}}, DiffEqBase.DEStats}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.Rosenbrock23Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, OrdinaryDiffEq.Rosenbrock23Tableau{Float64}, SciMLBase.TimeGradientWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Float64}, Vector{Float64}}, SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, DefaultLinSolve, SparseDiffTools.ForwardColorJacCache{Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, Float64}, Float64, 6}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.UJacobianWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Vector{Float64}}, Float64}, Float64, 6}}, Vector{Float64}, Vector{Vector{NTuple{6, Float64}}}, UnitRange{Int64}, Nothing}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{SciMLBase.TimeGradientWrapper{ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Float64}, Vector{Float64}}, Float64}, Float64, 1}}}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(opnorm), Nothing, CallbackSet{Tuple{}, Tuple{}}, typeof(DiffEqBase.ODE_DEFAULT_ISOUTOFDOMAIN), typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE), typeof(DiffEqBase.ODE_DEFAULT_UNSTABLE_CHECK), DataStructures.BinaryHeap{Float64, DataStructures.FasterForward}, DataStructures.BinaryHeap{Float64, DataStructures.FasterForward}, Nothing, Nothing, Int64, Tuple{}, Tuple{}, Tuple{}}, Vector{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::BrownFullBasicInit{Float64}, #unused#::Val{true})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/initialize_dae.jl:285
[6] _initialize_dae!
@ ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/initialize_dae.jl:40 [inlined]
[7] initialize_dae! (repeats 2 times)
@ ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/initialize_dae.jl:31 [inlined]
[8] __init(prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::Rosenbrock23{6, true, DefaultLinSolve, Val{:forward}}, timeseries_init::Tuple{}, ts_init::Tuple{}, ks_init::Tuple{}, recompile::Type{Val{true}}; saveat::Tuple{}, tstops::Tuple{}, d_discontinuities::Tuple{}, save_idxs::Nothing, save_everystep::Bool, save_on::Bool, save_start::Bool, save_end::Nothing, callback::Nothing, dense::Bool, calck::Bool, dt::Float64, dtmin::Nothing, dtmax::Float64, force_dtmin::Bool, adaptive::Bool, gamma::Rational{Int64}, abstol::Nothing, reltol::Nothing, qmin::Rational{Int64}, qmax::Int64, qsteady_min::Int64, qsteady_max::Rational{Int64}, beta1::Nothing, beta2::Nothing, qoldinit::Rational{Int64}, controller::Nothing, fullnormalize::Bool, failfactor::Int64, maxiters::Int64, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), internalopnorm::typeof(opnorm), isoutofdomain::typeof(DiffEqBase.ODE_DEFAULT_ISOUTOFDOMAIN), unstable_check::typeof(DiffEqBase.ODE_DEFAULT_UNSTABLE_CHECK), verbose::Bool, timeseries_errors::Bool, dense_errors::Bool, advance_to_tstop::Bool, stop_at_next_tstop::Bool, initialize_save::Bool, progress::Bool, progress_steps::Int64, progress_name::String, progress_message::typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE), userdata::Nothing, allow_extrapolation::Bool, initialize_integrator::Bool, alias_u0::Bool, alias_du0::Bool, initializealg::OrdinaryDiffEq.DefaultInit, kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/solve.jl:436
[9] __init(prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::Rosenbrock23{6, true, DefaultLinSolve, Val{:forward}}, timeseries_init::Tuple{}, ts_init::Tuple{}, ks_init::Tuple{}, recompile::Type{Val{true}}) (repeats 5 times)
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/solve.jl:67
[10] #__solve#471
@ ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/solve.jl:4 [inlined]
[11] __solve
@ ~/.julia/packages/OrdinaryDiffEq/Zi9Zh/src/solve.jl:4 [inlined]
[12] #solve_call#42
@ ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:61 [inlined]
[13] solve_call(_prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, args::Rosenbrock23{6, true, DefaultLinSolve, Val{:forward}})
@ DiffEqBase ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:48
[14] solve_up(prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, sensealg::Nothing, u0::Vector{Float64}, p::Vector{Float64}, args::Rosenbrock23{0, true, DefaultLinSolve, Val{:forward}}; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ DiffEqBase ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:87
[15] solve_up
@ ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:78 [inlined]
[16] #solve#43
@ ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:73 [inlined]
[17] solve(prob::ODEProblem{Vector{Float64}, Tuple{Float64, Float64}, true, Vector{Float64}, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#7#9", var"#8#10"}, Matrix{Float64}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, args::Rosenbrock23{0, true, DefaultLinSolve, Val{:forward}})
@ DiffEqBase ~/.julia/packages/DiffEqBase/OPDgm/src/solve.jl:68
[18] top-level scope
@ REPL[5]:1