<div>Looks good but it only gives 4 eigenvalues.  How do I get the 50 smallest eigenvalues?<br></div><div><br></div><div>---John</div><br><div class="gmail_quote">On Tue, Aug 2, 2011 at 3:40 AM, Jose E. Roman <span dir="ltr">&lt;<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
El 02/08/2011, a las 09:24, John Chludzinski escribió:<br>
<div><div class="h5"><br>
&gt; Jose,<br>
&gt;<br>
&gt; I used:<br>
&gt;<br>
&gt; MatCreateSeqAIJ(PETSC_COMM_SELF, n, n, 0, cnt, &amp;A);<br>
&gt; MatSetValue(A, i, j, *(a+i*SIZE+j), INSERT_VALUES);<br>
&gt;<br>
&gt; to create the input matrices.  Their sizes collapsed dramatically: from 192208072 to 239944. But again, using: ./ex7.exe -f1 k.dat -f2 m.dat -eps_gen_hermitian -eps_smallest_real &gt; x.out 2&gt;&amp;1, I get:<br>
&gt;<br>
&gt; Generalized eigenproblem stored in file.<br>
&gt;<br>
&gt; Reading REAL matrices from binary files...<br>
&gt; Number of iterations of the method: 500<br>
&gt; Number of linear iterations of the method: 4009<br>
&gt; Solution method: krylovschur<br>
&gt;<br>
&gt; Number of requested eigenvalues: 1<br>
&gt; Stopping condition: tol=1e-07, maxit=500<br>
&gt; Number of converged approximate eigenpairs: 0<br>
&gt;<br>
&gt; I&#39;m really only interested in the least 50 eigenvalues.<br>
&gt;<br>
&gt; Smallest 10 eigenvalues<br>
&gt; -0.000019<br>
&gt; -0.000004<br>
&gt; -0.000000<br>
&gt; 0.000005<br>
&gt; 0.000012<br>
&gt; 0.000016<br>
&gt; 2.795284<br>
&gt; 2.795307<br>
&gt; 21.235339<br>
&gt; 21.235340<br>
&gt; 81.582017<br>
&gt;<br>
&gt; Largest 10 eigenvalues<br>
&gt; 176431487319.625000<br>
&gt; 176431532012.467468<br>
&gt; 176431562378.361359<br>
&gt; 176435480628.136292<br>
&gt; 176435488209.944031<br>
&gt; 176435555689.747253<br>
&gt; 176435563270.922424<br>
&gt; 663312473823.916260<br>
&gt; 663312473823.917969<br>
&gt; 663312666285.928589<br>
&gt; 663312666285.929810<br>
&gt;<br>
<br>
</div></div>You will never get convergence to the small eigenvalues in this way, since the largest ones are huge. Try something like:<br>
./ex7.exe -f1 k.dat -f2 m.dat -eps_gen_hermitian -st_type sinvert -eps_target 1.0<br>
<font color="#888888"><br>
Jose<br>
<br>
</font></blockquote></div><br>