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
- Unusual Feast results (dfeast_scsrgv)

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

Chris_Hickey

Beginner

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

08-06-2020
09:48 AM

195 Views

Unusual Feast results (dfeast_scsrgv)

Hi,

I've been trying to solve the generalised symmetric eigenproblem Bx = eAx for symmetric A and positive definite symmetric B, and so to work with the eigensolvers, we solve Ax = (1/e)Bx.

The 99 expected eigenvalues are 1,2,...,99.

I'm trying three solvers, Spectra, Intel MKL extremal (mkl_sparse_d_gv) and interval (dfeast_scsrgv).

I've attached an excel document with the time taken and precision.

Unexpectedly, Spectra is the fastest by far. MKL extremal seems very slow.

For MKL interval using FEAST, I get some strange results.

The interval solver doesn't seem to like it when there are eigenvalues on the boundary of the interval, but acts very strangely in test 4, for the interval [1,10].

Furthermore, the larger M0, the faster the code seems to run, which seems unusual, how should M0 be

Link Copied

3 Replies

Gennady_F_Intel

Moderator

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

08-06-2020
08:58 PM

188 Views

Chris_Hickey

Beginner

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

08-06-2020
11:57 PM

182 Views

I've only tried it on my CPU so far, which is a Intel(R) Xeon(R) CPU E3-1505M v6 @ 3.00GHz, 3000 Mhz, 4 Core(s), 8 Logical Processor(s).

The version of MKL I'm using is 2020.0.1. build 2020208 for intel 64 architecture.

Chris_Hickey

Beginner

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

08-12-2020
09:40 AM

149 Views

As requested elsewhere, I've removed the Spectra calls, although the dependency on Eigen remains for reading in the matrix market files.

Attached is the new source file.

Thanks,

Chris

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.