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"
// ============================================================================

namespace AIDA { class IFunction ; };

namespace Genfun 
{
  namespace GaudiMathImplementation
  { 
    // ========================================================================
    // =========================================================================
    
    class 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 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 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();
      // assigmenet operator is disabled 
      Adapter3DoubleFunction& operator = ( const Adapter3DoubleFunction&);
    private:
      Function                 m_func  ;
    };
    
    // ========================================================================
    // ========================================================================
    class 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 assignement 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
// ============================================================================