Fix plotting 3d polars

Author: 1-Bart-1

bin/run_julia

@@ -19,7 +19,7 @@ done
 if [[ $# -gt 0 ]]; then
     JULIA_ARGS=("$@")
 else
-    JULIA_ARGS=(-i -e 'using VortexStepMethod; function menu(); include("examples/menu.jl"); end; function menu_cp(); include("examples_cp/menu_cp.jl"); end')
+    JULIA_ARGS=()
 fi
 
 julia --project "${JULIA_ARGS[@]}"

ext/VortexStepMethodControlPlotsExt.jl

@@ -793,8 +793,9 @@ function VortexStepMethod.plot_polar_data(body_aero::BodyAerodynamics,
             ax = fig.add_subplot(1, 3, idx, projection="3d")
 
             # Create interpolation matrix
-            interp_matrix = zeros(length(alphas), length(delta_tes))
-            interp_matrix .= [interp(alpha, delta_te) for alpha in alphas, delta_te in delta_tes]
+            interp_matrix = zeros(length(delta_tes), length(alphas))
+            interp_matrix .= [interp(alpha, delta_te)
+                              for delta_te in delta_tes, alpha in alphas]
             X = collect(delta_tes) .+ zeros(length(alphas))'
             Y = collect(alphas)' .+ zeros(length(delta_tes))
 

ext/VortexStepMethodMakieExt.jl

@@ -920,7 +920,7 @@ function VortexStepMethod.plot_polar_data(body_aero::BodyAerodynamics,
 
             # Create interpolation matrix
             interp_matrix = [interp(alpha, delta_te)
-                             for alpha in alphas, delta_te in delta_tes]
+                             for delta_te in delta_tes, alpha in alphas]
 
             # Create wireframe
             wireframe!(ax, delta_tes, alphas, interp_matrix;

test/plotting/test_plotting.jl

@@ -245,6 +245,16 @@ end
         @test fig !== nothing
     end
 
+    fig_rect = plot_polar_data(body_aero;
+        alphas=collect(deg2rad.(-5:1.0:15)),
+        delta_tes=collect(deg2rad.(-3:1.0:5)),
+        is_show=false)
+    if backend == "Makie"
+        @test fig_rect isa Figure
+    else
+        @test fig_rect !== nothing
+    end
+
     # Tests for both backends
     body_aero_empty = create_body_aero()
     empty!(body_aero_empty.panels)