Automatic switching between Natural and PALC

This is one of the continuation methods implemented in the package. It is set by the option AutoSwitch() in continuation. See also AutoSwitch for more information.

This algorithm uses a Natural continuation when the branch is "horizontal" and switches to PALC (or other which cab be specified) otherwise.

Example

We provide an example of use. We define a BifurcationProblem as usual and pass the continuation algorithm AutoSwitch.

using Plotsusing BifurcationKitconst BK = BifurcationKitfunction F(x, p)	(;α) = p	f = similar(x)	f[1] = (-2x[1]+x[2]) + α * exp(x[1])	f[2] = ( x[1]-2x[2]) + α * exp(x[2])	return fendsol0 = zeros(2)par = (α = 0.0, )prob = ODEBifProblem(F, sol0, par, (@optic _.α))

┌─ Bifurcation problem with uType Vector{Float64}
├─ Inplace: false
├─ Dimension: 2
├─ Symmetric: false
└─ Parameter: α

opts = ContinuationPar(dsmax = 0.15, max_steps = 40, )br = continuation(prob, AutoSwitch(), opts, normC = norminf)

 ┌─ Curve type: EquilibriumCont
 ├─ Number of points: 41
 ├─ Type of vectors: Vector{Float64}
 ├─ Parameter α starts at 0.0, ends at 0.006224449009229535
 ├─ Algo: AutoSwitch
 └─ Special points:

┌────┬────────────┬─────────────┬────────────────────────────┬───────────┬─────────────┬──────────┬──────┐
│ Nb │    type    │      α      │          interval          │ precision │   status    │    δ     │ step │
├────┼────────────┼─────────────┼────────────────────────────┼───────────┼─────────────┼──────────┼──────┤
│  1 │         bp │ +0.36787092 │ (+0.36747428, +0.36787092) │     4e-04 │   converged │ ( 1,  0) │   11 │
│  2 │         bp │ +0.14907908 │ (+0.14907908, +0.14973607) │     7e-04 │   converged │ ( 1,  0) │   21 │
│  3 │   endpoint │ +0.00518664 │                            │           │             │          │   41 │
└────┴────────────┴─────────────┴────────────────────────────┴───────────┴─────────────┴──────────┴──────┘

You can plot the result as usual:

plot(br)