[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