correct steinmetz

bifurcationkit · Jun 15, 2024 · 5c21adc · 5c21adc
1 parent 540dc0a
commit 5c21adc
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 8 deletions.
diff --git a/docs/src/tutorials/ode/steinmetz.md b/docs/src/tutorials/ode/steinmetz.md
@@ -141,6 +141,7 @@ opts_po_cont = ContinuationPar(p_min = 0., p_max = 2.0,
 	nev = 4,
 	max_steps = 50, 
 	tol_stability = 1e-5)
+
 br_coll = continuation(probcoll, cicoll, PALC(tangent = Bordered()), opts_po_cont;
     # verbosity = 3, plot = true,
     callback_newton = BK.cbMaxNorm(10),
@@ -153,7 +154,7 @@ br_coll = continuation(probcoll, cicoll, PALC(tangent = Bordered()), opts_po_con
 opts_pocl_fold = ContinuationPar(br_coll.contparams, detect_bifurcation = 1, plot_every_step = 10, dsmax = 4e-2)
 fold_po_cl = @time continuation(br_coll, 2, (@lens _.k7), opts_pocl_fold;
         # verbosity = 3, plot = true,
-        detect_codim2_bifurcation = 2,
+        detect_codim2_bifurcation = 1,
         update_minaug_every_step = 1,
         start_with_eigen = false,
         usehessian = true,

diff --git a/docs/src/tutorials/ode/tutorialsCodim2PO.md b/docs/src/tutorials/ode/tutorialsCodim2PO.md
@@ -201,25 +201,23 @@ probcoll, ci = generate_ci_problem(PeriodicOrbitOCollProblem(30, 3), re_make(pro
 
 prob2 = @set probcoll.prob_vf.lens = @lens _.ϵ
 brpo_pd = continuation(prob2, ci, PALC(), ContinuationPar(opts_po_cont, dsmax = 5e-3);
-	# verbosity = 3, plot = true,
 	argspo...,
 	bothside = true,
 	)
 
 opts_pocoll_pd = ContinuationPar(brpo_pd.contparams, max_steps = 40, p_min = 1.e-2, dsmax = 1e-2, ds = -1e-3)
-@set! opts_pocoll_pd.newton_options.tol = 1e-12
-pd_po_coll2 = continuation(brpo_pd, 2, (@lens _.b0), opts_pocoll_pd;
-		# verbosity = 3,
-		detect_codim2_bifurcation = 2,
+@reset opts_pocoll_pd.newton_options.tol = 1e-12
+pd_po_coll2 = continuation(brpo_pd, 2, (@optic _.b0), opts_pocoll_pd;
+		detect_codim2_bifurcation = 1,
 		start_with_eigen = false,
 		jacobian_ma = :minaug,
 		normC = norminf,
 		callback_newton = BK.cbMaxNorm(10),
 		bothside = true,
 		)
 
-plot(fold_po_coll1, fold_po_coll2, branchlabel = ["FOLD", "FOLD"])
-plot!(pd_po_coll2, vars = (:ϵ, :b0), branchlabel = "PD")
+plot(pd_po_coll2, vars = (:ϵ, :b0), branchlabel = "PD")
+plot!(fold_po_coll1, fold_po_coll2, branchlabel = ["FOLD", "FOLD"])
 ```
 
 ## References