2424 * License along with SU2. If not, see <http://www.gnu.org/licenses/>.
2525 */
2626
27- #include " ../../include/solvers/CScalarSolver.hpp"
28-
2927#include " ../../../Common/include/parallelization/omp_structure.hpp"
3028#include " ../../../Common/include/toolboxes/geometry_toolbox.hpp"
29+ #include " ../../include/solvers/CScalarSolver.hpp"
3130
32- template <class TVariable >
31+ template <class TVariable >
3332CScalarSolver<TVariable>::CScalarSolver(bool conservative) : CSolver(), Conservative(conservative) {}
3433
35- template <class TVariable >
34+ template <class TVariable >
3635CScalarSolver<TVariable>::CScalarSolver(CGeometry* geometry, CConfig* config, bool conservative)
3736 : CSolver(), Conservative(conservative) {
3837 nMarker = config->GetnMarker_All ();
@@ -75,12 +74,14 @@ CScalarSolver<TVariable>::CScalarSolver(CGeometry* geometry, CConfig* config, bo
7574 }
7675}
7776
78- template <class TVariable >
79- CScalarSolver<TVariable>::~CScalarSolver () { delete nodes; }
77+ template <class TVariable >
78+ CScalarSolver<TVariable>::~CScalarSolver () {
79+ delete nodes;
80+ }
8081
81- template <class TVariable >
82- void CScalarSolver<TVariable>::Upwind_Residual(CGeometry* geometry, CSolver** solver_container, CNumerics** numerics_container,
83- CConfig* config, unsigned short iMesh) {
82+ template <class TVariable >
83+ void CScalarSolver<TVariable>::Upwind_Residual(CGeometry* geometry, CSolver** solver_container,
84+ CNumerics** numerics_container, CConfig* config, unsigned short iMesh) {
8485 const bool implicit = (config->GetKind_TimeIntScheme () == EULER_IMPLICIT);
8586 const bool muscl = config->GetMUSCL_Turb ();
8687 const bool limiter = (config->GetKind_SlopeLimit_Turb () != NO_LIMITER);
@@ -239,7 +240,7 @@ void CScalarSolver<TVariable>::Upwind_Residual(CGeometry* geometry, CSolver** so
239240 }
240241}
241242
242- template <class TVariable >
243+ template <class TVariable >
243244void CScalarSolver<TVariable>::SumEdgeFluxes(CGeometry* geometry) {
244245 SU2_OMP_FOR_STAT (omp_chunk_size)
245246 for (unsigned long iPoint = 0 ; iPoint < nPoint; ++iPoint) {
@@ -255,8 +256,9 @@ void CScalarSolver<TVariable>::SumEdgeFluxes(CGeometry* geometry) {
255256 END_SU2_OMP_FOR
256257}
257258
258- template <class TVariable >
259- void CScalarSolver<TVariable>::BC_Periodic(CGeometry* geometry, CSolver** solver_container, CNumerics* numerics, CConfig* config) {
259+ template <class TVariable >
260+ void CScalarSolver<TVariable>::BC_Periodic(CGeometry* geometry, CSolver** solver_container, CNumerics* numerics,
261+ CConfig* config) {
260262 /* --- Complete residuals for periodic boundary conditions. We loop over
261263 the periodic BCs in matching pairs so that, in the event that there are
262264 adjacent periodic markers, the repeated points will have their residuals
@@ -269,8 +271,9 @@ void CScalarSolver<TVariable>::BC_Periodic(CGeometry* geometry, CSolver** solver
269271 }
270272}
271273
272- template <class TVariable >
273- void CScalarSolver<TVariable>::PrepareImplicitIteration(CGeometry* geometry, CSolver** solver_container, CConfig* config) {
274+ template <class TVariable >
275+ void CScalarSolver<TVariable>::PrepareImplicitIteration(CGeometry* geometry, CSolver** solver_container,
276+ CConfig* config) {
274277 const auto flowNodes = solver_container[FLOW_SOL]->GetNodes ();
275278
276279 /* --- Set shared residual variables to 0 and declare
@@ -329,8 +332,9 @@ void CScalarSolver<TVariable>::PrepareImplicitIteration(CGeometry* geometry, CSo
329332 SetResidual_RMS (geometry, config);
330333}
331334
332- template <class TVariable >
333- void CScalarSolver<TVariable>::CompleteImplicitIteration(CGeometry* geometry, CSolver** solver_container, CConfig* config) {
335+ template <class TVariable >
336+ void CScalarSolver<TVariable>::CompleteImplicitIteration(CGeometry* geometry, CSolver** solver_container,
337+ CConfig* config) {
334338 const bool compressible = (config->GetKind_Regime () == ENUM_REGIME::COMPRESSIBLE);
335339
336340 const auto flowNodes = solver_container[FLOW_SOL]->GetNodes ();
@@ -376,8 +380,9 @@ void CScalarSolver<TVariable>::CompleteImplicitIteration(CGeometry* geometry, CS
376380 CompleteComms (geometry, config, SOLUTION_EDDY);
377381}
378382
379- template <class TVariable >
380- void CScalarSolver<TVariable>::ImplicitEuler_Iteration(CGeometry* geometry, CSolver** solver_container, CConfig* config) {
383+ template <class TVariable >
384+ void CScalarSolver<TVariable>::ImplicitEuler_Iteration(CGeometry* geometry, CSolver** solver_container,
385+ CConfig* config) {
381386 PrepareImplicitIteration (geometry, solver_container, config);
382387
383388 /* --- Solve or smooth the linear system. ---*/
@@ -401,10 +406,10 @@ void CScalarSolver<TVariable>::ImplicitEuler_Iteration(CGeometry* geometry, CSol
401406 CompleteImplicitIteration (geometry, solver_container, config);
402407}
403408
404- template <class TVariable >
409+ template <class TVariable >
405410void CScalarSolver<TVariable>::SetResidual_DualTime(CGeometry* geometry, CSolver** solver_container, CConfig* config,
406- unsigned short iRKStep, unsigned short iMesh,
407- unsigned short RunTime_EqSystem) {
411+ unsigned short iRKStep, unsigned short iMesh,
412+ unsigned short RunTime_EqSystem) {
408413 const bool implicit = (config->GetKind_TimeIntScheme () == EULER_IMPLICIT);
409414 const bool first_order = (config->GetTime_Marching () == TIME_MARCHING::DT_STEPPING_1ST);
410415 const bool second_order = (config->GetTime_Marching () == TIME_MARCHING::DT_STEPPING_2ND);
@@ -641,9 +646,9 @@ void CScalarSolver<TVariable>::SetResidual_DualTime(CGeometry* geometry, CSolver
641646 } // end dynamic grid
642647}
643648
644- template <class TVariable >
649+ template <class TVariable >
645650void CScalarSolver<TVariable>::LoadRestart(CGeometry** geometry, CSolver*** solver, CConfig* config, int val_iter,
646- bool val_update_geo) {
651+ bool val_update_geo) {
647652 /* --- Restart the solution from file information ---*/
648653
649654 unsigned short iVar, iMesh;
0 commit comments