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
- Which routine should be used for A*(R^{-1})?

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

utab

Beginner

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

02-05-2013
12:30 AM

55 Views

Which routine should be used for A*(R^{-1})?

Dear all,

For a transformation operation, I needed to compute A*(R^{-1}) where A is a rectangular matrix. A has m rows, which is much larger than the column size, n. In general, the column size is between 2 and 10 and R is an upper triangular square matrix of size n.

A is meant to be a block of iteration vectors in my code. However, I have to do the inversion from right, an operation that is possible in MATLAB like 'A/R', is there way to achieve this directly or should I use the inverse of R which is less likely I suppose?

Could you please direct me on this matter?

Best regards,

Umut

Link Copied

3 Replies

mecej4

Black Belt

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

02-05-2013
08:35 AM

55 Views

utab

Beginner

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

02-06-2013
06:57 AM

55 Views

mecej4 wrote:

well that might be but that is not exactly what I am looking for...

most probably, what I want to do should be be efficiently accomplished by inverting the upper triangular matrix R first and them multiply that from the right. Since R is a small matrix doing something like

R^{-1} I = inverse_of_R

then

A*inverse_of_R

I am suspecting whether if i will experience some accuracy issues or not in this case with the R^{-1}I linear system solutions...

SergeyKostrov

Valued Contributor II

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

02-06-2013
07:43 AM

55 Views

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.