constrained.hh

Go to the documentation of this file.

 
#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::unique_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