mumpsoverlapFabric.hh

Go to the documentation of this file.

 
#ifndef mumpsoverlapfabric_hh
#define mumpsoverlapfabric_hh
 
#include "mumpsoverlap.hh"
#include "solverFabric.hh"
#include "sparseMatrix.hh"
#include "CRS.hh"
 
namespace concepts {
 
  // *********************************************************** MumpsFabric **
 
  template<class F = Real>
  class MumpsOverlapFabric : public SolverFabric<F> {
  public:
    MumpsOverlapFabric(const int part) : part_(part) {};
    virtual Operator<F>* operator()() {return new Operator<F>(1,1);};
    virtual MumpsOverlap<F>* operator()(Operator<F>& matrix) {
      if (dynamic_cast<concepts::CRSConvertable<F>*>(&matrix))
        return new MumpsOverlap<F>(matrix, this->part_);
      else {
        concepts::SparseMatrix<F> tmpMatrix(matrix);
        return new MumpsOverlap<F>(tmpMatrix, this->part_);
      }
    }
  protected:
    int part_;
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this)<<"()";
    }
  };
 
} // concepts
 
#endif // mumpsoverlapfabric_hh