On Tue, Dec 13, 2011 at 9:36 AM, Uwe Schlifkowitz <span dir="ltr"><<a href="mailto:uwe.schlifkowitz@uibk.ac.at">uwe.schlifkowitz@uibk.ac.at</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">
On 12.12.2011, at 20:05, Barry Smith wrote:<br>
><br>
> On Dec 12, 2011, at 11:49 AM, Uwe Schlifkowitz wrote:<br>
><br>
>> I'm continuing Clemens Domanig's work. Currently I'm stuck at his problem as described here: <a href="http://lists.mcs.anl.gov/pipermail/petsc-users/2011-August/009571.html" target="_blank">http://lists.mcs.anl.gov/pipermail/petsc-users/2011-August/009571.html</a><br>
>><br>
>> As far as I understand, getting the diagonal matrix from MUMPS is not possible, so Clemens (and I) resort to MATSOLVERPETSC to obtain D.<br>
><br>
> We could add support for MatGetDiagonal() for PETSc factored matrices if you really need it, but what are you using it for? Do you want the inertia? We could give that for PETSc factored matrices.<br>
><br>
><br>
> Barry<br>
><br>
<br>
K is an element stiffness matrix. K is related to displacement d and force f by K d = f. Since K is not invertible, a solution can be found by means of LDL^T decomposition where K = LDL^T .<br>
<br>
Here L is the lower triangular matrix and D is the diagonal matrix, which is what I am looking for. The number of negative values in D correspond to the number of negative eigenvalues of K and serve as an indicator of the occurrence of branching points.<br>
<br>
I am not sure if this has to do with inertia, but I am fairly new to all of this so suggestions are very welcome.</blockquote><div><br></div><div>Yes, it looks like what you want</div><div><br></div><div> <a href="http://en.wikipedia.org/wiki/Sylvester's_law_of_inertia">http://en.wikipedia.org/wiki/Sylvester's_law_of_inertia</a></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>
Uwe</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<br>