Add P to VSMSolution and bench2.jl

Author: ufechner7

src/solver.jl

@@ -14,7 +14,7 @@ Struct for storing the solution of the [solve!](@ref) function. Must contain all
 - moment_coefficient_distribution::Vector{Float64}: Pitching moment coefficient around the spanwise vector of each panel. [-]
 - solver_status::SolverStatus: enum, see [SolverStatus](@ref)
 """
-mutable struct VSMSolution
+mutable struct VSMSolution{P}
     gamma_distribution::Union{Nothing, Vector{Float64}}
     aero_force::MVec3          
     aero_moments::MVec3       
@@ -25,8 +25,8 @@ mutable struct VSMSolution
     solver_status::SolverStatus 
 end
 
-function VSMSolution()
-    VSMSolution(nothing, zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(3), zeros(3), FAILURE)
+function VSMSolution(P)
+    VSMSolution{P}(nothing, zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(3), zeros(3), FAILURE)
 end
 
 """
@@ -57,7 +57,7 @@ Main solver structure for the Vortex Step Method.See also: [solve](@ref)
 ## Solution
 sol::VSMSolution = VSMSolution(): The result of calling [solve!](@ref) 
 """
-@with_kw struct Solver
+@with_kw struct Solver{P}
     # General settings
     aerodynamic_model_type::Model = VSM
     density::Float64 = 1.225
@@ -77,7 +77,7 @@ sol::VSMSolution = VSMSolution(): The result of calling [solve!](@ref)
     is_only_f_and_gamma_output::Bool = false
 
     # Solution
-    sol::VSMSolution = VSMSolution()
+    sol::VSMSolution{P} = VSMSolution(P)
 end
 
 """

test/bench2.jl

@@ -40,8 +40,8 @@ vel_app = [cos(alpha), 0.0, sin(alpha)] .* v_a
 set_va!(wa, vel_app)
 
 # Step 4: Initialize solvers for both LLT and VSM methods
-llt_solver = Solver(aerodynamic_model_type=LLT)
-vsm_solver = Solver(aerodynamic_model_type=VSM)
+P = length(wa.panels)
+vsm_solver = Solver{P}(aerodynamic_model_type=VSM)
 
 # Step 5: Solve using both methods
 println("Rectangular wing, solve_base:")
@@ -50,5 +50,5 @@ println("Rectangular wing, solve_base:")
 # time Julia:   0.6 ms Ryzen 7950x
 #               0.47 ms laptop, performance mode, grid
 println("Rectangular wing, solve!:")
-@btime sol = solve!($vsm_solver, $wa) # 608 allocations
+@btime sol = solve!($vsm_solver, $wa)                  # 608 allocations
 nothing