Foam::labelList Foam::structuredRenumber::renumber ( const polyMesh& mesh, const pointField& points ) const { if (points.size() != mesh.nCells()) { FatalErrorInFunction << "Number of points " << points.size() << " should equal the number of cells " << mesh.nCells() << exit(FatalError); } const polyBoundaryMesh& pbm = mesh.boundaryMesh(); const labelHashSet patchIDs(pbm.patchSet(patches_)); label nFaces = 0; forAllConstIter(labelHashSet, patchIDs, iter) { nFaces += pbm[iter.key()].size(); } // Extract a submesh. labelHashSet patchCells(2*nFaces); forAllConstIter(labelHashSet, patchIDs, iter) { const labelUList& fc = pbm[iter.key()].faceCells(); forAll(fc, i) { patchCells.insert(fc[i]); } }
Foam::labelList Foam::structuredDecomp::decompose ( const polyMesh& mesh, const pointField& cc, const scalarField& cWeights ) { const polyBoundaryMesh& pbm = mesh.boundaryMesh(); const labelHashSet patchIDs(pbm.patchSet(patches_)); label nFaces = 0; forAllConstIter(labelHashSet, patchIDs, iter) { nFaces += pbm[iter.key()].size(); } // Extract a submesh. labelHashSet patchCells(2*nFaces); forAllConstIter(labelHashSet, patchIDs, iter) { const labelUList& fc = pbm[iter.key()].faceCells(); forAll(fc, i) { patchCells.insert(fc[i]); } }