eigen3.h

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * Copyright (C) 2012 Sergey Lisitsyn
 */

#ifndef EIGEN3_H_
#define EIGEN3_H_

#include <shogun/lib/config.h>

    //#define EIGEN_RUNTIME_NO_MALLOC
    #include <Eigen/Eigen>
    #include <Eigen/Dense>
    #if EIGEN_VERSION_AT_LEAST(3,0,93)
        #include <Eigen/Sparse>
    #else
        #define EIGEN_YES_I_KNOW_SPARSE_MODULE_IS_NOT_STABLE_YET
        #include <unsupported/Eigen/SparseExtra>

        #ifndef DOXYGEN_SHOULD_SKIP_THIS
        // Triplet definition for Eigen3 backword compatibility
        namespace Eigen {
        template <typename T> struct Triplet
        {
            Triplet(index_t colIndex, index_t rowIndex, T valueT) :
            ecol(colIndex), erow(rowIndex), evalue(valueT)
            {
            }
            index_t col() const { return ecol; };
            index_t row() const { return erow; };
            T value() const { return evalue; };
            index_t ecol;
            index_t erow;
            T evalue;
        };

        // SimplicialLLT definition for Eigen3 backword compatibility
        template <typename T> class SimplicialLLT
        : public SimplicialCholesky<T,Lower>
        {
        public:
            SimplicialLLT()
            {
                SimplicialCholesky<T>::setMode(SimplicialCholeskyLLt);
            }
            inline const T matrixL()
            {
                return SimplicialCholesky<T>::m_matrix;
            }
            inline const T matrixU()
            {
                return SimplicialCholesky<T>::m_matrix.transpose();
            }
        };
        }
        #endif //DOXYGEN_SHOULD_SKIP_THIS

    #endif  //EIGEN_VERSION_AT_LEAST(3,0,93)

#if ((EIGEN_WORLD_VERSION == 3) && (EIGEN_MAJOR_VERSION == 2) && \
    ((EIGEN_MINOR_VERSION == 91) || (EIGEN_MINOR_VERSION == 92)))
    // Regression has been introduced to eigen develop (3.3alpha1+):
    // http://eigen.tuxfamily.org/bz/show_bug.cgi?id=1229
    // until this is not fixed we need to copy the matrix and calculate the log
    #define EIGEN_WITH_LOG_BUG_1229 1
#endif
namespace shogun
{
template<class T> class SGSparseMatrix;

template<typename T> class EigenSparseUtil
{
    public:
    static Eigen::SparseMatrix<T> toEigenSparse(SGSparseMatrix<T> sg_matrix);
};

}


#endif