Skip to content
Snippets Groups Projects
Commit eb1080c9 authored by Henry Weller's avatar Henry Weller
Browse files

checkMesh: Provide the number of geometric and solution directions.

parent 79f20242
No related branches found
No related tags found
No related merge requests found
...@@ -493,12 +493,15 @@ Foam::label Foam::checkGeometry(const polyMesh& mesh, const bool allGeometry) ...@@ -493,12 +493,15 @@ Foam::label Foam::checkGeometry(const polyMesh& mesh, const bool allGeometry)
// Min length // Min length
scalar minDistSqr = magSqr(1e-6 * globalBb.span()); scalar minDistSqr = magSqr(1e-6 * globalBb.span());
// Non-empty directions // Geometric directions
const Vector<label> validDirs = (mesh.geometricD() + Vector<label>::one)/2; const Vector<label> validDirs = (mesh.geometricD() + Vector<label>::one)/2;
Info<< " Mesh (non-empty, non-wedge) directions " << validDirs << endl; Info<< " Mesh has " << mesh.nGeometricD()
<< " geometric (non-empty/wedge) directions " << validDirs << endl;
// Solution directions
const Vector<label> solDirs = (mesh.solutionD() + Vector<label>::one)/2; const Vector<label> solDirs = (mesh.solutionD() + Vector<label>::one)/2;
Info<< " Mesh (non-empty) directions " << solDirs << endl; Info<< " Mesh has " << mesh.nSolutionD()
<< " solution (non-empty) directions " << solDirs << endl;
if (mesh.nGeometricD() < 3) if (mesh.nGeometricD() < 3)
{ {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment