Make `numerical` return for the $BB$ subspace

In its current form, the function numerical returns H_tildes that return errors for the BB subspace for orders >1. Since U and U^\dagger both return scipy.sparse.linalg.LinearOperators in these settings, it makes sense for H_tilde to produce a similar output for consistancy.