Fortran90: Subroutine DSYEV and associating eigenvalues and eigenvectors.

Click For Summary
SUMMARY

The discussion focuses on using the LAPACK (Linear Algebra Package) subroutine DSYEV to compute eigenvalues and eigenvectors of a large symmetrical real matrix. It confirms that DSYEV outputs eigenvalues in ascending order and that the corresponding eigenvectors are stored in the matrix A, with each column representing an eigenvector associated with the respective eigenvalue. Specifically, setting JOBV = "V" ensures that the eigenvectors are returned in the matrix A.

PREREQUISITES
  • Understanding of LAPACK (Linear Algebra Package)
  • Familiarity with the DSYEV subroutine
  • Knowledge of eigenvalues and eigenvectors
  • Basic programming skills in Fortran90
NEXT STEPS
  • Review the LAPACK documentation for DSYEV at Netlib
  • Explore the NAG documentation for further insights on LAPACK functions
  • Learn how to visualize eigenvalues and eigenvectors using graphing libraries
  • Investigate other LAPACK subroutines for different matrix types and operations
USEFUL FOR

This discussion is beneficial for numerical analysts, software developers working with linear algebra, and researchers needing to compute and visualize eigenvalues and eigenvectors using Fortran90 and LAPACK.

Animastryfe
Messages
80
Reaction score
0
Greetings. I am using the LAPACK (Linear Algebra Package) software package to find the eigenvalues and eigenvectors of a large symmetrical real matrix. Specifically, I calculate a scalar from each eigenvector, and I want to graph it against its associated eigenvalue.

I am using the subroutine DSYEV of LAPACK to do this. However, DSYEV outputs the eigenvalues in ascending order, and I'm not sure how it orders the eigenvectors. Is there a way to associate each eigenvector with its eigenvalue?

Edit: The official page for DSYEV is here: http://www.netlib.org/lapack/double/dsyev.f
Here is another page about it: http://www.nag.co.uk/numeric/fl/nagdoc_fl22/xhtml/F08/f08faf.xml
 
Technology news on Phys.org
The eigenvectors are returned in the A matrix in the same order as the eigenvalues.

A(1:N, 1) is the vector for W(1), A(1:N, 2) is the vector for W(2), etc.
 
Thank you very much for the reply.

Just to be completely clear: the columns of the matrix A are the eigenvectors?
 
If you set JOBV = "V", then yes. Otherwise, no!
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 18 ·
Replies
18
Views
7K
  • · Replies 1 ·
Replies
1
Views
6K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 5 ·
Replies
5
Views
6K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 13 ·
Replies
13
Views
4K
  • · Replies 4 ·
Replies
4
Views
7K
  • · Replies 3 ·
Replies
3
Views
4K