11 #ifndef EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H
12 #define EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H
14 #include "./Tridiagonalization.h"
47 template<
typename _MatrixType>
53 typedef typename Base::Index Index;
54 typedef _MatrixType MatrixType;
162 template<
typename MatrixType>
164 compute(
const MatrixType& matA,
const MatrixType& matB,
int options)
166 eigen_assert(matA.cols()==matA.rows() && matB.rows()==matA.rows() && matB.cols()==matB.rows());
167 eigen_assert((options&~(EigVecMask|GenEigMask))==0
168 && (options&EigVecMask)!=EigVecMask
169 && ((options&GenEigMask)==0 || (options&GenEigMask)==
Ax_lBx
170 || (options&GenEigMask)==
ABx_lx || (options&GenEigMask)==
BAx_lx)
171 &&
"invalid option parameter");
173 bool computeEigVecs = ((options&EigVecMask)==0) || ((options&EigVecMask)==
ComputeEigenvectors);
178 int type = (options&GenEigMask);
185 MatrixType matC = matA.template selfadjointView<Lower>();
186 cholB.
matrixL().template solveInPlace<OnTheLeft>(matC);
187 cholB.
matrixU().template solveInPlace<OnTheRight>(matC);
193 cholB.
matrixU().solveInPlace(Base::m_eivec);
198 MatrixType matC = matA.template selfadjointView<Lower>();
206 cholB.
matrixU().solveInPlace(Base::m_eivec);
211 MatrixType matC = matA.template selfadjointView<Lower>();
219 Base::m_eivec = cholB.
matrixL() * Base::m_eivec;
227 #endif // EIGEN_GENERALIZEDSELFADJOINTEIGENSOLVER_H