Use ==, not === for enums

Author: ufechner7

test/test_body_aerodynamics.jl

@@ -317,7 +317,7 @@ end
                 index_reversed = length(body_aero.panels) - i + 1
                 panel = body_aero.panels[index_reversed]
 
-                evaluation_point = if model === VSM
+                evaluation_point = if model == VSM
                     panel.control_point
                 else  # LLT
                     panel.aero_center
@@ -333,7 +333,7 @@ end
                     atol=1e-4
                 )
 
-                if model === VSM
+                if model == VSM
                     @test isapprox(
                         panel.aero_center,
                         expected_controlpoints[i]["coordinates_aoa"],

test/test_wing_geometry.jl

@@ -15,7 +15,7 @@ Points are compared with approximate equality to handle floating point differenc
 function ==(a::Section, b::Section)
     return (isapprox(a.LE_point, b.LE_point; rtol=1e-5, atol=1e-5) &&
             isapprox(a.TE_point, b.TE_point; rtol=1e-5, atol=1e-5) &&
-            a.aero_model === b.aero_model &&
+            a.aero_model == b.aero_model &&
             all(a.aero_data .== b.aero_data))
 end
 
@@ -329,7 +329,7 @@ end
         @test -2.0 < new_sections[5].LE_point[2] < -1.0
 
         for section in new_sections
-            @test section.aero_model === INVISCID
+            @test section.aero_model == INVISCID
         end
     end
 end

test/thesis_oriol_cayon.jl

@@ -205,7 +205,7 @@ function create_geometry_general(coordinates, Uinf, N, ring_geo, model)
         normal = x_airf
         tangential = y_airf
 
-        if model === VSM
+        if model == VSM
             cp = Dict(
                 "coordinates" => VSMpoint,
                 "chord" => chord,
@@ -215,7 +215,7 @@ function create_geometry_general(coordinates, Uinf, N, ring_geo, model)
                 "coordinates_aoa" => LLpoint
             )
             push!(controlpoints, cp)
-        elseif model === LLT
+        elseif model  LLT
             cp = Dict(
                 "coordinates" => LLpoint,
                 "chord" => chord,
@@ -320,7 +320,7 @@ function thesis_induction_matrix_creation(ringvec, controlpoints, rings, Uinf, G
     airf_coord = [controlpoints[icp]["airf_coord"] for icp in 1:N]
 
     for icp in 1:N
-        if model === VSM
+        if model == VSM
             # Velocity induced by infinite bound vortex with Gamma = 1
             U_2D[icp,:] = velocity_induced_bound_2D(ringvec[icp])
         end
@@ -347,7 +347,7 @@ function thesis_induction_matrix_creation(ringvec, controlpoints, rings, Uinf, G
 
             # Different from old thesis code as this was considered wrong
             if icp == jring
-                if model === VSM
+                if model == VSM
                     MatrixU[icp,jring] -= U_2D[icp,1]
                     MatrixV[icp,jring] -= U_2D[icp,2]
                     MatrixW[icp,jring] -= U_2D[icp,3]
@@ -501,7 +501,7 @@ function velocity_induced_single_ring_semiinfinite(ring, controlpoint, model, Ui
                 XV1, Vf, XVP, -GAMMA, Uinf
             )
         elseif filament["id"] == "bound"
-            if model === VSM
+            if model == VSM
                 XV2 = filament["x2"]
                 tempvel = velocity_3D_bound_vortex(XV1, XV2, XVP, GAMMA)
             else