GSLBisection.h

// -*- C++ -*-
//
// GSLBisection.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_GSLBisection_H
#define HERWIG_GSLBisection_H
//
// This is the declaration of the GSLBisection class.
//
 
#include "ThePEG/Pointer/ReferenceCounted.h"
#include "Herwig/Utilities/GSLHelper.h"
#include <gsl/gsl_errno.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_roots.h>
 
namespace Herwig {
 
using namespace ThePEG;
 
class GSLBisection : public Pointer::ReferenceCounted {
 
public:
 
  struct GSLerror {};
 
  struct IntervalError {};
 
  GSLBisection() : abserr_(0), relerr_(1.E-8), maxPoints_(100) {}
 
  inline GSLBisection(double abserr, double relerr, int max) :
    abserr_(abserr), relerr_(relerr), maxPoints_(max) {}
 
 
  static void GSLsubstHandler(const char *, const char *, 
                              int, int){
    throw GSLerror();
  }
 
  template <class T>
  inline typename T::ArgType value(const T & function, 
                                   const typename T::ArgType lower,
                                   const typename T::ArgType upper) const;
 
private:
 
  GSLBisection & operator=(const GSLBisection &) = delete;
 
private:
 
  double abserr_;
 
  double relerr_;
 
  int maxPoints_;
};
 
}
 
#include "GSLBisection.tcc"
 
#endif /* HERWIG_GSLBisection_H */