testEigenSolver.cpp

Go to the documentation of this file.

#include <iostream>
#include <vector>
#include <list>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/math/linalg/CLinearAlgebraSolver.h"
#include "DGtal/math/linalg/CSparseMatrix.h"
#include "DGtal/math/linalg/CDynamicMatrix.h"
#include "DGtal/math/linalg/CDenseVector.h"
#include "DGtal/math/linalg/CDynamicVector.h"
#include "DGtal/math/linalg/EigenSupport.h"

using namespace std;
using namespace DGtal;

// Functions for testing Eigen solvers classes concepts

bool testEigenSolverConcepts()
{
    typedef EigenLinearAlgebraBackend LAB;
    typedef LAB::DenseVector Vector;
    typedef LAB::SparseMatrix Matrix;
    typedef LAB::Triplet Triplet;
    typedef std::vector<Triplet> TripletsVector;
    typedef std::list<Triplet> TripletsList;
    BOOST_CONCEPT_ASSERT(( concepts::CDynamicVector<Vector> ));
    BOOST_CONCEPT_ASSERT(( concepts::CDenseVector<Vector> ));
    BOOST_CONCEPT_ASSERT(( concepts::CDynamicMatrix<Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CSparseMatrix<Matrix, TripletsVector::const_iterator> ));
    BOOST_CONCEPT_ASSERT(( concepts::CSparseMatrix<Matrix, TripletsList::const_iterator> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverSimplicialLLT, Vector, Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverSimplicialLDLT, Vector, Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverConjugateGradient, Vector, Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverBiCGSTAB, Vector, Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverSparseLU, Vector, Matrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebraSolver<LAB::SolverSparseQR, Vector, Matrix> ));

    return true;
}

// Standard services - public :

int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class EigenSolver" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;

  bool res = testEigenSolverConcepts();
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //