On Thu, May 17, 2012 at 6:36 PM, Andrew Spott <span dir="ltr"><<a href="mailto:andrew.spott@gmail.com" target="_blank">andrew.spott@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Is there a way to explicitly say that a sparse matrix is symmetric or hermitian and only save/load half of it?</blockquote><div><br></div><div>PETSc distinguishes two kinds of information. Information about the operator (I am symmetric) can be set with</div>
<div>options. This can be queried by other code in order to optimize algorithms. It can also be checked by functions</div><div>like MatIsSymmetric().</div><div><br></div><div>Information about an implementation is in the type, like MATSBAIJ which says I only store half my values. There</div>
<div>is no need for a symmetric matrix to be stored this way, but you can choose this if you want.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
-Andrew</font></span></blockquote></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<br>