Further shared reading optimizations.

Author: jblueh

SU2_CFD/include/solvers/CFVMFlowSolverBase.inl

@@ -1621,6 +1621,8 @@ void CFVMFlowSolverBase<V, FlowRegime>::SetResidual_DualTime(CGeometry *geometry
 
     /*--- Loop over all nodes (excluding halos) ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -1656,6 +1658,8 @@ void CFVMFlowSolverBase<V, FlowRegime>::SetResidual_DualTime(CGeometry *geometry
     }
     END_SU2_OMP_FOR
 
+    AD::EndNoSharedReading();
+
   }
 
   else {
@@ -1733,6 +1737,8 @@ void CFVMFlowSolverBase<V, FlowRegime>::SetResidual_DualTime(CGeometry *geometry
     /*--- Loop over all nodes (excluding halos) to compute the remainder
      of the dual time-stepping source term. ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -1770,6 +1776,8 @@ void CFVMFlowSolverBase<V, FlowRegime>::SetResidual_DualTime(CGeometry *geometry
       }
     }
     END_SU2_OMP_FOR
+
+    AD::EndNoSharedReading();
   }
 
 }

SU2_CFD/src/integration/CMultiGridIntegration.cpp

@@ -547,10 +547,12 @@ void CMultiGridIntegration::SetForcing_Term(CSolver *sol_fine, CSolver *sol_coar
 
 void CMultiGridIntegration::SetResidual_Term(CGeometry *geometry, CSolver *solver) {
 
+  AD::StartNoSharedReading();
   SU2_OMP_FOR_STAT(roundUpDiv(geometry->GetnPointDomain(), omp_get_num_threads()))
   for (unsigned long iPoint = 0; iPoint < geometry->GetnPointDomain(); iPoint++)
     solver->LinSysRes.AddBlock(iPoint, solver->GetNodes()->GetResTruncError(iPoint));
   END_SU2_OMP_FOR
+  AD::EndNoSharedReading();
 
 }
 

SU2_CFD/src/solvers/CIncEulerSolver.cpp

@@ -2555,6 +2555,8 @@ void CIncEulerSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver
 
     /*--- Loop over all nodes (excluding halos) ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -2607,6 +2609,8 @@ void CIncEulerSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver
       }
     }
     END_SU2_OMP_FOR
+
+    AD::EndNoSharedReading();
   }
 
   else {
@@ -2694,6 +2698,8 @@ void CIncEulerSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver
     /*--- Loop over all nodes (excluding halos) to compute the remainder
      of the dual time-stepping source term. ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -2749,6 +2755,8 @@ void CIncEulerSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver
       }
     }
     END_SU2_OMP_FOR
+
+    AD::EndNoSharedReading();
   }
 
 }

SU2_CFD/src/solvers/CSolver.cpp

@@ -1929,6 +1929,8 @@ void CSolver::SetResidual_RMS(const CGeometry *geometry, const CConfig *config)
 
   if (geometry->GetMGLevel() != MESH_0) return;
 
+  AD::StartNoSharedReading();
+
   SU2_OMP_MASTER {
 
   /*--- Set the L2 Norm residual in all the processors. ---*/
@@ -1985,12 +1987,16 @@ void CSolver::SetResidual_RMS(const CGeometry *geometry, const CConfig *config)
   }
   END_SU2_OMP_MASTER
   SU2_OMP_BARRIER
+
+  AD::EndNoSharedReading();
 }
 
 void CSolver::SetResidual_BGS(const CGeometry *geometry, const CConfig *config) {
 
   if (geometry->GetMGLevel() != MESH_0) return;
 
+  AD::StartNoSharedReading();
+
   SU2_OMP_MASTER {
 
   /*--- Set the L2 Norm residual in all the processors. ---*/
@@ -2028,6 +2034,8 @@ void CSolver::SetResidual_BGS(const CGeometry *geometry, const CConfig *config)
   }
   END_SU2_OMP_MASTER
   SU2_OMP_BARRIER
+
+  AD::EndNoSharedReading();
 }
 
 void CSolver::SetRotatingFrame_GCL(CGeometry *geometry, const CConfig *config) {
@@ -4111,6 +4119,8 @@ void CSolver::ComputeResidual_Multizone(const CGeometry *geometry, const CConfig
 
   SU2_OMP_PARALLEL {
 
+  AD::StartNoSharedReading();
+
   /*--- Set Residuals to zero ---*/
   SU2_OMP_MASTER
   for (unsigned short iVar = 0; iVar < nVar; iVar++){
@@ -4152,6 +4162,8 @@ void CSolver::ComputeResidual_Multizone(const CGeometry *geometry, const CConfig
   END_SU2_OMP_CRITICAL
   SU2_OMP_BARRIER
 
+  AD::EndNoSharedReading();
+
   SetResidual_BGS(geometry, config);
 
   }

SU2_CFD/src/solvers/CTurbSolver.cpp

@@ -789,6 +789,8 @@ void CTurbSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver_con
 
     /*--- Loop over all nodes (excluding halos) ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -855,6 +857,8 @@ void CTurbSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver_con
     }
     END_SU2_OMP_FOR
 
+    AD::EndNoSharedReading();
+
   } else {
 
     /*--- For unsteady flows on dynamic meshes (rigidly transforming or
@@ -955,6 +959,8 @@ void CTurbSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver_con
     /*--- Loop over all nodes (excluding halos) to compute the remainder
      of the dual time-stepping source term. ---*/
 
+    AD::StartNoSharedReading();
+
     SU2_OMP_FOR_STAT(omp_chunk_size)
     for (iPoint = 0; iPoint < nPointDomain; iPoint++) {
 
@@ -1023,6 +1029,8 @@ void CTurbSolver::SetResidual_DualTime(CGeometry *geometry, CSolver **solver_con
     }
     END_SU2_OMP_FOR
 
+    AD::EndNoSharedReading();
+
   } // end dynamic grid
 
 }