FunAdapters.h

Go to the documentation of this file.

#ifndef GAUDIMATH_FUNADAPTERS_H
#define GAUDIMATH_FUNADAPTERS_H 1
// ============================================================================
// Include files
// ============================================================================
#include <memory>
// ============================================================================
// from CLHEP
// ============================================================================
#include "CLHEP/GenericFunctions/GenericFunctions.hh"
#include "CLHEP/GenericFunctions/Argument.hh"
#include "CLHEP/GenericFunctions/AbsFunction.hh"
// ============================================================================
#include "GaudiKernel/Kernel.h"

#if defined(__clang__) || defined(__CLING__)
#pragma clang diagnostic push
#pragma clang diagnostic ignored "-Winconsistent-missing-override"
#elif defined(__GNUC__) && __GNUC__ >= 5
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wsuggest-override"
#endif

namespace AIDA { class IFunction ; }

namespace Genfun
{
  namespace GaudiMathImplementation
  {
    // ========================================================================
    // =========================================================================

    class GAUDI_API AdapterIFunction : public AbsFunction
    {
    public:
      FUNCTION_OBJECT_DEF( AdapterIFunction )
    public:
      AdapterIFunction ( const AIDA::IFunction& fun) ;
      AdapterIFunction ( const AdapterIFunction& );
      ~AdapterIFunction() override;

      double operator() ( double a          ) const override;

      double operator() ( const Argument& x ) const override;

      unsigned int dimensionality () const override { return m_dim ; }

      bool  hasAnalyticDerivative() const override { return true ; }
      Genfun::Derivative partial( unsigned int i  ) const override;

    private:
      AdapterIFunction();
      AdapterIFunction& operator=(const AdapterIFunction&);
    private:
      const AIDA::IFunction*      m_fun ;
      size_t                      m_dim ;
      mutable std::vector<double> m_arg ;
    };
    FUNCTION_OBJECT_IMP( AdapterIFunction )

    // ========================================================================
    // ========================================================================
    class GAUDI_API Adapter2DoubleFunction : public AbsFunction
    {
    public:
      typedef double (*Function)( const double ,
                                  const double ) ;
    public:
      FUNCTION_OBJECT_DEF( Adapter2DoubleFunction )
    public:
      Adapter2DoubleFunction ( Function func );
      Adapter2DoubleFunction ( const Adapter2DoubleFunction& );
      ~Adapter2DoubleFunction () override;

      double operator() (       double    x ) const override;

      double operator() ( const Argument& x ) const override;

      unsigned int dimensionality() const override { return 2    ; }
      bool  hasAnalyticDerivative() const override { return true ; }
      Genfun::Derivative partial( unsigned int i  ) const override;
    public:
      double operator() ( const double x , const double y ) const ;
    private:
      // default constructor is disabled
      Adapter2DoubleFunction();
      // assigmenet operator is disabled
      Adapter2DoubleFunction& operator = ( const Adapter2DoubleFunction&);
    private:
      Function                 m_func  ;
    };
    FUNCTION_OBJECT_IMP( Adapter2DoubleFunction )

    // ========================================================================
    // ========================================================================
    class GAUDI_API Adapter3DoubleFunction : public AbsFunction
    {
    public:
      typedef double (*Function)( const double ,
                                  const double ,
                                  const double ) ;
    public:
      FUNCTION_OBJECT_DEF( Adapter3DoubleFunction )
    public:
      Adapter3DoubleFunction ( Function func );
      Adapter3DoubleFunction ( const Adapter3DoubleFunction& );
      ~Adapter3DoubleFunction () override ;

      double operator() (       double    x ) const override ;

      double operator() ( const Argument& x ) const override ;

      unsigned int dimensionality() const override { return 3    ; }
      bool  hasAnalyticDerivative() const override { return true ; }
      Genfun::Derivative partial( unsigned int i  ) const override ;
    public:
      double operator() ( const double x ,
                          const double y ,
                          const double z ) const ;
    private:
      // default constructor is disabled
      Adapter3DoubleFunction();
      // assignment operator is disabled
      Adapter3DoubleFunction& operator = ( const Adapter3DoubleFunction&);
    private:
      Function                 m_func  ;
    };
    FUNCTION_OBJECT_IMP( Adapter3DoubleFunction )

    // ========================================================================
    // ========================================================================
    class GAUDI_API SimpleFunction : public AbsFunction
    {
    public:
      typedef double (*Function1)( const double               ) ;
      typedef double (*Function2)( const double*              ) ;
      typedef double (*Function3)( const std::vector<double>& ) ;
    protected:
      enum Case { TrivialArg , ArrayArg , VectorArg } ;
    public:
      FUNCTION_OBJECT_DEF( SimpleFunction )

    public:

      SimpleFunction ( Function1 func ) ;

      SimpleFunction ( Function2    func ,
                       const size_t dim ) ;

      SimpleFunction ( Function3    func ,
                       const size_t dim ) ;

      SimpleFunction ( const SimpleFunction& ) ;

      // destructor
      ~SimpleFunction() override = default;
    public:

      unsigned int dimensionality         () const override { return m_DIM   ; }
      bool  hasAnalyticDerivative         () const override { return true    ; }
      double             operator()  ( double          ) const override ;
      double             operator()  ( const Argument& ) const override ;
      Genfun::Derivative partial     ( unsigned int i  ) const override ;

    private:

      // the default constructor  is disabled
      SimpleFunction() ;
      // the assignment operator is disabled
      SimpleFunction& operator=(const SimpleFunction& );

    private:

      Case                         m_case     ;
      size_t                       m_DIM      ;

      Function1                    m_func1    ;

      Function2                    m_func2    ;
      std::unique_ptr<double[]>    m_arg2     ;

      Function3                    m_func3    ;
      mutable std::vector<double>  m_arg3     ;
    };
    FUNCTION_OBJECT_IMP( SimpleFunction )

  } // end of namespace GaudiMathImeplementation
} // end of namespace Genfun

#if defined(__clang__) || defined(__CLING__)
#pragma clang diagnostic pop
#elif defined(__GNUC__) && __GNUC__ >= 5
#pragma GCC diagnostic pop
#endif

#endif // GAUDIMATH_FUNADAPTERS_H
// ============================================================================