[patchsploration]

OpenFOAM has the ability to define physically discontinuous patches. For example, I grabbed a screenshot from the programmers guide.

I'd like to find subsets of the physically continuous faces of such patches -- ideally by shared patch face vertices, or similar. For example, "Patch 1" in picture would have two of these subsets.

I'm thinking:

+ Use patch.faceCells() and get vertices ? How ?
+ Use a PatchPrimitive feature like edgeLoops() ?

Currently, I'm accomplishing this task by physical proximity of patch face centers and storing face IDs in HashTables, and it doesn't satisfy. Rest assured I have poked around on the forum for a better answer, but don't see it.

I appreciate any suggestion of a critical forward path!