[petsc-users] DMPlex box mesh periodicity bug (?)

[Thibault Bridel-Bertomeu]

Dear all,

I have noticed what I think is a bug with a 3D DMPlex box mesh with
periodic boundaries.
When I project a function onto it, it behaves as if the last row of cells
in X and in Y direction do not have the right coordinates.

I attach to this email a minimal example that reproduces the bug (files
mwe_periodic_3d.F90, wrapper_petsc.c, wrapper_petsc.h90, makefile), as well
as the output of this code (initmesh.vtu for an output of the DM,
solution.vtu for an output of the data projected onto the mesh). There is
also a screenshot of what's going on.
If one considers the function I project onto the mesh, what should normally
happen is that there is a "hole" is the density field around the x=0, y=0
region, the rest being equal to one.

I hope it is just a mishandling from my end !!

Thank you in advance for your help,

Cheers,

Thibault
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrapper_petsc.c
Type: application/octet-stream
Size: 2474 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0006.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mwe_periodic_3d.F90
Type: application/octet-stream
Size: 3790 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0007.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: makefile
Type: application/octet-stream
Size: 403 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0008.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrapper_petsc.h90
Type: application/octet-stream
Size: 288 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0009.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: initmesh.vtu
Type: application/octet-stream
Size: 955335 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0010.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solution.vtu
Type: application/octet-stream
Size: 988234 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0011.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Capture d?e?cran 2021-07-31 a? 11.53.06.png
Type: image/png
Size: 1082614 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210731/25463d14/attachment-0001.png>