<div dir="ltr">Solver details:<div><br></div><div><div>EPS Object: 1 MPI processes</div><div> type: jd</div><div> search subspace is orthogonalized</div><div> block size=1</div><div> type of the initial subspace: non-Krylov</div><div> size of the subspace after restarting: 6</div><div> number of vectors after restarting from the previous iteration: 1</div><div> problem type: generalized non-symmetric eigenvalue problem</div><div> extraction type: harmonic Ritz</div><div> selected portion of the spectrum: smallest eigenvalues in magnitude</div><div> number of eigenvalues (nev): 1</div><div> number of column vectors (ncv): 18</div><div> maximum dimension of projected problem (mpd): 18</div><div> maximum number of iterations: 10000</div><div> tolerance: 0.0001</div><div> convergence test: relative to the eigenvalue</div><div>BV Object: 1 MPI processes</div><div> type: svec</div><div> 18 columns of global length 4225</div><div> vector orthogonalization method: classical Gram-Schmidt</div><div> orthogonalization refinement: if needed (eta: 0.7071)</div><div> block orthogonalization method: Gram-Schmidt</div><div> doing matmult as a single matrix-matrix product</div><div>DS Object: 1 MPI processes</div><div> type: gnhep</div><div>ST Object: 1 MPI processes</div><div> type: precond</div><div> shift: 0.</div><div> number of matrices: 2</div><div> all matrices have different nonzero pattern</div><div> KSP Object: (st_) 1 MPI processes</div><div> type: bcgsl</div><div> Ell = 2</div><div> Delta = 0</div><div> maximum iterations=10000, initial guess is zero</div><div> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.</div><div> left preconditioning</div><div> using PRECONDITIONED norm type for convergence test</div><div> PC Object: (st_) 1 MPI processes</div><div> type: jacobi</div><div> linear system matrix = precond matrix:</div><div> Mat Object: () 1 MPI processes</div><div> type: seqaij</div><div> rows=4225, cols=4225</div><div> total: nonzeros=37249, allocated nonzeros=37249</div><div> total number of mallocs used during MatSetValues calls =0</div><div> not using I-node routines</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 16, 2017 at 2:12 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"><div class="HOEnZb"><div class="h5"><br>
> El 16 ago 2017, a las 21:36, Kong, Fande <<a href="mailto:fande.kong@inl.gov">fande.kong@inl.gov</a>> escribió:<br>
><br>
> Hi All,<br>
><br>
> How to understand the following messages:<br>
><br>
> 1 EPS nconv=0 first unconverged value (error) 2.06312 (3.29164033e-01)<br>
> 2 EPS nconv=0 first unconverged value (error) 2.03951 (1.76223074e-01)<br>
> 3 EPS nconv=0 first unconverged value (error) 2.01177 (5.71109559e-02)<br>
> 4 EPS nconv=0 first unconverged value (error) 2.01042 (4.84609300e-02)<br>
> 5 EPS nconv=0 first unconverged value (error) 2.00708 (3.19917457e-02)<br>
> 6 EPS nconv=0 first unconverged value (error) 2.00595 (2.62792109e-02)<br>
> 7 EPS nconv=0 first unconverged value (error) 2.00504 (2.13766150e-02)<br>
> 8 EPS nconv=0 first unconverged value (error) 2.00441 (1.85066774e-02)<br>
> 9 EPS nconv=0 first unconverged value (error) 2.00397 (1.73188449e-02)<br>
> 10 EPS nconv=0 first unconverged value (error) 2.00366 (1.54528517e-02)<br>
> 11 EPS nconv=0 first unconverged value (error) 2.00339 (1.32215899e-02)<br>
> 12 EPS nconv=0 first unconverged value (error) 2.00316 (1.32215899e-02)<br>
> 13 EPS nconv=0 first unconverged value (error) 2.00316 (1.17928920e-02)<br>
> 14 EPS nconv=0 first unconverged value (error) 2.00297 (1.04964387e-02)<br>
> 15 EPS nconv=0 first unconverged value (error) 2.0028 (9.58244972e-03)<br>
> 16 EPS nconv=0 first unconverged value (error) 2.00268 (9.06634973e-03)<br>
> 17 EPS nconv=0 first unconverged value (error) 2.00198 (3.43444441e-04)<br>
> 18 EPS nconv=1 first unconverged value (error) 2.25718 (1.79769313e+308)<br>
> 18 EPS converged value (error) #0 2.00197 (5.69451918e-09)<br>
><br>
><br>
> When the solver converged, the wrong eigenvalue and the wrong residual are printed out. Do we design like this way?<br>
><br>
> Fande,<br>
<br>
</div></div>Is this the POWER solver? Most solvers in EPS approximate several eigenvalues simultaneously, but this is not the case in POWER - when one eigenvalue converges there is no approximation available for the next one.<br>
<br>
I will think about a simple fix.<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>