operator/testsuite/compositions.hh

Go to the documentation of this file.

// testing linear combination and composition

#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::auto_ptr<sparseqr::Driver> qr;
    concepts::Permutation<Real> * Pr;
    std::auto_ptr<concepts::Permutation<Real> > Prt;
    sparseqr::GivensRotations<Real> * Qgt, * Qg;
    concepts::TrivExtendRestrict<Real>* ext, * restr;
    std::auto_ptr<concepts::Compose<Real> > PrQgt, QgPrt, restrPrQgt,
      QgPrtext, restrPrQgtA, restrPrQgtAQgPrtext, restrPrQgtM,
      restrPrQgtMQgPrtext;
  };

}

#endif // testCompositions_hh