precon_wrapper.hpp

// FEAT3: Finite Element Analysis Toolbox, Version 3
// Copyright (C) 2010 by Stefan Turek & the FEAT group
// FEAT3 is released under the GNU General Public License version 3,
// see the file 'copyright.txt' in the top level directory for details.
 
#pragma once
 
// includes, FEAT
#include <kernel/solver/base.hpp>
#include <kernel/solver/legacy_preconditioners.hpp>
 
namespace FEAT
{
  namespace Solver
  {
    template<
      typename Matrix_,
      typename Filter_,
      template<typename,typename> class Precon_>
    class PreconWrapper :
      public SolverBase<typename Matrix_::VectorTypeR>
    {
    public:
      typedef Matrix_ MatrixType;
      typedef typename MatrixType::VectorTypeR VectorType;
 
    protected:
      const Filter_& _filter;
      Precon_<MatrixType, VectorType> _precond;
 
    public:
      template<typename... Args_>
      explicit PreconWrapper(const Filter_& filter, Args_&&... args) :
        _filter(filter),
        _precond(std::forward<Args_>(args)...)
      {
      }
 
      virtual String name() const override
      {
        return _precond.name();
      }
 
      virtual Status apply(VectorType& vec_cor, const VectorType& vec_def) override
      {
        _precond.apply(vec_cor, vec_def);
        this->_filter.filter_cor(vec_cor);
        return Status::success;
      }
    }; // class PreconWrapper<...>
  } // namespace Solver
} // namespace FEAT