FunAdapters.h

Go to the documentation of this file.

// $Id: FunAdapters.h,v 1.2 2005/11/25 10:27:03 mato Exp $
// ============================================================================
#ifndef GAUDIMATH_FUNADAPTERS_H
#define GAUDIMATH_FUNADAPTERS_H 1
// ============================================================================
// Include files
// ============================================================================
// from CLHEP
// ============================================================================
#include "CLHEP/GenericFunctions/GenericFunctions.hh"
#include "CLHEP/GenericFunctions/Argument.hh"
#include "CLHEP/GenericFunctions/AbsFunction.hh"
// ============================================================================
#include "GaudiKernel/Kernel.h"

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& );
      virtual ~AdapterIFunction() ;

      virtual double operator() ( double a          ) const ;

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

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

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

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

    // ========================================================================
    // ========================================================================
    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& );
      virtual ~Adapter2DoubleFunction () ;

      virtual double operator() (       double    x ) const ;

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

      virtual unsigned int dimensionality() const { return 2    ; }
      virtual bool  hasAnalyticDerivative() const { return true ; }
      virtual Genfun::Derivative partial( unsigned int i  ) const ;
    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  ;
    };

    // ========================================================================
    // ========================================================================
    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& );
      virtual ~Adapter3DoubleFunction () ;

      virtual double operator() (       double    x ) const ;

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

      virtual unsigned int dimensionality() const { return 3    ; }
      virtual bool  hasAnalyticDerivative() const { return true ; }
      virtual Genfun::Derivative partial( unsigned int i  ) const ;
    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  ;
    };

    // ========================================================================
    // ========================================================================
    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
      virtual ~SimpleFunction();
    public:

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

    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    ;
      double*                      m_arg2     ;

      Function3                    m_func3    ;
      mutable std::vector<double>  m_arg3     ;
    };

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


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