eigenoperator.hh

#ifndef DUNE_FEM_EIGENOPERATOR_HH
#define DUNE_FEM_EIGENOPERATOR_HH
 
#ifdef HAVE_EIGEN
 
// system includes
#include <string>
 
// local includes
#include <dune/fem/operator/matrix/eigenmatrix.hh>
 
namespace Dune
{
 
  namespace Fem
  {
 
    template< class DomainFunction, class RangeFunction >
    struct EigenLinearOperator
    : public EigenMatrixObject< typename DomainFunction::DiscreteFunctionSpaceType, typename RangeFunction::DiscreteFunctionSpaceType >,
      public Fem::AssembledOperator< DomainFunction, RangeFunction >
    {
      typedef typename DomainFunction::DiscreteFunctionSpaceType DomainSpaceType;
      typedef typename RangeFunction::DiscreteFunctionSpaceType RangeSpaceType;
      typedef EigenLinearOperator< DomainFunction, RangeFunction > ThisType;
      typedef EigenMatrixObject< DomainSpaceType, RangeSpaceType > BaseType;
 
      static constexpr bool assembled = true ;
 
      using BaseType::apply;
      using BaseType::exportMatrix;
 
      EigenLinearOperator( const std::string & ,
                           const DomainSpaceType &domainSpace,
                           const RangeSpaceType &rangeSpace,
                           const SolverParameter& param = SolverParameter() ) :
        BaseType( domainSpace, rangeSpace, param )
      {}
 
      virtual void clear() { BaseType::clear(); }
 
      virtual void operator()( const DomainFunction &arg, RangeFunction &dest ) const
      {
        apply( arg, dest );
      }
 
    };
  } // namespace Fem
 
} // namespace Dune
 
#endif
 
#endif // #ifndef DUNE_FEM_SPLINEAR_HH