differentiablefunctionfromcallables.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_FUNCTIONS_COMMON_DIFFEREENTIONABEFUNCTIONFROMCALLABLES_HH
#define DUNE_FUNCTIONS_COMMON_DIFFEREENTIONABEFUNCTIONFROMCALLABLES_HH
 
 
#include <dune/common/typeutilities.hh>
#include <dune/common/hybridutilities.hh>
 
#include <dune/functions/common/signature.hh>
 
#include <dune/functions/common/differentiablefunction.hh>
#include <dune/functions/common/functionconcepts.hh>
 
 
 
namespace Dune {
namespace Functions {
 
 
 
template<class Signature, template<class> class DerivativeTraits, class... Callables>
class DifferentiableFunctionFromCallables;
 
 
 
template<class Range, class Domain, template<class> class DerivativeTraits, class F>
class DifferentiableFunctionFromCallables<Range(Domain), DerivativeTraits, F>
{
public:
 
  using Signature = Range(Domain);
 
  using RawSignature = typename SignatureTraits<Signature>::RawSignature;
 
  using DerivativeSignature = typename DerivativeTraits<RawSignature>::Range(Domain);
 
  using Derivative = DifferentiableFunction<DerivativeSignature, DerivativeTraits>;
 
  template<class FF, disableCopyMove<DifferentiableFunctionFromCallables, FF> = 0>
  DifferentiableFunctionFromCallables(FF&& f) :
    f_(std::forward<FF>(f))
  {}
 
  Range operator() (const Domain& x) const
  {
    return f_(x);
  }
 
  friend Derivative derivative(const DifferentiableFunctionFromCallables& t)
  {
    DUNE_THROW(Dune::NotImplemented, "Derivative not implemented");
  }
 
private:
  F f_;
};
 
 
 
template<class Range, class Domain, template<class> class DerivativeTraits, class F, class DF, class... Derivatives>
class DifferentiableFunctionFromCallables<Range(Domain), DerivativeTraits, F, DF, Derivatives...>
{
public:
 
  using Signature = Range(Domain);
  using RawSignature = typename SignatureTraits<Signature>::RawSignature;
  using DerivativeSignature = typename DerivativeTraits<RawSignature>::Range(Domain);
 
  using Derivative = DifferentiableFunctionFromCallables<DerivativeSignature, DerivativeTraits, DF, Derivatives...>;
 
  template<class FF, class DFF, class... DDFF>
  DifferentiableFunctionFromCallables(FF&& f, DFF&& df, DDFF&&... ddf) :
    f_(std::forward<FF>(f)),
    df_(std::forward<DFF>(df), std::forward<DDFF>(ddf)...)
  {}
 
  Range operator() (const Domain& x) const
  {
    return f_(x);
  }
 
  friend Derivative derivative(const DifferentiableFunctionFromCallables& t)
  {
    return t.df_;
  }
 
private:
  F f_;
  Derivative df_;
};
 
 
template<class Signature, template<class> class DerivativeTraits, class... F>
DifferentiableFunctionFromCallables<Signature, DerivativeTraits, F...>
  makeDifferentiableFunctionFromCallables(const SignatureTag<Signature, DerivativeTraits>& signatureTag, F&&... f)
{
  return DifferentiableFunctionFromCallables<Signature, DerivativeTraits, F...>(f...);
}
 
 
 
} // namespace Functions
} // namespace Dune
 
#endif //DUNE_FUNCTIONS_COMMON_DIFFEREENTIONABEFUNCTIONFROMCALLABLES_HH