33 #ifndef SHARK_LINALG_EIGENVALUES_H 34 #define SHARK_LINALG_EIGENVALUES_H 39 namespace shark{
namespace blas{
61 template<
class MatrixT,
class MatrixU,
class VectorT>
70 std::size_t n = A().size1();
72 eigenVectors().resize(n,n);
73 eigenVectors().clear();
74 eigenValues().resize(n);
75 eigenValues().clear();
78 eigenVectors()( 0 , 0 ) = 1;
79 eigenValues()( 0 ) = A()( 0 , 0 );
84 for (std::size_t i = 0; i < n; i++) {
85 for (std::size_t j = 0; j <= i; j++) {
86 eigenVectors()(i, j) = A()(i, j);