Alexander :<div>I need a simple code from you to reproduce the error for debugging.</div><div><br></div><div>If mumps' solver is activated, then you should get error (sequential and parallel)</div><div><div>[0]PETSC ERROR: No support for this operation for this object type!</div>
<div>[0]PETSC ERROR: MatMatSolve_MUMPS() is not implemented yet!</div><div><br></div><div>If petsc's direct solver is activated, then you should not be able to run it in parallel. </div><div><br></div><div>Please check ~petsc/src/mat/examples/tests/ex125.c</div>
<div>to see if you used same calling procedure as this example.</div><div>Note: I fixed a bug in this example and patched petsc-3.3.</div><div>The update copy of ex125.c is attached for your convenience.</div><div><br></div>
<div>Hong</div><br><div class="gmail_quote">On Mon, Aug 13, 2012 at 10:56 AM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@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">
Alexander :<div>petsc-mumps interface does not support MatMatSolve()</div><div>see ~petsc/src/mat/impls/aij/mpi/mumps/mumps.c</div><div>Likely MatMatSolve_Basic() is called which calls MatSolve()</div><div>repeatedly, thus slow.</div>

<div><br></div><div>In addition, MUMPS </div><div>only supports centralized rhs. Petsc interface must scatter </div><div>distributed rhs into a seqential rhs vector for MUMPS.</div><div>If num of rhs is large, would be slow here as well.</div>

<div><br></div><div>The error in your report likely is a bug in petsc. I'm investigating it</div><div>and will let you know once it is fixed.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Hong</div>
</font></span><div class="HOEnZb"><div class="h5"><div><br></div><div><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>On 13.08.2012 17:23, Hong Zhang wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Alexander :<br>
Hmm, isrow is a matrix ordering, likely a bug :-(<br>
Can you send us a stand-alone code that reproduces this error?<br>
</blockquote>
<br></div>
Hong,<br>
<br>
I can try to pull it out, but this may be tricky.<div><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Do you use an external package for MatMatSolve()?<br>
</blockquote>
<br></div>
Yes, MUMPS.<br>
I use ksp with preonly.<br>
<br>
Originally I thought it will be faster to solve multiple problems using MatMatSolve.<br>
This turned out to be not true however, but the code is still like that.<span><font color="#888888"><br>
<br>
-- <br>
Regards,<br>
Alexander<br>
<br>
</font></span></blockquote></div><br></div>
</div></div></blockquote></div><br></div>