research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

viewFactor: moved some of the "noise" to print only on debug

Browse source
This commit is contained in:
Henry 2014-12-29 11:11:43 +00:00
parent 9550a248cf
commit 5f6e909f72
1 changed files with 34 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

43
src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C
View file

@ -69,16 +69,19 @@ void Foam::radiation::viewFactor::initialise()
if (debug)
{
Pout<< "Selected patches:" << selectedPatches_ << endl;
Pout<< "Number of coarse faces:" << nLocalCoarseFaces_ << endl;
Pout<< "radiation::viewFactor::initialise() Selected patches:"
<< selectedPatches_ << endl;
Pout<< "radiation::viewFactor::initialise() Number of coarse faces:"
<< nLocalCoarseFaces_ << endl;
}
totalNCoarseFaces_ = nLocalCoarseFaces_;
reduce(totalNCoarseFaces_, sumOp<label>());
if (Pstream::master())
if (debug && Pstream::master())
{
Info<< "Total number of clusters : " << totalNCoarseFaces_ << endl;
InfoIn("radiation::viewFactor::initialise()")
<< "Total number of clusters : " << totalNCoarseFaces_ << endl;
}
labelListIOList subMap
@ -173,7 +176,11 @@ void Foam::radiation::viewFactor::initialise()
new scalarSquareMatrix(totalNCoarseFaces_, totalNCoarseFaces_, 0.0)
);
Info<< "Insert elements in the matrix..." << endl;
if (debug)
{
InfoIn("radiation::viewFactor::initialise()")
<< "Insert elements in the matrix..." << endl;
}
for (label procI = 0; procI < Pstream::nProcs(); procI++)
{
@ -191,7 +198,11 @@ void Foam::radiation::viewFactor::initialise()
bool smoothing = readBool(coeffs_.lookup("smoothing"));
if (smoothing)
{
Info<< "Smoothing the matrix..." << endl;
if (debug)
{
InfoIn("radiation::viewFactor::initialise()")
<< "Smoothing the matrix..." << endl;
}
for (label i=0; i<totalNCoarseFaces_; i++)
{
@ -543,6 +554,7 @@ void Foam::radiation::viewFactor::calculate()
}
Info<< "\nSolving view factor equations..." << endl;
// Negative coming into the fluid
LUsolve(C, q);
}
@ -566,7 +578,13 @@ void Foam::radiation::viewFactor::calculate()
}
}
}
Info<< "\nDecomposing C matrix..." << endl;
if (debug)
{
InfoIn("radiation::viewFactor::initialise()")
<< "\nDecomposing C matrix..." << endl;
}
LUDecompose(CLU_(), pivotIndices_);
}
@ -589,7 +607,12 @@ void Foam::radiation::viewFactor::calculate()
}
}
Info<< "\nLU Back substitute C matrix.." << endl;
if (debug)
{
InfoIn("radiation::viewFactor::initialise()")
<< "\nLU Back substitute C matrix.." << endl;
}
LUBacksubstitute(CLU_(), pivotIndices_, q);
iterCounter_ ++;
}
@ -643,7 +666,9 @@ void Foam::radiation::viewFactor::calculate()
const scalarField& Qrp = Qr_.boundaryField()[patchID];
const scalarField& magSf = mesh_.magSf().boundaryField()[patchID];
scalar heatFlux = gSum(Qrp*magSf);
Info<< "Total heat transfer rate at patch: "
InfoIn("radiation::viewFactor::initialise()")
<< "Total heat transfer rate at patch: "
<< patchID << " "
<< heatFlux << endl;
}