inexactInv.hh

Go to the documentation of this file.

 
#ifndef inexactInv_hh
#define inexactInv_hh
 
#include "eigens.hh"
 
namespace concepts {
  // forward declaration
  template<typename F>
  class Operator;
} // namespace concepts
 
namespace eigensolver {
 
  using concepts::Real;
 
  // ************************************************************ InexactInv **
 
  template<typename F>
  class InexactInv : public EigenSolver<F> {
  public:
    InexactInv(concepts::Operator<F>& A, concepts::Operator<F>& B,
         const concepts::Vector<F>& EFguess,
         const Real innerTol = 0.5, const uint innerMax = 100,
         const Real outerTol = 1e-6, const uint outerMax = 100);
 
    virtual ~InexactInv();
    virtual const concepts::Array<F>& getEV() {
      compute_();
      return returnValue1_; }
    virtual const concepts::Array<concepts::Vector<F>*>& getEF() {
      compute_();
      return returnValue2_;
    }
    virtual uint iterations() const { return iterations_; }
    virtual uint converged() const {
      return (epsilon_ < outerTol_) ? 1 : 0;
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Operator<F>& A_, & B_;
    concepts::Array<F> returnValue1_;
    F& eigenvalue_;
    concepts::Vector<F> eigenfunction_;
    concepts::Array<concepts::Vector<F>*> returnValue2_;
    Real innerTol_, outerTol_;
    uint innerMax_, outerMax_;
 
    uint iterations_, innerIter_;
    Real epsilon_;
 
    bool computed_;
 
    void compute_();
  };
 
} // namespace eigensolver
 
#endif // inexactInv_hh