Adjust by sqrt(chi) so as to work over the reals, and use GeneralizedSelfAdjointEigenSolver