<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Mar 5, 2018 at 9:55 AM, Åsmund Ervik <span dir="ltr"><<a href="mailto:Asmund.Ervik@sintef.no" target="_blank">Asmund.Ervik@sintef.no</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hmm, there's many elements to consider here. Learning to use Plex might be a better use of my time than setting up and maintaining a custom DMDA repartitioning. But DMDAs are a more natural fit to what I want to do...<br>
<br>
However, I'm not in a huge rush to do this immediately. Can I ping you back some time next week, Matt, to ask more about details on how to use Plex for this case? (e.g. what is the Plex equivalent of DMDAVecGetArrayF90?) Enjoy Japan.<br></blockquote><div><br></div><div>Yep. Jed is correct that the kind of adjacency in DA is quite unnatural in Plex. Things could be done to make this better, but they are not being done now.</div><div>However, this is not the case in 1D, where you can easily get neighbors and 2nd neighbors in the right order.</div><div><br></div><div>We can talk about it after SIAM PP.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Regards,<br>
Åsmund<br>
<div class="HOEnZb"><div class="h5"><br>
> -----Original Message-----<br>
> From: Jed Brown [mailto:<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>]<br>
> Sent: Monday, March 5, 2018 3:48 PM<br>
> To: Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>>; Tobin Isaac<br>
> <<a href="mailto:tisaac@cc.gatech.edu">tisaac@cc.gatech.edu</a>><br>
> Cc: Dave May <<a href="mailto:dave.mayhem23@gmail.com">dave.mayhem23@gmail.com</a>>; Åsmund Ervik<br>
> <<a href="mailto:Asmund.Ervik@sintef.no">Asmund.Ervik@sintef.no</a>>; <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><br>
> Subject: Re: [petsc-users] Load balancing / redistributing a 1D DM<br>
><br>
> Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
><br>
> > On Mon, Mar 5, 2018 at 9:01 AM, Tobin Isaac <<a href="mailto:tisaac@cc.gatech.edu">tisaac@cc.gatech.edu</a>><br>
> wrote:<br>
> ><br>
> >> This is a somewhat incomplete description of the steps in linear<br>
> >> partitioning.  The rest can be accomplished with PetscSF calls, but I<br>
> >> should wrap it up in a PetscPartitioner because it's a mistake-prone<br>
> >> operation.<br>
> >><br>
> ><br>
> > Jed likes to do everything by hand because it is transparent, but then<br>
> > you become the maintainer.<br>
> > I think this is easy to do in Plex, and we maintain the code. It is<br>
> > less transparent, which is the tradeoff.<br>
><br>
> DMDA gives you direct indexing of neighbors, second neighbors, etc.<br>
> DMPlex doesn't ensure that, though you could configure a DMPlex to sort<br>
> vertices according to a particular convention.  I don't know whether that<br>
> convention would be preserved on repartitioning/redistribution.<br>
><br>
> With your very expensive material models, indirect indexing wouldn't be a<br>
> performance issue, but still may impact readability.<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>