<div dir="ltr">Dear all, <div><br></div><div>I have noticed what I think is a bug with a 3D DMPlex box mesh with periodic boundaries.</div><div>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.</div><div><br></div><div>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.</div><div>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.</div><div><br></div><div>I hope it is just a mishandling from my end !!</div><div><br></div><div>Thank you in advance for your help, </div><div><br></div><div>Cheers, </div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div><div><div><div><div>Thibault </div></div></div></div></div><div><br></div></div></div></div></div></div></div></div></div>