- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Did anybody meet a striking disagreement between the output data from
LAPACK routines and other library routines when calculating
eigenvectors of a symmetric matrix?
I took several symmetric 3x3 matrices and used 2+3 subroutines from
LAPACK in various combinations to calculte eigenproblem for a matrix.
They are DSTEQR, DSYTD2 and DSTEDC, DSTEV, DSTEVD.
(1) Running any combination gives the same output data.
Which is great!
(2) But calculated eigenvectors are wrong!
A simple check of their orthoganality properties confirms this.
(3) Moreover, I processed the same matrices with NAG-77 library
and with MUPAD. In both cases I obtained the same correct results
without any problems.
LAPACK I used:
version 3.0 dated by 1999
Some matrix I used:
BRR=RESHAPE( (/10., 2., 3., 2., 10., 4., 3., 4., 10. /), (/ 3,3 /) )
Is anybody able to comment my post?
Link Copied
0 Replies
![](/skins/images/3344F5B3B76C91485ED0E980FD0CA95E/responsive_peak/images/icon_anonymous_message.png)
Reply
Topic Options
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page