bramblePasciakFabric.hh

Go to the documentation of this file.

 
#ifndef BramblePasciakfabric_hh
#define BramblePasciakfabric_hh
 
#if __GNUC__ == 2
#  include <float.h>
#  define EPS DBL_EPSILON
#else
#  include <limits>
#  define EPS std::numeric_limits<double>::epsilon()
#endif
#include "bramblePasciak.hh"
#include "operator/solverFabric.hh"
 
namespace vectorial {
 
  using concepts::Real;
 
  // ************************************************** BramblePasciakFabric **
 
  class BramblePasciakFabric : public concepts::SolverFabric<Real> {
  public:
    BramblePasciakFabric(concepts::Operator<Real>& A,
                         concepts::Operator<Real>& B,
                         concepts::Operator<Real>& Bt,
                         concepts::Operator<Real>& C,
                         concepts::Operator<Real>& W,
                         Real maxeps = EPS,
                         int maxit = 0, uint relres = 0)
      : A_(A), B_(B), Bt_(Bt), C_(C), W_(W),
        maxeps_(maxeps), maxit_(maxit), relres_(relres) { }
    virtual concepts::Operator<Real>* operator()
      (concepts::Operator<Real>& matrix) {
      return new BramblePasciak(A_, B_, Bt_, C_, W_, 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>& W_;
    Real maxeps_;
    int maxit_;
    uint relres_;
  };
 
} // namespace vectorial
 
#endif // BramblePasciakfabric_hh