zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: LUscalarMatrix: incorrect debug control

Browse Source
This commit is contained in:
mattijs
2013-04-16 15:38:33 +01:00
parent dd4a5564a4
commit cf4e70612e
1 changed files with 24 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

43
src/OpenFOAM/matrices/LUscalarMatrix/LUscalarMatrix.C
View File

@ -142,38 +142,43 @@ Foam::LUscalarMatrix::LUscalarMatrix
convert(ldum, interfaceCoeffs, interfaces); convert(ldum, interfaceCoeffs, interfaces);
} }
if (debug && Pstream::master(comm_)) if (Pstream::master(comm_))
{ {
label nRows = n(); label nRows = n();
label nColumns = m(); label nColumns = m();
Pout<< "LUscalarMatrix : size:" << nRows << endl; if (debug)
for (label rowI = 0; rowI < nRows; rowI++)
{ {
const scalar* row = operator[](rowI); Pout<< "LUscalarMatrix : size:" << nRows << endl;
for (label rowI = 0; rowI < nRows; rowI++)
Pout<< "cell:" << rowI << " diagCoeff:" << row[rowI] << endl;
Pout<< " connects to upper cells :";
for (label columnI = rowI+1; columnI < nColumns; columnI++)
{ {
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<< endl; Pout<< " connects to lower cells :";
Pout<< " connects to lower cells :"; for (label columnI = 0; columnI < rowI; columnI++)
for (label columnI = 0; columnI < rowI; columnI++)
{
if (mag(row[columnI]) > SMALL)
{ {
Pout<< ' ' << columnI << " (coeff:" << row[columnI] << ")"; if (mag(row[columnI]) > SMALL)
{
Pout<< ' ' << columnI << " (coeff:" << row[columnI]
<< ")";
}
} }
Pout<< endl;
} }
Pout<< endl; Pout<< endl;
} }
Pout<< endl;
pivotIndices_.setSize(n()); pivotIndices_.setSize(n());
LUDecompose(*this, pivotIndices_); LUDecompose(*this, pivotIndices_);