correction constructor of BifProblem

src/BifurcationPoints.jl

@@ -66,7 +66,7 @@ $(TYPEDFIELDS)
     "Continuation step at which the special occurs"
     step::Int64 = 0
 
-    "`status ∈ {:converged, :guess, :guessL}` indicates whether the bisection algorithm was successful in detecting the special (bifurcation) point. If `status == :guess`, the bissection algorithm failed to meet the requirements given in `::ContinuationPar`. Same for `status == :guessL` but the bissection algorithm stopped on the left of the bifurcation point."
+    "`status ∈ {:converged, :guess, :guessL}` indicates whether the bisection algorithm was successful in detecting the special (bifurcation) point. If `status == :guess`, the bisection algorithm failed to meet the requirements given in `::ContinuationPar`. Same for `status == :guessL` but the bisection algorithm stopped on the left of the bifurcation point."
     status::Symbol = :guess
 
     "`δ = (δr, δi)` where δr indicates the change in the number of unstable eigenvalues and δi indicates the change in the number of unstable eigenvalues with nonzero imaginary part. `abs(δr)` is thus an estimate of the dimension of the kernel of the Jacobian at the special (bifurcation) point."

src/Problems.jl

@@ -282,6 +282,9 @@ for (op, at) in (
                     iip = _isinplace(_F)
                     F = iip ? (x, p) -> _F(similar(x), x, p) : _F
                 end
+
+                Finp = nothing
+
                 J = isnothing(J) ? (x, p) -> ForwardDiff.jacobian(z -> F(z, p), x) : J
                 jvp = isnothing(jvp) ?
                       (x, p, dx) -> ForwardDiff.derivative(t -> F(x .+ t .* dx, p), zero(eltype(dx))) : dF

src/periodicorbit/PeriodicOrbitTrapeze.jl

@@ -376,7 +376,7 @@ function jacobian_potrap_block(pb::PeriodicOrbitTrapProblem, u0::AbstractVector,
     cylic_potrap_block!(pb, u0, par, Aγ)
 
     In = spdiagm( 0 => ones(N))
-    Aγ[Block(M, 1)] = -γ * In
+    Aγ[Block(M, 1)] = (-γ) * In
     Aγ[Block(M, M)] = In
     return Aγ
 end
@@ -535,7 +535,7 @@ end
 
         @. Jn = In -  tmpJ * (h/2)
         # setblock!(Jc, Jn, ii, ii)
-        J0.nzval[indx[ii,ii]] .= Jn.nzval
+        J0.nzval[indx[ii, ii]] .= Jn.nzval
     end
 
     # setblock!(Aγ, -γ * In, M, 1)

test/lorenz84.jl

@@ -145,7 +145,6 @@ testEV(hp_codim2_test)
 @reset opts_br.newton_options.verbose = false
 sn_codim2 = nothing
 for _jac in (:autodiff, :minaug, :finiteDifferences, :MinAugMatrixBased)
-    @error "" _jac
     # be careful here, Bordered predictor not good for Fold continuation
     sn_codim2 = @time continuation((@set br.alg.tangent = Secant()), 5, (@optic _.T), ContinuationPar(opts_br, p_max = 3.2, p_min = -0.1, detect_bifurcation = 1, dsmin=1e-5, ds = -0.001, dsmax = 0.015, n_inversion = 10, save_sol_every_step = 1, max_steps = 30, max_bisection_steps = 55) ; verbosity = 0,
         normC = norminf,

test/test_linear.jl

@@ -60,7 +60,7 @@ typeof(z_sim) == typeof(z_sim2)
 # test _axpy_op
 J0 = rand(100, 100)
 dx = rand(size(J0, 1))
-_o = dx = rand(size(J0, 1))
+_o = rand(ComplexF64, size(J0, 1))
 a₀ = rand(ComplexF64)
 a₁ = -1.432
 BK._axpy(J0, 0, a₁)