Lipo deform

Author: 1-Bart-1

src/obj_geometry.jl

@@ -309,7 +309,10 @@ function refine_obj_wing!(wing::AbstractWing{T}; recompute_mapping=true) where {
     end
 
     # 5. Compute panel mapping and update non_deformed_sections
-    recompute_mapping && VortexStepMethod.compute_refined_panel_mapping!(wing)
+    if recompute_mapping
+        VortexStepMethod.compute_refined_panel_mapping!(wing)
+        VortexStepMethod.compute_refined_section_interpolation!(wing)
+    end
     VortexStepMethod.update_non_deformed_sections!(wing)
 
     return nothing
@@ -584,6 +587,7 @@ function ObjWing(
             MVector{3, Float64}(spanwise_direction),
             sections, Section{Float64}[], remove_nan, use_prior_polar, 0.0,  # billowing_percentage
             Int16[],  # refined_panel_mapping empty
+            Int16[], Float64[],  # refined_section interpolation cache empty
             Section{Float64}[], zeros(Float64, n_panels), zeros(Float64, n_panels),
             Float64(mass), Float64(gamma_tip), Matrix{Float64}(inertia_tensor),
             MVector{3, Float64}(T_cad_body), MMatrix{3, 3, Float64, 9}(R_cad_body),

src/solver.jl

@@ -1015,6 +1015,8 @@ function _wing_with_eltype(wing::Wing{P, Float64}, ::Type{TD}) where {P, TD}
         wing.use_prior_polar,
         wing.billowing_percentage,
         copy(wing.refined_panel_mapping),
+        copy(wing.refined_section_left_idx),
+        Vector{TD}(wing.refined_section_weight),
         Section{TD}[_section_with_eltype(s, TD) for s in wing.non_deformed_sections],
         Vector{TD}(wing.theta_dist),
         Vector{TD}(wing.delta_dist),