DirPowIt.hh

Go to the documentation of this file.

 
#ifndef DirPowIt_hh
#define DirPowIt_hh
 
//#include <cfloat>
#include "eigensolver/eigens.hh"
#include "basics/vectorsMatrices.hh"
 
namespace concepts {
  // forward declaration
  template<typename F>
  class Operator;
} // namespace concepts
 
namespace eigensolver {
 
using concepts::Real;
using concepts::Cmplx;
 
// ********************************************************* DirPowIt**
 
  template <typename F, typename G>
  class DirPowIt : public EigenSolver<G> {
  public:
    DirPowIt(concepts::Operator<F>& A,
             const int nev = 3, const int m = 3, const int k = 0,
             const int maxiter = 300,
             concepts::Array<concepts::Vector<G>*>* start = 0,
             const Real tol = 0
      );
    virtual ~DirPowIt();
    virtual const concepts::Array<G>& getEV();
    virtual const concepts::Array<concepts::Vector<G>*>& getEF();
    virtual uint iterations() const { return (uint)iter_; }
    virtual uint converged() const {return 0;}
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Operator<F>& A_;
    int nev_;
    int m_;
    int k_;
    int iter_;
    int maxiter_;
    Real tol_;
    concepts::Array<G> eigenvalues_;
    concepts::Array<concepts::Vector<G>*> eigenvectors_;
    concepts::Array<concepts::Vector<G>*> temp_eigenvectors_;
    void sort_(concepts::Array<Real>&);
    concepts::Array<uint> permutation_;
    bool computed_;
    void compute_();
  };
 
} // namespace eigensolver
 
#endif // DirPowIt_hh