vectorial/bramblePasciakFabric.hh

Go to the documentation of this file.

// fabric class for BramblePasciak

#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 << "BramblePasciakFabric(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