operator/PETScFabric.hh

Go to the documentation of this file.

// fabric class for PETSc

#ifndef petscFabric_hh
#define petscFabric_hh

#include <list>
#include "operator/PETSc.hh"
#include "operator/solverFabric.hh"

namespace concepts {

  // *********************************************************** PETScFabric **

  class PETScFabric : public SolverFabric<Real> {
  public:
    PETScFabric(const Real maxeps, const std::string ksptype, 
                const std::string pctype = std::string("none"))
      : maxeps_(maxeps), ksptype_(ksptype), pctype_(pctype) {}
    virtual Operator<Real>* operator()(Operator<Real>& matrix) {
      if (pctype_ == "none" || dynamic_cast<concepts::PETScMat*>(&matrix)) {
        return new PETSc(matrix, maxeps_, ksptype_, pctype_);
      } else {
        concepts::SparseMatrix<Real>* sparsemat =
          dynamic_cast<concepts::SparseMatrix<Real>*>(&matrix);
        if (sparsemat) {
          petscmat_.push_front(new concepts::PETScMat(*sparsemat));
        } else {
          sparsemat = new concepts::SparseMatrix<Real>(matrix);
          petscmat_.push_front(new concepts::PETScMat(*sparsemat));
          delete sparsemat;
        }
        return new PETSc(*(petscmat_.front()), maxeps_, ksptype_, pctype_);
      }
    }
    ~PETScFabric() {
      while (!petscmat_.empty()) {
        delete petscmat_.front();
        petscmat_.pop_front();
      }
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "PETScFabric(maxeps = " << maxeps_ << ", ksptype = "
                << ksptype_ << ", pctype = " << pctype_ << ")";
    }
  private:
    const Real maxeps_;
    const std::string ksptype_;
    const std::string pctype_;
    std::list<concepts::PETScMat*> petscmat_;
  };

} // concepts

#endif // petscFabric_hh