bicgstabFabric.hh

Go to the documentation of this file.

 
#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()() {
      return NULL;
    }
    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 << concepts::typeOf(*this) << "(maxeps = " << maxeps_ << ", maxit = "
                << maxit_ << ", relres = " << relres_ << ')';
    }
  private:
    Real maxeps_;
    int maxit_;
    bool relres_;
    bool throwing_;
  };
 
} // concepts
 
#endif // BiCGStabfabric_hh