Lines Matching refs:MatrixType
17 template<typename MatrixType> void selfadjointeigensolver(const MatrixType& m)
25 typedef typename MatrixType::Scalar Scalar;
27 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
28 typedef Matrix<RealScalar, MatrixType::RowsAtCompileTime, 1> RealVectorType;
29 typedef typename std::complex<typename NumTraits<typename MatrixType::Scalar>::Real> Complex;
33 MatrixType a = MatrixType::Random(rows,cols);
34 MatrixType a1 = MatrixType::Random(rows,cols);
35 MatrixType symmA = a.adjoint() * a + a1.adjoint() * a1;
37 MatrixType b = MatrixType::Random(rows,cols);
38 MatrixType b1 = MatrixType::Random(rows,cols);
39 MatrixType symmB = b.adjoint() * b + b1.adjoint() * b1;
41 SelfAdjointEigenSolver<MatrixType> eiSymm(symmA);
43 SelfAdjointEigenSolver<MatrixType> eiSymmGen(symmA, symmB);
52 MatrixType _evec;
53 convert<MatrixType>(symmA, gSymmA);
54 convert<MatrixType>(symmB, gSymmB);
55 convert<MatrixType>(symmA, gEvec);
77 MatrixType normalized_eivec = eiSymmGen.eigenvectors()*eiSymmGen.eigenvectors().colwise().norm().asDiagonal().inverse();
95 MatrixType sqrtSymmA = eiSymm.operatorSqrt();
100 template<typename MatrixType> void eigensolver(const MatrixType& m)
108 typedef typename MatrixType::Scalar Scalar;
110 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
111 typedef Matrix<RealScalar, MatrixType::RowsAtCompileTime, 1> RealVectorType;
112 typedef typename std::complex<typename NumTraits<typename MatrixType::Scalar>::Real> Complex;
116 MatrixType a = MatrixType::Random(rows,cols);
117 MatrixType a1 = MatrixType::Random(rows,cols);
118 MatrixType symmA = a.adjoint() * a + a1.adjoint() * a1;
120 EigenSolver<MatrixType> ei0(symmA);
125 EigenSolver<MatrixType> ei1(a);