/*
* This file is part of MXE.
* See index.html for further information.
*/
// From http://arma.sourceforge.net/docs.html#eig_sym
#include
using namespace arma;
int main(int argc, char** argv)
{
mat A = randu(50,50);
mat B = trans(A)*A; // generate a symmetric matrix
vec eigval;
mat eigvec;
// use standard algorithm by default
eig_sym(eigval, eigvec, B);
// use divide & conquer algorithm
eig_sym(eigval, eigvec, B, "dc");
return 0;
}