tutoSymEigen.cpp File Reference

#include <STKpp.h>

Include dependency graph for tutoSymEigen.cpp:

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 4 of file tutoSymEigen.cpp.

{
  CSquareX a(5);
  a << 0, 1, 2, 3, 4,
       2, 3, 4, 5, 6,
       2, 1, 6, 7, 8,
       0, 3,-1, 2, 3,
       3,-1, 1, 1, 0;
  stk_cout << "STK++ Eigen decomposition:\n";
  stk_cout << "--------------------------\n";
  SymEigen<CSquareX> s(a.upperSymmetrize()); s.run();
  stk_cout << "D matrix:\n";
  stk_cout << s.eigenValues();
  CSquareX QDQt = s.eigenVectors() * s.eigenValues().diagonalize() * s.eigenVectors().transpose();
  stk_cout << "QDQ^T matrix:\n";
  stk_cout << QDQt;
  stk_cout << "\nlapack Eigen decomposition:\n";
  stk_cout << "--------------------------\n";
  lapack::SymEigen<CSquareX> ls(a); ls.run();
  stk_cout << "D matrix:\n";
  stk_cout << ls.eigenValues();
  QDQt = ls.eigenVectors() * ls.eigenValues().diagonalize() * ls.eigenVectors().transpose();
  stk_cout << "QDQ^T matrix:\n";
  stk_cout << QDQt;
  return 0;
}

References STK::ArrayBase< Derived >::diagonalize(), STK::ISymEigen< Derived >::eigenValues(), STK::ISymEigen< Derived >::eigenVectors(), STK::ISymEigen< Derived >::run(), stk_cout, STK::ArrayBase< Derived >::transpose(), and STK::ArrayBase< Derived >::upperSymmetrize().