Hi Matt
In fact, I created a routine for my specific mesh generation. I firstly, create
the DMPlex object at first, then setting the chart, then doing symmetrize and
interpolate. Next, I created coordinate section and giving the coordinates.
After finishing giving the coordinates, I declare the distribution in routine,
with this two lines:
DMPlexDistribute(*dm, 0, NULL, &dmDist);
if (dmDist) {DMDestroy(dm); *dm = dmDist;}
Thanks
Amir
From: Matthew Knepley [mailto:[email protected]]
Sent: Friday, February 16, 2018 6:34 AM
To: Mohammad Hassan Baghaei <[email protected]>
Cc: PETSc <[email protected]>
Subject: Re: [petsc-users] Dealing with DMPlexDistribute()
On Thu, Feb 15, 2018 at 4:40 PM, Mohammad Hassan Baghaei
<[email protected] <mailto:[email protected]> > wrote:
Hi
I am using DMPlex as interface for mesh generation. On single core, I got
around 30000 mesh cells. Whenever, I run on multiple core, say 3, in the output
file for DM, VTK, I got 3 times mesh cell and point numbers. Does it mean that
DMPlexDistribute() does not work properly!
Hi Amir,
It sounds like you are generating a mesh on every process. How are you
generating the mesh?
Thanks,
Matt
Thanks
Amir
--
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.caam.rice.edu/~mk51/>
