cgUzawaFabric.hh

Go to the documentation of this file.

 
#ifndef CGUzawafabric_hh
#define CGUzawafabric_hh
 
#if __GNUC__ == 2
#  include <float.h>
#  define EPS DBL_EPSILON
#else
#  include <limits>
#  define EPS std::numeric_limits<double>::epsilon()
#endif
#include "cgUzawa.hh"
#include "operator/solverFabric.hh"
 
using concepts::Real;
 
namespace vectorial {
 
  // ********************************************************* CGUzawaFabric **
 
  class CGUzawaFabric : public concepts::SolverFabric<Real> {
  public:
    CGUzawaFabric(concepts::Operator<Real>& A, concepts::Operator<Real>& B,
                  concepts::Operator<Real>& Bt, concepts::Operator<Real>& C,
                  concepts::Operator<Real>& Ai, 
                  Real maxeps = EPS,
                  int maxit = 0, uint relres = 0)
      : A_(A), B_(B), Bt_(Bt), C_(C), Ai_(Ai), 
        maxeps_(maxeps), maxit_(maxit), relres_(relres) { }
    virtual concepts::Operator<Real>* operator()
      (concepts::Operator<Real>& matrix) {
      return new CGUzawa(A_, B_, Bt_, C_, Ai_, maxeps_, maxit_, relres_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this)<<"(maxeps = " << maxeps_ << ", maxit = "
                << maxit_ << ", relres = " << relres_ << ")";
    }
  private:
    concepts::Operator<Real>& A_;
    concepts::Operator<Real>& B_;
    concepts::Operator<Real>& Bt_;
    concepts::Operator<Real>& C_;
    concepts::Operator<Real>& Ai_;
    Real maxeps_;
    int maxit_;
    uint relres_;
  };
 
} // namespace vectorial
 
#endif // CGUzawafabric_hh