arpackFabric.hh

Go to the documentation of this file.

 
#ifndef arpackFabric_hh
#define arpackFabric_hh
 
#include <memory>
#include "operator/compositions.hh"
#include "ARPACKsymm.hh"
#include "solverFabric.hh"
 
namespace eigensolver {
 
  // ********************************************************** ArPackFabric **
 
  class ArPackFabric : public SolverFabric<Real> {
  public:
    ArPackFabric(int kmax = 1, Real tol = 0.0, const int maxiter = 300,
     enum ArPackSymm::which target = ArPackSymm::SM,
     enum ArPackSymm::modus mode = ArPackSymm::REGINV,
     const Real sigma = 0.0,
     const concepts::Vector<Real>* start = 0) :
      solver_(nullptr), OP_(nullptr), B_(nullptr), tmp1_(nullptr), tmp2_(nullptr), kmax_(kmax),
      maxiter_(maxiter), tol_(tol), target_(target), mode_(mode),
      sigma_(sigma), start_(start) {}
    virtual ~ArPackFabric();
    virtual ArPackSymm& operator()(concepts::Operator<Real>& stiff,
             concepts::Operator<Real>& mass);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<ArPackSymm> solver_;
    std::unique_ptr<concepts::Operator<Real> > OP_, B_, tmp1_, tmp2_;
    int kmax_;
    int maxiter_;
    Real tol_;
    enum ArPackSymm::which target_;
    enum ArPackSymm::modus mode_;
    const Real sigma_;
    const concepts::Vector<Real>* start_;
    ArPackFabric(const ArPackFabric&);
 
    void normalMode_(concepts::Operator<Real>& stiff,
        concepts::Operator<Real>& mass);
    void reginvMode_(concepts::Operator<Real>& stiff,
        concepts::Operator<Real>& mass);
    void shiftinvMode_(concepts::Operator<Real>& stiff,
          concepts::Operator<Real>& mass);
    void bucklingMode_(concepts::Operator<Real>& stiff,
          concepts::Operator<Real>& mass);
    void cayleyMode_(concepts::Operator<Real>& stiff,
        concepts::Operator<Real>& mass);
  };
 
} // namespace eigensolver
 
#endif // arpackFabric_hh