[petsc-users] DMPlex boundary condition from array

[Pierre Seize]

Hi
I would like to solve div(u) = 0 with a FE method in a spherical mesh made
of tetrahedron.
My problem concern the boundary condition : a Dirichlet boundary condition
on the surface of the sphere where the value is given by another program,
in the form of a list of indices corresponding to the vertices on the
boundary and a list of values to be applied on those vertices.

I am able to use DMPlexCreateFromCellListPetsc to create the DM on process
0, but the moment I distribute the DM the vertices are renumbered and I
loose the correspondence with the boundary describing arrays.

I have two questions:
    - is there a way to get the mapping between old indices and new global
indices after the distribute operation?
    - how can I apply a boundary condition such as mine, i.e. set given
values at given vertices?

Thanks in advance
Pierre Seize
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240925/5e37e88a/attachment.html>