Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- is there any MKL function can outputs L matrix of cholesky decomposition

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

lin__ziije

Beginner

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

12-02-2019
10:53 PM

145 Views

is there any MKL function can outputs L matrix of cholesky decomposition

I want to solve some linear equations（Ax=b） with same matrix A and different vector b. So i decided to decompose matrix A with cholesky decomposition to get matrix L (A= LT*L).

Is there any MKL function can outputs L matrix of cholesky decomposition, i have tried some other library, but most of them can not directly output matrix L.

Link Copied

5 Replies

Kirill_V_Intel

Employee

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

12-03-2019
04:15 PM

145 Views

Hello,

Is the matrix sparse or dense? Intel MKL provides a way in both cases to solve one system with multiple right-hand sides but functionality is different. If you just want to solve with multiple right-hand sides, you don't actually need the matrix L. So, is there something else which you didn't mention in your description?

Thanks,

Kirill

lin__ziije

Beginner

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

12-03-2019
08:06 PM

145 Views

Voronin, Kirill (Intel) wrote:Hello,

Is the matrix sparse or dense? Intel MKL provides a way in both cases to solve one system with multiple right-hand sides but functionality is different. If you just want to solve with multiple right-hand sides, you don't actually need the matrix L. So, is there something else which you didn't mention in your description?

Thanks,

Kirill

thanks for your answer. Matrix A is a 5000x5000 matrix with about 60,000 (0.25%) none zeros elements. I actually need the matrix L for follow-up work.

so， MKL can not output matrix L driectly?

Kirill_V_Intel

Employee

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

12-03-2019
08:25 PM

145 Views

I see. For such a tiny matrix and such sparsity I believe there is no need to use sparse matrices. Have a look at https://software.intel.com/en-us/onemkl-developer-reference-c-getrf from LAPACK.

Best,

Kirill

lin__ziije

Beginner

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

12-03-2019
10:41 PM

145 Views

Voronin, Kirill (Intel) wrote:I see. For such a tiny matrix and such sparsity I believe there is no need to use sparse matrices. Have a look at https://software.intel.com/en-us/onemkl-developer-reference-c-getrf from LAPACK.

Best,

Kirill

I will try it, by the way ,have you ever try to solve problem like mine, do you know how long will take to do cholesky or LU decomposition. my program is run on a server with 2 Intel(R) Xeon(R) CPU E5-2609 v3 @ 1.90GHz.

Gennady_F_Intel

Moderator

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

12-05-2019
08:15 PM

145 Views

- 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.