diff --git a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C index 33922b5019..fd78acb546 100644 --- a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C +++ b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C @@ -587,7 +587,41 @@ void Foam::fvMatrix::relax(const scalar alpha) } } + + if (debug) + { + // Calculate amount of non-dominance. + label nNon = 0; + scalar maxNon = 0.0; + scalar sumNon = 0.0; + forAll(D, celli) + { + scalar d = (sumOff[celli] - D[celli])/D[celli]; + + if (d > 0) + { + nNon++; + maxNon = max(maxNon, d); + sumNon += d; + } + } + + reduce(nNon, sumOp