Refactoring (make private fields private)

Author: ufechner7

src/solver.jl

@@ -25,7 +25,7 @@ Struct for storing the solution of the [solve!](@ref) function. Must contain all
 """
 @with_kw mutable struct VSMSolution{P}
     ### private vectors of solve_base!
-    x_airf_array::Matrix{Float64} = zeros(P, 3)
+    _x_airf_array::Matrix{Float64} = zeros(P, 3)
     y_airf_array::Matrix{Float64} = zeros(P, 3)
     z_airf_array::Matrix{Float64} = zeros(P, 3)
     va_array::Matrix{Float64} = zeros(P, 3)
@@ -211,7 +211,7 @@ function solve!(solver::Solver, body_aero::BodyAerodynamics, gamma_distribution=
             gamma_new,
             solver.core_radius_fraction,
             solver.sol.z_airf_array,
-            solver.sol.x_airf_array,
+            solver.sol._x_airf_array,
             solver.sol.va_array,
             solver.br.va_norm_array,
             solver.br.va_unit_array
@@ -338,7 +338,7 @@ function solve(solver::Solver, body_aero::BodyAerodynamics, gamma_distribution=n
         solver.lr.alpha_array,
         solver.lr.v_a_array,
         solver.sol.chord_array,
-        solver.sol.x_airf_array,
+        solver.sol._x_airf_array,
         solver.sol.y_airf_array,
         solver.sol.z_airf_array,
         solver.sol.va_array,
@@ -368,15 +368,15 @@ function solve_base!(solver::Solver, body_aero::BodyAerodynamics, gamma_distribu
     relaxation_factor = solver.relaxation_factor
 
     # Clear arrays
-    solver.sol.x_airf_array .= 0
+    solver.sol._x_airf_array .= 0
     solver.sol.y_airf_array .= 0
     solver.sol.z_airf_array .= 0
     solver.sol.va_array .= 0
     solver.sol.chord_array .= 0
 
     # Fill arrays from panels
     for (i, panel) in enumerate(panels)
-        solver.sol.x_airf_array[i, :] .= panel.x_airf
+        solver.sol._x_airf_array[i, :] .= panel.x_airf
         solver.sol.y_airf_array[i, :] .= panel.y_airf
         solver.sol.z_airf_array[i, :] .= panel.z_airf
         solver.sol.va_array[i, :] .= panel.va
@@ -435,7 +435,7 @@ function gamma_loop!(
 )
     va_array = solver.sol.va_array
     chord_array = solver.sol.chord_array
-    x_airf_array = solver.sol.x_airf_array
+    x_airf_array = solver.sol._x_airf_array
     y_airf_array = solver.sol.y_airf_array
     z_airf_array = solver.sol.z_airf_array
     solver.lr.converged   = false

test/bench.jl

@@ -130,7 +130,7 @@ using LinearAlgebra
                 )
                 solver.sol.va_array .= va_array
                 solver.sol.chord_array .= chord_array
-                solver.sol.x_airf_array .= x_airf_array
+                solver.sol._x_airf_array .= x_airf_array
                 solver.sol.y_airf_array .= y_airf_array
                 solver.sol.z_airf_array .= z_airf_array
                 result = @benchmark gamma_loop!(