26 #ifndef EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H
27 #define EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H
62 template<
typename _MatrixType>
177 template<
typename MatrixType>
181 eigen_assert(matA.cols()==matA.rows() && matB.rows()==matA.rows() && matB.cols()==matB.rows());
185 || (options&GenEigMask)==
ABx_lx || (options&GenEigMask)==
BAx_lx)
186 &&
"invalid option parameter");
200 MatrixType matC = matA.template selfadjointView<Lower>();
201 cholB.
matrixL().template solveInPlace<OnTheLeft>(matC);
202 cholB.
matrixU().template solveInPlace<OnTheRight>(matC);
208 cholB.
matrixU().solveInPlace(Base::m_eivec);
213 MatrixType matC = matA.template selfadjointView<Lower>();
221 cholB.
matrixU().solveInPlace(Base::m_eivec);
226 MatrixType matC = matA.template selfadjointView<Lower>();
234 Base::m_eivec = cholB.
matrixL() * Base::m_eivec;
242 #endif // EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H