Subroutine Question: DSYEV; Eigensolver for Symmetric Matrix

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

Subroutine Question: DSYEV; Eigensolver for Symmetric Matrix

Postby DavidB » Wed Jan 11, 2006 1:58 am

Hello.

Is this forum the place for posting questions regarding specific routines, or should questions be sent to a specific individual?

My question:

I am looking over the code for DSYEV, a driver routine for solving for the eigenvalues and eigenvectors of a real, symmetric matrix. One of the parameters passed into DSYEV is UPLO. I am not clear on the purpose of this variable. The comments state the following:

“If UPLO = ‘U’, the leading NxN upper triangular part of A contains the upper triangular part of the matrix A. If UPLO = ‘L’, the leading NxN lower triangular part of A contains the lower triangular part of the matrix A. “

So the upper part of the matrix contains the upper part of the matrix, and the lower part of the matrix contains the lower part of the matrix?!?! That would seem obvious, so maybe I am just misunderstanding this variable. Somebody please clarify for me.
DavidB
 
Posts: 34
Joined: Thu Dec 15, 2005 1:49 pm
Location: Vancouver, B.C. Canada

Postby sven » Wed Jan 11, 2006 5:55 am

The point is that you do not need to supply anything in the strictly lower triangular part of A when UPLO = 'U', or in the strictly upper triangular part of A when UPLO = 'L'. See Section 5.3.1 of the LAPACK Users' Guide:

http://www.netlib.org/lapack/lug/node122.html

Sven Hammarling.
sven
 
Posts: 144
Joined: Wed Dec 22, 2004 4:28 am

Postby DavidB » Wed Jan 11, 2006 1:26 pm

Thank-you.

That clears things up. (I had downloaded DSYEV and copied the code to a Word document--which was over 130 pages long! So I printed off ONLY the code for the DSYEV sub-routine itself.) The User Guide has now been bookmarked and I will definitely be making use of it in the future.
DavidB
 
Posts: 34
Joined: Thu Dec 15, 2005 1:49 pm
Location: Vancouver, B.C. Canada


Return to User Discussion

Who is online

Users browsing this forum: Bing [Bot], Google [Bot] and 1 guest