[petsc-users] eigensolution error with slepc

[Manav Bhatia]

Following is the output of -eps_view from ex7 and from my code. I have included only what is different. 
There are differences in nev, ncv, mpd, etc.


-Manav



output from ex7:


EPS Object: 1 MPI processes
  type: arnoldi
  problem type: generalized symmetric eigenvalue problem
  selected portion of the spectrum: largest eigenvalues in magnitude
  postprocessing eigenvectors with purification
  number of eigenvalues (nev): 1
  number of column vectors (ncv): 16
  maximum dimension of projected problem (mpd): 16
  maximum number of iterations: 768
  tolerance: 1e-08
  convergence test: relative to the eigenvalue
BV Object: 1 MPI processes
  type: svec
  17 columns of global length 6150
  vector orthogonalization method: classical Gram-Schmidt
  orthogonalization refinement: if needed (eta: 0.7071)
  block orthogonalization method: Gram-Schmidt
  non-standard inner product
  Mat Object:   1 MPI processes
    type: seqaij
    rows=6150, cols=6150
    total: nonzeros=314892, allocated nonzeros=314892
    total number of mallocs used during MatSetValues calls =0
      using I-node routines: found 2050 nodes, limit used is 5
  doing matmult as a single matrix-matrix product


output from my code: 

EPS Object: 1 MPI processes
  type: arnoldi
  problem type: generalized symmetric eigenvalue problem
  selected portion of the spectrum: largest eigenvalues in magnitude
  postprocessing eigenvectors with purification
  number of eigenvalues (nev): 20
  number of column vectors (ncv): 100
  maximum dimension of projected problem (mpd): 100
  maximum number of iterations: 10000
  tolerance: 1e-05
  convergence test: relative to the eigenvalue
BV Object: 1 MPI processes
  type: svec
  101 columns of global length 6150
  vector orthogonalization method: classical Gram-Schmidt
  orthogonalization refinement: if needed (eta: 0.7071)
  block orthogonalization method: Gram-Schmidt
  non-standard inner product
  Mat Object:   1 MPI processes
    type: seqaij
    rows=6150, cols=6150
    total: nonzeros=314892, allocated nonzeros=314892
    total number of mallocs used during MatSetValues calls =0
      using I-node routines: found 2050 nodes, limit used is 5
  doing matmult as a single matrix-matrix product






> 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.
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160404/df7d7dc5/attachment.html>