The test for symmetry is not implemented for all matrix types. Looking at the code, it seems to only be SeqAIJ, but MatIsTranspose(A,A,...) would also work for MPIAIJ.<div class="gmail_extra"><br><br><div class="gmail_quote">
On Sun, Dec 2, 2012 at 8:45 AM, w_ang_temp <span dir="ltr"><<a href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="line-height:1.7;font-size:14px;font-family:arial"><p>Hello,</p>
<p>    I use MatIsSymmetric to know if the matrix A is symmetric.</p>
<p>According to my model, it should be symmetric due to the theory.</p>
<p>But I always get the result 'PetscBool  *flg = 0', although I </p>
<p>set 'tol' a large value(0.001).</p>
<p>    Because the matrix is of 20000 dimension, I can not output the</p>
<p>matrix to the txt. So I want to konw if there is something to be paid attention to</p>
<p>about the function 'MatIsSymmetric' in version 3.2. Or do I have some other ways</p>
<p>to determine the symmetry.I think symmetry is one of the most important thing</p>
<p>in my analysis.</p>
<p>    Thanks.</p>
<p>                                               Jim</p>
<div></div></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br></div>