Ensuring differentiability of composed system #90
Description
A standard use case would be gradient based parameter calibration. Thanks to DiffEqSensitivity it should be possible differentiate through most ODEProblems as long as the functions are setup correctly.
On the master branch following the same Cyber-Physical example in the gist we run into following error
Click to expand Error!
ERROR: LoadError: MethodError: no method matching readout(::ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, ::AlgebraicDynamics.DWDDynam.var"#56#58"{Int64, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}, Matrix{Float64}}, ::Tuple{Int64, Int64, Int64}, ::Float64)
Closest candidates are:
readout(::AbstractMachine, ::AbstractVector{T} where T, ::Any, ::Any) at /home/jagupt/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:192
readout(::AbstractMachine, ::Catlab.CategoricalAlgebra.Sets.SetFunction{Dom, Codom} where {S, Dom<:(Catlab.CategoricalAlgebra.FinSets.FinSet{S, T} where T), Codom<:Catlab.CategoricalAlgebra.Sets.SetOb}, ::Any...) at /home/jagupt/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:193
readout(::DelayMachine{T, I} where {T, I}, ::AbstractVector{T} where T, ::Any, ::Any) at /home/jagupt/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:191
...
Stacktrace:
[1] (::AlgebraicDynamics.DWDDynam.var"#59#60"{Vector{AlgebraicDynamics.DWDDynam.var"#56#58"{Int64, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}, Matrix{Float64}}}, Tuple{Tuple{Int64, Int64, Int64}, Float64}})(::Tuple{Int64, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}})
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:425
[2] iterate
@ ./generator.jl:47 [inlined]
[3] collect
@ ./array.jl:681 [inlined]
[4] map
@ ./abstractarray.jl:2323 [inlined]
[5] get_readouts
@ ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:424 [inlined]
[6] (::AlgebraicDynamics.DWDDynam.var"#v#44"{AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}, WiringDiagram{Any, Any, Any, Any}, Vector{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}}, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}})(u::Matrix{Float64}, xs::Matrix{Float64}, p::Tuple{Int64, Int64, Int64}, t::Float64)
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:352
[7] eval_dynamics(f::ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, u::Matrix{Float64}, xs::Matrix{Float64}, p::Tuple{Int64, Int64, Int64}, t::Float64)
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:208
[8] #24
@ ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:250 [inlined]
[9] ODEFunction
@ ~/.julia/packages/SciMLBase/x3z0g/src/scimlfunctions.jl:334 [inlined]
[10] initialize!(integrator::OrdinaryDiffEq.ODEIntegrator{CompositeAlgorithm{Tuple{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, false, Matrix{Float64}, Nothing, Float64, Tuple{Int64, Int64, Int64}, Float64, Float64, Float64, Float64, Vector{Matrix{Float64}}, OrdinaryDiffEq.ODECompositeSolution{Float64, 3, Vector{Vector{Float64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{Float64}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, CompositeAlgorithm{Tuple{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.CompositeInterpolationData{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64}, OrdinaryDiffEq.Rosenbrock23ConstantCache{Float64, SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{Float64}, Tuple{Int64, Int64, Int64}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{Int64, Int64, Int64}}, Matrix{Float64}, LinearAlgebra.LU{Float64, Matrix{Float64}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}}, DiffEqBase.DEStats}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64}, OrdinaryDiffEq.Rosenbrock23ConstantCache{Float64, SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{Float64}, Tuple{Int64, Int64, Int64}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{Int64, Int64, Int64}}, Matrix{Float64}, LinearAlgebra.LU{Float64, Matrix{Float64}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.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, UnitRange{Int64}, Int64, Tuple{}, Vector{Float64}, Tuple{}}, Matrix{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/perform_step/low_order_rk_perform_step.jl:569
[11] initialize!(integrator::OrdinaryDiffEq.ODEIntegrator{CompositeAlgorithm{Tuple{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, false, Matrix{Float64}, Nothing, Float64, Tuple{Int64, Int64, Int64}, Float64, Float64, Float64, Float64, Vector{Matrix{Float64}}, OrdinaryDiffEq.ODECompositeSolution{Float64, 3, Vector{Vector{Float64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{Float64}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, CompositeAlgorithm{Tuple{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.CompositeInterpolationData{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64}, OrdinaryDiffEq.Rosenbrock23ConstantCache{Float64, SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{Float64}, Tuple{Int64, Int64, Int64}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{Int64, Int64, Int64}}, Matrix{Float64}, LinearAlgebra.LU{Float64, Matrix{Float64}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}}, DiffEqBase.DEStats}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64}, OrdinaryDiffEq.Rosenbrock23ConstantCache{Float64, SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{Float64}, Tuple{Int64, Int64, Int64}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{Int64, Int64, Int64}}, Matrix{Float64}, LinearAlgebra.LU{Float64, Matrix{Float64}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.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, UnitRange{Int64}, Int64, Tuple{}, Vector{Float64}, Tuple{}}, Matrix{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Tsit5ConstantCache{Float64, Float64}, OrdinaryDiffEq.Rosenbrock23ConstantCache{Float64, SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{Float64}, Tuple{Int64, Int64, Int64}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{Int64, Int64, Int64}}, Matrix{Float64}, LinearAlgebra.LU{Float64, Matrix{Float64}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/perform_step/composite_perform_step.jl:39
[12] __init(prob::ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::CompositeAlgorithm{Tuple{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{5, false, DefaultLinSolve, Val{:forward}}}, AutoSwitch{Tsit5{typeof(OrdinaryDiffEq.trivial_limiter!), typeof(OrdinaryDiffEq.trivial_limiter!), Static.False}, Rosenbrock23{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, timeseries_init::Tuple{}, ts_init::Tuple{}, ks_init::Tuple{}, recompile::Type{Val{true}}; saveat::Vector{Float64}, tstops::Tuple{}, d_discontinuities::Tuple{}, save_idxs::UnitRange{Int64}, 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::Int64, beta1::Nothing, beta2::Nothing, qoldinit::Rational{Int64}, controller::Nothing, fullnormalize::Bool, failfactor::Int64, maxiters::Int64, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), internalopnorm::typeof(LinearAlgebra.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{Symbol, Any, Tuple{Symbol, Symbol}, NamedTuple{(:default_set, :sense), Tuple{Bool, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/solve.jl:456
[13] #__solve#493
@ ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/solve.jl:4 [inlined]
[14] __solve(::ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, ::Nothing; default_set::Bool, kwargs::Base.Iterators.Pairs{Symbol, Any, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:saveat, :save_idxs, :sense), Tuple{Vector{Float64}, UnitRange{Int64}, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ DifferentialEquations ~/.julia/packages/DifferentialEquations/7WPQg/src/default_solve.jl:8
[15] #solve_call#42
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:61 [inlined]
[16] #solve_up#44
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:93 [inlined]
[17] #solve#43
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:73 [inlined]
[18] _concrete_solve_adjoint(::ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, ::Nothing, ::ForwardDiffSensitivity{0, nothing}, ::Matrix{Float64}, ::Tuple{Int64, Int64, Int64}; saveat::Vector{Float64}, save_idxs::Nothing, kwargs::Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ DiffEqSensitivity ~/.julia/packages/DiffEqSensitivity/uakCr/src/concrete_solve.jl:346
[19] #_concrete_solve_adjoint#219
@ ~/.julia/packages/DiffEqSensitivity/uakCr/src/concrete_solve.jl:66 [inlined]
[20] #_solve_adjoint#61
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:345 [inlined]
[21] #rrule#59
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:310 [inlined]
[22] chain_rrule_kw
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/chainrules.jl:228 [inlined]
[23] macro expansion
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0 [inlined]
[24] _pullback
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:9 [inlined]
[25] _apply
@ ./boot.jl:804 [inlined]
[26] adjoint
@ ~/.julia/packages/Zygote/AlLTp/src/lib/lib.jl:200 [inlined]
[27] _pullback
@ ~/.julia/packages/ZygoteRules/AIbCs/src/adjoint.jl:65 [inlined]
[28] _pullback
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:73 [inlined]
[29] _pullback(::Zygote.Context, ::DiffEqBase.var"##solve#43", ::Nothing, ::Nothing, ::Nothing, ::Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, ::typeof(solve), ::ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[30] _apply(::Function, ::Vararg{Any, N} where N)
@ Core ./boot.jl:804
[31] adjoint
@ ~/.julia/packages/Zygote/AlLTp/src/lib/lib.jl:200 [inlined]
[32] _pullback
@ ~/.julia/packages/ZygoteRules/AIbCs/src/adjoint.jl:65 [inlined]
[33] _pullback
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:68 [inlined]
[34] _pullback(::Zygote.Context, ::CommonSolve.var"#solve##kw", ::NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, ::typeof(solve), ::ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Int64, Int64, Int64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[35] _pullback
@ /mnt/tools/julia/dev/AlgebraicDynamics/examples/Cyber-Physical_param_calibration.jl:116 [inlined]
[36] _pullback(ctx::Zygote.Context, f::var"#105#106", args::Tuple{Float64, Float64, Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[37] _pullback(f::Function, args::Tuple{Float64, Float64, Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface.jl:34
[38] pullback(f::Function, args::Tuple{Float64, Float64, Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface.jl:40
[39] gradient(f::Function, args::Tuple{Float64, Float64, Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface.jl:75
[40] top-level scope
@ /mnt/tools/julia/dev/AlgebraicDynamics/examples/Cyber-Physical_param_calibration.jl:113
Same is true in the batched case on fix_Problem branch with this gist, but we see a more useful error:
Click to expand Error!
[ Info: Precompiling AlgebraicDynamics [5fd6ff03-a254-427e-8840-ba658f502e32]
ERROR: LoadError: TypeError: in typeassert, expected Float64, got a value of type ForwardDiff.Dual{Nothing, Float64, 5}
Stacktrace:
[1] setindex!
@ ./array.jl:845 [inlined]
[2] macro expansion
@ ./multidimensional.jl:903 [inlined]
[3] macro expansion
@ ./cartesian.jl:64 [inlined]
[4] _unsafe_setindex!(::IndexLinear, ::Matrix{Float64}, ::Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ::Int64, ::Base.Slice{Base.OneTo{Int64}})
@ Base ./multidimensional.jl:898
[5] _setindex!(::IndexLinear, ::Matrix{Float64}, ::Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ::Int64, ::Base.Slice{Base.OneTo{Int64}})
@ Base ./multidimensional.jl:887
[6] setindex!
@ ./abstractarray.jl:1267 [inlined]
[7] fill_readins!(readins::Matrix{Float64}, d::WiringDiagram{Any, Any, Any, Any}, Inputs::Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}, readouts::Vector{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, xs::Matrix{Float64})
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:435
[8] (::AlgebraicDynamics.DWDDynam.var"#v#44"{AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}, WiringDiagram{Any, Any, Any, Any}, Vector{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}}, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}, Catlab.CategoricalAlgebra.Limits.Colimit{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FreeDiagrams.DiscreteDiagram{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Any, Vector{Catlab.CategoricalAlgebra.FinSets.FinSetInt}}, Catlab.CategoricalAlgebra.FreeDiagrams.Multicospan{Catlab.CategoricalAlgebra.FinSets.FinSetInt, Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}, Vector{Catlab.CategoricalAlgebra.FinSets.FinDomFunctionVector{Int64, UnitRange{Int64}, Catlab.CategoricalAlgebra.FinSets.FinSetInt}}}}})(u::Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, xs::Matrix{Float64}, p::Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, t::Float64)
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:355
[9] eval_dynamics(f::ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, u::Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, xs::Matrix{Float64}, p::Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, t::Float64)
@ AlgebraicDynamics.DWDDynam ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:208
[10] #24
@ ~/.julia/dev/AlgebraicDynamics/src/dwd_dynam.jl:250 [inlined]
[11] ODEFunction
@ ~/.julia/packages/SciMLBase/x3z0g/src/scimlfunctions.jl:334 [inlined]
[12] ode_determine_initdt(u0::Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, t::Float64, tdir::Float64, dtmax::Float64, abstol::ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, reltol::ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), prob::ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, integrator::OrdinaryDiffEq.ODEIntegrator{CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, false, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Nothing, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Float64, Float64, Float64, Float64, Vector{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.ODECompositeSolution{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, 3, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.CompositeInterpolationData{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}}, DiffEqBase.DEStats}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.DEOptions{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.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, UnitRange{Int64}, Int64, Tuple{}, Vector{Float64}, Tuple{}}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Nothing, OrdinaryDiffEq.DefaultInit})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/initdt.jl:201
[13] auto_dt_reset!(integrator::OrdinaryDiffEq.ODEIntegrator{CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, false, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Nothing, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Float64, Float64, Float64, Float64, Vector{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.ODECompositeSolution{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, 3, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.CompositeInterpolationData{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}}, DiffEqBase.DEStats}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.DEOptions{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.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, UnitRange{Int64}, Int64, Tuple{}, Vector{Float64}, Tuple{}}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Nothing, OrdinaryDiffEq.DefaultInit})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/integrators/integrator_interface.jl:329
[14] handle_dt!(integrator::OrdinaryDiffEq.ODEIntegrator{CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, false, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Nothing, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Float64, Float64, Float64, Float64, Vector{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.ODECompositeSolution{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, 3, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.CompositeInterpolationData{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, Vector{Float64}, Vector{Vector{Vector{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}}, DiffEqBase.DEStats}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.CompositeCache{Tuple{OrdinaryDiffEq.Vern9ConstantCache{OrdinaryDiffEq.Vern9Tableau{Float64, Float64}}, OrdinaryDiffEq.Rosenbrock5ConstantCache{SciMLBase.TimeDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, SciMLBase.UDerivativeWrapper{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, OrdinaryDiffEq.Rodas5Tableau{Float64, Float64}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, LinearAlgebra.LU{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}}, DefaultLinSolve}}, OrdinaryDiffEq.AutoSwitchCache{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, OrdinaryDiffEq.DEOptions{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.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, UnitRange{Int64}, Int64, Tuple{}, Vector{Float64}, Tuple{}}, Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, Nothing, OrdinaryDiffEq.DefaultInit})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/solve.jl:504
[15] __init(prob::ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, AutoSwitch{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}, timeseries_init::Tuple{}, ts_init::Tuple{}, ks_init::Tuple{}, recompile::Type{Val{true}}; saveat::Vector{Float64}, tstops::Tuple{}, d_discontinuities::Tuple{}, save_idxs::UnitRange{Int64}, 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::Int64, beta1::Nothing, beta2::Nothing, qoldinit::Rational{Int64}, controller::Nothing, fullnormalize::Bool, failfactor::Int64, maxiters::Int64, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), internalopnorm::typeof(LinearAlgebra.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{Symbol, Any, Tuple{Symbol, Symbol}, NamedTuple{(:default_set, :sense), Tuple{Bool, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/solve.jl:466
[16] __solve(::ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, ::CompositeAlgorithm{Tuple{Vern9, Rodas5{5, false, DefaultLinSolve, Val{:forward}}}, AutoSwitch{Vern9, Rodas5{0, false, DefaultLinSolve, Val{:forward}}, Rational{Int64}, Int64}}; kwargs::Base.Iterators.Pairs{Symbol, Any, NTuple{4, Symbol}, NamedTuple{(:default_set, :saveat, :save_idxs, :sense), Tuple{Bool, Vector{Float64}, UnitRange{Int64}, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/4BjSL/src/solve.jl:4
[17] __solve(::ODEProblem{Matrix{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, Tuple{Float64, Float64}, false, Tuple{ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}, ForwardDiff.Dual{ForwardDiff.Tag{ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Float64}, Float64, 5}}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, ::Nothing; default_set::Bool, kwargs::Base.Iterators.Pairs{Symbol, Any, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:saveat, :save_idxs, :sense), Tuple{Vector{Float64}, UnitRange{Int64}, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}})
@ DifferentialEquations ~/.julia/packages/DifferentialEquations/7WPQg/src/default_solve.jl:8
[18] #solve_call#42
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:61 [inlined]
[19] #solve_up#44
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:93 [inlined]
[20] #solve#43
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:73 [inlined]
[21] (::DiffEqSensitivity.var"#252#258"{0, Matrix{Float64}, Vector{Float64}, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, Nothing, ForwardDiffSensitivity{0, nothing}, Matrix{Float64}, Tuple{Float64, Float64, Float64}, Tuple{}, UnitRange{Int64}})()
@ DiffEqSensitivity ~/.julia/packages/DiffEqSensitivity/uakCr/src/concrete_solve.jl:521
[22] unthunk
@ ~/.julia/packages/ChainRulesCore/7ZiwT/src/tangent_types/thunks.jl:194 [inlined]
[23] wrap_chainrules_output
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/chainrules.jl:104 [inlined]
[24] map
@ ./tuple.jl:215 [inlined]
[25] map (repeats 3 times)
@ ./tuple.jl:216 [inlined]
[26] wrap_chainrules_output
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/chainrules.jl:105 [inlined]
[27] ZBack
@ ~/.julia/packages/Zygote/AlLTp/src/compiler/chainrules.jl:204 [inlined]
[28] (::Zygote.var"#kw_zpullback#40"{DiffEqSensitivity.var"#forward_sensitivity_backpass#255"{0, Vector{Float64}, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, Nothing, ForwardDiffSensitivity{0, nothing}, Matrix{Float64}, Tuple{Float64, Float64, Float64}, Tuple{}, UnitRange{Int64}}})(dy::Matrix{Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/chainrules.jl:230
[29] (::Zygote.var"#208#209"{Tuple{NTuple{6, Nothing}, Tuple{}}, Zygote.var"#kw_zpullback#40"{DiffEqSensitivity.var"#forward_sensitivity_backpass#255"{0, Vector{Float64}, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, Nothing, ForwardDiffSensitivity{0, nothing}, Matrix{Float64}, Tuple{Float64, Float64, Float64}, Tuple{}, UnitRange{Int64}}}})(Δ::Matrix{Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/lib/lib.jl:203
[30] (::Zygote.var"#1733#back#210"{Zygote.var"#208#209"{Tuple{NTuple{6, Nothing}, Tuple{}}, Zygote.var"#kw_zpullback#40"{DiffEqSensitivity.var"#forward_sensitivity_backpass#255"{0, Vector{Float64}, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, Nothing, ForwardDiffSensitivity{0, nothing}, Matrix{Float64}, Tuple{Float64, Float64, Float64}, Tuple{}, UnitRange{Int64}}}}})(Δ::Matrix{Float64})
@ Zygote ~/.julia/packages/ZygoteRules/AIbCs/src/adjoint.jl:67
[31] Pullback
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:73 [inlined]
[32] (::Zygote.Pullback{Tuple{DiffEqBase.var"##solve#43", Nothing, Nothing, Nothing, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, typeof(solve), ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}}, Any})(Δ::Matrix{Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[33] #208
@ ~/.julia/packages/Zygote/AlLTp/src/lib/lib.jl:203 [inlined]
[34] (::Zygote.var"#1733#back#210"{Zygote.var"#208#209"{Tuple{NTuple{6, Nothing}, Tuple{}}, Zygote.Pullback{Tuple{DiffEqBase.var"##solve#43", Nothing, Nothing, Nothing, Base.Iterators.Pairs{Symbol, InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}, Tuple{Symbol}, NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, typeof(solve), ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}}, Any}}})(Δ::Matrix{Float64})
@ Zygote ~/.julia/packages/ZygoteRules/AIbCs/src/adjoint.jl:67
[35] Pullback
@ ~/.julia/packages/DiffEqBase/YIT73/src/solve.jl:68 [inlined]
[36] (::Zygote.Pullback{Tuple{CommonSolve.var"#solve##kw", NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, typeof(solve), ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}}, Any})(Δ::Matrix{Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[37] Pullback
@ /mnt/tools/julia/dev/AlgebraicDynamics/examples/Cyber-Physical.jl:112 [inlined]
[38] (::Zygote.Pullback{Tuple{var"#11#12", Tuple{Float64, Float64, Float64}}, Tuple{Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Matrix{Float64}}}, Zygote.var"#1613#back#155"{typeof(identity)}, Zygote.var"#2163#back#380"{Zygote.var"#378#379"}, Zygote.Pullback{Tuple{Type{ODEProblem}, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}, Matrix{Float64}, Tuple{Int64, Float64}, Tuple{Float64, Float64, Float64}}, Tuple{Zygote.Pullback{Tuple{AlgebraicDynamics.DWDDynam.var"##ODEProblem#23", Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Type{ODEProblem}, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}, Matrix{Float64}, Tuple{Int64, Float64}, Tuple{Float64, Float64, Float64}}, Any}, Zygote.Pullback{Tuple{Type{NamedTuple}}, Tuple{}}, Zygote.var"#2076#back#357"{Zygote.var"#pairs_namedtuple_pullback#356"{(), NamedTuple{(), Tuple{}}}}}}, Zygote.Pullback{Tuple{CommonSolve.var"#solve##kw", NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, typeof(solve), ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}}, Any}, Zygote.Pullback{Tuple{Type{NamedTuple{(:sense,), T} where T<:Tuple}, Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Tuple{Zygote.Pullback{Tuple{Type{NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Tuple{Zygote.var"#1777#back#226"{Zygote.Jnew{NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Nothing, true}}}}}}, Zygote.ZBack{ChainRules.var"#sum_abs2_pullback#1395"{Array{Float64, 3}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Matrix{Float64}}}, Zygote.var"#1613#back#155"{typeof(identity)}, Zygote.Pullback{Tuple{Type{ZygoteVJP}}, Tuple{}}, Zygote.Pullback{Tuple{Core.var"#Type##kw", NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}, Type{InterpolatingAdjoint}}, Any}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Tuple{Int64, Float64}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Array{Float64, 3}}}, Zygote.Pullback{Tuple{typeof(Base.Broadcast.materialize), Array{Float64, 3}}, Tuple{}}, Zygote.var"#3766#back#1024"{Zygote.var"#1022#1023"{Array{Float64, 3}, Matrix{Float64}}}, Zygote.Pullback{Tuple{Type{NamedTuple{(:autojacvec,), T} where T<:Tuple}, Tuple{ZygoteVJP}}, Tuple{Zygote.Pullback{Tuple{Type{NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}}, Tuple{ZygoteVJP}}, Tuple{Zygote.var"#1777#back#226"{Zygote.Jnew{NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}, Nothing, true}}}}}}}})(Δ::Float64)
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface2.jl:0
[39] (::Zygote.var"#55#56"{Zygote.Pullback{Tuple{var"#11#12", Tuple{Float64, Float64, Float64}}, Tuple{Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Matrix{Float64}}}, Zygote.var"#1613#back#155"{typeof(identity)}, Zygote.var"#2163#back#380"{Zygote.var"#378#379"}, Zygote.Pullback{Tuple{Type{ODEProblem}, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}, Matrix{Float64}, Tuple{Int64, Float64}, Tuple{Float64, Float64, Float64}}, Tuple{Zygote.Pullback{Tuple{AlgebraicDynamics.DWDDynam.var"##ODEProblem#23", Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Type{ODEProblem}, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}, Matrix{Float64}, Tuple{Int64, Float64}, Tuple{Float64, Float64, Float64}}, Any}, Zygote.Pullback{Tuple{Type{NamedTuple}}, Tuple{}}, Zygote.var"#2076#back#357"{Zygote.var"#pairs_namedtuple_pullback#356"{(), NamedTuple{(), Tuple{}}}}}}, Zygote.Pullback{Tuple{CommonSolve.var"#solve##kw", NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, typeof(solve), ODEProblem{Matrix{Float64}, Tuple{Float64, Float64}, false, Tuple{Float64, Float64, Float64}, ODEFunction{false, AlgebraicDynamics.DWDDynam.var"#24#25"{ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}, Matrix{Float64}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}}, Any}, Zygote.Pullback{Tuple{Type{NamedTuple{(:sense,), T} where T<:Tuple}, Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Tuple{Zygote.Pullback{Tuple{Type{NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}}, Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Tuple{Zygote.var"#1777#back#226"{Zygote.Jnew{NamedTuple{(:sense,), Tuple{InterpolatingAdjoint{0, true, Val{:central}, ZygoteVJP, Bool}}}, Nothing, true}}}}}}, Zygote.ZBack{ChainRules.var"#sum_abs2_pullback#1395"{Array{Float64, 3}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Matrix{Float64}}}, Zygote.var"#1613#back#155"{typeof(identity)}, Zygote.Pullback{Tuple{Type{ZygoteVJP}}, Tuple{}}, Zygote.Pullback{Tuple{Core.var"#Type##kw", NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}, Type{InterpolatingAdjoint}}, Any}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, ContinuousMachine{Float64, AlgebraicDynamics.DWDDynam.DirectedInterface{Float64}}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Tuple{Int64, Float64}}}, Zygote.var"#1590#back#149"{Zygote.var"#147#148"{Zygote.Context, GlobalRef, Array{Float64, 3}}}, Zygote.Pullback{Tuple{typeof(Base.Broadcast.materialize), Array{Float64, 3}}, Tuple{}}, Zygote.var"#3766#back#1024"{Zygote.var"#1022#1023"{Array{Float64, 3}, Matrix{Float64}}}, Zygote.Pullback{Tuple{Type{NamedTuple{(:autojacvec,), T} where T<:Tuple}, Tuple{ZygoteVJP}}, Tuple{Zygote.Pullback{Tuple{Type{NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}}, Tuple{ZygoteVJP}}, Tuple{Zygote.var"#1777#back#226"{Zygote.Jnew{NamedTuple{(:autojacvec,), Tuple{ZygoteVJP}}, Nothing, true}}}}}}}}})(Δ::Float64)
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface.jl:41
[40] gradient(f::Function, args::Tuple{Float64, Float64, Float64})
@ Zygote ~/.julia/packages/Zygote/AlLTp/src/compiler/interface.jl:76
[41] top-level scope
@ /mnt/tools/julia/dev/AlgebraicDynamics/examples/Cyber-Physical.jl:109