operator/bicgstabFabric.hh

Go to the documentation of this file.

// fabric class for BiCGStab

#ifndef BiCGStabfabric_hh
#define BiCGStabfabric_hh

#if __GNUC__ == 2
#  include <float.h>
#  define EPS DBL_EPSILON
#else
#  include <limits>
#  define EPS std::numeric_limits<double>::epsilon()
#endif
#include "bicgstab.hh"
#include "solverFabric.hh"

namespace concepts {

  // ************************************************************** BiCGStabFabric **

  template<typename F, typename G = F>
  class BiCGStabFabric : public PrecondSolverFabric<F,G> {
  public:
    BiCGStabFabric(Real maxeps = EPS,
                   int maxit = 0, bool relres = false, bool throwing = true)
      : maxeps_(maxeps), maxit_(maxit), relres_(relres),
        throwing_(throwing) {}
    virtual Operator<F>* operator()(Operator<F>& matrix) {
      return new BiCGStab<F>(matrix, maxeps_, maxit_, relres_, throwing_);
    }
    virtual Operator<F>* operator()(Operator<F>& matrix, Operator<G>& precond)
    {
      return new BiCGStab<F>(matrix, precond, maxeps_, maxit_, relres_,
                             throwing_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "BiCGStabFabric(maxeps = " << maxeps_ << ", maxit = "
                << maxit_ << ", relres = " << relres_ << ')';
    }
  private:
    Real maxeps_;
    int maxit_;
    bool relres_;
    bool throwing_;
  };

} // concepts

#endif // BiCGStabfabric_hh