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 <[email protected]
<mailto:[email protected]>> wrote:
Michael Lange <[email protected]
<mailto:[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
