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 .