Remove errors from #1281. No auto for su2double with rhs math operation and MAXNVAR fix for flow.

TobiKattmann · TobiKattmann · commit 19d3329da5b5 · 2021-05-27T12:05:32.000+02:00
diff --git a/SU2_CFD/src/solvers/CHeatSolver.cpp b/SU2_CFD/src/solvers/CHeatSolver.cpp
@@ -363,8 +363,8 @@ void CHeatSolver::Upwind_Residual(CGeometry *geometry, CSolver **solver_containe
     nVarFlow = solver_container[FLOW_SOL]->GetnVar();
 
     /*--- Define some auxiliary vectors related to the primitive flow solution ---*/
-    su2double Primitive_Flow_i[MAXNVAR] = {0};
-    su2double Primitive_Flow_j[MAXNVAR] = {0};
+    vector<su2double> Primitive_Flow_i(nVarFlow, 0.0);
+    vector<su2double> Primitive_Flow_j(nVarFlow, 0.0);
 
     for (auto iEdge = 0ul; iEdge < geometry->GetnEdge(); iEdge++) {
 
@@ -421,7 +421,7 @@ void CHeatSolver::Upwind_Residual(CGeometry *geometry, CSolver **solver_containe
         Temp_i_Corrected = Temp_i + Project_Temp_i_Grad;
         Temp_j_Corrected = Temp_j + Project_Temp_j_Grad;
 
-        numerics->SetPrimitive(Primitive_Flow_i, Primitive_Flow_j);
+        numerics->SetPrimitive(Primitive_Flow_i.data(), Primitive_Flow_j.data());
         numerics->SetTemperature(Temp_i_Corrected, Temp_j_Corrected);
       }
 
@@ -571,9 +571,9 @@ void CHeatSolver::BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_conta
   su2double *Normal, *Coord_i, *Coord_j, Area, dist_ij, Twall, dTdn;
   const bool implicit = (config->GetKind_TimeIntScheme_Flow() == EULER_IMPLICIT);
 
-  const auto laminar_viscosity = config->GetMu_ConstantND();
-  const auto Prandtl_Lam = config->GetPrandtl_Lam();
-  const auto thermal_diffusivity = flow ? laminar_viscosity/Prandtl_Lam : config->GetThermalDiffusivity_Solid();
+  const su2double laminar_viscosity = config->GetMu_ConstantND();
+  const su2double Prandtl_Lam = config->GetPrandtl_Lam();
+  const su2double thermal_diffusivity = flow ? laminar_viscosity/Prandtl_Lam : config->GetThermalDiffusivity_Solid();
 
   //su2double Prandtl_Turb = config->GetPrandtl_Turb();
   //laminar_viscosity = config->GetViscosity_FreeStreamND(); // TDE check for consistency for CHT
@@ -832,8 +832,8 @@ void CHeatSolver::BC_ConjugateHeat_Interface(CGeometry *geometry, CSolver **solv
 
   const bool implicit = (config->GetKind_TimeIntScheme_Flow() == EULER_IMPLICIT);
 
-  const auto Temperature_Ref = config->GetTemperature_Ref();
-  const auto rho_cp_solid = config->GetDensity_Solid()*config->GetSpecific_Heat_Cp();
+  const su2double Temperature_Ref = config->GetTemperature_Ref();
+  const su2double rho_cp_solid = config->GetDensity_Solid()*config->GetSpecific_Heat_Cp();
 
   if (flow) {
 
