correct mistake in aBS

bifurcationkit · Nov 16, 2024 · 705a236 · 705a236
1 parent c9b7e97
commit 705a236
Show file tree

Hide file tree

Showing 3 changed files with 15 additions and 12 deletions.
diff --git a/src/bifdiagram/BranchSwitching.jl b/src/bifdiagram/BranchSwitching.jl
@@ -219,7 +219,7 @@ function multicontinuation(br::AbstractBranchResult,
                         bpnf::NdBranchPoint,
                         options_cont::ContinuationPar = br.contparams;
                         δp = nothing,
-                        ampfactor = getvectoreltype(br)(1),
+                        ampfactor = _getvectoreltype(br)(1),
                         perturb = identity,
                         plot_solution = plot_solution(br.prob),
                         kwargs...)

diff --git a/src/codim2/codim2.jl b/src/codim2/codim2.jl
@@ -234,17 +234,17 @@ function continuation(br::AbstractBranchResult,
     _options_cont = detect_codim2_parameters(detect_codim2_bifurcation, options_cont; update_minaug_every_step, kwargs...)
 
     if br.specialpoint[ind_bif].type == :hopf
-    return continuation_hopf(prob, br, ind_bif, lens2, _options_cont;
-        start_with_eigen,
-        compute_eigen_elements,
-        update_minaug_every_step,
-        kwargs...)
+        return continuation_hopf(prob, br, ind_bif, lens2, _options_cont;
+            start_with_eigen,
+            compute_eigen_elements,
+            update_minaug_every_step,
+            kwargs...)
     else
-    return continuation_fold(prob, br, ind_bif, lens2, _options_cont;
-        start_with_eigen,
-        compute_eigen_elements,
-        update_minaug_every_step,
-        kwargs...)
+        return continuation_fold(prob, br, ind_bif, lens2, _options_cont;
+            start_with_eigen,
+            compute_eigen_elements,
+            update_minaug_every_step,
+            kwargs...)
     end
 end
 ####################################################################################################

diff --git a/src/periodicorbit/Floquet.jl b/src/periodicorbit/Floquet.jl
@@ -450,9 +450,12 @@ end
 
 @views function (eig::FloquetColl)(JacColl, nev; kwargs...)
     pbcoll = JacColl.pb
-    𝒯 = eltype(pbcoll)
     J = _get_matrix(JacColl)
     n, m, Ntst = size(pbcoll)
+    _eig_floquet_col(J, n, m, Ntst, nev)
+end
+
+function _eig_floquet_col(J::AbstractMatrix{𝒯}, n, m, Ntst, nev) where 𝒯
     nbcoll = n * m
     N = n
     In = LinearAlgebra.I(N)