Issue when executing example

[felix-russo]

Hi!

I freshly installed ChaosTools (v3.5.1) and wanted to test it using one of the examples that you provide:

using ChaosTools

function lorenz_rule(u, p, t)
       σ = p[1]; ρ = p[2]; β = p[3]
       du1 = σ * (u[2]-u[1])
       du2 = u[1] * (ρ-u[3]) - u[2]
       du3 = u[1] * u[2] - β * u[3]
       return SVector{3}(du1, du2, du3)
end

lor = CoupledODEs(lorenz_rule, fill(10.0, 3), [10, 32, 8/3])
λλ = lyapunovspectrum(lor, 10000; Δt = 0.1)

However, I get

ERROR: MethodError: no method matching reinit!(::OrdinaryDiffEqCore.ODEIntegrator{…}, ::SMatrix{…}; reinit_dae::Bool, reset_dt::Bool, t0::Float64)
This error has been manually thrown, explicitly, so the method may exist but be intentionally marked as unimplemented.

Closest candidates are:
  reinit!(::OrdinaryDiffEqCore.ODEIntegrator, ::Any; t0, tf, erase_sol, tstops, saveat, d_discontinuities, reset_dt, reinit_callbacks, initialize_save, reinit_cache, reinit_retcode) got unsupported keyword argument "reinit_dae"
   @ OrdinaryDiffEqCore ~/.julia/packages/OrdinaryDiffEqCore/gALQU/src/integrators/integrator_interface.jl:324
  reinit!(::SciMLBase.DEIntegrator, ::Any...; kwargs...)
   @ SciMLBase ~/.julia/packages/SciMLBase/VAClc/src/integrator_interface.jl:281
  reinit!(::DiffEqBase.NonlinearTerminationModeCache{dep_retcode}, ::Any, ::Any, ::Any...; abstol, reltol, kwargs...) where dep_retcode
   @ DiffEqBase ~/.julia/packages/DiffEqBase/HW4ge/src/termination_conditions.jl:258
  ...

Stacktrace:
 [1] kwerr(::@NamedTuple{…}, ::Function, ::OrdinaryDiffEqCore.ODEIntegrator{…}, ::SMatrix{…})
   @ Base ./error.jl:165
 [2] reinit!(ds::CoupledODEs{…}, u::SMatrix{…}; p::Vector{…}, t0::Float64, kw::@Kwargs{})
   @ DynamicalSystemsBase ~/.julia/packages/DynamicalSystemsBase/GpQyA/src/core_systems/continuous_time_ode.jl:148
 [3] reinit!(tands::TangentDynamicalSystem{…}, u::SVector{…}; p::Vector{…}, t0::Float64, Q0::Matrix{…})
   @ DynamicalSystemsBase ~/.julia/packages/DynamicalSystemsBase/GpQyA/src/derived_systems/tangent_space.jl:246
 [4] reinit!
   @ ~/.julia/packages/DynamicalSystemsBase/GpQyA/src/derived_systems/tangent_space.jl:233 [inlined]
 [5] lyapunovspectrum(tands::TangentDynamicalSystem{…}, N::Int64; Δt::Float64, Ttr::Int64, show_progress::Bool, u0::SVector{…})
   @ ChaosTools ~/.julia/packages/ChaosTools/PXjUl/src/chaosdetection/lyapunovs/lyapunovspectrum.jl:71
 [6] lyapunovspectrum(ds::CoupledODEs{…}, N::Int64, k::Int64; kwargs::@Kwargs{…})
   @ ChaosTools ~/.julia/packages/ChaosTools/PXjUl/src/chaosdetection/lyapunovs/lyapunovspectrum.jl:54
 [7] lyapunovspectrum
   @ ~/.julia/packages/ChaosTools/PXjUl/src/chaosdetection/lyapunovs/lyapunovspectrum.jl:52 [inlined]
 [8] top-level scope
   @ ~/Code/play.jl:12
Some type information was truncated. Use `show(err)` to see complete types.

Not sure if this is a bug or I am doing something wrong!

Thanks for your help!

[Datseris]

Hi your code works perfectly fine for me:

julia> function lorenz_rule(u, p, t)
              σ = p[1]; ρ = p[2]; β = p[3]
              du1 = σ * (u[2]-u[1])
                     du2 = u[1] * (ρ-u[3]) - u[2]
              du3 = u[1] * u[2] - β * u[3]
              return SVector{3}(du1, du2, du3)
       end
lorenz_rule (generic function with 1 method)

julia> 

julia> lor = CoupledODEs(lorenz_rule, fill(10.0, 3), [10, 32, 8/3])      
3-dimensional CoupledODEs
 deterministic: true
 discrete time: false
 in-place:      false
 dynamic rule:  lorenz_rule
 ODE solver:    Tsit5
 ODE kwargs:    (abstol = 1.0e-6, reltol = 1.0e-6)
 parameters:    [10.0, 32.0, 2.6666666666666665]
 time:          0.0
 state:         [10.0, 10.0, 10.0]


julia> λλ = lyapunovspectrum(lor, 10000; Δt = 0.1)
3-element Vector{Float64}:
   0.9905076825780865
   0.0037076648444580836
 -14.660818767837927

can you post other package versions you use? Manifest in particular? e.g., pkg> st -m DynamicalSystemsBase OrdinaryDiffEqCore.

[felix-russo]

Update OrdinaryDiffEqCore to Version 1.36, all is resolved.

Thanks for the swift reply and the help!