diff --git a/src/OpenFOAM/matrices/LUscalarMatrix/LUscalarMatrix.C b/src/OpenFOAM/matrices/LUscalarMatrix/LUscalarMatrix.C index ceb0661489..76cb4f9576 100644 --- a/src/OpenFOAM/matrices/LUscalarMatrix/LUscalarMatrix.C +++ b/src/OpenFOAM/matrices/LUscalarMatrix/LUscalarMatrix.C @@ -142,38 +142,43 @@ Foam::LUscalarMatrix::LUscalarMatrix convert(ldum, interfaceCoeffs, interfaces); } - if (debug && Pstream::master(comm_)) + if (Pstream::master(comm_)) { label nRows = n(); label nColumns = m(); - Pout<< "LUscalarMatrix : size:" << nRows << endl; - for (label rowI = 0; rowI < nRows; rowI++) + if (debug) { - const scalar* row = operator[](rowI); - - Pout<< "cell:" << rowI << " diagCoeff:" << row[rowI] << endl; - - Pout<< " connects to upper cells :"; - for (label columnI = rowI+1; columnI < nColumns; columnI++) + Pout<< "LUscalarMatrix : size:" << nRows << endl; + for (label rowI = 0; rowI < nRows; rowI++) { - if (mag(row[columnI]) > SMALL) + const scalar* row = operator[](rowI); + + Pout<< "cell:" << rowI << " diagCoeff:" << row[rowI] << endl; + + Pout<< " connects to upper cells :"; + for (label columnI = rowI+1; columnI < nColumns; columnI++) { - Pout<< ' ' << columnI << " (coeff:" << row[columnI] << ")"; + if (mag(row[columnI]) > SMALL) + { + Pout<< ' ' << columnI << " (coeff:" << row[columnI] + << ")"; + } } - } - Pout<< endl; - Pout<< " connects to lower cells :"; - for (label columnI = 0; columnI < rowI; columnI++) - { - if (mag(row[columnI]) > SMALL) + Pout<< endl; + Pout<< " connects to lower cells :"; + for (label columnI = 0; columnI < rowI; columnI++) { - Pout<< ' ' << columnI << " (coeff:" << row[columnI] << ")"; + if (mag(row[columnI]) > SMALL) + { + Pout<< ' ' << columnI << " (coeff:" << row[columnI] + << ")"; + } } + Pout<< endl; } Pout<< endl; } - Pout<< endl; pivotIndices_.setSize(n()); LUDecompose(*this, pivotIndices_);