Re: [petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Berend van Wachem
Dear Matt, Thanks for your quick response. I have a DMPlex with a polyhedral mesh, and have defined a number of vectors with data at the cell center. I have generated data for a number of timesteps, and I write the data for each point to a file together with the (x,y,z) co-ordinate of the cell

Re: [petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Berend van Wachem
Dear Matt, The problem is that I haven't figured out how to write a polyhedral DMplex in parallel. So, currently, I can write the Vec data in parallel, but the cones for the cells/faces/edges/nodes for the mesh from just one process to a file (after gathering the DMplex to a single process).

Re: [petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Matthew Knepley
On Mon, Jan 22, 2024 at 10:49 AM Berend van Wachem wrote: > Dear Petsc-Team, > > Is there a good way to define a unique integer number in each element > (e.g. a cell) of a DMPlex mesh, which is in the same location, > regardless of the number of processors or the distribution of the mesh > over

Re: [petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Matthew Knepley
On Mon, Jan 22, 2024 at 1:57 PM Berend van Wachem wrote: > Dear Matt, > > Thanks for your quick response. > I have a DMPlex with a polyhedral mesh, and have defined a number of > vectors with data at the cell center. I have generated data > for a number of timesteps, and I write the data for

Re: [petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Matthew Knepley
On Mon, Jan 22, 2024 at 2:26 PM Berend van Wachem wrote: > Dear Matt, > > The problem is that I haven't figured out how to write a polyhedral DMplex > in parallel. So, currently, I can write the Vec data > in parallel, but the cones for the cells/faces/edges/nodes for the mesh > from just one

[petsc-users] Unique number in each element of a DMPlex mesh

2024-01-22 Thread Berend van Wachem
Dear Petsc-Team, Is there a good way to define a unique integer number in each element (e.g. a cell) of a DMPlex mesh, which is in the same location, regardless of the number of processors or the distribution of the mesh over the processors? So, for instance, if I have a DMPlex box mesh,