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
 
// 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::SparseMatrix<DenseVector::Scalar, Eigen::ColMajor, DenseVector::Index> SparseMatrix;
    typedef Eigen::Triplet<double, SparseMatrix::StorageIndex> Triplet;
    
    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
 
#endif // !defined EigenSupport_h
 
#undef EigenSupport_RECURSES
#endif // else defined(EigenSupport_RECURSES)