Addressed Comments

Author: fmpmorgado

SU2_CFD/src/solvers/CNEMOEulerSolver.cpp

@@ -256,7 +256,7 @@ void CNEMOEulerSolver::CommonPreprocessing(CGeometry *geometry, CSolver **solver
       config->SetNonphysical_Points(ErrorCounter);
 
       if ((rank == MASTER_NODE) && (ErrorCounter != 0))
-        cout << "Warning. The original solution contains "<< ErrorCounter << " points that are not physical." << endl;
+        cout << "Warning. The initial solution contains "<< ErrorCounter << " points that are not physical." << endl;
   }
 
   /*--- Artificial dissipation ---*/

SU2_CFD/src/solvers/CSolverFactory.cpp

@@ -429,7 +429,6 @@ CSolver* CSolverFactory::CreateFlowSolver(SUB_SOLVER_TYPE kindFlowSolver, CSolve
       break;
     case SUB_SOLVER_TYPE::NEMO_EULER:
       flowSolver = new CNEMOEulerSolver(geometry, config, iMGLevel);
-      //flowSolver->Preprocessing(geometry, solver, config, iMGLevel, NO_RK_ITER, RUNTIME_FLOW_SYS, false);
       break;
     case SUB_SOLVER_TYPE::NEMO_NAVIER_STOKES:
       flowSolver = new CNEMONSSolver(geometry, config, iMGLevel);

SU2_CFD/src/variables/CNEMOEulerVariable.cpp

@@ -45,10 +45,7 @@ CNEMOEulerVariable::CNEMOEulerVariable(su2double val_pressure,
                                                                          config ),
                                        Gradient_Reconstruction(config->GetReconstructionGradientRequired() ? Gradient_Aux : Gradient_Primitive) {
 
-  vector<su2double> energies;
   unsigned short iDim, iSpecies;
-  su2double soundspeed, rho;
-  su2double sqvel = 0.0;
 
   /*--- Setting variable amounts ---*/
   nDim            = ndim;
@@ -148,12 +145,10 @@ CNEMOEulerVariable::CNEMOEulerVariable(su2double val_pressure,
   fluidmodel->SetTDStatePTTv(val_pressure, val_massfrac, val_temperature, val_temperature_ve);
 
   /*--- Compute necessary quantities ---*/
-  rho = fluidmodel->GetDensity();
-  soundspeed = fluidmodel->ComputeSoundSpeed();
-  for (iDim = 0; iDim < nDim; iDim++){
-    sqvel += val_mach[iDim]*soundspeed * val_mach[iDim]*soundspeed;
-  }
-  energies = fluidmodel->ComputeMixtureEnergies();
+  const su2double rho = fluidmodel->GetDensity();
+  const su2double soundspeed = fluidmodel->ComputeSoundSpeed();
+  const su2double sqvel = GeometryToolbox::SquaredNorm(nDim, val_mach) * pow(soundspeed,2);
+  const auto& energies = fluidmodel->ComputeMixtureEnergies();
 
   /*--- Loop over all points --*/
   for(unsigned long iPoint = 0; iPoint < nPoint; ++iPoint){