fabs around solution in ComputeUnderRelaxation()

bmunguia · bmunguia · commit 6207f6046afc · 2020-08-03T10:02:22.000-07:00
diff --git a/SU2_CFD/include/solvers/CFVMFlowSolverBase.inl b/SU2_CFD/include/solvers/CFVMFlowSolverBase.inl
@@ -597,7 +597,7 @@ void CFVMFlowSolverBase<V, R>::ComputeUnderRelaxationFactor(CSolver** solver_con
 
       if ((iVar == 0) || (iVar == nVar - 1)) {
         const unsigned long index = iPoint * nVar + iVar;
-        su2double ratio = fabs(LinSysSol[index]) / (nodes->GetSolution(iPoint, iVar) + EPS);
+        su2double ratio = fabs(LinSysSol[index]) / (fabs(nodes->GetSolution(iPoint, iVar)) + EPS);
         if (ratio > allowableRatio) {
           localUnderRelaxation = min(allowableRatio / ratio, localUnderRelaxation);
         }
diff --git a/SU2_CFD/src/solvers/CTurbSolver.cpp b/SU2_CFD/src/solvers/CTurbSolver.cpp
@@ -674,7 +674,7 @@ void CTurbSolver::ComputeUnderRelaxationFactor(CSolver **solver_container, const
          turbulence variables can change over a nonlinear iteration. */
 
         const unsigned long index = iPoint * nVar + iVar;
-        su2double ratio = fabs(LinSysSol[index]) / (nodes->GetSolution(iPoint, iVar) + EPS);
+        su2double ratio = fabs(LinSysSol[index]) / (fabs(nodes->GetSolution(iPoint, iVar)) + EPS);
         if (ratio > allowableRatio) {
           localUnderRelaxation = min(allowableRatio / ratio, localUnderRelaxation);
         }

Original file line number	Diff line number	Diff line change
`@@ -597,7 +597,7 @@ void CFVMFlowSolverBase<V, R>::ComputeUnderRelaxationFactor(CSolver** solver_con`
`597`	`597`
`598`	`598`	`if ((iVar == 0) \|\| (iVar == nVar - 1)) {`
`599`	`599`	`const unsigned long index = iPoint * nVar + iVar;`
`600`		`- su2double ratio = fabs(LinSysSol[index]) / (nodes->GetSolution(iPoint, iVar) + EPS);`
	`600`	`+ su2double ratio = fabs(LinSysSol[index]) / (fabs(nodes->GetSolution(iPoint, iVar)) + EPS);`
`601`	`601`	`if (ratio > allowableRatio) {`
`602`	`602`	`localUnderRelaxation = min(allowableRatio / ratio, localUnderRelaxation);`
`603`	`603`	`}`
Original file line number	Diff line number	Diff line change
`@@ -674,7 +674,7 @@ void CTurbSolver::ComputeUnderRelaxationFactor(CSolver **solver_container, const`
`674`	`674`	`turbulence variables can change over a nonlinear iteration. */`
`675`	`675`
`676`	`676`	`const unsigned long index = iPoint * nVar + iVar;`
`677`		`- su2double ratio = fabs(LinSysSol[index]) / (nodes->GetSolution(iPoint, iVar) + EPS);`
	`677`	`+ su2double ratio = fabs(LinSysSol[index]) / (fabs(nodes->GetSolution(iPoint, iVar)) + EPS);`
`678`	`678`	`if (ratio > allowableRatio) {`
`679`	`679`	`localUnderRelaxation = min(allowableRatio / ratio, localUnderRelaxation);`
`680`	`680`	`}`