32 void Foam::blockMesh::check(
const polyMesh& bm,
const dictionary&
dict)
const
41 for (label cej=cei+1; cej<edges_.size(); cej++)
43 if (edges_[cei].compare(edges_[cej]) != 0)
45 Info<<
" Curved edge ";
47 Info<<
" is a duplicate of curved edge " << edges_[cei]
80 Info<<
" Curved edge ";
82 Info<<
" does not correspond to a block edge." <<
endl;
92 for (label cfj=cfi+1; cfj<faces_.size(); cfj++)
94 if (faces_[cfi].compare(faces_[cfj]) != 0)
96 Info<<
" Curved face ";
98 Info<<
" is a duplicate of curved face ";
120 Info<<
" Curved face ";
122 Info<<
" does not correspond to a block face." <<
endl;
131 label nBoundaryFaces = 0;
134 nBoundaryFaces +=
cells[celli].nFaces();
137 nBoundaryFaces -= 2*bm.nInternalFaces();
139 label nDefinedBoundaryFaces = 0;
142 nDefinedBoundaryFaces +=
patches[patchi].size();
149 <<
tab <<
tab <<
"Number of internal faces : "
150 << bm.nInternalFaces() <<
nl
151 <<
tab <<
tab <<
"Number of boundary faces : "
152 << nBoundaryFaces <<
nl
153 <<
tab <<
tab <<
"Number of defined boundary faces : "
154 << nDefinedBoundaryFaces <<
nl
155 <<
tab <<
tab <<
"Number of undefined boundary faces : "
156 << nBoundaryFaces - nDefinedBoundaryFaces <<
nl;
158 if ((nBoundaryFaces - nDefinedBoundaryFaces) > 0)
161 <<
"(Warning : only leave undefined the front and back planes "
162 <<
"of 2D planar geometries!)" <<
endl;
165 Info<<
tab <<
"Checking patch -> block consistency" <<
endl;
175 const face& patchFace = Patch[patchFacei];
176 bool patchFaceOK =
false;
182 forAll(cellFaces, cellFacei)
184 if (patchFace ==
faces[cellFaces[cellFacei]])
191 patchFace.areaNormal(
points)
197 <<
"Face " << patchFacei
198 <<
" of patch " << patchi
199 <<
" (" <<
patches[patchi].name() <<
")"
212 <<
"Face " << patchFacei
213 <<
" of patch " << patchi
214 <<
" (" <<
patches[patchi].name() <<
")"
215 <<
" does not match any block faces" <<
endl;
230 <<
"Block mesh topology incorrect, stopping mesh generation!"