EigenSupport.h

#pragma once

#if defined(EigenSupport_RECURSES)
#error Recursive header files inclusion detected in EigenSupport.h
#else // defined(EigenSupport_RECURSES)

#define EigenSupport_RECURSES

#if !defined EigenSupport_h

#define EigenSupport_h

#if defined(WITH_EIGEN)

// Eigen API extensions
#define EIGEN_DENSEBASE_PLUGIN "DGtal/math/linalg/EigenDenseBaseAddons.h"
#define EIGEN_SPARSEMATRIX_PLUGIN "DGtal/math/linalg/EigenSparseMatrixAddons.h"

// Inclusions
#include <iostream>
#if defined(__GNUG__) && !defined(__clang__)
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wpragmas"
#pragma GCC diagnostic ignored "-Wdeprecated-declarations"
#pragma GCC diagnostic ignored "-Wignored-attributes"
#pragma GCC diagnostic ignored "-Wshadow"
#endif
#if defined(__clang__)
#pragma clang diagnostic push
#pragma clang diagnostic ignored "-Wdocumentation"
#pragma clang diagnostic ignored "-Wignored-attributes"
#endif
#include <Eigen/Dense>
#include <Eigen/Sparse>
#include <Eigen/SparseCholesky>
#include <Eigen/IterativeLinearSolvers>
#include <Eigen/SparseLU>
#include <Eigen/SparseQR>
#if defined(__clang__)
#pragma clang diagnostic pop
#endif
#if defined(__GNUG__) && !defined(__clang__)
#pragma GCC diagnostic pop
#endif

namespace DGtal
{

// struct EigenLinearAlgebraBackend
struct EigenLinearAlgebraBackend
{
    typedef Eigen::VectorXd DenseVector;
    typedef Eigen::MatrixXd DenseMatrix;
    typedef Eigen::VectorXi IntegerVector;

    typedef Eigen::Triplet<double> Triplet;
    typedef Eigen::SparseMatrix<DenseVector::Scalar, Eigen::ColMajor, DenseVector::Index> SparseMatrix;

    typedef Eigen::SimplicialLLT<SparseMatrix> SolverSimplicialLLT;
    typedef Eigen::SimplicialLDLT<SparseMatrix> SolverSimplicialLDLT;
    typedef Eigen::ConjugateGradient<SparseMatrix> SolverConjugateGradient;
    typedef Eigen::BiCGSTAB<SparseMatrix> SolverBiCGSTAB;
    typedef Eigen::SparseLU<SparseMatrix> SolverSparseLU;
    typedef Eigen::SparseQR<SparseMatrix, Eigen::COLAMDOrdering<SparseMatrix::Index> > SolverSparseQR;
};


inline
std::ostream&
operator<<(std::ostream & os, const Eigen::ComputationInfo& info)
{
    switch (info)
    {
        case Eigen::Success:
            os << "success";
            break;
        case Eigen::NumericalIssue:
            os << "numerical_issue";
            break;
        case Eigen::NoConvergence:
            os << "no_convergence";
            break;
        case Eigen::InvalidInput:
            os << "invalid_input";
            break;
    }

    return os;
}

} // namespace DGtal

#else // defined WITH_EIGEN
#error "DGtal Eigen support hasn't been enabled"
#endif // defined WITH_EIGEN
#endif // !defined EigenSupport_h

#undef EigenSupport_RECURSES
#endif // else defined(EigenSupport_RECURSES)