MatMatMult() in petsc is not well-tested for complex - could be buggy.<div>Can you send us the matrices A and B in petsc binary format for investigation?</div><div><br></div><div>Hong<br><br><div class="gmail_quote">On Mon, Feb 6, 2012 at 5:55 AM, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear PETSc team,<br>
<br>
I try to use:<br>
call MatMatMult(A,B,MAT_INITIAL_<u></u>MATRIX,PETSC_DEFAULT_DOUBLE_<u></u>PRECISION,C,ierr);CHKERRQ(<u></u>ierr)<br>
<br>
Where both A and B are rectangular, but A is sparse and B is dense. Both are double complex and distributed.<br>
The product PETSc gives me contains some errors in some part of the matrix.<br>
I output A, B and C then computed product in matlab.<br>
<br>
Attached you see figure plotted as:<br>
imagesc(log10(abs(C-Cm)))<br>
<br>
Where Cm -- product computed in matlab.<br>
<br>
The pattern and amplitude vary depending on the number of cores I use. This picture is obtained for 48 cores (I've tried 12, 64 cores as well).<br>
<br>
Where should I look for possible explanation?<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Regards,<br>
Alexander<br>
</font></span></blockquote></div><br></div>