eigensolver/inexactInv.hh

Go to the documentation of this file.

/* inexact inverse iteration
 */

#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