We recommend to people that they just call MatSetNullSpace() and solve the singular system. What if reductions are very expensive? Are there Krylov methods that accommodate a pipelined or deferred null space removal? Shouldn't be hard to do, but I haven't seen papers on it.