[petsc-dev] DMPlex with vertical numbering?

[Matthew Knepley]

On Tue, Nov 12, 2013 at 9:59 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> Michael Lange <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/20131112/ca3ce72e/attachment.html>