Add function to explicitly compute the Schur complement
Currently, the function asSchurComplement provides a way to create a matrix that, when multiplied with a vector, behaves like the Schur complement. A function that wraps PETSc's MatSchurComplementComputeExplicitOperator would allow to explicitly compute that matrix e.g. for testing.