[petsc-users] DMPlex PETSCViewer for Surface Component

Mike Michell mi.mike1021 at gmail.com
Fri Aug 5 08:09:32 CDT 2022

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,

> 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 <jed at jedbrown.org> wrote:
>> I would create a sub-DM containing only the part you want to view.
>> Mike Michell <mi.mike1021 at gmail.com> 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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220805/d4e9ada1/attachment.html>

More information about the petsc-users mailing list