<div dir="ltr">Hi Jose,<div><br></div><div>This is only a very simple test. 2 outer iterations perform 24 MatMult operation, which means I have to solve a Finite element system with millions of DOFs for 24 times. This is only for a polymer system with two beads. For a polymer chain with hundreds of beads and springs, I need more outer iterations to get its max/min eigenvalues, which means I have to perform more computations of the finite element system. Those solves happen only for one time step. When polymer beads move, I need to redo everything. </div><div><br></div><div>In other words, the eigenvalues (max and min) have to be solved every time step or every several time steps. Therefore, it is very expensive! My question is that if It is possible to use the eigenvalues obtained from last time step as an initial guess to accelerate the convergence. Thank you for your help.</div><div><br></div><div>Xujun</div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 16, 2015 at 4:18 PM, 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 class=""><br>
> El 16/9/2015, a las 21:44, Xujun Zhao <<a href="mailto:xzhao99@gmail.com">xzhao99@gmail.com</a>> escribió:<br>
><br>
> The following is another test, which shows the outer iteration is 2, and MatMult is performed 24 times. So it means the first outer iteration performs 16, and the subsequent iterations performs 8. (24=16+8)<br>
><br>
> What parameters can I tune to accelerate the convergence? Will changing ncv work better?<br>
> The following nconv=1 and unconverged value, what do they mean?<br>
><br>
> 2 EPS nconv=1 first unconverged value (error) 0.0140231 (1.32666590e-02)<br>
><br>
> --->test: EPS solve results:<br>
>           the eigenvalue: 0.009301<br>
>           Solution method: krylovschur<br>
>           Number of iterations of the method: 2<br>
>           EPS converged reason : 2<br>
<br>
</span>2 outer iterations is very fast, your problem seems very easy. You can try changing ncv but I would bet you won’t get the result with much less than 24 MatMult’s. Just try.<br>
<br>
The monitor you are using shows the number of converged eigenvalues and the first non-converged Ritz value.<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>