Turn on suggestions

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

Showing results for

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

pennisi__giovanni

Beginner

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

12-15-2016
09:17 AM

74 Views

Feast RCI

Hello,

I have to solve the generalized eigen problem `Ax` = `λ` `Bx` , where A is symmetric, positive-definite, and B is symmetric and not positive definite. I read that the FEAST algorithm can solve only the eigen system where B is positive definite. So I tried to solve the reverse problem, swapping A and B, and solving for 1/λ, but then the range of eigenvalues gets inverted, and I need to find the largest eigenvalues instead of the smallest. Unfortunately the matrix A is big and the algorithm takes too much time to find the eigenvalues and the eigenvectors, for the reverse problem! I understood that is possible to use the feast_rci to solve the direct problem where B can be not positive definite, am I correct? If yes where I can find a simple example that uses PARDISO for factorization and for the other operations required by the feast_rci? I've already read this page https://software.intel.com/en-us/node/521741#24EF00BB-984C-442C-889B-58B7C042C86C but if you can provide a simple example that uses PARDISO or other mkl function, it could be easier for me to understand how to solve my problem!

Thanks in advance,

Giovanni.

Link Copied

7 Replies

Jing_X_Intel

Employee

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

12-15-2016
04:50 PM

74 Views

Hi

sample codes of Feast and PARDISO can be found in the example folder of mkl ($(MKLROOT)\examples\examples_core_c.zip\solvers_eec\source).

pennisi__giovanni

Beginner

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

12-28-2016
02:43 AM

74 Views

Hi Jing,

Thanks for your reply, I studied the example code you suggested and I adapted it to my code, but due to the fact tha Feast RCI algorithm still depends on the interval emin-emax, if the upper bound (emax) is too far from the grater eigeinvalue in that interval, the value of this eigenvalue is very different from the expected one . Instead if a choose a narrower interval where the upper bound is very close to the last searched eigenvalue the result is correct! Is there a way to estimate the upper bound of the search interval?

Thanks in advance,

Giovanni.

pennisi__giovanni

Beginner

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

12-28-2016
07:26 AM

74 Views

Doing some more test, it seems that the FEAST RCI algorithm finds the greatest eigenvalues instead of the lowest one for every interval passed. If the interval is narrower enough to contain only the searched eigenvalues, obviously the values returned are correct!

Do you have any suggestions?

Gennady_F_Intel

Moderator

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

12-28-2016
07:33 PM

74 Views

pennisi__giovanni

Beginner

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

12-29-2016
12:19 AM

74 Views

Hi Gennady,

Thanks for your reply, I'm using Intel MKL Math library 11.3 and I think that includes only FEAST 2.1. If I am correct, where can I download the experimental package?

Thanks in advance,

Giovanni.

Gennady_F_Intel

Moderator

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

12-29-2016
06:00 AM

74 Views

yes, it will work with version 11.3 ( the current is MKL 2017) of Intel MKL.

pennisi__giovanni

Beginner

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

01-04-2017
07:53 AM

74 Views

Hi Gennady F.,

Where can I find the documentation that describes the use of FEAST 3.0? Or Can you give me a small example that explains the estimation of eigenvalues in a search interval?

Thanks in advance,

Giovanni.

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.