<div dir="ltr">On Fri, Aug 16, 2013 at 10:09 AM, Bill McGrory <span dir="ltr"><<a href="mailto:mcgrory@aerosoftinc.com" target="_blank">mcgrory@aerosoftinc.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">More information pertaining to the referenced post.<br>
<br>
I am using default parameters for my KSP solver, so my routines are pretty simple I don't pass any options in through the command line, so this is all default stuff here.<br>
<br>
        I create, and fill A and b, and then call<br>
MatAssemblyBegin(A,MAT_FINAL_<u></u>ASSEMBLY);<br>
<br>
                MatAssemblyEnd(A,MAT_FINAL_<u></u>ASSEMBLY);   <br>
<br>
        VecAssemblyBegin(b);<br>
        VecAssemblyEnd(b);<br>
<br>
        VecSet(x,0.);<br>
        <br>
        // Solve A x = b<br>
<br>
        KSPCreate(PETSC_COMM_WORLD,&<u></u>ksp);<br>
<br>
        KSPSetOperators(ksp,A,A,<u></u>DIFFERENT_NONZERO_PATTERN);<br>
<br>
        KSPSetFromOptions(ksp);<br>
<br>
        KSPSolve(ksp,b,x);<br>
<br>
The missing diagonals, I see are in the preconditioner, not the original Matrix, so I thought I would check my original.<br></blockquote><div><br></div><div>From the code above, this makes no sense. You do not have a separate preconditioner matrix.</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
When I make a call to MatMissingDiagonal, after assembling my matrix, I get the error telling me that MatMissingDiagonal is not supported for a mpibaij matrix.<br>
<br>
Do I have any alternative for querying my assembled matrix?<br></blockquote><div><br></div><div>You can always use MatGetValues().</div><div><br class="">   Matt</div><div>  </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

Thanks again<span class=""><font color="#888888"><br>
Bill<br>
<br>
<br>
</font></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>