SmallESTest.hh

 
#ifndef SmallESTest_hh
#define SmallESTest_hh
 
#include <iostream>
#include <cstdlib>
#include <unistd.h> // for command line parsing
#include "basics/testsuite.hh"
#include "basics/typedefs.hh"
 
#include "eigensolver/SmallES.hh"
#include "operator/superLU.hh"
#include "basics/exceptions.hh"
#include "toolbox/array.hh"
 
 
namespace test {
 
  class SmallESTest : public TestCase {
  public:
    virtual ~SmallESTest();
    virtual void run();
    void SmallESTestReal();
    void SmallESTestCmplx();
  private:
    /*
      IMPORTANT: The num Test doen't work for Real valued input with complex eigenvalues!!!!
                 The structure of the output will be changed soon (Issue 178)
    */
    template <typename F>
    void numTestArray_(const concepts::Array<F>& Id,
                       const concepts::Array<F>& ev1,
                       const concepts::Array<concepts::Vector<F>* >& ef1);
    template <typename F>
    void numTestArray_(const concepts::Matrix<F>& Id,
                       const concepts::Array<F>& ev1,
                       const concepts::Array<concepts::Vector<F>* >& ef1);
    template <typename F>
    void numTestArray_(const concepts::ElementMatrix<F>& Id,
                       const concepts::Array<F>& ev1,
                       const concepts::Array<concepts::Vector<F>* >& ef1);
    template <typename F>
    void numTestArray_(const concepts::DenseMatrix<F>& Id,
                       const concepts::Array<F>& ev1,
                       const concepts::Array<concepts::Vector<F>* >& ef1){
      numTestArray_(reinterpret_cast<const concepts::Matrix<F> &>(Id), ev1, ef1);
    };
 
 
 
  };
 
} // namespace test
 
# endif // SmallESTest_hh