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.
Sw^-1 Sb
We should always be using pinv. There should be no other option.
pinv
References:
Start: None
Due: None