# BifurcationKit MethodError

Hi all. Currently, I have a reaction network defined in Catalyst, and I’ve been following the tutorial at Bifurcation Diagrams · Catalyst.jl. However, I’ve been running into an error and am not entirely sure how to resolve it.

``````begin
odefun = ODEFunction(convert(ODESystem,rn), jac=true)
F = (u,p) -> odefun(u,p,0)
J = (u,p) -> odefun.jac(u,p,0)
end
``````
``````begin
params = [92.331732,0.86466471,79.9512906,97.932525,1,3.4134082,0.61409879,0.61409879,3.4134082,0.98168436,0.98168436,4.7267833,0.17182818,0.68292191,1,3.2654672,0.61699064,0.61699064,37.913879,0.86466471,0.86466471,0.99326205,0.99326205,1,5.9744464,1.7182818,1.7182818,1.7182818,1.7182818,0.55950727,1.0117639]
p_idx = 25
p_span = (0.1,20.)
plot_var_idx = 19
end;
``````
``````begin
opts = ContinuationPar( dsmax = 0.05,        # Maximum arclength value of the pseudo-arc length continuation method.
dsmin = 1e-4,        # Minimum arclength value of the pseudo-arc length continuation method.
ds=0.001,            # Initial arclength value of the pseudo-arc length continuation method (should be positive).
maxSteps = 100000,   # The maximum number of steps.
pMin = p_span[1],    # Minimum p-vale (if hit, the method stops).
pMax = p_span[2],    # Maximum p-vale (if hit, the method stops).
detectBifurcation=3, # Value in {0,1,2,3} determening to what extent bofurcation points are detected (0 means nothing is done, 3 both them and there localisation are detected).
newtonOptions = NewtonPar(tol = 1e-9, verbose = false, maxIter = 15)) #Parameters to the newton solver (when finding fixed points) see BifurcationKit documentation.

DO = DeflationOperator( 2,      # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
dot,    # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
1.,     # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
[fill(0.,length(rn.states))]); # Guess(es) of the fixed point for the initial parameter set. Do not need to be exact.
end
``````
``````begin
params_input = setindex!(copy(params),p_span[1],p_idx)                                # The input parameter values have to start at the first index of our parameter span.
branches, = continuation(F, J, params_input, (@lens _[p_idx]), opts , DO,             # Gives our input.
verbosity = 0, plot=false,                                                    # We do not want to display, or plot, intermediary results.
recordFromSolution = (x, p) -> x[plot_var_idx],                                   # How we wish to print the output in the diagram. Here we simply want the value of the target varriable.
perturbSolution = (x,p,id) -> (x  .+ 0.8 .* rand(length(x))),                     # Parameter for the continuation method, see BifurcationKit documentation.
callbackN = BifurcationKit.cbMaxNorm(1e7))                      # Parameter for the continuation method, see BifurcationKit documentation.
end
``````
``````MethodError: no method matching continuation(::var"#81#83"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#462"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x022aeffd, 0xe8b3b74e, 0xce2bcce9, 0x51b2ee73, 0x5b2ef1bd)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#466"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x9ca978cf, 0x116effc5, 0x2d7e2d5e, 0x86d966f8, 0x3b1dd272)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#487#generated_observed#469"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::var"#82#84"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#462"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x022aeffd, 0xe8b3b74e, 0xce2bcce9, 0x51b2ee73, 0x5b2ef1bd)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#466"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x9ca978cf, 0x116effc5, 0x2d7e2d5e, 0x86d966f8, 0x3b1dd272)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#487#generated_observed#469"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::Vector{Float64}, ::Setfield.IndexLens{Tuple{Int64}}, ::BifurcationKit.ContinuationPar{Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}}, ::BifurcationKit.DeflationOperator{Int64, typeof(LinearAlgebra.dot), Float64, Vector{Float64}}; verbosity=0, plot=false, recordFromSolution=var"#77#79"{Int64}(19), perturbSolution=var"#78#80"{typeof(+), typeof(rand), typeof(*), typeof(length)}(+, rand, *, length), callbackN=BifurcationKit.cbMaxNorm{Float64}(1.0e7))

Closest candidates are:

continuation(!Matched::BifurcationKit.AbstractBifurcationProblem, ::Tv, ::Any, !Matched::Tv, !Matched::Real, ::Any, !Matched::Setfield.Lens, !Matched::BifurcationKit.ContinuationPar; kwargs...) where Tv at C:\Users\Dragolonth\.julia\packages\BifurcationKit\fnCOr\src\bifdiagram\BranchSwitching.jl:8

continuation(!Matched::BifurcationKit.ContIterable, ::Any, !Matched::Real, ::Any, !Matched::Real) at C:\Users\Dragolonth\.julia\packages\BifurcationKit\fnCOr\src\bifdiagram\BranchSwitching.jl:25 got unsupported keyword arguments "verbosity", "plot", "recordFromSolution", "perturbSolution", "callbackN"

continuation(!Matched::BifurcationKit.TWProblem, ::Any, !Matched::BifurcationKit.AbstractContinuationAlgorithm, !Matched::BifurcationKit.ContinuationPar; kwargs...) at C:\Users\Dragolonth\.julia\packages\BifurcationKit\fnCOr\src\wave\WaveProblem.jl:212

``````

Hi,

It seems you have istalled BifurcationKit@0.2.0 on your system but your code looks like it was written for BifurcationKit@0.1.12. So either do

`] add BifurcationKit@0.1.12`

or modify your code for the new version 0.2.0 for which you need to define a `BifurcationProblem`

Thank you for your help @rveltz. I added v0.1.12 of BifurcationKit, and after running the same code, I was met with this error:

``````LinearAlgebra.SingularException(1)

checknonsingular@factorization.jl:19[inlined]
checknonsingular@factorization.jl:21[inlined]
#lu!#146@lu.jl:82[inlined]
#lu#153@lu.jl:279[inlined]
lu@lu.jl:278[inlined]
\(::Matrix{Float64}, ::Vector{Float64})@generic.jl:1142
#_#5@LinearSolver.jl:81[inlined]
DefaultLS@LinearSolver.jl:81[inlined]
var"#newton#59"(::typeof(LinearAlgebra.norm), ::BifurcationKit.cbMaxNorm{Float64}, ::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol}, NamedTuple{(:iterationC, :p), Tuple{Int64, Float64}}}, ::typeof(BifurcationKit.newton), ::var"#1#3"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::var"#2#4"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::Vector{Float64}, ::Vector{Float64}, ::BifurcationKit.NewtonPar{Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}})@Newton.jl:121
var"#iterate#141"(::Int64, ::typeof(iterate), ::BifurcationKit.ContIterable{var"#1#3"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, var"#2#4"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, Vector{Float64}, Vector{Float64}, Setfield.IndexLens{Tuple{Int64}}, Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}, BifurcationKit.SecantPred, BifurcationKit.BorderingBLS{BifurcationKit.DefaultLS, Float64}, BifurcationKit.var"#586#602"{BifurcationKit.var"#586#587#603"}, var"#5#7"{Int64}, typeof(LinearAlgebra.norm), BifurcationKit.DotTheta{BifurcationKit.var"#168#170", BifurcationKit.var"#169#171"}, typeof(BifurcationKit.finaliseDefault), BifurcationKit.cbMaxNorm{Float64}, Nothing, Nothing})@Continuation.jl:276
iterate@Continuation.jl:259[inlined]
getStatesContResults(::BifurcationKit.DefContIterable{BifurcationKit.ContIterable{var"#1#3"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, var"#2#4"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, Vector{Float64}, Vector{Float64}, Setfield.IndexLens{Tuple{Int64}}, Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}, BifurcationKit.SecantPred, BifurcationKit.BorderingBLS{BifurcationKit.DefaultLS, Float64}, BifurcationKit.var"#586#602"{BifurcationKit.var"#586#587#603"}, var"#5#7"{Int64}, typeof(LinearAlgebra.norm), BifurcationKit.DotTheta{BifurcationKit.var"#168#170", BifurcationKit.var"#169#171"}, typeof(BifurcationKit.finaliseDefault), BifurcationKit.cbMaxNorm{Float64}, Nothing, Nothing}, var"#6#8"{typeof(+), typeof(rand), typeof(*), typeof(length)}, BifurcationKit.var"#591#607", BifurcationKit.var"#592#608"}, ::Vector{Vector{Float64}})@DeflatedContinuation.jl:73
deflatedContinuation(::BifurcationKit.DefContIterable{BifurcationKit.ContIterable{var"#1#3"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, var"#2#4"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, Vector{Float64}, Vector{Float64}, Setfield.IndexLens{Tuple{Int64}}, Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}, BifurcationKit.SecantPred, BifurcationKit.BorderingBLS{BifurcationKit.DefaultLS, Float64}, BifurcationKit.var"#586#602"{BifurcationKit.var"#586#587#603"}, var"#5#7"{Int64}, typeof(LinearAlgebra.norm), BifurcationKit.DotTheta{BifurcationKit.var"#168#170", BifurcationKit.var"#169#171"}, typeof(BifurcationKit.finaliseDefault), BifurcationKit.cbMaxNorm{Float64}, Nothing, Nothing}, var"#6#8"{typeof(+), typeof(rand), typeof(*), typeof(length)}, BifurcationKit.var"#591#607", BifurcationKit.var"#592#608"}, ::BifurcationKit.DeflationOperator{Int64, typeof(LinearAlgebra.dot), Float64, Vector{Float64}}, ::BifurcationKit.ContinuationPar{Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}}, ::Int64, ::Bool)@DeflatedContinuation.jl:170
var"#continuation#576"(::Int64, ::Int64, ::Int64, ::Int64, ::Bool, ::BifurcationKit.SecantPred, ::BifurcationKit.BorderingBLS{BifurcationKit.DefaultLS, Float64}, ::BifurcationKit.DotTheta{BifurcationKit.var"#168#170", BifurcationKit.var"#169#171"}, ::var"#5#7"{Int64}, ::BifurcationKit.var"#586#602"{BifurcationKit.var"#586#587#603"}, ::var"#6#8"{typeof(+), typeof(rand), typeof(*), typeof(length)}, ::BifurcationKit.cbMaxNorm{Float64}, ::BifurcationKit.var"#591#607", ::BifurcationKit.var"#592#608", ::typeof(LinearAlgebra.norm), ::typeof(BifurcationKit.continuation), ::var"#1#3"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::var"#2#4"{SciMLBase.ODEFunction{true, ModelingToolkit.var"#f#312"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xd6f5a074, 0x205d69d5, 0xe5c13464, 0xf31a028a, 0xc5bbe85a)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xa8c000fc, 0xa8412c32, 0x5a4f87dc, 0x7000a9b4, 0x60d8d03c)}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, ModelingToolkit.var"#_jac#316"{RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0xec3151cf, 0x34c553f9, 0xd483c496, 0xd11d0a74, 0x4553ee4f)}, RuntimeGeneratedFunctions.RuntimeGeneratedFunction{(:ˍ₋out, :ˍ₋arg1, :ˍ₋arg2, :t), ModelingToolkit.var"#_RGF_ModTag", ModelingToolkit.var"#_RGF_ModTag", (0x2e0c14b2, 0x6fc06165, 0x41e5edb9, 0x84d3450f, 0x2916623f)}}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, ModelingToolkit.var"#298#generated_observed#319"{Bool, ModelingToolkit.ODESystem, Dict{Any, Any}}, Nothing}}, ::Vector{Float64}, ::Setfield.IndexLens{Tuple{Int64}}, ::BifurcationKit.ContinuationPar{Float64, BifurcationKit.DefaultLS, BifurcationKit.DefaultEig{typeof(real)}}, ::BifurcationKit.DeflationOperator{Int64, typeof(LinearAlgebra.dot), Float64, Vector{Float64}})@DeflatedContinuation.jl:165
top-level scope@Local: 3[inlined]
``````

Do you have any ideas as to what might be causing this issue?

It is a mistake in BK. I have not pushed the fix, for now you can use detectBifurcation=1 or use the previous version.i apologize because I am away from my computer for a whike