Severe bug in LDA (default algorithm)
The default algorithm in this package is wrong. We cannot use symmetric to evaluate the eigenvalue decomposition of Sw^-1 Sb, because that multiplication is not guaranteed symmetric.
We should always be using pinv. There should be no other option.
References: