Hello

In the attached test, I read a small grid made in gmsh with periodic bc.

This is a 2d mesh.

The cell numbers are shown in the figure.

All faces have length = 2.5

But using PetscFVFaceGeom I am getting length of 7.5 for some faces. E.g.,

face: 59, centroid = 3.750000, 2.500000, normal = 0.000000, -7.500000

===> Face length incorrect = 7.500000, should be 2.5

support[0] = 11, cent = 8.750000, 3.750000, area = 6.250000

support[1] = 15, cent = 8.750000, 1.250000, area = 6.250000


There are also errors in the orientation of normal.

If we disable periodicity in geo file, this error goes away.

Thanks
praveen

Attachment: dmplex.c
Description: Binary data

Attachment: ug_periodic.geo
Description: Binary data

Reply via email to