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 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 ?