On Fri, Feb 24, 2012 at 2:39 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
On Feb 24, 2012, at 2:23 PM, Jed Brown wrote:<br>
<br>
> We need a richer field interface. Some way for a client of DM to discover vector and tensor structure.<br>
<br>
</div>   No idea what a tensor is, but DMGetSubDMS(dm,PetscInt *cnt,IS *ises,DM *dms);   and DMPushSubDMType(DM,"name of a type of splitting like u,v,p or velocity,p"); DMPopSubDMType(); and DMGetAvailableSubDMTypes(DM,char **);<br>

<br>
  Do you need more than that?<br></blockquote><div><br></div><div>Yuck Yuck Yuck. Why are we persisting in shoving all of the index stuff into the DM interface. I think this is so</div><div>unwieldy and complicated. PetscSection is small, powerful, and does everything suggested and more. I see</div>
<div>no argument in support of stuff like the above.</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">
    Barry<br>
<br>
Note that current DMDA's default subdms served up would be a DM for each field consisting of one of the degrees of freedom per node and it could also serve up  DMs that represent any collections of those for example dof 0 and 1 per node, then dof 2 per node. We'll need to come up with some naming convention for these.<br>

<div class="HOEnZb"><div class="h5"><br>
<br>
> Possibly also (eventually) a way to solve equations of state so that we can formulate reduced problems in non-conservative variables.<br>
><br>
> On Feb 24, 2012 2:19 PM, "Matthew Knepley" <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
> On Fri, Feb 24, 2012 at 1:58 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
> <a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/518ff70e8a0a" target="_blank">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/518ff70e8a0a</a><br>
><br>
> Matt, I want to get rid of these conditionals, not add more. We should have a DM base interface for getting fields on which to split, that common interface should _replace_ the DMComposite specialization.<br>
><br>
> I agree. However, I won't put anything in until I do it by hand once.<br>
><br>
> Also, I was tempted to just promote DMGetGlobalISes(), but its not quite right. I am planning<br>
> to put the IS method in PetscSection. I am hoping eventually DMDA uses PetscSection for<br>
> layout.<br>
><br>
>    Matt<br>
><br>
> --<br>
> 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<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>