applier.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef APPLIER_H_
#define APPLIER_H_

#include <dune/istl/solver.hh>

#include <opm/elasticity/matrixops.hpp>

namespace Opm {
  namespace Elasticity {

  template<class T>
struct OperatorApplier
{
  explicit OperatorApplier(T& t) : A(t)
  {
  }

  void apply(Vector& v, Vector& d);

  void pre(Vector& x, Vector& b);

  void post(Vector& x);

  T& A;
};

typedef OperatorApplier<Dune::InverseOperator<Vector, Vector> > InverseApplier;
typedef OperatorApplier<Dune::Preconditioner<Vector, Vector> > PreApplier;

  template<>
void InverseApplier::apply(Vector& v, Vector& d)
{
  Dune::InverseOperatorResult r;
  A.apply(v, d, r);
}

  template<>
void PreApplier::apply(Vector& v, Vector& d)
{
  A.apply(v, d);
}

  template<>
void InverseApplier::pre(Vector& /* x */, Vector& /* b */)
{
}

  template<>
void PreApplier::pre(Vector& x, Vector& b)
{
  A.pre(x,b);
}

  template<>
void InverseApplier::post(Vector& /* x */)
{
}

  template<>
void PreApplier::post(Vector& x)
{
  A.post(x);
}

}
}

#endif