spaces and comments

Author: FlorianDm

SU2_CFD/include/numerics/turbulent/turb_sources.hpp

@@ -336,39 +336,39 @@ class CSourcePieceWise_TurbSST final : public CNumerics {
   inline void ResidualAxisymmetric(su2double alfa_blended, su2double zeta){
 
     if (Coord_i[1] < EPS) return;
-     
+
     su2double yinv, rhov, k, w;
     su2double sigma_k_i, sigma_w_i;
     su2double pk_axi, pw_axi, cdk_axi, cdw_axi;
 
     AD::SetPreaccIn(Coord_i[1]);
-    
+
     yinv = 1.0/Coord_i[1];
     rhov = Density_i*V_i[2];
     k = TurbVar_i[0];
     w = TurbVar_i[1];
-    
+
     /*--- Compute blended constants ---*/
     sigma_k_i = F1_i*sigma_k_1+(1.0-F1_i)*sigma_k_2;
     sigma_w_i = F1_i*sigma_w_1+(1.0-F1_i)*sigma_w_2;
 
     /*--- Production ---*/
     pk_axi = max(0.0,2.0/3.0*rhov*k*(2.0/zeta*(yinv*V_i[2]-PrimVar_Grad_i[2][1]-PrimVar_Grad_i[1][0])-1.0));
     pw_axi = alfa_blended*zeta/k*pk_axi;
-    
+
     /*--- Convection-Diffusion ---*/
     cdk_axi = rhov*k-(Laminar_Viscosity_i+sigma_k_i*Eddy_Viscosity_i)*TurbVar_Grad_i[0][1];
     cdw_axi = rhov*w-(Laminar_Viscosity_i+sigma_w_i*Eddy_Viscosity_i)*TurbVar_Grad_i[1][1];
-    
+
     /*--- Add terms to the residuals ---*/
     Residual[0] += yinv*Volume*(pk_axi-cdk_axi);
     Residual[1] += yinv*Volume*(pw_axi-cdw_axi);
 
-    /*--- Add contribution to the jacobian for implicit time integration--- (ignore to conserve diagonal dominance)*/ 
-    //Jacobian_i[0][0] += yinv*Volume*(sigma_k_i*TurbVar_Grad_i[0][1]/zeta-V_i[2]);
-    //Jacobian_i[0][1] += yinv*Volume*(-sigma_k_i*k*TurbVar_Grad_i[0][1]/(zeta*zeta));
-    //Jacobian_i[1][0] += yinv*Volume*(sigma_w_i*TurbVar_Grad_i[1][1]/zeta);
-    //Jacobian_i[1][1] += yinv*Volume*(-sigma_w_i*k*TurbVar_Grad_i[1][1]/(zeta*zeta)-V_i[2]);
+    /*--- Add contribution to the jacobian for implicit time integration (ignore to conserve diagonal dominance) ---*/
+    Jacobian_i[0][0] += yinv*Volume*(sigma_k_i*TurbVar_Grad_i[0][1]/zeta-V_i[2]);
+    Jacobian_i[0][1] += yinv*Volume*(-sigma_k_i*k*TurbVar_Grad_i[0][1]/(zeta*zeta));
+    Jacobian_i[1][0] += yinv*Volume*(sigma_w_i*TurbVar_Grad_i[1][1]/zeta);
+    Jacobian_i[1][1] += yinv*Volume*(-sigma_w_i*k*TurbVar_Grad_i[1][1]/(zeta*zeta)-V_i[2]);
 
   }