compositions.hh

Go to the documentation of this file.

 
#ifndef testCompositions_hh
#define testCompositions_hh
 
#include <memory>
#include "basics/testcase.hh"
#include "operator/sparseMatrix.hh"
#include "operator/permutation.hh"
#include "sparseqr/driver.hh"
#include "sparseqr/givensRotations.hh"
 
namespace test {
 
  using concepts::Real;
 
  // *************************************************** CompositionTestBase **
 
  class CompositionTestBase : public TestCase {
  protected:
    template<typename F>
    void compare(const concepts::Matrix<F>& a, concepts::Operator<F>& b);
  };
 
  // ****************************************************** CompositionsTest **
 
  class CompositionsTest : public CompositionTestBase {
  public:
    CompositionsTest();
    virtual ~CompositionsTest() {}
    virtual void run();
 
 
    void test1();
    void test2();
    void test3();
    void test4();
    void test5();
  private:
    //concepts::DummySpace<Real> spc;
    uint dim;
    concepts::SparseMatrix<Real> a, b, c;
  };
 
  // ************************************************** MoreCompositionsTest **
 
  class MoreCompositionsTest : public CompositionTestBase {
  public:
    MoreCompositionsTest();
    virtual ~MoreCompositionsTest() {}
    virtual void run();
 
 
 
    void testPermutationConsistency();
    void testTrivExtendRestrictConsistency();
    void testCompositions();
  private:
//     concepts::DummySpace<Real> spc;
    uint dim;
    concepts::SparseMatrix<Real> a;
  };
 
  // ************************************************** DeepCompositionsTest **
 
  class DeepCompositionsTest : public CompositionTestBase {
  public:
    DeepCompositionsTest();
    virtual ~DeepCompositionsTest() {}
    virtual void run();
    void createSparseQR();
 
 
    void testCompositions();
  private:
    concepts::DummySpace<Real> spc, cspc;
    concepts::SparseMatrix<Real> a, b, c;
    std::unique_ptr<sparseqr::Driver> qr;
    concepts::Permutation<Real> * Pr;
    std::unique_ptr<concepts::Permutation<Real> > Prt;
    sparseqr::GivensRotations<Real> * Qgt, * Qg;
    concepts::TrivExtendRestrict<Real>* ext, * restr;
    std::unique_ptr<concepts::Compose<Real> > PrQgt, QgPrt, restrPrQgt,
      QgPrtext, restrPrQgtA, restrPrQgtAQgPrtext, restrPrQgtM,
      restrPrQgtMQgPrtext;
  };
 
}
 
#endif // testCompositions_hh