Do we really need a new routine for this? Why doesn't MatSetNullSpace() work?