Remove the shifting from the qr algorithm.

This commit is contained in:
Sébastien Crozet 2014-08-17 22:42:16 +02:00
parent 844092205e
commit ff95e2a9d8

View File

@ -123,9 +123,9 @@ pub fn eigen_qr<N: Float,
unsafe { shifter.unsafe_set((i, i), shift.clone()) }
}
let (q, r) = qr(&(eigenvalues - shifter));
let (q, r) = qr(&eigenvalues);// - shifter));
eigenvalues = r * q + shifter;
eigenvalues = r * q /*+ shifter*/;
eigenvectors = eigenvectors * q;
}