63 const scalar maxNonOrtho
71 const scalar maxConcave
79 const scalar maxIntSkew
83 const scalar maxBounSkew
87 const scalar minWeight
91 const scalar minVolRatio
99 const scalar minTriangleTwist
108 label nWrongFaces = 0;
110 Info<<
"Checking faces in error :" <<
endl;
113 if (maxNonOrtho < 180.0-SMALL)
115 polyMeshGeometry::checkFaceDotProduct
129 Info<<
" non-orthogonality > "
130 <<
setw(3) << maxNonOrtho
132 << nNewWrongFaces-nWrongFaces <<
endl;
134 nWrongFaces = nNewWrongFaces;
139 polyMeshGeometry::checkFacePyramids
153 Info<<
" faces with face pyramid volume < "
154 <<
setw(5) << minVol <<
" : "
155 << nNewWrongFaces-nWrongFaces <<
endl;
157 nWrongFaces = nNewWrongFaces;
160 if (maxConcave < 180.0-SMALL)
162 polyMeshGeometry::checkFaceAngles
175 Info<<
" faces with concavity > "
176 <<
setw(3) << maxConcave
178 << nNewWrongFaces-nWrongFaces <<
endl;
180 nWrongFaces = nNewWrongFaces;
183 if (minArea > -SMALL)
185 polyMeshGeometry::checkFaceArea
197 Info<<
" faces with area < "
198 <<
setw(5) << minArea
200 << nNewWrongFaces-nWrongFaces <<
endl;
202 nWrongFaces = nNewWrongFaces;
205 if (maxIntSkew > 0 || maxBounSkew > 0)
207 polyMeshGeometry::checkFaceSkewness
223 Info<<
" faces with skewness > "
224 <<
setw(3) << maxIntSkew
225 <<
" (internal) or " <<
setw(3) << maxBounSkew
226 <<
" (boundary) : " << nNewWrongFaces-nWrongFaces <<
endl;
228 nWrongFaces = nNewWrongFaces;
231 if (minWeight >= 0 && minWeight < 1)
233 polyMeshGeometry::checkFaceWeights
248 Info<<
" faces with interpolation weights (0..1) < "
249 <<
setw(5) << minWeight
251 << nNewWrongFaces-nWrongFaces <<
endl;
253 nWrongFaces = nNewWrongFaces;
256 if (minVolRatio >= 0)
258 polyMeshGeometry::checkVolRatio
271 Info<<
" faces with volume ratio of neighbour cells < "
272 <<
setw(5) << minVolRatio
274 << nNewWrongFaces-nWrongFaces <<
endl;
276 nWrongFaces = nNewWrongFaces;
283 polyMeshGeometry::checkFaceTwist
298 Info<<
" faces with face twist < "
299 <<
setw(5) << minTwist
301 << nNewWrongFaces-nWrongFaces <<
endl;
303 nWrongFaces = nNewWrongFaces;
306 if (minTriangleTwist > -1)
310 polyMeshGeometry::checkTriangleTwist
324 Info<<
" faces with triangle twist < "
325 <<
setw(5) << minTriangleTwist
327 << nNewWrongFaces-nWrongFaces <<
endl;
329 nWrongFaces = nNewWrongFaces;
334 polyMeshGeometry::checkCellDeterminant
341 polyMeshGeometry::affectedCells(mesh, checkFaces),
347 Info<<
" faces on cells with determinant < "
348 <<
setw(5) << minDet <<
" : "
349 << nNewWrongFaces-nWrongFaces <<
endl;
351 nWrongFaces = nNewWrongFaces;
356 return nWrongFaces > 0;
411 const scalar maxNonOrtho
419 const scalar maxConcave
427 const scalar maxIntSkew
431 const scalar maxBounSkew
435 const scalar minWeight
439 const scalar minVolRatio
443 const scalar minTwist
447 const scalar minTriangleTwist
456 label nWrongFaces = 0;
458 Info<<
"Checking faces in error :" <<
endl;
461 if (maxNonOrtho < 180.0-SMALL)
474 Info<<
" non-orthogonality > "
475 <<
setw(3) << maxNonOrtho
477 << nNewWrongFaces-nWrongFaces <<
endl;
479 nWrongFaces = nNewWrongFaces;
496 Info<<
" faces with face pyramid volume < "
497 <<
setw(5) << minVol <<
" : "
498 << nNewWrongFaces-nWrongFaces <<
endl;
500 nWrongFaces = nNewWrongFaces;
503 if (maxConcave < 180.0-SMALL)
516 Info<<
" faces with concavity > "
517 <<
setw(3) << maxConcave
519 << nNewWrongFaces-nWrongFaces <<
endl;
521 nWrongFaces = nNewWrongFaces;
524 if (minArea > -SMALL)
526 meshGeom.
checkFaceArea(report, minArea, checkFaces, &wrongFaces);
530 Info<<
" faces with area < "
531 <<
setw(5) << minArea
533 << nNewWrongFaces-nWrongFaces <<
endl;
535 nWrongFaces = nNewWrongFaces;
538 if (maxIntSkew > 0 || maxBounSkew > 0)
552 Info<<
" faces with skewness > "
553 <<
setw(3) << maxIntSkew
554 <<
" (internal) or " <<
setw(3) << maxBounSkew
555 <<
" (boundary) : " << nNewWrongFaces-nWrongFaces <<
endl;
557 nWrongFaces = nNewWrongFaces;
560 if (minWeight >= 0 && minWeight < 1)
573 Info<<
" faces with interpolation weights (0..1) < "
574 <<
setw(5) << minWeight
576 << nNewWrongFaces-nWrongFaces <<
endl;
578 nWrongFaces = nNewWrongFaces;
581 if (minVolRatio >= 0)
594 Info<<
" faces with volume ratio of neighbour cells < "
595 <<
setw(5) << minVolRatio
597 << nNewWrongFaces-nWrongFaces <<
endl;
599 nWrongFaces = nNewWrongFaces;
617 Info<<
" faces with face twist < "
618 <<
setw(5) << minTwist
620 << nNewWrongFaces-nWrongFaces <<
endl;
622 nWrongFaces = nNewWrongFaces;
625 if (minTriangleTwist > -1)
640 Info<<
" faces with triangle twist < "
641 <<
setw(5) << minTriangleTwist
643 << nNewWrongFaces-nWrongFaces <<
endl;
645 nWrongFaces = nNewWrongFaces;
661 Info<<
" faces on cells with determinant < "
662 <<
setw(5) << minDet <<
" : "
663 << nNewWrongFaces-nWrongFaces <<
endl;
665 nWrongFaces = nNewWrongFaces;
670 return nWrongFaces > 0;