jdbsym.hh

Go to the documentation of this file.

 
#ifndef jdbsym_hh
#define jdbsym_hh
 
#include "eigens.hh"
 
namespace concepts {
  // forward declaration
  template<typename F>
  class Operator;
  template<typename F>
  class SolverFabric;
} // namespace concepts
 
namespace eigensolver {
 
using concepts::Real;
 
  // **************************************************************** JdbSym **
 
  class JdbSym : public EigenSolver<Real> {
  public:
    JdbSym(concepts::Operator<Real>& stiff, concepts::Operator<Real>& mass,
     Real tol, uint maxit = 150, Real tau = 0.0, uint jdtype = 1,
     uint kmax = 1, concepts::SolverFabric<Real>* fabric = 0,
     const concepts::Array<concepts::Vector<Real>*>* start = 0);
    virtual ~JdbSym();
    virtual const concepts::Array<Real>& getEV();
    virtual const concepts::Array<concepts::Vector<Real>*>& getEF();
    virtual uint iterations() const { return (uint)it_; }
    virtual uint converged() const { return (uint)k_conv_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    void compute_();
    bool computed_;
    concepts::Operator<Real>& stiff_;
    concepts::Operator<Real>& mass_;
    Real tol_;
    uint maxit_;
    Real tau_;
    const uint jdtype_;
    uint kmax_;
    int it_;
    int k_conv_;
 
    concepts::Array<Real> lambda_;
    concepts::Array<Real> eigenvectors_;
    concepts::Array<concepts::Vector<Real>*> ev_;
 
    uint stiffCounter_, massCounter_;
 
    concepts::SolverFabric<Real>* fabric_;
    const concepts::Array<concepts::Vector<Real>*>* start_;
 
    static int jdbsymLock_;
  };
 
}
 
#endif // jdbsym_hh