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
- Linear System Solving For Tiny Matrix

- 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

04-20-2017
02:36 AM

36 Views

Linear System Solving For Tiny Matrix

We are looking for to solve Linear Systems of tiny complex matrix like 10x10 in a very fast way.

We would like to use iterative methods because we have to solve a sequence with about 100 linear systems and solution of each system is a very good initial conditions for next system.

We've verified the direct solution method in the last MKL 2017 Update 2 but the speed is not enough.

I believe iterative methods could be the right solution.

What do you suggest?

Thank you

Gianluca

Link Copied

3 Replies

Alexander_K_Intel2

Employee

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

04-20-2017
04:07 AM

36 Views

Hi,

Quick question - is your system sparse or dense?

Thanks,

Alex

Gianluca_G_1

Beginner

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

04-20-2017
06:30 AM

36 Views

System is dense up to size 4x4 and sparse for bigger size systems

Alexander_K_Intel2

Employee

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

04-20-2017
07:07 AM

36 Views

if your systems don't change significantly from system to system I suggest to factorize first system and further use it's factorization as preconditioner for further systems. For sparse system you just need to call pardiso solver with phase 33 and set new matrix as input parameter. For dense you can implement simple iterative solver and use dgetrs functionality for preconditioner handle

Thanks,

Alex

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