differentiableoperator.hh

#ifndef DUNE_FEM_DIFFERENTIABLEOPERATOR_HH
#define DUNE_FEM_DIFFERENTIABLEOPERATOR_HH
 
#include <dune/fem/operator/common/operator.hh>
 
namespace Dune
{
 
  namespace Fem
  {
 
    template< class JacobianOperator >
    class DifferentiableOperator
    : public virtual Dune::Fem::Operator< typename JacobianOperator::DomainFunctionType,
                                          typename JacobianOperator::RangeFunctionType >
    {
      typedef Dune::Fem::Operator< typename JacobianOperator::DomainFunctionType,
                                   typename JacobianOperator::RangeFunctionType > BaseType;
 
    public:
      typedef JacobianOperator JacobianOperatorType;
 
      typedef typename BaseType::DomainFunctionType DomainFunctionType;
      typedef typename BaseType::RangeFunctionType RangeFunctionType;
 
      typedef typename DomainFunctionType::RangeFieldType DomainFieldType;
      typedef typename RangeFunctionType::RangeFieldType RangeFieldType;
 
      virtual void jacobian ( const DomainFunctionType &u, JacobianOperatorType &jOp ) const = 0;
    };
 
  } // namespace Fem
 
} // namespace Dune
 
#endif // #ifndef DUNE_FEM_DIFFERENTIABLEOPERATOR_HH