The theorem on which this
paper is based is an easy generalization of the fact that the nullspace of an operator
is the orthocomplement of (the closure of) the range of its adjoint. Its significance,
here, is the observation that this may be applied to give a computationally feasible
algorithm for the problem of the title.