[petsc-users] Nodes coordinates in distributed dmplex

[Alejandro D Otero]

Hi, sorry if this is an obvious question, but I cannot figure out how to
recover finite element nodes coordinates once I have distributed a mesh
stored as a dmplex. I am using petsc4py as interface to petsc rutines.

I first created a dmplex using:
dm.createFromCellList()

In a sequential run I got the coordinates with:
Coords = dm.getCoordinates()

which gave a sequential vector with the coordinates of the mesh nodes.

When I distribute the mesh with:
dm.distribute()

each mpi process has it own dm but the indexing of the vector resulting
from getCoordinates() or getCoordinatesLocal() seems not consistent with
the local numbering of the cells and nodes.

Which is the correct way of doing this in PETSc philosophy?

Thanks in advance,
Alejandro
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150709/32fc9b9c/attachment-0001.html>