particle: Added some extra debugging statements

This commit is contained in:
Will Bainbridge
2018-07-11 19:21:48 +01:00
parent 23d523d9bf
commit 95cf9666e7
2 changed files with 50 additions and 0 deletions

View File

@ -152,6 +152,11 @@ void Foam::particle::rotate(const bool reverse)
void Foam::particle::changeTet(const label tetTriI)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
const bool isOwner = mesh_.faceOwner()[tetFacei_] == celli_;
const label firstTetPtI = 1;
@ -227,6 +232,11 @@ void Foam::particle::changeTet(const label tetTriI)
void Foam::particle::changeFace(const label tetTriI)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
// Get the old topology
const triFace triOldIs(currentTetIndices().faceTriIs(mesh_));
@ -342,6 +352,11 @@ void Foam::particle::changeFace(const label tetTriI)
void Foam::particle::changeCell()
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
// Set the cell to be the one on the other side of the face
const label ownerCellI = mesh_.faceOwner()[tetFacei_];
const bool isOwner = celli_ == ownerCellI;
@ -354,6 +369,11 @@ void Foam::particle::changeCell()
void Foam::particle::changeToMasterPatch()
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
label thisPatch = patch();
forAll(mesh_.cells()[celli_], cellFaceI)
@ -395,6 +415,11 @@ void Foam::particle::locate
const string boundaryMsg
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
// Find the cell, if it has not been given
if (celli < 0)
{
@ -591,6 +616,11 @@ Foam::scalar Foam::particle::track
const scalar fraction
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
scalar f = trackToFace(displacement, fraction);
while (onInternalFace())
@ -610,6 +640,11 @@ Foam::scalar Foam::particle::trackToCell
const scalar fraction
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
const scalar f = trackToFace(displacement, fraction);
if (onInternalFace())
@ -627,6 +662,11 @@ Foam::scalar Foam::particle::trackToFace
const scalar fraction
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
scalar f = 1;
label tetTriI = onFace() ? 0 : -1;

View File

@ -108,6 +108,11 @@ void Foam::particle::hitFace
trackingData& td
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
typename TrackCloudType::particleType& p =
static_cast<typename TrackCloudType::particleType&>(*this);
typename TrackCloudType::particleType::trackingData& ttd =
@ -183,6 +188,11 @@ void Foam::particle::trackToAndHitFace
trackingData& td
)
{
if (debug)
{
Info << "Particle " << origId() << nl << FUNCTION_NAME << nl << endl;
}
trackToFace(direction, fraction);
hitFace(direction, cloud, td);