[petsc-dev] ML Reuse redux

[Jed Brown]

On Mon, Aug 29, 2011 at 11:10, John Fettig <john.fettig at gmail.com> wrote:

> I think the way to do it is similar to the way it is done in the
> function ML_Gen_MultiLevelHierarchy_UsingSmoothedAggr_ReuseExistingAgg,
> although that function rebuilds R, A, and P.  So we just have to
> rebuild A.  Here's some code that does that, it seems to work as well
> as the previous code I sent that only worked in serial (this works in
> parallel too).
>

Would you be willing to make this into a Mercurial patch against petsc-dev?
Just commit your changes locally, hg export them, and send us the patch.


>  Could this be made to be the default if you set the
> flag "SAME_NONZERO_PATTERN"?  Then if you want to rebuild the
> preconditioner, call with "DIFFERENT_NONZERO_PATTERN".
>

Perhaps it can be made the default, but there should be a runtime option to
recompute even when SAME_NONZERO_PATTERN is set (because the nonzero pattern
of the operator, which is what that flag is supposed to indicate, did not
change).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110829/e1e3c21d/attachment.html>