Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
7220 ディスカッション

Disable the permutation during the schur complement computation

Simplice_D_
ビギナー
847件の閲覧回数

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 件の賞賛
2 返答(返信)
Gennady_F_Intel
モデレーター
847件の閲覧回数

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

Simplice_D_
ビギナー
847件の閲覧回数

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 ?

返信