eigensolver/ARPACKsymm.hh

Go to the documentation of this file.

// ArPackSymm eigensolver using dsaupd

#ifndef ArPackSymmdsymm_hh
#define ArPackSymmdsymm_hh

#include <cfloat>
#include "eigens.hh"

namespace concepts {
  // forward declaration
  template<typename F>
  class Operator;
} // namespace concepts

namespace eigensolver {

using concepts::Real;

  // **************************************************************** ArPackSymm **

  class ArPackSymm : public EigenSolver<Real> {
  public:
    enum which { LA, SA, LM, SM, BE };
    enum modus { NORMAL = 1, REGINV = 2, SHIFTINV = 3, BUCKLING = 4,
     CAYLEY = 5 };
    ArPackSymm(concepts::Operator<Real>& OP,
               concepts::Operator<Real>& A,
               concepts::Operator<Real>& B,
               const int kmax = 1, const Real tol = 0.0,
               const int maxiter = 300,
               enum which target = SM, enum modus mode = REGINV,
               const Real sigma = 0.0,
               const concepts::Vector<Real>* start = 0);
    virtual ~ArPackSymm();
    virtual const concepts::Array<Real>& getEV();
    virtual const concepts::Array<concepts::Vector<Real>*>& getEF();
    virtual uint iterations() const { return (uint)iter_; }
    virtual uint converged() const { return (uint)k_conv_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Operator<Real>& OP_;
    concepts::Operator<Real>& A_;
    concepts::Operator<Real>& B_;
    int kmax_;
    int maxiter_;
    Real tol_;
    int numop_, numopb_, numreo_, k_conv_, iter_;
    concepts::Array<Real> eigenvalues_;
    concepts::Array<Real> eigenvectors_;
    concepts::Array<concepts::Vector<Real>*> ev_;
    enum which target_;
    enum modus mode_;
    const Real sigma_;
    const concepts::Vector<Real>* start_;
    bool computed_;
    void compute_();
  };

} // namespace eigensolver

#endif // ArPackSymm_hh