<div dir="ltr">You may try the contour integral method (CISS). Obviously, it is computationally more expensive. I used the contour integral method on MATLAB and it worked quite well. I expected the same result with SLEPc.</div><br><div class="gmail_quote"><div dir="ltr">On Thu, 28 Jun 2018 at 23:22, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">This is a known issue of Krylov solvers. Sometimes you can workaround it by increasing the number of eigenvalues to compute (nev, for instance request twice as much as you need) and maybe reduce the mpd parameter to force more restarts. This approach is quite wasteful and still there is no guarantee that you get the full multiplicities right.<br>
<br>
To avoid the issue completely, use the lobpcg solver. However, SLEPc's implementation of lobpcg is not very well tuned, I have to take some time to improve it. You can also try the BLOPEX interface (configure --download-blopex in SLEPc).<br>
<br>
Jose<br>
<br>
<br>
> El 28 jun 2018, a las 22:53, Greg Meyer <<a href="mailto:gregory.meyer@gmail.com" target="_blank">gregory.meyer@gmail.com</a>> escribió:<br>
> <br>
> Hi all,<br>
> <br>
> I'm using SLEPc to solve for some eigenvalues. If I ask for the eigenvalues with the leftmost real part (using SMALLEST_REAL), and there are repeated eigenvalues, it appears that the solver does not always find all of the repeated values before jumping up to the next unique eigenvalue. For example if the lowest eigenvalues of my matrix are [-12, -10, -10, -10, -8], and I request 4 eigenvalues, I may get [-12, -10, -10, -8]. I'm using the default solver. Is there any way, or any solver, that ensures that no eigenvalues are skipped?<br>
> <br>
> Thanks in advance,<br>
> Greg<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Vinh Pham-Xuan, Ph.D<br>School of Electronic Engineering<br>Dublin City University<br>Ireland<br></div></div></div>
<br>
<pre style="background-color:rgb(255,255,255)"><i style="color:rgb(146,146,146);font-family:Verdana;white-space:normal;font-size:10px"><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal"><i style="color:rgb(146,146,146);font-family:verdana;font-size:10px"><p style="font-style:normal;font-family:verdana,arial,helvetica,sans-serif;display:inline"><strong>Séanadh Ríomhphoist/</strong></p></i><i style="color:rgb(146,146,146);font-family:verdana;font-size:10px"><p style="font-style:normal;font-family:verdana,arial,helvetica,sans-serif;display:inline"><strong>Email Disclaimer</strong></p></i></div><pre style="white-space:pre-wrap;color:rgb(34,34,34);font-style:normal"><i><div style="font-family:arial,sans-serif;font-style:normal;white-space:normal"><pre style="white-space:pre-wrap"><i><p style="display:inline"><span style="color:rgb(146,146,146);font-family:verdana,arial,helvetica,sans-serif;white-space:normal">Tá an ríomhphost seo agus aon chomhad a sheoltar leis faoi rún agus is lena úsáid ag an seolaí agus sin amháin é. <font color="#1155cc"><a href="https://www4.dcu.ie/iss/seanadh-riomhphoist.shtml" target="_blank">Is féidir tuilleadh a léamh anseo.</a></font></span><span style="color:rgb(146,146,146);font-family:verdana,arial,helvetica,sans-serif;font-style:normal;white-space:normal"><a href="https://www4.dcu.ie/iss/seanadh-riomhphoist.shtml" target="_blank"> </a></span></p></i></pre></div><div style="font-family:arial,sans-serif;font-style:normal;white-space:normal"><i style="white-space:pre-wrap"><p style="color:rgb(146,146,146);font-family:verdana,arial,helvetica,sans-serif;white-space:normal;font-style:normal;display:inline"><em>This e-mail and any files transmitted with it are confidential and are intended solely for use by the addressee. <font color="#1155cc"><a href="https://www4.dcu.ie/iss/email-disclaimer.shtml" target="_blank">Read more here.</a></font> </em></p></i></div></i></pre></i></pre>