[petsc-users] specifying vertex coordinates using DMPlexCreateFromCellListParallel

[Matthew Knepley]

On Tue, Apr 11, 2017 at 9:21 AM, Hassan Raiesi <
Hassan.Raiesi at aero.bombardier.com> wrote:

> Hello,
>
>
>
> I’m trying to use DMPlexCreateFromCellListParallel to create a DM from an
> already partitioned mesh,
>
> It requires an array of numVertices*spaceDim numbers, but how should one
> order the coordinates of the vertices?
>

Global order. Here is the idea. You must read the file in chunks so that
each proc can read its own chunk in parallel
without talking to anyone else.


> we only pass the global vertex numbers using ‘const int cells[]’ to define
> the cell-connectivity, so passing the vertex coordinates in local ordering
> wouldn’t make sense?
>

Yes.


> If it needs to be in global ordering, should I sort the global index of
> the node numbers owned by each rank (as they wont be continuous).
>

Nope.

  Thanks,

    Matt


>
>
> Thank you
>
>
>
> Hassan Raiesi,
>
> Bombardier Aerospace
>
> www.bombardier.com
>
>
>



-- 
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

http://www.caam.rice.edu/~mk51/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170625/78d8f9e2/attachment.html>