Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

Disable the permutation during the schur complement computation

Simplice_D_
Beginner
601 Views

Hi All,

I am using Pardiso from MKL to compute the Schur complement of a matrix, and I realize that it permutes the matrix during the reordering and numerical factorization.

Is there a way to disable that permutation ? I am doing a partial factorization of a submatrix involved in a more large project, and I just want to compute the Schur complement without any ordering or numerical pivoting.

Thanks

 

0 Kudos
2 Replies
Gennady_F_Intel
Moderator
601 Views

yes, you are right, that's not possible right now however we are planning to implement this feature into one of the future versions.

-- Gennady

0 Kudos
Simplice_D_
Beginner
601 Views

Thanks for your answer Gennady.

By assuming the permutation, in which format is the Schur complement matrix stored ? I have tried both dense row-major and column-major layout, and I do not get the same Schur complement as in my Matlab code. Some entries are incorrect. 

Is there any reported bug about that functionality ?

0 Kudos
Reply