Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- a generalized eigenvalue problem question

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

utab

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

03-05-2012
06:36 AM

60 Views

a generalized eigenvalue problem question

To solve a generalized eigenvalue problem of the kind

AX = lambda BX

, say I use dpotrf -> dsygst -> dsyevd path, then the eigenvectors should be transformed back, right?

Because this is the definition when it is converted to a standart problem. And the vectors are not directly the vectors of the original generalized problem. However, I have difficulty in locating the transformation routine. Could you help me with this?

Best,

Umut

Link Copied

1 Reply

mecej4

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

03-06-2012
04:22 AM

60 Views

If you want to recover the original eigenvectors, you may need to declare extra arrays for saving copies of the transformation matrices and write extra statements to copy arrays to be saved, since Lapack routines often overwrite input matrices with result matrices.

I have not performed these manipulations myself, and the example included with MKL (dsygstx.f) does not compute eigenvectors. Therefore, I cannot tell you in detail what to do. What one needs to do is to sit down with pencil and paper, read the documentation, and make notes as to how to manage and carry out the work. The coding is then usually quite straightforward.

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

For more complete information about compiler optimizations, see our Optimization Notice.