Thank you Jed and Matt. That is exactly what I wanted to do. I was wondering if there is an example of mapping between DMPLEXs. I have several requirements of the mapping for my code, such as - mapping solution vectors on cell-centroids into vertices - mapping from volume dmplex into sub-dm for wall surface region Any example will help. Many thanks, Mike
Yes, you can use DMPlexCreateSubmesh() (and friends depending on exactly what kind of submesh you want). This will allow you to create a vector over only this mesh, and map your volumetric solution to that subvector. Then you can view the subvector (which pulls in the submesh).
Thanks,
Matt
On Mon, Aug 1, 2022 at 10:59 PM Jed Brown <[email protected]> wrote:
I would create a sub-DM containing only the part you want to view.
Mike Michell <[email protected]> writes:
Hi,
I am a user of DMPlex object in 3D grid topology. Currently the solution field is printed out using viewer PETSCVIEWERVTK in .vtu format. By doing that the entire volume components are written in the file. I was wondering if there is an option that I can tell PETSc viewer to print out only surface component, instead of the entire volume.
Best, Mike
-- 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://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>