matrix_eigen_unittest.cpp

Go to the documentation of this file.

/* +------------------------------------------------------------------------+
   |                     Mobile Robot Programming Toolkit (MRPT)            |
   |                          http://www.mrpt.org/                          |
   |                                                                        |
   | Copyright (c) 2005-2018, Individual contributors, see AUTHORS file     |
   | See: http://www.mrpt.org/Authors - All rights reserved.                |
   | Released under BSD License. See details in http://www.mrpt.org/License |
   +------------------------------------------------------------------------+ */

#include <Eigen/Dense>
#include <gtest/gtest.h>

using namespace Eigen;
using namespace std;

/** A macro for obtaining the name of the current function:  */
#if defined(_MSC_VER) && (_MSC_VER >= 1300)
#define __CURRENT_FUNCTION_NAME__ __FUNCTION__
#elif defined(_MSC_VER) && (_MSC_VER < 1300)
// Visual C++ 6 HAS NOT A __FUNCTION__ equivalent.
#define __CURRENT_FUNCTION_NAME__ ::system::extractFileName(__FILE__).c_str()
#else
#define __CURRENT_FUNCTION_NAME__ __PRETTY_FUNCTION__
#endif

#if 0

template <int ColRowOrder>
void do_test_EigenVal4x4_sym_vs_generic_eigen()
{
    using Mat44 = Matrix<double,4,4,ColRowOrder>;

    const double   dat_C1[] = {
        13.737245,10.248641,-5.839599,11.108320,
        10.248641,14.966139,-5.259922,11.662222,
        -5.839599,-5.259922,9.608822,-4.342505,
        11.108320,11.662222,-4.342505,12.121940 };
    const Mat44 C1(dat_C1);  // It doesn't mind the row/col major order since data are symetric

    // Symetric --------------------
    // This solver returns the eigenvectors already sorted.
    Eigen::SelfAdjointEigenSolver<Mat44> eigensolver(C1);
//  MatrixXd eVecs_s = eigensolver.eigenvectors();
//  MatrixXd eVals_s = eigensolver.eigenvalues();

    cout << endl << __CURRENT_FUNCTION_NAME__ << endl
        << "SelfAdjointEigenSolver:\n"
        << "eigvecs: " << endl << eigensolver.eigenvectors() << endl
        << "eigvals: " << endl << eigensolver.eigenvalues() << endl;

    // Generic ---------------------
    Eigen::EigenSolver<Mat44> es(C1, true);
//  MatrixXd eVecs_g = es.eigenvectors().real();
//  MatrixXd eVals_g = es.eigenvalues().real();

    cout << endl
        << "EigenSolver:\n"
        << "eigvecs: " << endl << es.eigenvectors() << endl
        << "eigvals: " << endl << es.eigenvalues() << endl;
}

// Compare the two ways of computing matrix eigenvectors: generic & for symmetric matrices:
TEST(MatricesEigen,EigenVal4x4_sym_vs_generic)
{
    do_test_EigenVal4x4_sym_vs_generic_eigen<Eigen::ColMajor>();
    do_test_EigenVal4x4_sym_vs_generic_eigen<Eigen::RowMajor>();
}

#endif