cgFabric.hh

Go to the documentation of this file.

 
#ifndef CGfabric_hh
#define CGfabric_hh
 
#if __GNUC__ == 2
#  include <float.h>
#  define EPS DBL_EPSILON
#else
#  include <limits>
#  define EPS std::numeric_limits<double>::epsilon()
#endif
#include "cg.hh"
#include "solverFabric.hh"
 
namespace concepts {
 
  // ************************************************************** CGFabric **
 
  template<typename F>
  class CGFabric : public SolverFabric<F> {
  public:
    CGFabric(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 CG<F>(matrix, maxeps_, maxit_, relres_, throwing_);
    }
    Operator<F>* operator()() { return NULL; }
  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 // CGfabric_hh