PlayaAmesosSolver.hpp

/* @HEADER@ */
//   
/* @HEADER@ */

#ifndef PLAYA_AMESOSSOLVER_HPP
#define PLAYA_AMESOSSOLVER_HPP

#include "PlayaDefs.hpp"
#include "PlayaLinearSolverBaseDecl.hpp"
#include "PlayaHandleable.hpp"
#include "PlayaPrintable.hpp"
#include "Teuchos_Describable.hpp"

#include "Teuchos_Array.hpp"
#include "Teuchos_RefCountPtr.hpp"
#include "Teuchos_ParameterList.hpp"

namespace Playa
{
  using namespace Teuchos;

  class AmesosSolver : public LinearSolverBase<double>,
                       public Playa::Handleable<LinearSolverBase<double> >,
                       public Printable,
                       public Describable
  {
  public:
    AmesosSolver(const ParameterList& params);

    virtual ~AmesosSolver(){;}

    void print(std::ostream& os) const 
    {
      os << description() << std::endl;
    }
    
    std::string description() const {return "AmesosSolver[" + kernel_ + "]";}

    

    virtual SolverState<double> solve(const LinearOperator<double>& op,
                                      const Vector<double>& rhs,
                                      Vector<double>& soln) const ;

    virtual RCP<LinearSolverBase<double> > getRcp() 
    {return rcp(this);}


  protected:

  private:
    std::string kernel_;
  };
  
}

#endif