`LinearAlgebra.jl` is handled in a wrong way with kernel 1.11

[horvathadam314]

The following lines shall be added to the beginning of NormalForms.jl.

using LinearAlgebra
const LA = LinearAlgebra 

Without this, the following error occurs:

UndefVarError: `dot` not defined in `DDEBifurcationKit`
Suggestion: check for spelling errors or missing imports.
Hint: a global variable of this name also exists in LinearAlgebra.
Hint: a global variable of this name may be made accessible by importing ColorVectorSpace in the current active module Main

[rveltz]

What is kernel 1.11?

[horvathadam314]

What is kernel 1.11?

Sorry, I meant Julia v1.11. I haven't test it with v1.12, but with v1.11 the tutorial code that is presented for DDE-s returns with the error message I've mentioned above. Adding the lines mentioned in the issue solved the problem for me locally.

[rveltz]

tutorial code that is presented for DDE-s returns with the error message I've mentioned above

can you give the link to the tutorial?

Also, what is your BK version?

[horvathadam314]

The details are listed below. Let me know if there is anything else I can help you with.

can you give the link to the tutorial?

Here you can find the link for the tutorial that I've used: link .

Also, what is your BK version?

• BK version:

  [0f109fa4] BifurcationKit v0.5.6

• DDEBK version:

  [954e4062] DDEBifurcationKit v0.0.5 `~/.julia/dev/DDEBifurcationKit`

• julia version and other information that may be useful

julia> versioninfo()
Julia Version 1.11.7
Commit f2b3dbda30a (2025-09-08 12:10 UTC)
Build Info:
  Official https://julialang.org/ release
Platform Info:
  OS: macOS (arm64-apple-darwin24.0.0)
  CPU: 8 × Apple M2
  WORD_SIZE: 64
  LLVM: libLLVM-16.0.6 (ORCJIT, apple-m2)
Threads: 1 default, 0 interactive, 1 GC (on 4 virtual cores)

Code and output

This is the code I use right now. It is the verbatim copy of the tutorial that is linked above.

Code

using Revise, DDEBifurcationKit, Plots
using BifurcationKit
const BK = BifurcationKit

function neuronVF(x, xd, p)
   (; κ, β, a12, a21, τs, τ1, τ2) = p
   [
      -κ * x[1] + β * tanh(xd[3][1]) + a12 * tanh(xd[2][2]),
      -κ * x[2] + β * tanh(xd[3][2]) + a21 * tanh(xd[1][1])
   ]
end

delaysF(par) = [par.τ1, par.τ2, par.τs]

pars = (κ = 0.5, β = -1, a12 = 1, a21 = 0.5, τ1 = 0.2, τ2 = 0.2, τs = 1.5)
x0 = [0.01, 0.001]

prob = ConstantDDEBifProblem(neuronVF, delaysF, x0, pars, (@optic _.τs))

optn = NewtonPar(eigsolver = DDE_DefaultEig(maxit=200))
opts = ContinuationPar(p_max = 13., p_min = 0., newton_options = optn, ds = -0.01, detect_bifurcation = 3, nev = 5, dsmax = 0.2, n_inversion = 4)
br = continuation(prob, PALC(), opts; verbosity = 0, plot = true, bothside = true, normC = norminf)

scene = plot(br)

hopfpt = BK.get_normal_form(br, 2)

# continuation of the first Hopf point
brhopf = continuation(br, 3, (@optic _.a21),
         ContinuationPar(br.contparams, detect_bifurcation = 1, dsmax = 0.04, max_steps = 230, p_max = 15., p_min = -1.,ds = -0.02);
         detect_codim2_bifurcation = 2,
         # bothside = true,
         start_with_eigen = true)

# continuation of the second Hopf point
brhopf2 = continuation(br, 2, (@optic _.a21),
         ContinuationPar(br.contparams, detect_bifurcation = 1, dsmax = 0.1, max_steps = 56, p_max = 15., p_min = -1.,ds = -0.01, n_inversion = 4);
         detect_codim2_bifurcation = 2,
         start_with_eigen = true,
         bothside=true)

scene = plot(brhopf, brhopf2, legend = :top)

prob2 = ConstantDDEBifProblem(neuronVF, delaysF, x0, pars, (@optic _.a21))
br2 = BK.continuation(prob2, PALC(), ContinuationPar(opts, ds = 0.1, p_max = 3., n_inversion = 4); verbosity = 0, plot = false, normC = norminf)

# continuation parameters
opts_po_cont = ContinuationPar(dsmax = 0.1, ds= -0.0001, dsmin = 1e-4, p_max = 10., p_min=-0., max_steps = 120, detect_bifurcation = 0, save_sol_every_step=1)
@reset opts_po_cont.newton_options.tol = 1e-8
@reset opts_po_cont.newton_options.verbose = false

# arguments for periodic orbits
args_po = (	record_from_solution = (x, p; k...) -> begin
			xtt = BK.get_periodic_orbit(p.prob, x, nothing)
			return (max = maximum(xtt[1,:]),
					min = minimum(xtt[1,:]),
					period = getperiod(p.prob, x, nothing))
		end,
		plot_solution = (x, p; k...) -> begin
			xtt = BK.get_periodic_orbit(p.prob, x, nothing)
			plot!(xtt.t, xtt[1,:]; label = "V1", k...)
			plot!(xtt.t, xtt[2,:]; label = "V2", k...)
			plot!(br2; subplot = 1, putspecialptlegend = false)
			end,
		normC = norminf)

probpo = PeriodicOrbitOCollProblem(60, 4; N = 2, jacobian = BK.AutoDiffDense())
br_pocoll = @time continuation(
	br2, 1, opts_po_cont,
	probpo;
	verbosity = 0,	plot = false,
	args_po...,
	δp = 0.003,
	normC = norminf,
	)
scene = plot(br2, br_pocoll)

scene = plot(layout = 2)
for ii = 1:10:110
	solpo = BK.get_periodic_orbit(br_pocoll.γ.prob.prob, br_pocoll.sol[ii].x, nothing)
	plot!(scene, solpo.t ./ solpo.t[end], solpo.u[1,:], label = "", subplot = 1)
end
xlabel!(scene, "t / period", subplot = 1)
plot!(scene, br_pocoll, vars = (:param, :period), subplot = 2, xlims=(2.2,2.4))
scene

Error

UndefVarError: `dot` not defined in `DDEBifurcationKit`
Suggestion: check for spelling errors or missing imports.
Hint: a global variable of this name may be made accessible by importing LinearAlgebra in the current active module Main
Hint: a global variable of this name may be made accessible by importing ColorVectorSpace in the current active module Main

Stacktrace:
  [1] macro expansion
    @ ./logging/logging.jl:385 [inlined]
  [2] hopf_normal_form(prob::ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, pt::BifurcationKit.Hopf{Vector{Float64}, Nothing, Float64, Float64, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, Vector{ComplexF64}, Vector{ComplexF64}, @NamedTuple{a::ComplexF64, b::ComplexF64}}, ls::DefaultLS; autodiff::Bool, verbose::Bool)
    @ DDEBifurcationKit ~/.julia/dev/DDEBifurcationKit/src/NormalForms.jl:52
  [3] (::DDEBifurcationKit.var"#testBT_GH#132"{typeof(LinearAlgebra.norm), @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, NewtonPar{Float64, DefaultLS, DDE_DefaultEig{Float64, typeof(real), Nothing}}})(iter::ContIterable{BifurcationKit.HopfCont, BifurcationKit.HopfMAProblem{DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, BifurcationKit.AutoDiff, RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:a21}, typeof(BifurcationKit.plot_default), DDEBifurcationKit.var"#125#133"{DDEBifurcationKit.var"#125#126#134"{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Tuple{Symbol, Symbol}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}}}}, PALC{Secant, BorderingBLS{DefaultLS, Float64, typeof(VectorInterface.inner), typeof(LinearAlgebra.norm)}, Float64, BifurcationKit.DotTheta{BifurcationKit.var"#344#346", BifurcationKit.var"#345#347"}}, Float64, DefaultLS, DDEBifurcationKit.HopfDDEEig{DDE_DefaultEig{Float64, typeof(real), Nothing}}, typeof(LinearAlgebra.norm), DDEBifurcationKit.var"#updateMinAugHopf#129"{DDEBifurcationKit.var"#updateMinAugHopf#123#130"{Int64, @Kwargs{bothside::Bool}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, Float64}}, typeof(BifurcationKit.cb_default), PairOfEvents{ContinuousEvent{ComposedFunction{typeof(BifurcationKit.convert_to_tuple_eve), DDEBifurcationKit.var"#testBT_GH#132"{typeof(LinearAlgebra.norm), @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, NewtonPar{Float64, DefaultLS, DDE_DefaultEig{Float64, typeof(real), Nothing}}}}, Tuple{String, String}, Float64, typeof(BifurcationKit.default_finalise_event!), Nothing}, BifurcationKit.BifEvent{typeof(BifurcationKit.detect_bifurcation_event)}}}, state::ContState{BorderedArray{RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, Float64}, Float64, Vector{ComplexF64}, Matrix{ComplexF64}, Tuple{Tuple{Float64, Float64, Int64}, Tuple{Float64, Float64, Int64}}})
    @ DDEBifurcationKit ~/.julia/dev/DDEBifurcationKit/src/codim2/Hopf.jl:229
  [4] call_composed
    @ ./operators.jl:1054 [inlined]
  [5] call_composed
    @ ./operators.jl:1053 [inlined]
  [6] ComposedFunction
    @ ./operators.jl:1050 [inlined]
  [7] AbstractEvent
    @ ~/.julia/packages/BifurcationKit/Y3DPO/src/events/Event.jl:6 [inlined]
  [8] PairOfEvents
    @ ~/.julia/packages/BifurcationKit/Y3DPO/src/events/Event.jl:222 [inlined]
  [9] update_event!(it::ContIterable{BifurcationKit.HopfCont, BifurcationKit.HopfMAProblem{DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, BifurcationKit.AutoDiff, RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:a21}, typeof(BifurcationKit.plot_default), DDEBifurcationKit.var"#125#133"{DDEBifurcationKit.var"#125#126#134"{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Tuple{Symbol, Symbol}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}}}}, PALC{Secant, BorderingBLS{DefaultLS, Float64, typeof(VectorInterface.inner), typeof(LinearAlgebra.norm)}, Float64, BifurcationKit.DotTheta{BifurcationKit.var"#344#346", BifurcationKit.var"#345#347"}}, Float64, DefaultLS, DDEBifurcationKit.HopfDDEEig{DDE_DefaultEig{Float64, typeof(real), Nothing}}, typeof(LinearAlgebra.norm), DDEBifurcationKit.var"#updateMinAugHopf#129"{DDEBifurcationKit.var"#updateMinAugHopf#123#130"{Int64, @Kwargs{bothside::Bool}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, Float64}}, typeof(BifurcationKit.cb_default), PairOfEvents{ContinuousEvent{ComposedFunction{typeof(BifurcationKit.convert_to_tuple_eve), DDEBifurcationKit.var"#testBT_GH#132"{typeof(LinearAlgebra.norm), @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, NewtonPar{Float64, DefaultLS, DDE_DefaultEig{Float64, typeof(real), Nothing}}}}, Tuple{String, String}, Float64, typeof(BifurcationKit.default_finalise_event!), Nothing}, BifurcationKit.BifEvent{typeof(BifurcationKit.detect_bifurcation_event)}}}, state::ContState{BorderedArray{RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, Float64}, Float64, Vector{ComplexF64}, Matrix{ComplexF64}, Tuple{Tuple{Float64, Float64, Int64}, Tuple{Float64, Float64, Int64}}})
    @ BifurcationKit ~/.julia/packages/BifurcationKit/Y3DPO/src/Continuation.jl:317
 [10] iterate_from_two_points(it::ContIterable{BifurcationKit.HopfCont, BifurcationKit.HopfMAProblem{DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, BifurcationKit.AutoDiff, RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:a21}, typeof(BifurcationKit.plot_default), DDEBifurcationKit.var"#125#133"{DDEBifurcationKit.var"#125#126#134"{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Tuple{Symbol, Symbol}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}}}}, PALC{Secant, BorderingBLS{DefaultLS, Float64, typeof(VectorInterface.inner), typeof(LinearAlgebra.norm)}, Float64, BifurcationKit.DotTheta{BifurcationKit.var"#344#346", BifurcationKit.var"#345#347"}}, Float64, DefaultLS, DDEBifurcationKit.HopfDDEEig{DDE_DefaultEig{Float64, typeof(real), Nothing}}, typeof(LinearAlgebra.norm), DDEBifurcationKit.var"#updateMinAugHopf#129"{DDEBifurcationKit.var"#updateMinAugHopf#123#130"{Int64, @Kwargs{bothside::Bool}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, Float64}}, typeof(BifurcationKit.cb_default), PairOfEvents{ContinuousEvent{ComposedFunction{typeof(BifurcationKit.convert_to_tuple_eve), DDEBifurcationKit.var"#testBT_GH#132"{typeof(LinearAlgebra.norm), @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, PropertyLens{:a21}, DDEBifurcationKit.HopfDDEProblem{ConstantDDEBifProblem{BifFunction{typeof(neuronVF), Nothing, DDEBifurcationKit.var"#4#18"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}, Nothing, Nothing, Nothing, Nothing, DDEBifurcationKit.var"#7#22"{DDEBifurcationKit.var"#d1Fad#20"{DDEBifurcationKit.var"#3#17"{typeof(neuronVF)}}}, DDEBifurcationKit.var"#9#24", DDEBifurcationKit.var"#11#26", DDEBifurcationKit.var"#13#28", Float64, BifurcationKit.Jet{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}, typeof(delaysF), Vector{Float64}, Vector{Float64}, @NamedTuple{κ::Float64, β::Int64, a12::Int64, a21::Float64, τ1::Float64, τ2::Float64, τs::Float64}, PropertyLens{:τs}, typeof(BifurcationKit.plot_default), typeof(BifurcationKit.record_sol_default), typeof(BifurcationKit.save_solution_default), Float64}, Vector{ComplexF64}, Float64, DefaultLS, DefaultLS, MatrixBLS{DefaultLS}, MatrixBLS{DefaultLS}, LinearAlgebra.UniformScaling{Bool}}, NewtonPar{Float64, DefaultLS, DDE_DefaultEig{Float64, typeof(real), Nothing}}}}, Tuple{String, String}, Float64, typeof(BifurcationKit.default_finalise_event!), Nothing}, BifurcationKit.BifEvent{typeof(BifurcationKit.detect_bifurcation_event)}}}, u₀::RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, p₀::Float64, u₁::RecursiveArrayTools.ArrayPartition{Float64, NTuple{4, Vector{Float64}}}, p₁::Float64; _verbosity::UInt8)
...
    @ ~/.julia/dev/DDEBifurcationKit/src/codim2/Hopf.jl:278 [inlined]
 [20] #continuation#422
    @ ~/.julia/packages/BifurcationKit/Y3DPO/src/codim2/codim2.jl:308 [inlined]
 [21] top-level scope
    @ ~/1_work/GIT_repositories/WheelDynamics_Delay_PhD/notebooks/jl_notebook_cell_df34fa98e69747e1a8f8a730347b8e2f_W0sZmlsZQ==.jl:36
Output is truncated. View as a [scrollable element](command:cellOutput.enableScrolling?70f256ce-7c44-4d6b-83c6-0b08b4b90af0) or open in a [text editor](command:workbench.action.openLargeOutput?70f256ce-7c44-4d6b-83c6-0b08b4b90af0). Adjust cell output [settings](command:workbench.action.openSettings?%5B%22%40tag%3AnotebookOutputLayout%22%5D)...

[rveltz]

Let me know if there is anything else I can help you with.

Can you try it on julia 1.10 please, I am quite busy today. I think it needs import LinearAlgebra as LA in DDEBifurcationKit.jl and make sure that the norms are called LA.norm in the package. This stems from BK now importing packages instead of using.

[rveltz]

seems good on master. Re-open if problematic.