mxe/src/armadillo-test.cpp

24 lines
424 B
C++

/*
* This file is part of MXE. See LICENSE.md for licensing information.
*/
#include <armadillo>
using namespace arma;
int main()
{
mat A = randu<mat>(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;
}