<div class="gmail_quote">On Tue, Mar 15, 2011 at 23:43, 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;">
I'm just asking that all matrix classes be able to propagate this information NOT just MatNest().</blockquote></div><br><div>Okay, great. My point was just that MatNest pretty much already does this [1], with MatAIJ, you need some new API to declare part of the matrix as unchanged (otherwise I don't know how to efficiently determine it).</div>
<div><br></div><div>[1] MatAssemblyEnd_Nest() still calls MatAssemblyEnd() on submatrices. It could keep track of whether submatrices have been "gotten" (if not, they can't have been modified) or (better) submatrices could increment state in MatAssemblyEnd() only when they have been modified on some process.</div>