<div class="gmail_quote">On Mon, Jun 6, 2011 at 13:44, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov">karpeev@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=":2on">A related matter: should there be a way to reset the MatSetValuesLocal method?<br></div></blockquote><div><br></div><div>PCShellSetOperation() can reset any operation.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":2on">
The setting of local values frequently depends on factors other than<br>
the type of the matrix.<br></div></blockquote><div><br></div><div>That's what the ISLocalToGlobalMapping is for.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":2on">
The default implementation depends only on the class (Mat) rather than<br>
the type (e.g., MATMPIAIJ)<br>
of the matrix.</div></blockquote></div><br><div>The definition of the "local space" is something that should not depend on the implementation.</div><div><br></div><div>It calls through mat->ops->setvalueslocal, so I don't understand your concern.</div>