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

Thanks, Tracy

Link Copied

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

**?gemv**to compute A X, followed by a call to

**?dot**.

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

The result of:

X' * A * X

For symmetric A has such a nice structure (Result is symmetric) that it is a pity MKL doesn't have an optimized function for it.

Quadratic Forms are very common in Machine Learning, Optimization, etc... It would benefit many users.

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

In case of dense, there is no single call solution as mecej4 mentioned.

However, if X and A are sparse matrices, then we added exactly this type of product to our Inspector-Executor Sparse BLAS routines. The structure is indeed quite nice. This type of product shows up a lot in multi-scale finite element methods as well where X could be a projection or elongation matrix. We call it the symmetric product with api -- mkl_sparse_sypr().

See reference documentation mkl_sparse_sypr for more details on how to use it.

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

Blas Level 2 is a tricky thing. For example, my dsptrd faster dsptrd Intel MKL twice.

See https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/288316

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