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><div><br></div><div>Hong</div><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 class="im">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 class="im"><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 class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Regards,<br>
Alexander<br>
<br>
</font></span></blockquote></div><br></div>