<div dir="ltr">Does shell matrix work?</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 3, 2015 at 11:04 AM, Xujun Zhao <span dir="ltr"><<a href="mailto:xzhao99@gmail.com" target="_blank">xzhao99@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">One problem is that I don't have the explicit form of matrix D, but only u = D*v which can be obtained from my PETSc solver. How should I set up my EPSSolver in this case? Are there examples in SLEPc? Thanks.<span class="HOEnZb"><font color="#888888"><div><br></div><div>Xujun</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 3, 2015 at 10:59 AM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><br>
El 03/06/2015, a las 17:47, Xujun Zhao escribió:<br>
<br>
> Hi Jose,<br>
><br>
> Thank you for your reply. How about the computational cost compared to one EPSSolve() with all eigenvalues? what methods does SLEPc use for each solve? Because it may be cheaper for largest eigenvalue if the power method is used, but I don't if it is still so for smallest eigenvalue?<br>
><br>
> Xujun<br>
><br>
<br>
</span>Don't compute all eigenvalues.<br>
<br>
For the largest eigenvalue, don't use the power iteration. The default solver (Krylov-Schur) will be very fast for that. For the smallest eigenvalue, convergence may be slow if eigenvalues are small and poorly separated - it may be necessary to do shift-and-invert, in which case the cost may blow up.<br>
<span><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>
</div></div></blockquote></div><br></div>