On Wed, 10 Nov 2010, Jed Brown wrote:
On Tue, Nov 9, 2010 at 15:43, Tim Kroeger <[email protected]> wrote: 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.
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).
Yes, this fixes it. Thank you very much again. Best Regards, Tim -- Dr. Tim Kroeger CeVis -- Center of Complex Systems and Visualization University of Bremen [email protected] Universitaetsallee 29 [email protected] D-28359 Bremen Phone +49-421-218-7710 Germany Fax +49-421-218-4236