mortar_schur.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MORTAR_SCHUR_HPP_
#define MORTAR_SCHUR_HPP_

#include <dune/istl/matrixmatrix.hh>
#include <opm/elasticity/applier.hpp>
#include <opm/elasticity/matrixops.hpp>

namespace Opm {
namespace Elasticity {

  template<class T>
class MortarBlockEvaluator : public Dune::LinearOperator<Vector, Vector> {
  public:
    // define the category
    enum {
      category=Dune::SolverCategory::sequential
    };

    MortarBlockEvaluator(T& Ai_,
                         const Matrix& B_) :
      Ai(Ai_), B(B_), op(Ai)
    {
    }

    void apply(const Vector& x, Vector& y) const
    {
      y = 0;
      applyscaleadd(1.0, x, y);
    }

    void applyscaleadd(field_type alpha, const Vector& x, Vector& y) const
    {
      Vector temp1(B.N());
      B.mv(x, temp1);
      Vector temp(temp1.size());
      op.pre(temp, temp1);
      Dune::InverseOperatorResult res;
      temp = 0;
      op.pre(temp, temp1);
      op.apply(temp, temp1);
      op.post(temp);
      B.usmtv(alpha, temp, y);
    }
  protected:
    T& Ai;            
    const Matrix& B;  
    mutable OperatorApplier<T> op; 
};

}
}

#endif