<div dir="ltr"><div><div><div><div>Yes, you are probably right: My code is not yet bug free (by all means!). However, I have been working very hard on it and I will see if I find something. B should not ever be really singular!<br>

<br></div>I do, however, still have some questions:<br><br></div>- Does slepc detect whether my matrix is tridiagonal and does it apply a direct method to solve it? Or do you talk about it in general: that it would be better for me to use a direct method that is not part of slepc?<br>

</div>- How do I make use of the information that both A and B have the same non-zero structure? Is there an easy way?<br></div>- What is the call "MatSetOption(mat,MAT_HERMITIAN,PETSC_TRUE,ierr)" do? Is it important? It does not seem to influence anything...<br>

<br>Thank you for your time and pacience!<br><br>Toon<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 19 August 2014 15:40, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">Toon Weyens <<a href="mailto:tweyens@fis.uc3m.es">tweyens@fis.uc3m.es</a>> writes:<br>
<br>
> Okay I reread again your answers and it seems that I missed the fact that I<br>
> do not need shift invert to find the EV with the largest real part?<br>
<br>
</div>Correct, inversion is way more expensive and used to find internal<br>
eigenvalues.<br>
<div class=""><br>
> Then the "trick" of ex. 13 does not apply (they look for something<br>
> close to 0).  So what could I do to avoid the LU?<br>
<br>
</div>Is your matrix B singular?  What are its boundary conditions?  You have<br>
to solve with B for the generalized eigenvalue problem.  That could be<br>
done iteratively, but a direct solve makes sense for a tridiagonal<br>
system and the zero pivot likely indicates a deeper issue with the<br>
formulation that won't vanish simply by avoiding LU.<br>
</blockquote></div><br></div>