<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Apr 4, 2016 at 9:42 AM, Manav Bhatia <span dir="ltr"><<a href="mailto:bhatiamanav@gmail.com" target="_blank">bhatiamanav@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ok. So, I ran ex7 with the same command-line options in your email, and got a result. This is on my Mac OS X, without any changes to the lapack/blas/slepc library.<br>
<br>
I also ran my code on a linux machine with lapack/blas build from source (obtained from netlib), and got the same error as on my mac.<br>
<br>
Not sure how to go about debugging this.<br></blockquote><div><br></div><div>Sounds like a memory error. Run with valgrind.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
-Manav<br>
<br>
<br>
> On Apr 4, 2016, at 3:46 AM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>> wrote:<br>
><br>
><br>
>> El 3 abr 2016, a las 22:17, Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com">bhatiamanav@gmail.com</a>> escribió:<br>
>><br>
>> I just send you the matrices.<br>
>><br>
>> Thanks,<br>
>> Manav<br>
><br>
> I cannot reproduce the problem. I was able to compute the eigenvalues without problems with Krylov-Schur (note that I scaled matrix A by 1e7 because the eigenvalues are tiny):<br>
><br>
> $ ./ex7 -f1 A.petsc -f2 B.petsc -eps_gen_hermitian<br>
><br>
> Generalized eigenproblem stored in file.<br>
><br>
> Reading REAL matrices from binary files...<br>
> Number of iterations of the method: 1<br>
> Number of linear iterations of the method: 17<br>
> Solution method: krylovschur<br>
><br>
> Number of requested eigenvalues: 1<br>
> Stopping condition: tol=1e-08, maxit=768<br>
> Linear eigensolve converged (6 eigenpairs) due to CONVERGED_TOL; iterations 1<br>
> ---------------------- --------------------<br>
>            k             ||Ax-kBx||/||kx||<br>
> ---------------------- --------------------<br>
>       11.410228            3.18851e-09<br>
>        1.323631            2.33005e-09<br>
>        1.323631            3.98884e-09<br>
>        0.224497            1.24334e-08<br>
>        0.168344            9.81685e-08<br>
>        0.129229            4.85372e-08<br>
> ---------------------- --------------------<br>
><br>
> Note that I am solving the problem as a GHEP, and this gives far better residuals than solving as a GNHEP. The Arnoldi solver did not fail either.<br>
><br>
> I don't know what caused your error. You can try configuring with the source Lapack: --download-f2cblaslapack<br>
><br>
> Jose<br>
><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>