jdbsymFabric.hh

Go to the documentation of this file.

 
#ifndef jdbsymFabric_hh
#define jdbsymFabric_hh
 
#include <memory>
#include "jdbsym.hh"
#include "solverFabric.hh"
 
namespace eigensolver {
 
  // ********************************************************** JdbSymFabric **
 
  class JdbSymFabric : public SolverFabric<Real> {
  public:
    JdbSymFabric(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)
      : solver_(nullptr), tol_(tol), tau_(tau), maxit_(maxit), jdtype_(jdtype)
      , kmax_(kmax), fabric_(fabric), start_(start) {}
    virtual ~JdbSymFabric();
    virtual JdbSym& operator()(concepts::Operator<Real>& stiff,
                               concepts::Operator<Real>& mass);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<JdbSym> solver_;
    Real tol_, tau_;
    uint maxit_, jdtype_, kmax_;
    concepts::SolverFabric<Real>* fabric_;
    const concepts::Array<concepts::Vector<Real>*>* start_;
  };
 
} // eigensolver
 
#endif // jdbsymFabric_hh