SmallES.hh

Go to the documentation of this file.

 
#ifndef SmallES_hh
#define SmallES_hh
 
#include <cfloat>
#include "eigensolver/eigens.hh"
#include "operator/matrix.hh"
 
namespace eigensolver {
 
using concepts::Real;
using concepts::Cmplx;
 
  // ********************************************************* SmallES **
 
  template <typename F = Cmplx>
  class SmallES : public EigenSolver<F> {
  public:
 
    SmallES(concepts::Array<F>& A);
 
    template <class G>
    SmallES(concepts::Matrix<G>& A);
 
    template <class G>
    SmallES(concepts::ElementMatrix<G>& A);
 
    template <class G>
    SmallES(concepts::DenseMatrix<G>& A) : SmallES(reinterpret_cast<concepts::Matrix<G> &>(A)) {};
 
 
    virtual ~SmallES();
    virtual const concepts::Array<F>& getEV();
    virtual const concepts::Array<concepts::Vector<F>*>& getEF();
    virtual uint iterations() const { return 0; }
    virtual uint converged() const { return 0; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    bool flag_;
    concepts::Array<F>* A_;
    concepts::Array<F> eigenvalues_;
    concepts::Array<concepts::Vector<F>*> eigenvectors_;
    bool computed_;
    void compute_();
  };
 
} // namespace eigensolver
 
#endif // SmallES_hh