🟢 Sphere

• 🟢 Sphere

In this tutorial, we show how to cover a sphere as solution of

\[F(u) := \|u\|^2-1 = 0\]

We use this model as a mean to introduce the basics of MultiParamContinuation.jl.

It is easy to encode the manifold as follows

using CairoMakie, MultiParamContinuation
const MPC = MultiParamContinuation

F(u,p) = [u[1]^2 + u[2]^2 + u[3]^2 - 1]

prob = ManifoldProblem(F,
                    [1,0.,0.],
                    nothing
                        )

show(prob)

2-d Manifold Problem
    ├─ n = 3
    └─ m = 1

We now compute the covering of the manifold

S = MPC.continuation(prob,
            Henderson(),
            CoveringPar(max_charts = 20000,
                    max_steps = 250,
                    verbose = 0,
                    newton_options = NonLinearSolveSpec(;maxiters = 5),
                    R0 = .2,
                    )
            )
show(S)

Surface 2-d
   ├─ # charts = 234
   └─ problem =
          2-d Manifold Problem
              ├─ n = 3
              └─ m = 1

You plot the result as follows

MPC.plotd(S;
    draw_tangent = true,
    plot_center = false,
    draw_edges = true,
    ind_plot = [1,3])