[petsc-users] eigensolution error with slepc

[Manav Bhatia]

It turns out that the matrices that I sent you was from a load step before the one that failed. 

I am calling eps_solve once for each load step, and I get the error at the 14th load step. So, the matrices that I sent were from step 13. 

The eigenvalues you got match the 1/lambda that I get from my code at load step 13. 

I will try to output the offending matrices and share those with you. 

Regards,
Manav


> On Apr 4, 2016, at 9:47 AM, Jose E. Roman <jroman at dsic.upv.es> wrote:
> 
> 
>> El 4 abr 2016, a las 16:42, Manav Bhatia <bhatiamanav at gmail.com> escribió:
>> 
>> 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. 
>> 
>> 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.
>> 
>> Not sure how to go about debugging this. 
>> 
>> -Manav
> 
> What is different in your code with respect to ex7? Run with -eps_view and compare both.
>