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.