From 3e348fdcf9f5360cbc4bc8bb558d11c82f4be752 Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Mon, 11 Jul 2016 15:24:52 +0200 Subject: Workaround MSVC bug --- Eigen/src/Eigenvalues/GeneralizedEigenSolver.h | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) (limited to 'Eigen/src/Eigenvalues') diff --git a/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h b/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h index 650617ca7..19730f9e4 100644 --- a/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h +++ b/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h @@ -310,6 +310,7 @@ GeneralizedEigenSolver::compute(const MatrixType& A, const MatrixTyp if (m_realQZ.info() == Success) { m_matS = m_realQZ.matrixS(); + const MatrixType &matT = m_realQZ.matrixT(); if (computeEigenvectors) m_eivec = m_realQZ.matrixZ().transpose(); @@ -322,7 +323,7 @@ GeneralizedEigenSolver::compute(const MatrixType& A, const MatrixTyp if (i == A.cols() - 1 || m_matS.coeff(i+1, i) == Scalar(0)) { m_alphas.coeffRef(i) = m_matS.coeff(i, i); - m_betas.coeffRef(i) = m_realQZ.matrixT().coeff(i,i); + m_betas.coeffRef(i) = matT.coeff(i,i); ++i; } else @@ -332,7 +333,8 @@ GeneralizedEigenSolver::compute(const MatrixType& A, const MatrixTyp // T = [a 0] // [0 b] - RealScalar a = m_realQZ.matrixT().coeff(i, i), b = m_realQZ.matrixT().coeff(i+1, i+1); + RealScalar a = matT.coeff(i, i), + b = matT(i+1, i+1); // NOTE: using operator() instead of coeff() workarounds a MSVC bug. Matrix S2 = m_matS.template block<2,2>(i,i) * Matrix(b,a).asDiagonal(); Scalar p = Scalar(0.5) * (S2.coeff(0,0) - S2.coeff(1,1)); -- cgit v1.2.3