topic Thanks for your answer in Intel® oneAPI Math Kernel Library

Disable the permutation during the schur complement computation

Simplice_D_ — Wed, 28 Jun 2017 13:34:29 GMT

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

yes, you are right, that's

Gennady_F_Intel — Wed, 28 Jun 2017 18:06:42 GMT

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

Thanks for your answer

Simplice_D_ — Wed, 05 Jul 2017 17:38:10 GMT

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 ?