aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/QR/SelfAdjointEigenSolver.h
diff options
context:
space:
mode:
Diffstat (limited to 'Eigen/src/QR/SelfAdjointEigenSolver.h')
-rw-r--r--Eigen/src/QR/SelfAdjointEigenSolver.h6
1 files changed, 3 insertions, 3 deletions
diff --git a/Eigen/src/QR/SelfAdjointEigenSolver.h b/Eigen/src/QR/SelfAdjointEigenSolver.h
index 2d833bf13..d366bafc7 100644
--- a/Eigen/src/QR/SelfAdjointEigenSolver.h
+++ b/Eigen/src/QR/SelfAdjointEigenSolver.h
@@ -230,17 +230,17 @@ compute(const MatrixType& matA, const MatrixType& matB, bool computeEigenvectors
Cholesky<MatrixType> cholB(matB);
// compute C = inv(U') A inv(U)
- MatrixType matC = cholB.matrixL().inverseProduct(matA);
+ MatrixType matC = cholB.matrixL().solveTriangular(matA);
// FIXME since we currently do not support A * inv(U),
// let's do (inv(U') A')' :
- matC = (cholB.matrixL().inverseProduct(matC.adjoint())).adjoint();
+ matC = (cholB.matrixL().solveTriangular(matC.adjoint())).adjoint();
compute(matC, computeEigenvectors);
if (computeEigenvectors)
{
// transform back the eigen vectors: evecs = inv(U) * evecs
- m_eivec = cholB.matrixL().adjoint().template marked<Upper>().inverseProduct(m_eivec);
+ m_eivec = cholB.matrixL().adjoint().template marked<Upper>().solveTriangular(m_eivec);
}
}