GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-22 17:49:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

fix a bug -- bad read found by valgrind

Browse Source
This commit is contained in:
Benoit Jacob 2009-01-08 14:23:30 +00:00
parent af27fb7590
commit e2d2a7d222
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/QR/SelfAdjointEigenSolver.h
View File

@ -357,10 +357,12 @@ static void ei_tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, int st
subdiag[k - 1] = c * subdiag[k-1] - s * z; subdiag[k - 1] = c * subdiag[k-1] - s * z;
x = subdiag[k]; x = subdiag[k];
z = -s * subdiag[k+1];
if (k < end - 1) if (k < end - 1)
{
z = -s * subdiag[k+1];
subdiag[k + 1] = c * subdiag[k+1]; subdiag[k + 1] = c * subdiag[k+1];
}
// apply the givens rotation to the unit matrix Q = Q * G // apply the givens rotation to the unit matrix Q = Q * G
// G only modifies the two columns k and k+1 // G only modifies the two columns k and k+1