From 6373c3cd0014ad8e0a6952915c3c6342b57c1863 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Fri, 14 Aug 2009 13:49:29 +0200
Subject: [PATCH] oops bis, I forgot that SelfAdjointEigneSolver directly
 called the selector...

---
 Eigen/src/QR/SelfAdjointEigenSolver.h | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/Eigen/src/QR/SelfAdjointEigenSolver.h b/Eigen/src/QR/SelfAdjointEigenSolver.h
index a8e89ba1a..d07d28424 100644
--- a/Eigen/src/QR/SelfAdjointEigenSolver.h
+++ b/Eigen/src/QR/SelfAdjointEigenSolver.h
@@ -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
     if (matrixQ)
     {
-      #ifdef EIGEN_DEFAULT_TO_ROW_MAJOR
-      ei_apply_rotation_in_the_plane_selector<Scalar,Dynamic>::run(matrixQ+k, matrixQ+k+1, n, c, s, n, n);
-      #else
-      ei_apply_rotation_in_the_plane_selector<Scalar,1>::run(matrixQ+k*n, matrixQ+k*n+n, n, c, s, 1, 1);
-      #endif
+      Map<Matrix<Scalar,Dynamic,Dynamic> > q(matrixQ,n,n);
+      q.applyJacobiOnTheRight(k,k+1,c,s);
     }
   }
 }