Once you have called this routine you can call KSPSetOperators() on the resulting ksp to provide the operator for the Galerkin problem,
an alternative is to use PCGalerkinSetComputeSubmatrix() to provide a routine that computes the submatrix as needed.