<div dir="ltr">To close this thread let me point out the problem.<div><br></div><div>The problem occured due to the distribution of the number of processors over the total number of nodes used (as Jose E. Roman pointed out.) </div><div>The problem was solved and all the 300 (good) eigenvalues, properly converged,</div><div>were obtained when the correct distribution of the number of processors over the nodes was ascertained using <b>numactl </b>.</div><div><br></div><div>TLDR, incorrect distribution of processors using numactl may lead to bogus eigenvalues.</div><div><br></div><div>Thanks Jose E. Roman !</div><div><br></div><div>Problem is now solved.</div><div><br></div><div>With regards,</div><div> Vijay</div><div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 15, 2015 at 3:00 PM, Vijay Gopal Chilkuri <span dir="ltr"><<a href="mailto:vijay.gopal.c@gmail.com" target="_blank">vijay.gopal.c@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">Yes, those seem to be the right eigenvalues.<div>Ok so the solution is to recompile PETSc/SLEPc with a basic configuration and test with --with-debugging=1 </div><div><br></div>Would it make a difference if I use <b>-esp_type lanczos </b>or some other diagonalization procedure ?<div><br></div><div>I'll run the test with a new version of PETSc/SLEPc and report back.</div><div><br></div><div>Thanks a  lot,</div><div> Vijay</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 15, 2015 at 2:55 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><br>
El 14/05/2015, a las 19:13, Vijay Gopal Chilkuri escribió:<br>
<br>
> oups sorry, I send you a smaller one (540540) this should finish in a few minutes.<br>
><br>
> It requires the same makefile and irpf90.a library.<br>
> so just replace the old problem.c file with this and it should compile.<br>
><br>
> Thanks again,<br>
>  Vijay<br>
><br>
<br>
</span>I was able to compute 300 eigenvalues of this matrix of size 540540. All eigenvalues are in the range -4.70811 .. -4.613807, and the associated residual is always below 1e-9.<br>
<br>
There must be something in your software configuration that is causing problems. I would suggest trying with a basic PETSc/SLEPc configuration, with no openmp flags, using --download-fblaslapack (instead of MKL). Also, although it should not make any difference, you may want to try with a smaller number of MPI processes (rather than 741).<br>
<span><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div>