operator/cgFabric.hh

Go to the documentation of this file.

// fabric class for CG

#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_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "CGFabric(maxeps = " << maxeps_ << ", maxit = "
    << maxit_ << ", relres = " << relres_ << ')';
    }
  private:
    Real maxeps_;
    int maxit_;
    bool relres_;
    bool throwing_;
  };

} // concepts

#endif // CGfabric_hh