ARPACK.hh

Go to the documentation of this file.

 
#ifndef ARPACK_hh
#define ARPACK_hh
 
#include <cfloat>
#include "eigensolver/eigens.hh"
 
namespace concepts {
  // forward declaration
  template<typename F>
  class Operator;
} // namespace concepts
 
namespace eigensolver {
 
using concepts::Real;
using concepts::Cmplx;
 
  // ********************************************************* ArPack **
 
  template <typename F, typename G = Real, typename H = Real>
  class ArPack : public EigenSolver<F> {
  public:
    enum which { LM, SM, LR, SR, LI, SI };
    enum modus { NORMAL = 1, REGINV = 2, SHIFTINV = 3, SHIFTINVIMAG = 4 };
    ArPack(concepts::Operator<G>& OP,
           concepts::Operator<F>& A,
           concepts::Operator<H>& B, const int kmax = 1,
           const Real tol = 0.0, const int maxiter = 300,
           enum which target = SM, enum modus mode = REGINV,
           const F sigma = 0.0,
           const concepts::Vector<F>* start = 0,
           const concepts::Array<F>* resid = 0,
           const bool schur = false);
    virtual ~ArPack();
    virtual const concepts::Array<F>& getEV();
    virtual const concepts::Array<concepts::Vector<F>*>& getEF();
    virtual uint iterations() const { return (uint)iter_; }
    virtual uint converged() const { return (uint)k_conv_; }
    concepts::Array<F> getRESID() { return resid_out_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Operator<G>& OP_;
    concepts::Operator<F>& A_;
    concepts::Operator<H>& B_;
    int kmax_;
    int maxiter_;
    Real tol_;
    int numop_, numopb_, numreo_, k_conv_, iter_;
    concepts::Array<F> eigenvalues_;
    concepts::Array<F> eigenvectors_;
    concepts::Array<concepts::Vector<F>*> ev_;
    enum which target_;
    enum modus mode_;
    const F sigma_;
    const concepts::Vector<F>* start_;
    const concepts::Array<F>* resid_in_;
    const bool schur_;
    bool computed_;
    concepts::Array<F> resid_out_;
    void compute_();
 
    void computeCmplx_();
  };
 
} // namespace eigensolver
 
#endif // ARPACK_hh