improve examples

Author: rveltz

examples/ellipsoid.jl

@@ -12,15 +12,13 @@ using MultiParamContinuation
 using Test, LinearAlgebra
 const MPC = MultiParamContinuation
 
-F(u,p) = [u[1]^2 + u[2]^2 + p.R * u[3]^2 - 1] # sphere
+F(u,p) = [u[1]^2 + u[2]^2 + p.R * u[3]^2 - 1]
 
 prob = ManifoldProblem(F, 
                     [0. ,0, 1],
                     (R = 2.,);
-                    # get_tangent
                         )
 
-# problem quand j'ajoute une charte
 S = MPC.continuation(prob,
             Henderson(np0 = 4,
                         θmin = 0.001,

examples/plane2.jl

@@ -0,0 +1,37 @@
+using Revise
+cd(@__DIR__)
+using Pkg
+pkg"activate ."
+using MultiParamContinuation
+
+using GLMakie
+Makie.inline!(false)
+Makie.inline!(true)
+
+
+using Test, LinearAlgebra
+const MPC = MultiParamContinuation
+
+F(u,p) = [u[1] + u[2] - u[3]]
+
+prob = ManifoldProblem(F, zeros(3), nothing)
+
+S = continuation(prob,
+            Henderson(np0 = 4),
+            CoveringPar(max_charts = 20000, 
+                    max_steps = 100,
+                    R0 = .1,
+                    ϵ = Inf,
+                    delta_angle = Inf,
+                    )
+                )
+
+MPC.plotd(S; 
+    plot_circle = true, 
+    # draw_tangent = false, 
+    draw_edges = true,
+    plot_center = true,
+    put_ids = true,
+    ind_plot = [1,2])
+
+

examples/sphere_bifurcationKit.jl

@@ -19,7 +19,6 @@ F(u,p) = [u[1]^2 + u[2]^2 + u[3]^2 - 1] # sphere
 prob = ManifoldProblem_BK(F, 
                     [1.1,0.,0.],
                     nothing;
-                    # get_tangent
                         )
 
 # problem quand j'ajoute une charte

examples/sphere_event.jl

@@ -11,7 +11,7 @@ using MultiParamContinuation, LinearAlgebra
 
 const MPC = MultiParamContinuation
 
-F(u,p) = [u[1]^2 + u[2]^2 + u[3]^2 - 1] # sphere
+F(u,p) = [u[1]^2 + u[2]^2 + u[3]^2 - 1]
 
 function get_tangent(u, par)
     # compute the normal

examples/transcritical.jl

@@ -0,0 +1,51 @@
+using Revise
+cd(@__DIR__)
+using Pkg
+pkg"activate ."
+
+using GLMakie
+Makie.inline!(false)
+Makie.inline!(true)
+
+using MultiParamContinuation
+
+const MPC = MultiParamContinuation
+
+function F(u,p) 
+    x,y,z = u
+    [(z-x-.5*(y-1)*(y-1))*(x+y+z) + y/10]
+end
+
+prob = ManifoldProblem(F, [-0.6,-0.6,-0.6], nothing;
+            finalize_solution = Cube(1.1))
+
+S = continuation(prob,
+            Henderson(np0 = 4, 
+                      θmin = 0.1,
+                    #   use_curvature = true,
+                      ),
+            CoveringPar(max_charts = 20000,
+                    max_steps = 2000,
+                    verbose = 0,
+                    newton_options = NonLinearSolveSpec(;maxiters = 6, abstol = 1e-12, reltol = 1e-10),
+                    R0 = .1,
+                    ϵ = 0.005,
+                    delta_angle = 10.15,
+                    ))
+
+MPC.plotd(S; 
+    # draw_circle = true,
+    draw_tangent = true, 
+    # plot_center = true,
+    # put_ids = true,
+    ind_plot = [1,3])
+
+step!(S, 1000);fig = MPC.plotd(S; draw_circle = false)
+
+
+MPC.plot2d(S; 
+    # draw_circle = true, 
+    draw_tangent = true, 
+    plot_center = true,
+    # put_ids = true,
+    ind_plot = [1,2])