[petsc-dev] DMPlex with vertical numbering?

Michael Lange michael.lange at imperial.ac.uk
Wed Nov 13 10:21:39 CST 2013


Hi,

Yes, a PetscSFMoveRoots() wrapper would be greatly appreciated, 
especially since we want to ultimately wrap and access this via petsc4py.

Many thanks
Michael


On 12/11/13 17:02, Matthew Knepley wrote:
> On Tue, Nov 12, 2013 at 9:59 AM, Jed Brown <jedbrown at mcs.anl.gov 
> <mailto:jedbrown at mcs.anl.gov>> wrote:
>
>     Michael Lange <michael.lange at imperial.ac.uk
>     <mailto:michael.lange at imperial.ac.uk>> writes:
>
>     > Hi Matt,
>     >
>     > Thanks for the clarification, this does make sense. I have
>     achieved the
>     > local DoF ordering I want by only permuting the offsets in the
>     created
>     > section. However, I am now trying to extract the according halo
>     exchange
>     > information from the associated PetscSF object and, although the
>     local
>     > leafs have been renumbered according to my permutation, the
>     roots they
>     > are connected to have not been updated. Is there an easy way to
>     update
>     > the roots or re-initialise the PetscSF from the reordered section?
>
>     You can SFBcast the new root numbers and create a new SF.
>
>     Since this is a somewhat common operation that can be optimized by the
>     SF implementation, it may make sense to provide PetscSFMoveRoots() or
>     some such.
>
>
> I already do this in DMPlexPermute(). We can abstract the code out.
>
>    Matt
>
> -- 
> 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131113/c92f345b/attachment.html>


More information about the petsc-dev mailing list