<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jul 10, 2015 at 4:28 PM, Julian Andrej <span dir="ltr"><<a href="mailto:juan@tf.uni-kiel.de" target="_blank">juan@tf.uni-kiel.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
i'm trying to solve a generalized eigenvalue problem which is the<br>
stokes equations discretized by finite elements (with fenics) and<br>
producing a banded matrix (reordered structure) of the well known<br>
block form<br>
<br>
[N    Q]<br>
[QT  0]<br>
<br>
The domain is the unit square with dirichlet boundary condition<br>
evaluating to zero at all boundary nodes.<br>
<br>
A is the block matrix in banded reordered structure and M is the mass matrix.<br></blockquote><div><br></div><div>You have not written an equation here. Are you solving</div><div><br></div><div>  Q^T N^{-1} Q x = \lambda M x</div><div><br></div><div>If so, the Schur complement can be rank deficient by 1 if you have Neumann</div><div>conditions on the pressure. I don't know how SLEPc handles this.</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">
I'm using slepc4py for the eigenvalue calculation.<br>
<br>
E = SLEPc.EPS().create()<br>
E.setOperators(A, M)<br>
E.setDimensions(NEV, PETSc.DECIDE)<br>
E.setFromOptions()<br>
<br>
I always get the error<br>
[0] Zero pivot row 1 value 0 tolerance 2.22045e-14<br>
<br>
I cannot find a combination which solves the eigenvalues for this problem.<br>
<br>
The system itself solves fine with a KSP solver object from PETSc<br>
using tfmqr with the icc PC.<br>
<br>
I can assemble the matrix with a penalty term for the pressure and<br>
calculate the eigenvalues with a direct solver, but i try to avoid<br>
that.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>