42 Info<<
"bool polyMeshZipUpCells(polyMesh& mesh) const: "
43 <<
"zipping up topologically open cells" <<
endl;
64 label nChangedFacesInMesh = 0;
71 nChangedFacesInMesh = 0;
84 const edgeList cellEdges = Cells[celli].edges(oldFaces);
85 const labelList cellPoints = Cells[celli].labels(oldFaces);
93 edgeList curFaceEdges = oldFaces[curFaces[facei]].edges();
95 forAll(curFaceEdges, faceEdgeI)
97 const edge& curEdge = curFaceEdges[faceEdgeI];
99 forAll(cellEdges, cellEdgeI)
101 if (cellEdges[cellEdgeI] == curEdge)
103 edgeUsage[cellEdgeI]++;
110 edgeList singleEdges(cellEdges.size());
111 label nSingleEdges = 0;
115 if (edgeUsage[edgeI] == 1)
117 singleEdges[nSingleEdges] = cellEdges[edgeI];
120 else if (edgeUsage[edgeI] != 2)
123 <<
"edge " << cellEdges[edgeI] <<
" in cell " << celli
124 <<
" used " << edgeUsage[edgeI] <<
" times. " <<
nl
125 <<
"Should be 1 or 2 - serious error "
126 <<
"in mesh structure. " <<
endl;
131 Info<<
"face: " << oldFaces[curFaces[facei]]
135 Info<<
"Cell edges: " << cellEdges <<
nl
136 <<
"Edge usage: " << edgeUsage <<
nl
137 <<
"Cell points: " << cellPoints <<
endl;
141 Info<<
"vertex create \"" << cellPoints[cpI]
148 problemCells.
insert(celli);
153 if (nSingleEdges == 0)
continue;
155 singleEdges.setSize(nSingleEdges);
162 Info<<
"face: " << oldFaces[curFaces[facei]] <<
endl;
165 Info<<
"Cell edges: " << cellEdges <<
nl
166 <<
"Edge usage: " << edgeUsage <<
nl
167 <<
"Single edges: " << singleEdges <<
nl
168 <<
"Cell points: " << cellPoints <<
endl;
172 Info<<
"vertex create \"" << cellPoints[cpI]
184 forAll(singleEdges, edgeI)
186 const edge& curEdge = singleEdges[edgeI];
188 forAll(cellPoints, pointi)
192 cellPoints[pointi] == curEdge.
start()
193 || cellPoints[pointi] == curEdge.
end()
196 pointUsage[pointi]++;
201 boolList singleEdgeUsage(singleEdges.size(),
false);
205 forAll(singleEdges, edgeI)
207 bool blockedHead =
false;
208 bool blockedTail =
false;
210 label newEdgeStart = singleEdges[edgeI].start();
211 label newEdgeEnd = singleEdges[edgeI].end();
214 forAll(cellPoints, pointi)
216 if (cellPoints[pointi] == newEdgeStart)
218 if (pointUsage[pointi] > 2)
223 else if (cellPoints[pointi] == newEdgeEnd)
225 if (pointUsage[pointi] > 2)
232 if (blockedHead && blockedTail)
235 singleEdgeUsage[edgeI] =
true;
244 label nEdgesToInsert = 0;
247 forAll(singleEdges, edgeI)
251 bool blockHead =
false;
252 bool blockTail =
false;
254 if (!singleEdgeUsage[edgeI])
257 singleEdgeUsage[edgeI] =
true;
259 label newEdgeStart = singleEdges[edgeI].start();
260 label newEdgeEnd = singleEdges[edgeI].
end();
262 pointChain.
insert(newEdgeStart);
263 pointChain.
append(newEdgeEnd);
266 Info<<
"found edge to start with: "
267 << singleEdges[edgeI] <<
endl;
271 forAll(cellPoints, pointi)
273 if (cellPoints[pointi] == newEdgeStart)
275 if (pointUsage[pointi] > 2)
278 Info<<
"start head blocked" <<
endl;
284 else if (cellPoints[pointi] == newEdgeEnd)
286 if (pointUsage[pointi] > 2)
289 Info<<
"start tail blocked" <<
endl;
297 bool stopSearching =
false;
302 stopSearching =
false;
304 forAll(singleEdges, addEdgeI)
306 if (!singleEdgeUsage[addEdgeI])
310 singleEdges[addEdgeI].start();
313 singleEdges[addEdgeI].end();
316 Info<<
"Trying candidate "
317 << singleEdges[addEdgeI] <<
endl;
323 if (pointChain.
first() == addStart)
326 pointChain.
insert(addEnd);
328 singleEdgeUsage[addEdgeI] =
true;
330 else if (pointChain.
first() == addEnd)
332 pointChain.
insert(addStart);
334 singleEdgeUsage[addEdgeI] =
true;
340 if (!blockTail && !singleEdgeUsage[addEdgeI])
342 if (pointChain.
last() == addStart)
345 pointChain.
append(addEnd);
347 singleEdgeUsage[addEdgeI] =
true;
349 else if (pointChain.
last() == addEnd)
351 pointChain.
append(addStart);
353 singleEdgeUsage[addEdgeI] =
true;
358 label curEdgeStart = pointChain.
first();
359 label curEdgeEnd = pointChain.
last();
362 Info<<
"curEdgeStart: " << curEdgeStart
363 <<
" curEdgeEnd: " << curEdgeEnd <<
endl;
366 forAll(cellPoints, pointi)
368 if (cellPoints[pointi] == curEdgeStart)
370 if (pointUsage[pointi] > 2)
379 else if (cellPoints[pointi] == curEdgeEnd)
381 if (pointUsage[pointi] > 2)
393 if (curEdgeStart == curEdgeEnd)
404 stopSearching =
true;
408 Info<<
"current pointChain: " << pointChain
412 if (stopSearching)
break;
415 }
while (stopSearching);
419 Info<<
"completed patch chain: " << pointChain <<
endl;
422 if (pointChain.size() > 2)
424 edgesToInsert[nEdgesToInsert] = pointChain;
429 edgesToInsert.setSize(nEdgesToInsert);
432 Info<<
"edgesToInsert: " << edgesToInsert <<
endl;
436 forAll(edgesToInsert, edgeToInsertI)
444 const labelList& unorderedEdge = edgesToInsert[edgeToInsertI];
454 for (label i = 1; i < dist.size(); i++)
456 dist[i] = (
Points[unorderedEdge[i]] - startPoint) & dir;
460 labelList orderedEdge(unorderedEdge.size(), -1);
465 label nextPoint = -1;
466 scalar minDist = GREAT;
470 if (!
used[i] && dist[i] < minDist)
478 orderedEdge[epI] = unorderedEdge[nextPoint];
479 used[nextPoint] =
true;
483 Info<<
"unorderedEdge: " << unorderedEdge <<
nl
484 <<
"orderedEdge: " << orderedEdge <<
endl;
488 forAll(orderedEdge, checkI)
492 label checkJ = checkI + 1;
493 checkJ < orderedEdge.size();
497 if (orderedEdge[checkI] == orderedEdge[checkJ])
500 <<
"Duplicate point found in edge to insert. "
501 <<
nl <<
"Point: " << orderedEdge[checkI]
502 <<
" edge: " << orderedEdge <<
endl;
504 problemCells.
insert(celli);
520 labelList facesSharingEdge(startPF.size() + endPF.size());
525 facesSharingEdge[nfse++] = startPF[pfI];
529 facesSharingEdge[nfse++] = endPF[pfI];
532 forAll(facesSharingEdge, facei)
534 bool faceChanges =
false;
537 const label currentFaceIndex = facesSharingEdge[facei];
540 oldFaces[currentFaceIndex].edges();
542 forAll(curFaceEdges, cfeI)
544 if (curFaceEdges[cfeI] == testEdge)
553 nChangedFacesInMesh++;
566 face& newFace = newFaces[currentFaceIndex];
568 bool allPointsPresent =
true;
572 bool curPointFound =
false;
576 if (newFace[nfI] == orderedEdge[oeI])
578 curPointFound =
true;
584 allPointsPresent && curPointFound;
588 if (allPointsPresent)
590 Info<<
"All points present" <<
endl;
594 if (!allPointsPresent)
607 Info<<
"Not all points present." <<
endl;
610 label nNewFacePoints = 0;
612 bool edgeAdded =
false;
614 forAll(newFaceEdges, curFacEdgI)
617 if (newFaceEdges[curFacEdgI] == testEdge)
627 + orderedEdge.size() - 2
632 newFaceEdges[curFacEdgI].start()
640 i < orderedEdge.size() - 1;
644 newFace[nNewFacePoints] =
654 label i = orderedEdge.size() - 1;
659 newFace[nNewFacePoints] =
669 newFace[nNewFacePoints] =
670 newFaceEdges[curFacEdgI].start();
677 << oldFaces[currentFaceIndex] <<
nl
678 <<
"newFace: " << newFace <<
endl;
686 label checkJ = checkI + 1;
687 checkJ < newFace.size();
691 if (newFace[checkI] == newFace[checkJ])
694 <<
"Duplicate point found "
695 <<
"in the new face. " <<
nl
697 << orderedEdge[checkI]
701 problemCells.
insert(celli);
712 Info<<
"This edge modifies an already modified "
713 <<
"edge. Point insertions skipped."
722 if (problemCells.size())
726 <<
"Found " << problemCells.size() <<
" problem cells." <<
nl
727 <<
"Cells: " << problemCells.sortedToc()
731 Info<<
"Cycle " << ++nCycles
732 <<
" changed " << nChangedFacesInMesh <<
" faces." <<
endl;
747 patchSizes[patchi] =
bMesh[patchi].size();
748 patchStarts[patchi] =
bMesh[patchi].start();
765 mesh.faceZones().clearAddressing();
770 }
while (nChangedFacesInMesh > 0 || nCycles > 100);
773 mesh.setInstance(
mesh.time().timeName());
775 if (nChangedFacesInMesh > 0)
778 <<
"with the original mesh"