Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
This community is designed for sharing of public information. Please do not share Intel or third-party confidential information here.
6590 Discussions

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

Beginner
121 Views

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?

2 Replies
Employee
121 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

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