GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-13 20:26:03 +08:00
Code Issues Packages Projects Releases Wiki Activity

oops bis, I forgot that SelfAdjointEigneSolver directly called the selector...

Browse Source
This commit is contained in:
Gael Guennebaud 2009-08-14 13:49:29 +02:00
parent 8abec72259
commit 6373c3cd00
1 changed files with 2 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -377,11 +377,8 @@ static void ei_tridiagonal_qr_step(RealScalar* diag, RealScalar* subdiag, int st
// G only modifies the two columns k and k+1 // G only modifies the two columns k and k+1
if (matrixQ) if (matrixQ)
{ {
#ifdef EIGEN_DEFAULT_TO_ROW_MAJOR Map<Matrix<Scalar,Dynamic,Dynamic> > q(matrixQ,n,n);
ei_apply_rotation_in_the_plane_selector<Scalar,Dynamic>::run(matrixQ+k, matrixQ+k+1, n, c, s, n, n); q.applyJacobiOnTheRight(k,k+1,c,s);
#else
ei_apply_rotation_in_the_plane_selector<Scalar,1>::run(matrixQ+k*n, matrixQ+k*n+n, n, c, s, 1, 1);
#endif
} }
} }
} }