mumpsFabric.hh

Go to the documentation of this file.

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