<div dir="auto"><div style="font-size:12.8px" dir="auto"><div style="width:380.19px;margin:16px 0px"><div dir="auto">Hi<div dir="auto">I would like to solve div(u) = 0 with a FE method in a spherical mesh made of tetrahedron.</div><div dir="auto">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.</div><div dir="auto"><br></div><div dir="auto">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.</div><div dir="auto"><br></div><div dir="auto">I have two questions:</div><div dir="auto">    - is there a way to get the mapping between old indices and new global indices after the distribute operation?</div><div dir="auto">    - how can I apply a boundary condition such as mine, i.e. set given values at given vertices?</div><div dir="auto"><br></div><div dir="auto">Thanks in advance </div><div dir="auto">Pierre Seize </div></div></div></div></div>