Merge branch 'develop' into feature_preconditioning

Author: Cristopher-Morales

SU2_CFD/include/numerics/flow/flow_diffusion.hpp

@@ -349,9 +349,7 @@ class CAvgGradInc_Flow final : public CAvgGrad_Base {
  */
 class CGeneralAvgGrad_Flow final : public CAvgGrad_Base {
 private:
-  su2double Mean_SecVar[2],  /*!< \brief Mean secondary variables. */
-  Mean_Thermal_Conductivity, /*!< \brief Mean value of the thermal conductivity. */
-  Mean_Cp;                   /*!< \brief Mean value of the Cp. */
+  su2double Mean_SecVar[2];  /*!< \brief Mean secondary variables. */
 
   /*!
    * \brief Compute the Jacobian of the heat flux vector

SU2_CFD/src/numerics/flow/flow_sources.cpp

@@ -489,13 +489,9 @@ CNumerics::ResidualType<> CSourceGravity::ComputeResidual(const CConfig* config)
   return ResidualType<>(residual, jacobian, nullptr);
 }
 
-CSourceRotatingFrame_Flow::CSourceRotatingFrame_Flow(unsigned short val_nDim, unsigned short val_nVar, const CConfig* config) :
-                           CSourceBase_Flow(val_nDim, val_nVar, config) {
-
-  Gamma = config->GetGamma();
-  Gamma_Minus_One = Gamma - 1.0;
-
-}
+CSourceRotatingFrame_Flow::CSourceRotatingFrame_Flow(unsigned short val_nDim, unsigned short val_nVar,
+                                                     const CConfig* config)
+    : CSourceBase_Flow(val_nDim, val_nVar, config) {}
 
 CNumerics::ResidualType<> CSourceRotatingFrame_Flow::ComputeResidual(const CConfig* config) {
 
@@ -557,9 +553,6 @@ CSourceIncRotatingFrame_Flow::CSourceIncRotatingFrame_Flow(unsigned short val_nD
 
   implicit = (config->GetKind_TimeIntScheme_Flow() == EULER_IMPLICIT);
 
-  Gamma = config->GetGamma();
-  Gamma_Minus_One = Gamma - 1.0;
-
   /*--- Retrieve the angular velocity vector from config. ---*/
   for (unsigned short iDim = 0; iDim < 3; iDim++)
     Omega[iDim] = config->GetRotation_Rate(iDim)/config->GetOmega_Ref();
@@ -620,10 +613,7 @@ CNumerics::ResidualType<> CSourceIncRotatingFrame_Flow::ComputeResidual(const CC
 
 CSourceVorticityConfinement::CSourceVorticityConfinement(unsigned short val_nDim, unsigned short val_nVar,
                                                          const CConfig* config)
-    : CSourceBase_Flow(val_nDim, val_nVar, config) {
-  Gamma = config->GetGamma();
-  Gamma_Minus_One = Gamma - 1.0;
-}
+    : CSourceBase_Flow(val_nDim, val_nVar, config) {}
 
 CNumerics::ResidualType<> CSourceVorticityConfinement::ComputeResidual(const CConfig* config) {
   /*--- density, \rho ---*/