<div class="gmail_quote">On Tue, Nov 9, 2010 at 15:43, Tim Kroeger <span dir="ltr"><<a href="mailto:tim.kroeger@cevis.uni-bremen.de">tim.kroeger@cevis.uni-bremen.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Okay, I implemented a workaround on user side (using a temporary vector and MatMult()). This seems to resolve the problem now. I'll leave it like this in the libMesh part, so that libMesh remains compatible with PETSc-3.1-p4. Nevertheless, if you fix this within a few days and launch PETSc-3.1-p5 with this fix (or just send me a patch), let me know, I can then test it.</blockquote>
</div><br><div>I fixed this, and also a bug with nontrivial MatDiagonalScale after MatGetSubMatrix. You can get it from hg, or in the next patch level (I don't know when that will be).</div><div><br></div><div>Jed</div>