<div class="gmail_quote">On Wed, Dec 14, 2011 at 19:34, 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">
We could support a matrix-free SNESVI if extracting a submatrix from a MF matrix was implemented by pre-/post-multiplying it by a suitable P, no?</blockquote></div><br><div>You can MatGetSubMatrix() a MatMFFD. The implementation calls MatCreateSubMatrix() and is used by PCFieldSplit.</div>