LinearSolver.hpp

#pragma once

#include <memory>
#include <string>

#include <amdis/CreatorMap.hpp>
#include <amdis/CreatorInterface.hpp>
#include <amdis/linearalgebra/LinearSolverInterface.hpp>

namespace AMDiS
{
  template <class Mat, class VecX, class VecY = VecX>
  class LinearSolver
      : public LinearSolverInterface<Mat,VecX,VecY>
  {
    using M = typename Mat::BaseMatrix;
    using X = typename VecX::BaseVector;
    using Y = typename VecX::BaseVector;

    using Map = CreatorMap<LinearSolverInterface<M,X,Y>>;

  public:
    LinearSolver(std::string const& name, std::string const& prefix)
      : solver_(named(Map::get(name, prefix))->createWithString(prefix))
    {}

    void init(Mat const& A) override
    {
      solver_->init(A.matrix());
    }

    void finish() override
    {
      solver_->finish();
    }

    void apply(VecX& x, VecY const& b, Dune::InverseOperatorResult& stat) override
    {
      solver_->apply(x.vector(), b.vector(), stat);
    }

  private:
    std::shared_ptr<LinearSolverInterface<M,X,Y>> solver_;
  };

} // end namespace AMDiS