Add a Makie extension.

Author: kellertuer

Project.toml

@@ -1,7 +1,7 @@
 name = "ManifoldsBase"
 uuid = "3362f125-f0bb-47a3-aa74-596ffd7ef2fb"
+version = "2.4.0"
 authors = ["Seth Axen <seth.axen@gmail.com>", "Mateusz Baran <mateuszbaran89@gmail.com>", "Ronny Bergmann <manopt@ronnybergmann.net>", "Antoine Levitt <antoine.levitt@gmail.com>"]
-version = "2.3.5"
 
 [workspace]
 projects = ["test", "docs", "tutorials"]
@@ -13,19 +13,22 @@ Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 
 [weakdeps]
+Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Quaternions = "94ee1d12-ae83-5a48-8b1c-48b8ff168ae0"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [extensions]
+ManifoldsBaseMakieExt = "Makie"
 ManifoldsBasePlotsExt = "Plots"
 ManifoldsBaseQuaternionsExt = "Quaternions"
 ManifoldsBaseRecursiveArrayToolsExt = "RecursiveArrayTools"
 ManifoldsBaseStatisticsExt = "Statistics"
 
 [compat]
 LinearAlgebra = "1.6"
+Makie = "0.24"
 Markdown = "1.6"
 Plots = "1"
 Printf = "1.6"

ext/ManifoldsBaseMakieExt.jl

@@ -0,0 +1,61 @@
+module ManifoldsBaseMakieExt
+
+using ManifoldsBase
+using Makie
+using Printf: @sprintf
+
+function ManifoldsBase.plot_slope(
+        x, y; slope = 2, line_base = 0, a = 0, b = 2.0, i = 1, j = length(x),
+    )
+    fig = Makie.Figure()
+    # Setup the log log plot
+    ax = Makie.Axis(
+        fig[1, 1];
+        xscale = log10, yscale = log10, xlabel = "t", ylabel = L"E(t)", title = "Slope plot",
+    )
+    # Main error data
+    Makie.lines!(ax, x, y; label = L"E(t)", linewidth = 3, color = :lightblue)
+    # Reference slope
+    s_line = [exp10(line_base + t * slope) for t in log10.(x)]
+    Makie.lines!(
+        ax, x, s_line; label = "slope s=$slope",
+        linestyle = :dash, color = :black, linewidth = 2
+    )
+
+    # Best-fit slope segment if i and j are provided
+    if (i != 0) && (j != 0)
+        best_line = [exp10(a + t * b) for t in log10.(x[i:j])]
+        Makie.lines!(
+            ax, x[i:j], best_line;
+            label = "best slope $(@sprintf("%.4f", b))",
+            color = :blue, linestyle = :dot, linewidth = 2,
+        )
+    end
+    Makie.axislegend(ax; position = :lt)
+    return fig
+end
+function ManifoldsBase.plot_check_geodesic(T, N, e_norm, e_pt, e_alpha)
+    fig = Figure()
+    ax = Axis(
+        fig[1, 1];
+        xlabel = "t", ylabel = "error magnitude", title = "Geodesic checks",
+    )
+    # Speed deviation: |‖Xᵢ‖ − mean‖X‖|
+    Makie.lines!(
+        ax, T[1:(N - 1)], e_norm;
+        label = L"|\|X_i\| - \overline{\|X\|}|", linewidth = 3, color = :lightblue
+    )
+    # Parallel-transport mismatch (norm): ‖Xᵢ − PT_{i+1→i}(X_{i+1})‖
+    lines!(
+        ax, T[1:(N - 2)], e_pt;
+        label = "PT difference (norm)", linewidth = 3, color = :blue
+    )
+    # Parallel-transport mismatch (angle): |1 − cos∠(Xᵢ, PT_{i+1→i}X_{i+1})|
+    lines!(
+        ax, T[1:(N - 2)], e_alpha;
+        label = "PT difference (angle)", linewidth = 3, color = :green
+    )
+    axislegend(ax; position = :rt)
+    return fig
+end
+end

ext/ManifoldsBasePlotsExt.jl

@@ -3,48 +3,26 @@ module ManifoldsBasePlotsExt
 using ManifoldsBase
 using Plots
 using Printf: @sprintf
-import ManifoldsBase: plot_slope
 
 function ManifoldsBase.plot_slope(
-        x,
-        y;
-        slope = 2,
-        line_base = 0,
-        a = 0,
-        b = 2.0,
-        i = 1,
-        j = length(x),
+        x, y;
+        slope = 2, line_base = 0, a = 0, b = 2.0, i = 1, j = length(x),
     )
     fig = plot(
-        x,
-        y;
-        xaxis = :log,
-        yaxis = :log,
-        label = "\$E(t)\$",
-        linewidth = 3,
-        legend = :topleft,
-        color = :lightblue,
+        x, y; xaxis = :log, yaxis = :log, label = "\$E(t)\$",
+        linewidth = 3, legend = :topleft, color = :lightblue,
     )
     s_line = [exp10(line_base + t * slope) for t in log10.(x)]
     plot!(
-        fig,
-        x,
-        s_line;
-        label = "slope s=$slope",
-        linestyle = :dash,
-        color = :black,
-        linewidth = 2,
+        fig, x, s_line; label = "slope s=$slope",
+        linestyle = :dash, color = :black, linewidth = 2,
     )
     if (i != 0) && (j != 0)
         best_line = [exp10(a + t * b) for t in log10.(x[i:j])]
         plot!(
-            fig,
-            x[i:j],
-            best_line;
+            fig, x[i:j], best_line;
             label = "best slope $(@sprintf("%.4f", b))",
-            color = :blue,
-            linestyle = :dot,
-            linewidth = 2,
+            color = :blue, linestyle = :dot, linewidth = 2,
         )
     end
     return fig