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
- PARDISO-Solving a large matrix containing many identical nonzero sparsity structure submatrices

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

Long_N_2

Beginner

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

04-30-2015
01:26 AM

47 Views

PARDISO-Solving a large matrix containing many identical nonzero sparsity structure submatrices

Hi friends,

I am solving coupled problem analysis in FEM. I want to solve a matrix A.x = b, A could be

A = [A11, A12; A21, A22] and B has a size of [A11, A12]

A could contain 3x3, 4x4...nxn submatrices. The submatrices have the same nonzero sparsity structure (IA and JA). At the moment the solution may be reordering IA and JA from submatrices to a large matrix A to input into Pardiso Solver. However it will increase the storage memory for IA and JA.

It there any way to solve matrix A without reodering IA and JA? As I know multi right hand side method solves the only A matrix with different rhs, it can not apply for this case. Can anyone give me a hint?

Link Copied

2 Replies

Alexander_K_Intel2

Employee

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

04-30-2015
01:58 AM

47 Views

Hi,

Probably i didn't get your question clearly... You wrote that matrix 2x2 blocked but further said that it is not generally 2x2 blocked, but number of blocks could be bigger - am i correct?

Thanks,

Alex

Long_N_2

Beginner

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

04-30-2015
03:36 AM

47 Views

Hi Alex,

2x2 blocked matrix is an example for my case, let us solve the problem of matrix A has 2x2 block. If we solve this broblem, we can also solve nxn blocked matrix in the same way.

Thanks for your considering

Long

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.