LDA.cpp eigennized

I have changed the LAPACK codes into EIGEN3 for the binary classification through LDA.
Question: As in PCA, there are probably multiple ways to do this. Computing the covariance matrix is fine if data is not high-dimensional, but if it is, then computing matrix is impossible. Could you investigate a bit and come up with a fix for that case?