eigensolver/testsuite/constrained.hh

Go to the documentation of this file.

// testing constraint EVP

#ifndef testConstrained_hh
#define testConstrained_hh

#include "basics/testcase.hh"
#include "operator/sparseMatrix.hh"
#include "toolbox/array.hh"

using concepts::Real;

namespace test {

  // ******************************************************** ConstrainedEVP **

  class ConstrainedEVP : public TestCase {
  public:
    ConstrainedEVP(uint dim, uint cdim, uint gdim, uint kmax) :
      spc(dim), cspc(cdim), gspc(gdim),
      a(0), b(spc, spc), c(spc, cspc), kmax_(kmax) {}
    virtual void run();
#ifdef HAS_ARPACK

    void testARPACKregularInv();
    void testARPACKshiftInv();
#endif
#ifdef HAS_JDBSym

    void testJDBSym();
#endif

    void testInexactInv();
  private:
    void testResults(const concepts::Array<Real>& res);
  protected:
    concepts::DummySpace<Real> spc, cspc, gspc;
    std::auto_ptr<concepts::Operator<Real> > a;
    concepts::SparseMatrix<Real> b, c;
    const Real* exact_;
    const uint kmax_;
  };

  // ********************************************************** GolubExample **

  class GolubExample : public ConstrainedEVP {
  public:
    GolubExample();
    virtual ~GolubExample() {}
  private:
    static const Real exact[4];
  };

  // ******************************************************* GolubExampleSum **

  class GolubExampleSum : public ConstrainedEVP {
  public:
    GolubExampleSum();
    virtual ~GolubExampleSum() {}
  private:
    static const Real exact[4];
    concepts::SparseMatrix<Real> a1, a2;
  };

  // ************************************************ MaxwellTransmissionEVP **

  class MaxwellTransmissionEVP : public ConstrainedEVP {
  public:
    MaxwellTransmissionEVP();
    virtual ~MaxwellTransmissionEVP() {}
  private:
    static const Real exact[20];
    void matlab_matrix_entry(concepts::Matrix<Real>& m,
           uint i, uint j, Real entry);
  };

} // namespace test

#endif // testConstrained_hh