cgUzawa.hh

Go to the documentation of this file.

 
#ifndef cgUzawa_hh
#define cgUzawa_hh
 
#include <string>
 
#include "basics/typedefs.hh"
#include "operator/compositions.hh"
#include "vectorial/vectorial.hh"
 
namespace vectorial {
 
  using concepts::Real;
 
  // *************************************************************** CGUzawa **
 
  class CGUzawa : public concepts::Operator<Real> {
  public:
    CGUzawa(concepts::Operator<Real>& A, concepts::Operator<Real>& B,
            concepts::Operator<Real>& Bt, concepts::Operator<Real>& C,
            concepts::Operator<Real>& Ai, 
            Real maxeps, int maxit = 0, uint relres = false);
 
    virtual void operator()(const concepts::Function<Real>& fncY,
                            concepts::Function<Real>& fncX);
    void operator()(const concepts::Vector<Real>& fncY,
                    concepts::Vector<Real>& fncX);
 
    uint iterations() const { return it_; }
 
    Real epsilon() const { return eps_; }
  protected:
    std::ostream& info(std::ostream& os) const;
  private:
    uint bdim_[2];
    concepts::Operator<Real>& A_;
    concepts::Operator<Real>& B_;
    concepts::Operator<Real>& Bt_;
    concepts::Operator<Real>& C_;
    concepts::Operator<Real>& Ai_;
    Real maxeps_;
    uint maxit_;
    Real eps_;
    uint it_;
    bool relres_;
  };
 
} // namespace vectorial
 
#endif // cgUzawa_hh