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
- FGMRES for pentadiagonal 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

Elsayed_A_

Beginner

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

10-23-2016
10:29 AM

78 Views

FGMRES for pentadiagonal matrix

Hello,

I am new to Fortran intel MKL. I am trying to use the FGMRES with ILUT preconditioner for solving the Penta diagonal matrix (A*x=b.) resulting from discretization of 3D convection-diffusion equation. The coefficient matrix is naturally stored in a band format.

Since ILUT is using CSR storage, how can I convert the band storage format to CSR format?

Also, for matrix-vector multiplication step, for my case do I have to multiply inv(A)*b? I am confused about this.

Regards,

Elsayed

Link Copied

3 Replies

Gennady_F_Intel

Moderator

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

10-23-2016
08:50 PM

78 Views

Elsayed_A_

Beginner

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

10-23-2016
08:57 PM

78 Views

Gennady,

Thanks for your prompt reply. Actually, this is what I am doing. My question then is how to convert from band to dense representation.

Gennady_F_Intel

Moderator

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

10-24-2016
02:40 AM

78 Views

at this moment, MKL doesn't provide such sort of conversion routine

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.