src/armadillo-test.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

/*
 * This file is part of MXE.
 * See index.html for further 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;
}