<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 4, 2013 at 12:43 PM, Paul Mullowney <span dir="ltr"><<a href="mailto:paulm@txcorp.com" target="_blank">paulm@txcorp.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":27w">I'm guessing this question has been asked and answered before on this list. Apologies for the repeat. Anyway, suppose I have a symmetric matrix stored as only L in a Matrix Market file. Then, I read in the L part with MatSetValues() and then do:<br>

<br>
MatSetOption(A,MAT_SYMMETRIC,<u></u>PETSC_TRUE);<br>
<br>
Will PETSc then do a KSP solve (with say an ICC preconditioner) on the full matrix even though the upper part hasn't been explicitly created/inserted via MatSetValues()? The example below suggests not but I want to be absolutely sure.<br>

<br>
<a href="http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex2.c.html" target="_blank">http://www.mcs.anl.gov/petsc/<u></u>petsc-current/src/ksp/ksp/<u></u>examples/tutorials/ex2.c.html</a><br>
<br>
Is there any way to work with just L?</div></blockquote></div><br>Insert it as the upper triangular part in a MATSBAIJ, then everything you want will work.</div></div>