Problem with eigs()

I have a sparse symmetric positive-definite real matrix. I
am trying to find the largest and smallest eigenvalues using
eigs() routine. Here's the error I am getting:

eigs(A,1,'be')

??? Error using ==> eigs at 231
Error with ARPACK routine dsaupd: info = -13

Is this something solvable? Does it mean there's a problem
with the matrix?

I am able to estimate the condition number of this matrix
with condest(), but I would like to be able to use eigs(),
so that I can estimate the condition number of a matrix with
a preconditioner, by solving generalized eigenvalue problem,
eigs(A,M,1,'be'), but it fails the same way whether I use
specify the preconditioner or not...

Andrey
.

Relevant Pages

  • Re: Compute eigen values/vectors of a 22k x 22k matrix
    ... Any CAS program (for example Maple, Mathematica, or Maxima) should be ... and 1's (without exploiting the symmetry, ... compute the eigenvalues would be. ... I have a routine (to find eigenvalues/eigenvectors by Jacobi ...
    (sci.math)
  • Re: Compute eigen values/vectors of a 22k x 22k matrix
    ... Any CAS program (for example Maple, Mathematica, or Maxima) should be ... and 1's (without exploiting the symmetry, ... compute the eigenvalues would be. ... I have a routine (to find eigenvalues/eigenvectors by Jacobi ...
    (sci.math)
  • Re: Random orthogonal matrix
    ... I am already aware of that but I assume that some routine specialized ... computational cost to finding eigenvalues and eigenvectors would be ... are those of the Statistics Department or of Purdue University. ... Herman Rubin, Department of Statistics, Purdue University ...
    (sci.math.num-analysis)
  • Re: Compute eigen values/vectors of a 22k x 22k matrix
    ... Any CAS program (for example Maple, Mathematica, or Maxima) should be ... and 1's (without exploiting the symmetry, ... compute the eigenvalues would be. ... I have a routine (to find eigenvalues/eigenvectors by Jacobi ...
    (sci.math)
  • Re: eigenvalue solver
    ... eig handle any kind of matrics but it computes all eigenvalues. ... It doesn't have option to find only a few eigenvalues. ... How about using eigs on ... To use eigs subroutine, B must be positive definite,but this is not ...
    (comp.soft-sys.matlab)