On Tue, Nov 12, 2013 at 9:59 AM, Jed Brown <[email protected]> wrote:
> Michael Lange <[email protected]> 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
