inexactFabric.hh

Go to the documentation of this file.

 
#ifndef inexactFabric_hh
#define inexactFabric_hh
 
#include <memory>
#include "inexactInv.hh"
#include "solverFabric.hh"
 
namespace eigensolver {
 
  // ****************************************************** InexactInvFabric **
 
  class InexactInvFabric : public SolverFabric<Real> {
  public:
    InexactInvFabric(const concepts::Vector<Real>& EFguess,
         const Real innerTol = 0.5, const uint innerMax = 100,
         const Real outerTol = 1e-6, const uint outerMax = 100) :
      solver_(0), EFguess_(EFguess), innerTol_(innerTol), outerTol_(outerTol),
      innerMax_(innerMax), outerMax_(outerMax) {}
    virtual ~InexactInvFabric();
    virtual InexactInv<Real>& operator()(concepts::Operator<Real,Real>& stiff,
           concepts::Operator<Real,Real>& mass);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<InexactInv<Real> > solver_;
 
    concepts::Vector<Real> EFguess_;
    Real innerTol_, outerTol_;
    uint innerMax_, outerMax_;
 
    InexactInvFabric(const InexactInvFabric&);
  };
 
} // namespace eigensolver
 
#endif // inexactFabric_hh