On Fri, Jul 31, 2026 at 6:28 AM Matteo Semplice via petsc-users <
petsc-users(a)lists.mcs.anl.gov> wrote:
> Dear PETSc,
>
> I do not understand the behaviour of DMPlexGetCellCoordinates on a
> Plex generated to a periodic MSH mesh.
>
> I have in the code
>
> DM dmMesh;
> PetscCall( DMCreate(MPI_COMM_WORLD, &dmMesh) );
> PetscCall( DMSetType(dmMesh, DMPLEX) );
> PetscCall( DMSetFromOptions(dmMesh) );
>
> and pass the option
>
> -dm_plex_filename periodiccube.msh
>
> to load a periodic mesh of tetrahedra for the unit cube.
>
> When calling DMPlexGetCellCoordinates on 3d elements (cells), I get DG
> coordinates (independently on passing or not -dm_localize)
>
> But when calling DMPlexGetCellCoordinates on 1d or 2d elements (edges
> and faces), I DO NOT get DG coordinates (independently on passing or not
> -dm_localize).
>
> I can work around this, but is there a way to have DG coordinates also
> for intermediate elements in the Plex? I see an option in
> DMCreateBoxMesh, but no options for loaded meshes.
>
Sorry it took me so long to respond.
The intention is for -dm_localize_height to work. Does that work for you?
Thanks,
Matt
> Thanks
>
> Matteo
>
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
https://urldefense.us/v3/__https://www.cse.buffalo.edu/*knepley/__;fg!!G_uC… <https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCf… >
