A Bi-level Nonlinear Eigenvector Algorithm for Wasserstein Discriminant Analysis