The function computes the Schur complement of the block structured symmetric matrix with respect to the picking matrix as follows:

Let be a given picking matrix

and let be the complementary picking matrix (i.e., the orthogonal complement) such that . For a given symmetric matrix , the permutation operation then yields the matrix

Subsequently, the function computes the Schur complement as .