Merge pull request #1138 from su2code/fix_warnings

pcarruscag · web-flow · commit 4f11b35e0eb8 · 2020-12-18T12:31:33.000Z
Fix compiler warnings
diff --git a/Common/src/geometry/CMultiGridGeometry.cpp b/Common/src/geometry/CMultiGridGeometry.cpp
@@ -1084,7 +1084,7 @@ void CMultiGridGeometry::SetControlVolume(CConfig *config, CGeometry *fine_grid,
 
   unsigned long iFinePoint, iCoarsePoint, iEdge, iParent;
   long FineEdge, CoarseEdge;
-  unsigned short iChildren, iDim;
+  unsigned short iChildren;
   bool change_face_orientation;
   su2double Coarse_Volume, Area;
 
diff --git a/Common/src/toolboxes/CSquareMatrixCM.cpp b/Common/src/toolboxes/CSquareMatrixCM.cpp
@@ -87,7 +87,7 @@ void CSquareMatrixCM::MatMatMult(const char                          side,
     /*--- Left side: mat_out = this * mat_in. Set some sizes
           and allocate the memory for mat_out. ---*/
     const int M = Size(), N = mat_in.cols();
-    assert(M == mat_in.rows());
+    assert(M == static_cast<int>(mat_in.rows()));
 
     mat_out.resize(M,N);
 
@@ -117,7 +117,7 @@ void CSquareMatrixCM::MatMatMult(const char                          side,
     /*--- Right_side: mat_out = mat_in * this. Set some sizes
           and allocate the memory for mat_out. ---*/
     const int M = mat_in.rows(), N = Size();
-    assert(N == mat_in.cols());
+    assert(N == static_cast<int>(mat_in.cols()));
 
     mat_out.resize(M,N);
 
diff --git a/SU2_CFD/include/solvers/CFVMFlowSolverBase.inl b/SU2_CFD/include/solvers/CFVMFlowSolverBase.inl
@@ -1696,7 +1696,7 @@ template <class V, ENUM_REGIME FlowRegime>
 void CFVMFlowSolverBase<V, FlowRegime>::Momentum_Forces(const CGeometry* geometry, const CConfig* config) {
   unsigned long iVertex, iPoint;
   unsigned short iDim, iMarker, Boundary, Monitoring, iMarker_Monitoring;
-  su2double Area, factor, RefVel2 = 0.0, RefTemp, RefDensity = 0.0, Mach2Vel, Mach_Motion, MassFlow, Density;
+  su2double factor, RefVel2 = 0.0, RefTemp, RefDensity = 0.0, Mach2Vel, Mach_Motion, MassFlow, Density;
   const su2double *Normal = nullptr, *Coord = nullptr;
   string Marker_Tag, Monitoring_Tag;
   su2double AxiFactor;
@@ -1790,7 +1790,6 @@ void CFVMFlowSolverBase<V, FlowRegime>::Momentum_Forces(const CGeometry* geometr
           Normal = geometry->vertex[iMarker][iVertex]->GetNormal();
           Coord = geometry->nodes->GetCoord(iPoint);
           Density = nodes->GetDensity(iPoint);
-          Area = GeometryToolbox::Norm(nDim, Normal);
           MassFlow = 0.0;
           su2double Velocity[MAXNDIM] = {0.0}, MomentDist[MAXNDIM] = {0.0};
           for (iDim = 0; iDim < nDim; iDim++) {
diff --git a/SU2_CFD/src/interfaces/cht/CConjugateHeatInterface.cpp b/SU2_CFD/src/interfaces/cht/CConjugateHeatInterface.cpp
@@ -40,45 +40,28 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
                                                 const CConfig *donor_config, unsigned long Marker_Donor,
                                                 unsigned long Vertex_Donor, unsigned long Point_Donor) {
 
-  unsigned short nDim, iDim;
-  unsigned long iPoint, PointNormal;
+  const auto nDim = donor_geometry->GetnDim();
 
-  su2double *Coord, *Coord_Normal, *Normal, *Edge_Vector, dist, dist2, Area,
-      Twall, Tnormal, dTdn, rho_cp_solid, Prandtl_Lam, laminar_viscosity,
+  su2double Twall, Tnormal, dTdn, rho_cp_solid, Prandtl_Lam, laminar_viscosity,
       thermal_diffusivity, thermal_conductivity=0.0, thermal_conductivityND,
       heat_flux_density=0.0, conductivity_over_dist=0.0;
 
-  nDim = donor_geometry->GetnDim();
-
-  Edge_Vector = new su2double[nDim];
-
   /*--- Check whether the current zone is a solid zone or a fluid zone ---*/
 
-  bool compressible_flow = ((donor_config->GetKind_Solver() == NAVIER_STOKES)
-               || (donor_config->GetKind_Solver() == RANS)
-               || (donor_config->GetKind_Solver() == DISC_ADJ_NAVIER_STOKES)
-               || (donor_config->GetKind_Solver() == DISC_ADJ_RANS));
-  bool incompressible_flow = (((donor_config->GetKind_Solver() == INC_NAVIER_STOKES)
-               || (donor_config->GetKind_Solver() == INC_RANS)
-               || (donor_config->GetKind_Solver() == DISC_ADJ_INC_NAVIER_STOKES)
-               || (donor_config->GetKind_Solver() == DISC_ADJ_INC_RANS))
-               && (donor_config->GetEnergy_Equation()));
-  bool heat_equation = ((donor_config->GetKind_Solver() == HEAT_EQUATION)
-               || (donor_config->GetKind_Solver() == DISC_ADJ_HEAT));
+  const bool compressible_flow = (donor_config->GetKind_Regime() == COMPRESSIBLE);
+  const bool incompressible_flow = (donor_config->GetKind_Regime() == INCOMPRESSIBLE) && donor_config->GetEnergy_Equation();
+  const bool heat_equation = (donor_config->GetKind_Solver() == HEAT_EQUATION) ||
+                             (donor_config->GetKind_Solver() == DISC_ADJ_HEAT);
 
-  Coord         = donor_geometry->nodes->GetCoord(Point_Donor);
+  const auto Coord = donor_geometry->nodes->GetCoord(Point_Donor);
+  const auto PointNormal = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNormal_Neighbor();
+  const auto Coord_Normal = donor_geometry->nodes->GetCoord(PointNormal);
 
-  Normal        = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNormal();
-  PointNormal   = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNormal_Neighbor();
-  Coord_Normal  = donor_geometry->nodes->GetCoord(PointNormal);
+  Twall = 0.0; Tnormal = 0.0; dTdn = 0.0;
 
-  Twall = 0.0; Tnormal = 0.0; dTdn = 0.0; dist2 = 0.0; Area = 0.0;
-
-  for (iDim = 0; iDim < nDim; iDim++) {
-    Edge_Vector[iDim] = Coord_Normal[iDim] - Coord[iDim];
-  }
-  dist = GeometryToolbox::Norm(nDim, Edge_Vector);
-  Area = GeometryToolbox::Norm(nDim, Normal);
+  su2double Edge_Vector[3] = {0.0};
+  GeometryToolbox::Distance(nDim, Coord_Normal, Coord, Edge_Vector);
+  su2double dist = GeometryToolbox::Norm(nDim, Edge_Vector);
 
   /*--- Retrieve temperature solution and its gradient ---*/
 
@@ -101,6 +84,8 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
     Tnormal = donor_solution->GetNodes()->GetSolution(PointNormal,0);
 
     // TODO: Check if these improve accuracy, if needed at all
+    //    const auto Normal = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNormal();
+    //    su2double Area = GeometryToolbox::Norm(nDim, Normal);
     //    for (iDim = 0; iDim < nDim; iDim++) {
     //      dTdn += (Twall - Tnormal)/dist * (Edge_Vector[iDim]/dist) * (Normal[iDim]/Area);
     //    }
@@ -136,7 +121,7 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
   }
   else if (incompressible_flow) {
 
-    iPoint = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNode();
+    const auto iPoint = donor_geometry->vertex[Marker_Donor][Vertex_Donor]->GetNode();
 
     thermal_conductivityND  = donor_solution->GetNodes()->GetThermalConductivity(iPoint);
     heat_flux_density       = thermal_conductivityND*dTdn;
@@ -188,10 +173,6 @@ void CConjugateHeatInterface::GetDonor_Variable(CSolver *donor_solution, CGeomet
     Donor_Variable[3] = Tnormal*donor_config->GetTemperature_Ref();
   }
 
-
-
-  delete [] Edge_Vector;
-
 }
 
 void CConjugateHeatInterface::SetTarget_Variable(CSolver *target_solution, CGeometry *target_geometry,
diff --git a/SU2_CFD/src/numerics/NEMO/convection/lax.cpp b/SU2_CFD/src/numerics/NEMO/convection/lax.cpp
@@ -32,8 +32,8 @@ CCentLax_NEMO::CCentLax_NEMO(unsigned short val_nDim,
                              unsigned short val_nVar,
                              unsigned short val_nPrimVar,
                              unsigned short val_nPrimVarGrad,
-                             CConfig *config) : CNEMONumerics(val_nDim, val_nVar, 
-                                                              val_nPrimVar, 
+                             CConfig *config) : CNEMONumerics(val_nDim, val_nVar,
+                                                              val_nPrimVar,
                                                               val_nPrimVarGrad,
                                                               config) {
 
@@ -47,14 +47,14 @@ CCentLax_NEMO::CCentLax_NEMO(unsigned short val_nDim,
   MeanV    = new su2double[nPrimVar];
   MeandPdU = new su2double[nVar];
   ProjFlux = new su2double[nVar];
-  Flux     = new su2double[nVar]; 
+  Flux     = new su2double[nVar];
 
   mean_eves.resize(nSpecies,0.0);
 
 }
 
 CCentLax_NEMO::~CCentLax_NEMO(void) {
-  
+
   delete [] Diff_U;
   delete [] MeanU;
   delete [] MeanV;
@@ -65,7 +65,7 @@ CCentLax_NEMO::~CCentLax_NEMO(void) {
 
 CNumerics::ResidualType<> CCentLax_NEMO::ComputeResidual(const CConfig *config) {
 
-  unsigned short iDim, iSpecies, iVar;
+  unsigned short iDim, iVar;
   su2double rho_i, rho_j, h_i, h_j, a_i, a_j;
   su2double ProjVel_i, ProjVel_j;
 
@@ -87,7 +87,7 @@ CNumerics::ResidualType<> CCentLax_NEMO::ComputeResidual(const CConfig *config)
     MeanV[iVar] = 0.5*(V_i[iVar]+V_j[iVar]);
 
   /*--- Compute NonEq specific variables ---*/
-  vector<su2double> mean_eves = fluidmodel->ComputeSpeciesEve(MeanV[TVE_INDEX]); 
+  vector<su2double> mean_eves = fluidmodel->ComputeSpeciesEve(MeanV[TVE_INDEX]);
   fluidmodel->ComputedPdU(MeanV, mean_eves, MeandPdU);
 
   /*--- Get projected flux tensor ---*/
@@ -104,12 +104,12 @@ CNumerics::ResidualType<> CCentLax_NEMO::ComputeResidual(const CConfig *config)
   /*--- Dissipation --*/
   Local_Lambda_i = (fabs(ProjVel_i)+a_i*Area);
   Local_Lambda_j = (fabs(ProjVel_j)+a_j*Area);
-  MeanLambda = 0.5*(Local_Lambda_i+Local_Lambda_j);  
+  MeanLambda = 0.5*(Local_Lambda_i+Local_Lambda_j);
 
   Phi_i = pow(Lambda_i/(4.0*MeanLambda+EPS),Param_p);
   Phi_j = pow(Lambda_j/(4.0*MeanLambda+EPS),Param_p);
   StretchingFactor = 4.0*Phi_i*Phi_j/(Phi_i+Phi_j+EPS);
-  
+
   /*--- Computes differences btw. conservative variables ---*/
   for (iVar = 0; iVar < nVar; iVar++)
     Diff_U[iVar] = U_i[iVar] - U_j[iVar];
diff --git a/SU2_CFD/src/output/CFlowOutput.cpp b/SU2_CFD/src/output/CFlowOutput.cpp
@@ -811,7 +811,7 @@ void CFlowOutput::Add_CpInverseDesignOutput(CConfig *config){
 
 void CFlowOutput::Set_CpInverseDesign(CSolver *solver, CGeometry *geometry, CConfig *config){
 
-  unsigned short iMarker, icommas, Boundary, iDim;
+  unsigned short iMarker, icommas, Boundary;
   unsigned long iVertex, iPoint, (*Point2Vertex)[2], nPointLocal = 0, nPointGlobal = 0;
   su2double XCoord, YCoord, ZCoord, Pressure, PressureCoeff = 0, Cp, CpTarget, *Normal = nullptr, Area, PressDiff = 0.0;
   bool *PointInDomain;
diff --git a/SU2_CFD/src/python_wrapper_structure.cpp b/SU2_CFD/src/python_wrapper_structure.cpp
@@ -634,7 +634,6 @@ passivedouble CDriver::GetThermalConductivity(unsigned short iMarker, unsigned l
 
 vector<passivedouble> CDriver::GetVertexUnitNormal(unsigned short iMarker, unsigned long iVertex){
 
-  unsigned short iDim;
   su2double *Normal;
   su2double Area;
   vector<su2double> ret_Normal(3, 0.0);
diff --git a/SU2_CFD/src/solvers/CAdjEulerSolver.cpp b/SU2_CFD/src/solvers/CAdjEulerSolver.cpp
@@ -56,7 +56,7 @@ CAdjEulerSolver::CAdjEulerSolver(CGeometry *geometry, CConfig *config, unsigned
   unsigned short iDim, iVar, iMarker, nLineLets;
   ifstream restart_file;
   string filename, AdjExt;
-  su2double myArea_Monitored, Area, *Normal;
+  su2double myArea_Monitored, *Normal;
 
   adjoint = true;
 
@@ -4385,7 +4385,7 @@ void CAdjEulerSolver::BC_Engine_Inflow(CGeometry *geometry, CSolver **solver_con
 void CAdjEulerSolver::BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker) {
 
   unsigned long iVertex, iPoint, Point_Normal;
-  su2double Area, *Normal, *V_domain, *V_exhaust, *Psi_domain, *Psi_exhaust;
+  su2double *Normal, *V_domain, *V_exhaust, *Psi_domain, *Psi_exhaust;
   unsigned short iVar, iDim;
 
   bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
@@ -4409,8 +4409,6 @@ void CAdjEulerSolver::BC_Engine_Exhaust(CGeometry *geometry, CSolver **solver_co
       for (iDim = 0; iDim < nDim; iDim++) Normal[iDim] = -Normal[iDim];
       conv_numerics->SetNormal(Normal);
 
-      Area = GeometryToolbox::Norm(nDim, Normal);
-
       /*--- Index of the closest interior node ---*/
 
       Point_Normal = geometry->vertex[val_marker][iVertex]->GetNormal_Neighbor();
diff --git a/SU2_CFD/src/solvers/CAdjNSSolver.cpp b/SU2_CFD/src/solvers/CAdjNSSolver.cpp
@@ -43,7 +43,7 @@ CAdjNSSolver::CAdjNSSolver(CGeometry *geometry, CConfig *config, unsigned short
   su2double RefDensity  = config->GetDensity_FreeStreamND();
   su2double Gas_Constant    = config->GetGas_ConstantND();
   su2double Mach_Motion     = config->GetMach_Motion();
-  su2double Area=0.0, *Normal = nullptr, myArea_Monitored;
+  su2double *Normal = nullptr, myArea_Monitored;
   su2double RefVel2, Mach2Vel, Weight_ObjFunc, factor;
   su2double *Velocity_Inf;
 
diff --git a/SU2_CFD/src/solvers/CFEM_DG_EulerSolver.cpp b/SU2_CFD/src/solvers/CFEM_DG_EulerSolver.cpp
@@ -6699,7 +6699,7 @@ void CFEM_DG_EulerSolver::MultiplyResidualByInverseMassMatrix(
 
       /* Multiply the residual with the inverse of the mass matrix.
          Use the array workArray as temporary storage. */
-      for(unsigned long mm=0; mm<nVar*volElem[l].nDOFsSol; ++mm)
+      for(unsigned short mm=0; mm<nVar*volElem[l].nDOFsSol; ++mm)
         workArray[mm] = res[mm];
 
       blasFunctions->gemm(volElem[l].nDOFsSol, nVar, volElem[l].nDOFsSol,
diff --git a/SU2_CFD/src/solvers/CHeatSolver.cpp b/SU2_CFD/src/solvers/CHeatSolver.cpp
@@ -709,7 +709,7 @@ void CHeatSolver::Viscous_Residual(CGeometry *geometry, CSolver **solver_contain
 
 void CHeatSolver::Set_Heatflux_Areas(CGeometry *geometry, CConfig *config) {
 
-  unsigned short iMarker, iMarker_HeatFlux, Monitoring, iDim;
+  unsigned short iMarker, iMarker_HeatFlux, Monitoring;
   unsigned long iPoint, iVertex;
   string HeatFlux_Tag, Marker_Tag;
 
@@ -831,7 +831,6 @@ void CHeatSolver::BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_conta
 void CHeatSolver::BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config,
                                                      unsigned short val_marker) {
 
-  unsigned short iDim;
   unsigned long iVertex, iPoint;
   su2double Wall_HeatFlux, Area, *Normal;
 
@@ -1067,7 +1066,6 @@ void CHeatSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container,
 void CHeatSolver::BC_ConjugateHeat_Interface(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics, CConfig *config, unsigned short val_marker) {
 
   unsigned long iVertex, iPoint;
-  unsigned short iDim;
 
   su2double thermal_diffusivity, rho_cp_solid, Temperature_Ref, T_Conjugate, Tinterface,
       Tnormal_Conjugate, HeatFluxDensity, HeatFlux, Area;
@@ -1283,7 +1281,7 @@ void CHeatSolver::Heat_Fluxes(CGeometry *geometry, CSolver **solver_container, C
 void CHeatSolver::SetTime_Step(CGeometry *geometry, CSolver **solver_container, CConfig *config,
                                unsigned short iMesh, unsigned long Iteration) {
 
-  unsigned short iDim, iMarker;
+  unsigned short iMarker;
   unsigned long iEdge, iVertex, iPoint = 0, jPoint = 0;
   su2double Area, Vol, laminar_viscosity, eddy_viscosity, thermal_diffusivity, Prandtl_Lam, Prandtl_Turb, Mean_ProjVel, Mean_BetaInc2, Mean_DensityInc, Mean_SoundSpeed, Lambda;
   su2double Global_Delta_Time = 0.0, Global_Delta_UnstTimeND = 0.0, Local_Delta_Time = 0.0, Local_Delta_Time_Inv, Local_Delta_Time_Visc, CFL_Reduction, K_v = 0.25;
diff --git a/SU2_CFD/src/solvers/CIncEulerSolver.cpp b/SU2_CFD/src/solvers/CIncEulerSolver.cpp
@@ -2742,7 +2742,6 @@ void CIncEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container,
                              CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker) {
   unsigned short iDim;
   unsigned long iVertex, iPoint, Point_Normal;
-  su2double Area;
   su2double *V_outlet, *V_domain, P_Outlet = 0.0, P_domain;
   su2double mDot_Target, mDot_Old, dP, Density_Avg, Area_Outlet;
   su2double Damping = config->GetInc_Outlet_Damping();
@@ -2779,8 +2778,6 @@ void CIncEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container,
       for (iDim = 0; iDim < nDim; iDim++) Normal[iDim] = -Normal[iDim];
       conv_numerics->SetNormal(Normal);
 
-      Area = GeometryToolbox::Norm(nDim, Normal);
-
       /*--- Current solution at this boundary node ---*/
 
       V_domain = nodes->GetPrimitive(iPoint);
diff --git a/SU2_CFD/src/solvers/CNEMOEulerSolver.cpp b/SU2_CFD/src/solvers/CNEMOEulerSolver.cpp
@@ -229,7 +229,7 @@ CNEMOEulerSolver::CNEMOEulerSolver(CGeometry *geometry, CConfig *config,
     Density_Inf     = FluidModel->GetDensity();
     Soundspeed_Inf  = FluidModel->GetSoundSpeed();
     Gamma           = FluidModel->ComputeGamma();
-    Gamma_Minus_One = Gamma - 1.0; 
+    Gamma_Minus_One = Gamma - 1.0;
 
     sqvel = 0.0;
     for (iDim = 0; iDim < nDim; iDim++){
@@ -289,7 +289,7 @@ CNEMOEulerSolver::~CNEMOEulerSolver(void) {
 void CNEMOEulerSolver::SetInitialCondition(CGeometry **geometry, CSolver ***solver_container, CConfig *config, unsigned long TimeIter) {
 
   unsigned long iPoint;
-  unsigned short iMesh, iDim;
+  unsigned short iMesh;
   const bool restart = (config->GetRestart() || config->GetRestart_Flow());
   const bool rans = false;
   const bool dual_time = ((config->GetTime_Marching() == DT_STEPPING_1ST) ||
@@ -341,7 +341,6 @@ void CNEMOEulerSolver::SetInitialCondition(CGeometry **geometry, CSolver ***solv
 void CNEMOEulerSolver::CommonPreprocessing(CGeometry *geometry, CSolver **solver_container, CConfig *config, unsigned short iMesh,
                                            unsigned short iRKStep, unsigned short RunTime_EqSystem, bool Output) {
 
-  bool disc_adjoint     = config->GetDiscrete_Adjoint();
   bool implicit         = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
   bool center           = (config->GetKind_ConvNumScheme_Flow() == SPACE_CENTERED);
   bool center_jst       = (config->GetKind_Centered_Flow() == JST) && (iMesh == MESH_0);
@@ -915,7 +914,7 @@ void CNEMOEulerSolver::Upwind_Residual(CGeometry *geometry, CSolver **solver_con
         numerics->SetdTvedU(dTvedU_i, dTvedU_j);
         numerics->SetEve   (Eve_i,    Eve_j   );
         numerics->SetCvve  (Cvve_i,   Cvve_j  );
-        numerics->SetGamma (Gamma_i,  Gamma_i );
+        numerics->SetGamma (Gamma_i,  Gamma_j );
       }
     } else {
       /*--- Set variables without reconstruction ---*/
@@ -2242,8 +2241,8 @@ void CNEMOEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solution_contain
 
       /*--- Compute Gamma using domain state ---*/
       Gamma = nodes->GetGamma(iPoint);
-      Gamma_Minus_One = Gamma - 1.0; 
-      
+      Gamma_Minus_One = Gamma - 1.0;
+
       /*--- Retrieve the specified back pressure for this outlet. ---*/
       if (gravity) P_Exit = config->GetOutlet_Pressure(Marker_Tag) - geometry->nodes->GetCoord(iPoint, nDim-1)*STANDARD_GRAVITY;
       else         P_Exit = config->GetOutlet_Pressure(Marker_Tag);
diff --git a/SU2_CFD/src/solvers/CNEMONSSolver.cpp b/SU2_CFD/src/solvers/CNEMONSSolver.cpp
diff --git a/SU2_CFD/src/solvers/CRadP1Solver.cpp b/SU2_CFD/src/solvers/CRadP1Solver.cpp
diff --git a/SU2_CFD/src/solvers/CTurbSASolver.cpp b/SU2_CFD/src/solvers/CTurbSASolver.cpp
diff --git a/SU2_CFD/src/solvers/CTurbSSTSolver.cpp b/SU2_CFD/src/solvers/CTurbSSTSolver.cpp
