diff --git a/src/OpenFOAM/matrices/lduMatrix/solvers/PBiCGStab/PBiCGStab.C b/src/OpenFOAM/matrices/lduMatrix/solvers/PBiCGStab/PBiCGStab.C index b7c1098afe88959ad833636704ad4b8b4df7b20b..e68a20760449020f72794a8d4b4f3f067f5e11ba 100644 --- a/src/OpenFOAM/matrices/lduMatrix/solvers/PBiCGStab/PBiCGStab.C +++ b/src/OpenFOAM/matrices/lduMatrix/solvers/PBiCGStab/PBiCGStab.C @@ -206,7 +206,11 @@ Foam::solverPerformance Foam::PBiCGStab::solve solverPerf.finalResidual() = gSumMag(sA, matrix().mesh().comm())/normFactor; - if (solverPerf.checkConvergence(tolerance_, relTol_)) + if + ( + solverPerf.nIterations() >= minIter_ + && solverPerf.checkConvergence(tolerance_, relTol_) + ) { for (label cell=0; cell<nCells; cell++) {