fix(lda): raise PartialEq tolerance from 1e-10 to 1e-6 for f32 safety

1e-10 rounds to 0.0 in f32 (min positive normal ~1.2e-7), making
(a - b).abs() > 0.0 almost always true for distinct values. The 1e-6
floor works correctly for both f32 and f64.

Addresses review feedback on #457.

Author: Mec-iS

src/decomposition/lda.rs

@@ -71,7 +71,7 @@ pub struct LDA<T: Number + RealNumber, X: Array2<T> + EVDDecomposable<T>> {
 
 impl<T: Number + RealNumber, X: Array2<T> + EVDDecomposable<T>> PartialEq for LDA<T, X> {
     fn eq(&self, other: &Self) -> bool {
-        let tol = T::from(1e-10).unwrap();
+        let tol = T::from(1e-6).unwrap();
         if self.n_features != other.n_features
             || self.eigenvalues.len() != other.eigenvalues.len()
             || self