<div class="gmail_quote">On Wed, Nov 30, 2011 at 14:04, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":1uq">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>
</div></blockquote></div><br><div>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.</div>