superLUfabric.hh

Go to the documentation of this file.

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