cleanup

Author: bigfooted

SU2_CFD/include/solvers/CFVMFlowSolverBase.inl

@@ -125,7 +125,7 @@ void CFVMFlowSolverBase<V, R>::Allocate(const CConfig& config) {
   AllocVectorOfMatrices(nVertex, nDim, Inlet_FlowDir);
   PointSource.resize(nPointDomain,nVar);
   PointSource.setConstant(0.0);
-  
+
   /*--- Force definition and coefficient arrays for all of the markers ---*/
 
   AllocVectorOfVectors(nVertex, CPressure);
@@ -1005,8 +1005,6 @@ void CFVMFlowSolverBase<V, R>::LoadRestart_impl(CGeometry **geometry, CSolver **
   END_SU2_OMP_SAFE_GLOBAL_ACCESS
 }
 
-
-
 template <class V, ENUM_REGIME R>
 void CFVMFlowSolverBase<V, R>::LoadRestart(CGeometry **geometry, CSolver ***solver,
                                            CConfig *config, int iter, bool update_geo) {

SU2_CFD/include/solvers/CIncEulerSolver.hpp

@@ -204,19 +204,6 @@ class CIncEulerSolver : public CFVMFlowSolverBase<CIncEulerVariable, ENUM_REGIME
                        CNumerics **numerics_container,
                        CConfig *config,
                        unsigned short iMesh) final;
-  /*!
-   * \brief Source term integration.
-   * \param[in] geometry - Geometrical definition of the problem.
-   * \param[in] solver_container - Container vector with all the solutions.
-   * \param[in] numerics_container - Description of the numerical method.
-   * \param[in] config - Definition of the particular problem.
-   * \param[in] iMesh - Index of the mesh in multigrid computations.
-   */
-  // void Custom_Source_Residual(CGeometry *geometry,
-  //                      CSolver **solver_container,
-  //                      CNumerics **numerics_container,
-  //                      CConfig *config,
-  //                      unsigned short iMesh) final;
 
   /*!
    * \brief Source term integration.

SU2_CFD/include/solvers/CSolver.hpp

@@ -1557,7 +1557,6 @@ class CSolver {
    * \param[in] geometry - Geometrical definition of the problem.
    * \param[in] solver_container - Container vector with all the solutions.
    * \param[in] numerics_container - Description of the numerical method.
-   * \param[in] second_numerics - Description of the second numerical method.
    * \param[in] config - Definition of the particular problem.
    * \param[in] iMesh - Index of the mesh in multigrid computations.
    */
@@ -1571,7 +1570,6 @@ class CSolver {
    * \param[in] geometry - Geometrical definition of the problem.
    * \param[in] solver_container - Container vector with all the solutions.
    * \param[in] numerics_container - Description of the numerical method.
-   * \param[in] second_numerics - Description of the second numerical method.
    * \param[in] config - Definition of the particular problem.
    * \param[in] iMesh - Index of the mesh in multigrid computations.
    */

SU2_CFD/include/solvers/CSpeciesSolver.hpp

@@ -40,7 +40,7 @@ class CSpeciesSolver : public CScalarSolver<CSpeciesVariable> {
  protected:
   unsigned short Inlet_Position;             /*!< \brief Column index for scalar variables in inlet files. */
   vector<su2activematrix> Inlet_SpeciesVars; /*!< \brief Species variables at inlet profiles. */
-  su2activematrix SpeciesPointSource;        /*!< \brief Value of the Flow Direction. */
+  su2activematrix SpeciesPointSource;        /*!< \brief User defined source term. */
 
  public:
   /*!

SU2_CFD/src/numerics/flow/flow_sources.cpp

@@ -73,20 +73,13 @@ CNumerics::ResidualType<> CSourceAxisymmetric_Flow::ComputeResidual(const CConfi
       sq_vel += Velocity_i *Velocity_i;
     }
 
-    // P = (gamma-1)*rho * e
-    // E = e+0.5*U^2
-    // do we have primitive variables available?
     Pressure_i = Gamma_Minus_One*U_i[0]*(U_i[nDim+1]/U_i[0]-0.5*sq_vel);
-    // (rho*e +p)/rho = e+p/rho
     Enthalpy_i = (U_i[nDim+1] + Pressure_i) / U_i[0];
-    //cout <<"P="<<V_i[nDim+1] << " rho=" << V_i[nDim+2] << ", "<<U_i[0]<<" " << Pressure_i<<" " << Enthalpy_i <<" "<< V_i[nDim+3]<< endl;
 
     residual[0] = yinv*Volume*U_i[2];
     residual[1] = yinv*Volume*U_i[1]*U_i[2]/U_i[0];
     residual[2] = yinv*Volume*(U_i[2]*U_i[2]/U_i[0]);
-    // we need: (rho*e+p)
     residual[3] = yinv*Volume*Enthalpy_i*U_i[2];
-    //residual[3] = yinv*Volume*(U_i[nDim+1]+Pressure_i)*U_i[2];
 
     /*--- Inviscid component of the source term. ---*/
 

SU2_CFD/src/solvers/CIncEulerSolver.cpp

@@ -1813,38 +1813,6 @@ void CIncEulerSolver::Source_Residual(CGeometry *geometry, CSolver **solver_cont
 
 }
 
-// void CIncEulerSolver::Custom_Source_Residual(CGeometry *geometry, CSolver **solver_container,
-//                                       CNumerics **numerics_container, CConfig *config, unsigned short iMesh) {
-
-//   /*--- Pick one numerics object per thread. ---*/
-//   CNumerics* numerics = numerics_container[SOURCE_SECOND_TERM + omp_get_thread_num()*MAX_TERMS];
-
-//   unsigned short iVar;
-//   unsigned long iPoint;
-//   AD::StartNoSharedReading();
-
-//   SU2_OMP_FOR_STAT(omp_chunk_size)
-//   for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
-
-//     /*--- Load the volume of the dual mesh cell ---*/
-
-//     numerics->SetVolume(geometry->nodes->GetVolume(iPoint));
-
-//     /*--- Get control volume size. ---*/
-//     su2double Volume = geometry->nodes->GetVolume(iPoint);
-
-//     /*--- Compute the residual for this control volume and subtract. ---*/
-//     for (iVar = 0; iVar < nVar; iVar++) {
-//       LinSysRes[iPoint*nVar+iVar] += PointSource[iPoint][iVar] * Volume;
-//     }
-//   }
-//   END_SU2_OMP_FOR
-
-//   AD::EndNoSharedReading();
-
-// }
-
-
 void CIncEulerSolver::Source_Template(CGeometry *geometry, CSolver **solver_container, CNumerics *numerics,
                                    CConfig *config, unsigned short iMesh) {