<div dir="ltr"><div>On Wed, Jan 16, 2013 at 12:09 PM, 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></div><div class="gmail_extra"><div class="gmail_quote">
<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"><div><br><br><br><br><br></div>
<div></div>
<div></div>
<div><br></div>>At 2012-12-02 23:10:51,"Jed Brown" <<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">>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"> </div>
<div class="gmail_extra"> </div>
<div class="gmail_extra">I made a test. I get all the elements of the matrix and compare the element A(i,j) with</div>
<div class="gmail_extra">A(j,i). The result shows that the matrix is symmetric. But the result of MatIsTranspose(A,A,0.0,matsysflg,ierr) </div>
<div class="gmail_extra">shows that it is unsymmetric(matsysflg=0). I think there is something that I donot know about using the function.</div></blockquote></div></blockquote><div><br></div><div style>Obviously your comparison has a tolerance.</div>
<div style><br></div><div style>  Matt</div><div> </div><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"><blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">
<div class="gmail_extra">
<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 style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div style="LINE-HEIGHT:1.7;FONT-FAMILY:arial;FONT-SIZE:14px">
<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></blockquote></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br><br clear="all"><div>
<br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>