Help with SSYEV

Open discussion regarding features, bugs, issues, vendors, etc.

Help with SSYEV

Postby vibhav » Fri Sep 08, 2006 10:09 pm

Hi All,

I am new to this forum and fortran. I am having a simple problem using lapack.

I am trying to solve for eigen-vector and value for a real symmetric matrix (positive definite). I am using SSYEV/DSYEV function.

For my matrix I am suppose to get all the positive eigen-values. But I some of my eigenvalues are negative.

Can anyone help me out with this or point me in the right direction. Is there other functions which is better than SSYEV?

Thanks

Vibhav
vibhav
 
Posts: 1
Joined: Fri Sep 08, 2006 10:00 pm

Postby Julien Langou » Sun Sep 10, 2006 8:35 pm

Hello

how much "negative" are your eigenvalues?

If the absolute value of your negative eigenvalue is smaller than machine epsilon times the norm of A, you should consider this eigenvalue to be zero. (Which make your matrix numerically not definite so you still might not be happy ...)

If you are looking for more accurate symmetric eigensolvers, you should give a try to xSYEVX, but this will not help that much. The main problem is that with the current LAPACK solvers you loose componentwise accuracy when we go to tridiagonal form. You are normwise accurate. Some new solvers are developped right now to remedy this problem.

Julien Langou
Julien Langou
 
Posts: 832
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 6 guests