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
- Avoid matrix inversion is it always possible?

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

Gianluca_G_1

Beginner

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

11-08-2017
08:13 AM

29 Views

Avoid matrix inversion is it always possible?

We need to solve a linear system with a dense matrix of complex numbers.

When the size of the matrix is lower than 4 GB (Microsoft limits for 64 bit systems) we can use MKL routines directly to solve:

Instead when the matrix is bigger than 4 GB we divide the matrix in blocks:

The linear system becomes:

From (1) we obtain:

**Solving equation (6) the solution x2 is obtained.**

**Using x2 in (3) it is easy to calculate x1.**

**This algorithm is interesting but introduces a matrix inversion.**

**This is a big problem when the linear system is ill conditioned.**

**Any suggestion to solve the problem avoiding the matrix inversion?**

Thankyou

Gianluca

1 Reply

Alexander_K_Intel2

Employee

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

11-08-2017
08:19 AM

29 Views

Answered on initial thread https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/748079#comment-1914725

Thanks,

Alex

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