[petsc-users] Question on DMPlex

Wed Feb 9 21:04:22 CST 2022

Dear developers,

I have a mesh with topological dimension 2 (with triangular faces) embedded in a 3D space stored in a hdf5 file. This file contains the coordinates (x,y,z) of each vertex and also, it contains a list of vertices that make each element. I was able to read the mesh into a DMPlex using DMPlexCreateFromFile (without interpolation, i.e. interpolate=PETSC_FALSE ). When I do this, the elements are stored at stratum height 0 and the vertices are stored at stratum height 1 (as checked via DMPlexGetHeightStratum). However, when I interpolate the mesh (interpolate=PETSC_TRUE), it seems that some entities are created in stratum height 2 and 3. From my understanding of DMPlex, this should not happen for 2D elements because there are edges but no 3D cells.

Does DMPlex assume that the mesh in 3D has a topological dimension 3 and is it possible to use 2D mesh embeded in 3D space?

Best regards

François Fillion-Gourdeau, PhD
INRS-EMT
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220210/ed65d885/attachment.html>