Now 112, 600 allocations

Author: ufechner7

src/solver.jl

@@ -16,6 +16,9 @@ Struct for storing the solution of the [solve!](@ref) function. Must contain all
 """
 mutable struct VSMSolution{P}
     panel_width_array::Vector{Float64}
+    cl_array::Vector{Float64}
+    cd_array::Vector{Float64}
+    cm_array::Vector{Float64}
     gamma_distribution::Union{Nothing, Vector{Float64}}
     aero_force::MVec3          
     aero_moments::MVec3       
@@ -27,7 +30,8 @@ mutable struct VSMSolution{P}
 end
 
 function VSMSolution(P)
-    VSMSolution{P}(zeros(P), nothing, zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(3), zeros(3), FAILURE)
+    VSMSolution{P}(zeros(P), zeros(P), zeros(P), zeros(P), nothing, 
+                   zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(MVec3), zeros(3), zeros(3), FAILURE)
 end
 
 """
@@ -118,9 +122,9 @@ function solve!(solver::Solver, body_aero::BodyAerodynamics, gamma_distribution=
 
     # Initialize arrays
     n_panels = length(panels)
-    cl_array = zeros(n_panels)
-    cd_array = zeros(n_panels)
-    cm_array = zeros(n_panels)
+    cl_array = solver.sol.cl_array
+    cd_array = solver.sol.cd_array
+    cm_array = solver.sol.cm_array
     panel_width_array = solver.sol.panel_width_array
     solver.sol.moment_distribution = zeros(n_panels)
     solver.sol.moment_coefficient_distribution = zeros(n_panels)