GSLFunAdapters.h

Go to the documentation of this file.

#ifndef GAUDIMATH_GSLFUNADAPTERS_H
#define GAUDIMATH_GSLFUNADAPTERS_H 1
// ============================================================================
// Include files
// ============================================================================
#include <memory>
// ============================================================================
// from CLHEP
// ============================================================================
#include "CLHEP/GenericFunctions/GenericFunctions.hh"
#include "CLHEP/GenericFunctions/Argument.hh"
#include "CLHEP/GenericFunctions/AbsFunction.hh"
// ============================================================================
// GaudiGSL/GaudiMath
// ============================================================================
#include "GaudiMath/FunAdapters.h"
// ============================================================================
// GSL
// ============================================================================
#include "gsl/gsl_mode.h"
#include "gsl/gsl_sf_result.h"
// ============================================================================
#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 Genfun
{
  namespace GaudiMathImplementation
  {
    typedef Genfun::GaudiMathImplementation::SimpleFunction GSLFunction ;

    class GAUDI_API GSLFunctionWithError : public AbsFunction
    {
    public:
      typedef int (*Function)( double , gsl_sf_result* ) ;
    public:
      FUNCTION_OBJECT_DEF( GSLFunctionWithError )
    public:
      GSLFunctionWithError ( Function function ) ;
      GSLFunctionWithError ( const GSLFunctionWithError& ) ;
      ~GSLFunctionWithError () override  = default;
      double operator() (       double    x ) const override;
      double operator() ( const Argument& x ) const override;
      bool         hasAnalyticDerivative () const override { return true ; }
      unsigned int dimensionality        () const override { return    1 ; }
      Genfun::Derivative partial ( unsigned int i ) const override;

      Function             function () const ;
      const gsl_sf_result& result   () const ;
      double               error    () const ;
    private:
      // default constructor is disabled
      GSLFunctionWithError () ;
      // assignement operator is disabled
      GSLFunctionWithError& operator=( const GSLFunctionWithError & ) ;
    private:
      Function        m_function ;
      std::unique_ptr<gsl_sf_result>  m_result   ;
    };
    FUNCTION_OBJECT_IMP( GSLFunctionWithError )

    class GAUDI_API GSLFunctionWithMode : public AbsFunction
    {
    public:
      typedef double (*Function)( double , gsl_mode_t ) ;
    public:
      FUNCTION_OBJECT_DEF( GSLFunctionWithMode )
    public:
      GSLFunctionWithMode ( Function          fun  ,
                            const gsl_mode_t& mode ) ;
      GSLFunctionWithMode ( const GSLFunctionWithMode& ) ;
      ~GSLFunctionWithMode () override = default;
      double operator() (       double    x ) const override;
      double operator() ( const Argument& x ) const override;
      bool         hasAnalyticDerivative () const override { return true ; }
      unsigned int dimensionality        () const override { return    1 ; }
      Genfun::Derivative partial ( unsigned int i ) const override ;
      Function             function () const ;
      const gsl_mode_t&    mode     () const ;
    private:
      // default constructor is disabled
      GSLFunctionWithMode () ;
      // assignement operator is disabled
      GSLFunctionWithMode& operator=( const GSLFunctionWithError & ) ;
    private:
      Function        m_function ;
      gsl_mode_t*     m_mode     ;
    };
    FUNCTION_OBJECT_IMP( GSLFunctionWithMode )

    class GAUDI_API GSLFunctionWithModeAndError : public AbsFunction
    {
    public:
      typedef int (*Function)( double          ,
                               gsl_mode_t      ,
                               gsl_sf_result*  ) ;
    public:
      FUNCTION_OBJECT_DEF( GSLFunctionWithModeAndError )
    public:
      GSLFunctionWithModeAndError ( Function          fun  ,
                                    const gsl_mode_t& mode ) ;
      GSLFunctionWithModeAndError ( const GSLFunctionWithModeAndError& ) ;
      ~GSLFunctionWithModeAndError () override = default;
      double operator() (       double    x ) const override;
      double operator() ( const Argument& x ) const override;
      bool         hasAnalyticDerivative () const override { return true ; }
      unsigned int dimensionality        () const override { return    1 ; }
      Genfun::Derivative partial ( unsigned int i ) const override;
      Function             function () const ;
      const gsl_mode_t&    mode     () const ;
      const gsl_sf_result& result   () const ;
      double               error    () const ;
    private:
      // default constructor is disabled
      GSLFunctionWithModeAndError () ;
      // assignement operator is disabled
      GSLFunctionWithModeAndError& operator=
      ( const GSLFunctionWithModeAndError & ) ;
    private:
      Function        m_function ;
      std::unique_ptr<gsl_mode_t>     m_mode     ;
      std::unique_ptr<gsl_sf_result>  m_result   ;
    };
    FUNCTION_OBJECT_IMP( GSLFunctionWithModeAndError )

  } // end of namespace GaudiMath Implemnentation
} // 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_GSLFUNADAPTERS_H
// ============================================================================