improve show method for DCResult

correct prediction in show method for NS normal form
remove figure size in Makie recipe

src/DeflatedContinuation.jl

@@ -75,7 +75,8 @@ Base.length(br::DCResult) = length(br.branches)
 
 function Base.show(io::IO, brdc::DCResult; comment = "", prefix = " ")
     printstyled(io, "Deflated continuation result, # branches = $(length(brdc.branches))", "\n", color=:cyan, bold = true)
-    for br in brdc.branches
+    for (ii, br) in pairs(brdc.branches)
+        printstyled("\nBranch #", ii, ":\n", color = :cyan, bold = true)
         show(io, br)
     end
 end

src/periodicorbit/BifurcationPoints.jl

@@ -129,7 +129,7 @@ function Base.show(io::IO, ns::NeimarkSackerPO)
     println(io, " bifurcation point of periodic orbit\n┌─ ", get_lens_symbol(ns.nf.lens)," ≈ $(ns.p).")
     println(io, "├─ Frequency θ ≈ ", ns.ω)
     println(io, "├─ Period at the periodic orbit T ≈ ", abs(ns.T))
-    println(io, "├─ Second period of the bifurcated torus ≈ ", abs(2pi/ns.ω))
+    println(io, "├─ Second period of the bifurcated torus ≈ ", abs(2pi*ns.ω*ns.T))
     if ns.prm
         println(io, "├─ Normal form z -> z⋅eⁱᶿ(1 + a⋅δp + b⋅|z|²)")
     else

src/periodicorbit/PeriodicOrbitCollocation.jl

@@ -364,16 +364,16 @@ function generate_ci_problem(pb::PeriodicOrbitOCollProblem,
     N = length(u0)
 
     n, m, Ntst = size(pb)
-    nunknows = N * (1 + m * Ntst)
+    n_unknows = N * (1 + m * Ntst)
 
     par = sol.prob.p
     prob_vf = re_make(bifprob, params = par)
 
     pbcoll = setproperties(pb,
                             N = N,
                             prob_vf = prob_vf,
-                            ϕ = zeros(nunknows),
-                            xπ = zeros(nunknows),
+                            ϕ = zeros(n_unknows),
+                            xπ = zeros(n_unknows),
                             cache = POCollCache(eltype(pb), N, m))
 
     # find best period candidate

src/periodicorbit/cop.jl

@@ -302,9 +302,9 @@ end
         _ipiv2perm!(p, F.ipiv, size(F, 1))
         _invperm!(pinv, p)
 
-        @assert p[nbcoll+1] == nbcoll+1 "Pivoting strategy failed!!"
+        @assert p[nbcoll+1] == nbcoll+1 "Pivoting strategy failed!! Please report to the website of BifurcationKit. You may try the default linear solver `defaultLS` as a backup."
         if dim > 0
-            @assert p[nbcoll+2] == nbcoll+2 "Pivoting strategy failed!!"
+            @assert p[nbcoll+2] == nbcoll+2 "Pivoting strategy failed!! Please report to the website of BifurcationKit. You may try the default linear solver `defaultLS` as a backup."
         end
 
         # Lₜ = LowerTriangular(F.L) # zero allocation?

src/plotting/RecipesMakie.jl

@@ -140,7 +140,7 @@ function plot(contres::AbstractBranchResult; kP...)
     ind1, ind2 = get_plot_vars(contres, nothing)
     xlab, ylab = get_axis_labels(ind1, ind2, contres)
 
-    fig = Figure(size = (1200, 700))
+    fig = Figure()
     ax1 = fig[1, 1] = Axis(fig, xlabel = String(xlab), ylabel = String(ylab), tellheight = true)
 
     plot!(ax1, contres; kP...)
@@ -153,7 +153,7 @@ function plot(brs::AbstractBranchResult...;
         branchlabel = fill("", length(brs)),
         kP...)
     if length(brs) == 0; return ;end
-    fig = Figure(size = (1200, 700))
+    fig = Figure()
     ax1 = fig[1, 1] = Axis(fig)
 
     for (id, contres) in pairs(brs)