Closed oameye closed 6 months ago
The problem is solved when one defines the Mathieu function out-of-place.
using DynamicalSystems
function mathieu(x,p,t)
dx1 = x[2]
dx2 = -p[:ω₀]^2*(1-p[:λ]*cos(2*p[:ω]*t))*x[1] - p[:γ]*x[2]
return SA[dx1, dx2]
end
fixed = Dict(:ω₀ => 1.0, :γ => 0.005, :ω => 1.0, :λ => 0.03)
ds = ContinuousDynamicalSystem(mathieu, [1.0, 0], fixed)
tr = trajectory(ds, 10.0)
local_growth_rates(ds, tr)
problem will be solved in DynamicalSystems.jl v3.0, coming out by the end of the month
Describe the bug When computing the
local_growth_rates
of the continuous Mathieu differential equation, I get the error that no method exists forrecursivecopy!(::Matrix{Float64}, ::Vector{SVector{2, Float64}})
.Minimal Working Example
Package versions
When running:
I get: