DefiniteIntegral.hh

Go to the documentation of this file.

// -*- C++ -*-
// $Id: DefiniteIntegral.hh,v 1.2 2003/09/06 14:04:13 boudreau Exp $
//-------------------------------------------------------------//
//                                                             //
// This functional performs Romberg integration on a function  //
// between lower bound and upper bound b.                      //
//                                                             //
// Two types:  OPEN: use open quadrature formula               //
//                   for improper integrals                    //
//             CLOSED (default) use closed quadrature          //
//                   formula.                                  //
//                                                             //
//-------------------------------------------------------------//
 
#ifndef _DefiniteIntegral_h_
#define _DefiniteIntegral_h_
#include "CLHEP/GenericFunctions/AbsFunctional.hh"
 
namespace Genfun {
 
  /**
   * @author
   * @ingroup genfun
   */
  class DefiniteIntegral:public AbsFunctional {
 
  public:
  
    // Type definition:
    typedef enum {CLOSED, OPEN} Type;
 
    // Constructor:
    DefiniteIntegral(double a, double b, Type=CLOSED);
 
    // Copy Constructor:
    DefiniteIntegral(const DefiniteIntegral &);
 
    // Assignment Operator:
    DefiniteIntegral & operator=(const DefiniteIntegral &) ;
 
    // Destructor:
    ~DefiniteIntegral();
 
    // Take the definite integral of a function between the bounds:
    virtual double operator [] (const AbsFunction & function) const;
 
    // Retrieve the number of function calls for the last operation:
    unsigned int numFunctionCalls() const;
 
    // Algorithmic parameters:
    
    // Desired precision (default 1.0E-06)
    void setEpsilon(double eps);
 
    // Maximum number of iterations (default 20(closed) 14 (open))
    void setMaxIter (unsigned int maxIter);
 
    // Minimum order:
    void setMinOrder (unsigned int order);
 
 
  private:
 
    class Clockwork;
    Clockwork *c;
 
  };
} // namespace Genfun
#endif