Fixed deform

Author: 1-Bart-1

src/wing_geometry.jl

@@ -456,14 +456,37 @@ function deform!(wing::Wing)
     chord = zeros(MVec3)
     normal = zeros(MVec3)
 
-    for i in 1:wing.n_panels
-        section1 = wing.non_deformed_sections[i]
-        section2 = wing.non_deformed_sections[i+1]
-        local_y .= normalize(section1.LE_point - section2.LE_point)
-        chord .= section1.TE_point .- section1.LE_point
+    # Process all sections (n_panels + 1)
+    # Each section gets angle(s) from adjacent panel(s)
+    for i in 1:(wing.n_panels + 1)
+        section = wing.non_deformed_sections[i]
+
+        # Determine the angle for this section
+        if i == 1
+            # First section: use angle from first panel
+            theta = wing.theta_dist[1]
+        elseif i == wing.n_panels + 1
+            # Last section: use angle from last panel
+            theta = wing.theta_dist[wing.n_panels]
+        else
+            # Middle sections: average of adjacent panels
+            theta = 0.5 * (wing.theta_dist[i-1] + wing.theta_dist[i])
+        end
+
+        # Compute local coordinate system
+        if i < wing.n_panels + 1
+            section2 = wing.non_deformed_sections[i+1]
+            local_y .= normalize(section.LE_point - section2.LE_point)
+        else
+            # For last section, use same local_y as previous
+            section_prev = wing.non_deformed_sections[i-1]
+            local_y .= normalize(section_prev.LE_point - section.LE_point)
+        end
+
+        chord .= section.TE_point .- section.LE_point
         normal .= chord × local_y
-        @. wing.refined_sections[i].TE_point = section1.LE_point +
-            cos(wing.theta_dist[i]) * chord - sin(wing.theta_dist[i]) * normal
+        @. wing.refined_sections[i].TE_point = section.LE_point +
+            cos(theta) * chord - sin(theta) * normal
     end
     return nothing
 end

test/ram_geometry/test_kite_geometry.jl

@@ -232,11 +232,6 @@ using Serialization
         original_te_points = [copy(wing.refined_sections[i].TE_point)
             for i in 1:n_sections]
 
-        @show wing.refined_sections[1].LE_point
-        @show wing.refined_sections[1].TE_point
-        @show wing.refined_sections[end].LE_point
-        @show wing.refined_sections[end].TE_point
-
         # Apply group_deform! with non-zero angles (2 groups, each controlling 2 panels)
         theta_angles = [deg2rad(15.0), deg2rad(20.0)]
         delta_angles = [deg2rad(5.0), deg2rad(10.0)]

test/runtests.jl

@@ -4,7 +4,25 @@ using Test, VortexStepMethod
 cd(@__DIR__)  # ensure we're in test/ no matter how tests are launched
 include("test_data_utils.jl")
 
+# Support selective test execution via ]test test_args=["pattern"]
+const test_patterns = isempty(ARGS) ? String[] : ARGS
+
 println("Running tests...")
+if !isempty(test_patterns)
+    println("Filtering tests matching: ", test_patterns)
+end
+
+# Helper to check if a test file matches any pattern
+function should_run_test(test_path::String)
+    isempty(test_patterns) && return true
+    for pattern in test_patterns
+        # Match directory (e.g., "solver") or specific file (e.g., "test_group_coefficients")
+        if occursin(pattern, test_path)
+            return true
+        end
+    end
+    return false
+end
 
 # keep your env check as-is...
 const build_is_production_build_env_name = "BUILD_IS_PRODUCTION_BUILD"
@@ -14,25 +32,25 @@ const build_is_production_build = let v = get(ENV, build_is_production_build_env
 end::Bool
 
 @testset verbose = true "Testing VortexStepMethod..." begin
-    if build_is_production_build
+    if build_is_production_build && should_run_test("bench")
         include("bench.jl")
     end
-    include("body_aerodynamics/test_body_aerodynamics.jl")
-    include("body_aerodynamics/test_results.jl")
-    include("filament/test_bound_filament.jl")
-    include("filament/test_semi_infinite_filament.jl")
-    include("panel/test_panel.jl")
-    include("plotting/test_plotting.jl")
-    include("polars/test_polars.jl")
-    include("ram_geometry/test_kite_geometry.jl")
-    include("settings/test_settings.jl")
-    include("solver/test_solver.jl")
-    include("solver/test_group_coefficients.jl")
-    include("VortexStepMethod/test_VortexStepMethod.jl")
-    include("wake/test_wake.jl")
-    include("wing_geometry/test_wing_geometry.jl")
-    include("yaml_geometry/test_yaml_geometry.jl")
-    include("Aqua.jl")
+    should_run_test("body_aerodynamics/test_body_aerodynamics.jl") && include("body_aerodynamics/test_body_aerodynamics.jl")
+    should_run_test("body_aerodynamics/test_results.jl") && include("body_aerodynamics/test_results.jl")
+    should_run_test("filament/test_bound_filament.jl") && include("filament/test_bound_filament.jl")
+    should_run_test("filament/test_semi_infinite_filament.jl") && include("filament/test_semi_infinite_filament.jl")
+    should_run_test("panel/test_panel.jl") && include("panel/test_panel.jl")
+    should_run_test("plotting/test_plotting.jl") && include("plotting/test_plotting.jl")
+    should_run_test("polars/test_polars.jl") && include("polars/test_polars.jl")
+    should_run_test("ram_geometry/test_kite_geometry.jl") && include("ram_geometry/test_kite_geometry.jl")
+    should_run_test("settings/test_settings.jl") && include("settings/test_settings.jl")
+    should_run_test("solver/test_solver.jl") && include("solver/test_solver.jl")
+    should_run_test("solver/test_group_coefficients.jl") && include("solver/test_group_coefficients.jl")
+    should_run_test("VortexStepMethod/test_VortexStepMethod.jl") && include("VortexStepMethod/test_VortexStepMethod.jl")
+    should_run_test("wake/test_wake.jl") && include("wake/test_wake.jl")
+    should_run_test("wing_geometry/test_wing_geometry.jl") && include("wing_geometry/test_wing_geometry.jl")
+    should_run_test("yaml_geometry/test_yaml_geometry.jl") && include("yaml_geometry/test_yaml_geometry.jl")
+    should_run_test("Aqua.jl") && include("Aqua.jl")
 end
 
 nothing