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 Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Basic question regarding matrix inversion vs solver

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

Highlighted

Mads_Trolle

Beginner

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

01-25-2016
01:56 AM

17 Views

Dear all,

Given a matrix A, I need to calculate element (1,1) of its inverse... i.e. if B = inv(A), i need B(1,1), and not the entire matrix B.

Hence, doing a complete inversion of matrix A seems a bit excessive... Does mkl contain any smart solver-based routines i might use...?

My first impression is to use

call getrf(A,i)

call getrs(A,I)

where I = [1,0,0,...], to get the first column of B...

Best,

Mads

Accepted Solutions

Highlighted

Ying_H_Intel

Employee

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

01-26-2016
08:04 PM

17 Views

Hi Mads,

You are right. If the matrix is dense, then I’d suggest to use ?gesv(A,I) as it is a bit simpler than using factorization and solve routines. If the matrix A is sparse, it is possible to use PARDISO solver with right-hand side “I” and partial solve to get just B(1,1).

Best

Ying

2 Replies

Highlighted

Ying_H_Intel

Employee

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

01-26-2016
08:04 PM

18 Views

Hi Mads,

You are right. If the matrix is dense, then I’d suggest to use ?gesv(A,I) as it is a bit simpler than using factorization and solve routines. If the matrix A is sparse, it is possible to use PARDISO solver with right-hand side “I” and partial solve to get just B(1,1).

Best

Ying

Highlighted

Mads_Trolle

Beginner

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

01-29-2016
02:03 AM

17 Views

Dear Ying,

Thank you! Exactly what I was looking for :)

Best,

Mads

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