On Wed, Nov 30, 2011 at 3:10 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;">
<br>
On Nov 30, 2011, at 2:07 PM, Jed Brown wrote:<br>
<br>
> On Wed, Nov 30, 2011 at 14:04, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
> Ok, so how can we resolve this? The problem is that the default ordering is set when the PCCreate_Cholesky() and is set to be natural in order "to be safe" if sbaij matrix is used. We could change it to default to ND and then when PCSetUp_Cholesky() is called simply set the ordering to natural if it is a sbaij matrix. Kind of a hack but what do you suggest that doesn't involve tracking some other variable indicating if the user has set the ordering type to something else ...<br>

><br>
> Why the hell set it so soon? This is like choosing which side/norm to use in a Krylov method. Make the decision lazily (when it is requested or needed to do something), otherwise leave it undefined.<br>
<br>
   There is currently no MatOrdering object nor MatColoring object one simply calls MatGetOrdering with an MatOrderingType and MatGetColoring with a MatColoringType.<br>
<br>
   These should both be totally refactored to match the paradigm of other PETSc objects; for example KSP  or MatPartitioning. Make a ticket in your cumbersome ticketing web interface :-)</blockquote><div><br></div><div>Piling on: Make a ticket to make the DMDA HDF5 function respect timesteps like the regular HDF5 output.</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;"><font color="#888888"><br>
   Barry<br>
<br>
<br>
<br>
</font></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>