[Phil_]

Hi all,

I'm trying to get all the faces sharing a point.

For the internal faces the following works:

labelList facesSharingPoint = mesh.pointFaces()[pointLabelI];

In order to get the boundary faces I tried:

labelList boundaryFacesSharingPoint = mesh.boundaryMesh()[patchI]pointFaces()[pointLabelI];

This fails, as the points are defined on the patches by local labels. How is it possible to return the boundary faces through a global point label?

Philip

[Jerryfan]

Hi,



Each of the fvPatch object maintains a polyPatch instance. There is a start() function defining in polyPatch which returns the start label of this patch in the polyMesh face list. This returned label is in the global level. The other thing about the patch label is that for each of the patch, its storage in the global polyMesh face list is continuous. Therefore, if you add the start label of that patch to the local labels you obtained, you will be able to get the global label.

forAll (mesh.BoundaryMesh(), PatchI)
{
label pstart=mesh.boundaryMesh()[patchI].start();
labelList boundaryFacesSharingPoint = mesh.boundaryMesh()[patchI]pointFaces();
forAll (boundaryFacesSharingPoint, pI)
{
boundaryFacesSharingPoint[pI] += pstart;
}
}

I didn't test the above code. But the idea should work for you. I hope it will help you.

[Phil_]

Hi Jerryfan,

thanks for your suggestion.

Maybe I'm not fully understanding your idea, but the main issue here is the relation of global vertice labels to local patch vertice labels.

mesh.pointFaces()

returns a labelListList containing the global vertice labels and the global labels of the internal faces sharing the respective vertice. I need to know the boundary faces sharing a vertice too.

mesh.boundaryMesh()[patchI]pointFaces()

returns a labelListList containing the local vertice labels of the patch and the local labels of the boundary faces sharing the respective vertice.

So the problem is, that I only know the global vertice label. Is there a function etc. to get the relation of global to local patch vertice labels?

Philip

[Phil_]

I found an much easier way to accomplish the original task, so finding all the faces sharing a point is no longer needed

[Jerryfan]

Would you mind sharing it?

[Phil_]

It has nothing to do with the original question, so I guess it's not worth sharing. I'm
avoiding any search for faces sharing a point and instead utilize

mesh.pointPoints()[pI];

to do stuff with the direct neighbor points.