Genfun::DefiniteIntegral Class Reference

#include <DefiniteIntegral.hh>

Inheritance diagram for Genfun::DefiniteIntegral:

Classes
class	Clockwork

Public Types
enum	Type { CLOSED , OPEN }

Public Member Functions
	DefiniteIntegral (double a, double b, Type=CLOSED)
	DefiniteIntegral (const DefiniteIntegral &)
DefiniteIntegral &	operator= (const DefiniteIntegral &)
	~DefiniteIntegral ()
virtual double	operator[] (const AbsFunction &function) const
unsigned int	numFunctionCalls () const
void	setEpsilon (double eps)
void	setMaxIter (unsigned int maxIter)
void	setMinOrder (unsigned int order)
Public Member Functions inherited from Genfun::AbsFunctional
	AbsFunctional ()
virtual	~AbsFunctional ()
virtual double	operator[] (const AbsFunction &function) const =0

Detailed Description

Author

Definition at line 25 of file DefiniteIntegral.hh.

Member Enumeration Documentation

Type

enum Genfun::DefiniteIntegral::Type

Enumerator
CLOSED
OPEN

Definition at line 30 of file DefiniteIntegral.hh.

{CLOSED, OPEN} Type;

Constructor & Destructor Documentation

DefiniteIntegral() [1/2]

Genfun::DefiniteIntegral::DefiniteIntegral	(	double	a,
		double	b,
		Type	type = CLOSED
	)

Definition at line 119 of file DefiniteIntegral.cc.

                                                                 :
    c(new Clockwork()) {
    c->a              = a;
    c->b              = b;
    c->type           = type;
    c->nFunctionCalls = 0;
    c->MAXITER        = type==OPEN ? 20 : 14; 
    c->EPS            = 1.0E-6;
    c->K              = 5;
  }

DefiniteIntegral() [2/2]

Genfun::DefiniteIntegral::DefiniteIntegral

(

const DefiniteIntegral &

right

)

Definition at line 134 of file DefiniteIntegral.cc.

  :AbsFunctional(right), c(new Clockwork(*right.c)) {
  }

~DefiniteIntegral()

Genfun::DefiniteIntegral::~DefiniteIntegral

(

)

Definition at line 130 of file DefiniteIntegral.cc.

                                      {
    delete c;
  }

Member Function Documentation

numFunctionCalls()

unsigned int Genfun::DefiniteIntegral::numFunctionCalls

(

)

const

Definition at line 219 of file DefiniteIntegral.cc.

                                                        {
    return c->nFunctionCalls;
  }

operator=()

DefiniteIntegral & Genfun::DefiniteIntegral::operator=

(

const DefiniteIntegral &

right

)

Definition at line 138 of file DefiniteIntegral.cc.

                                                                                 {
    if (this!=&right) {
      delete c;
      c = new Clockwork(*right.c);
    }
    return *this;
  }

operator[]()

double Genfun::DefiniteIntegral::operator[] ( const AbsFunction &  function ) const

virtual

Implements Genfun::AbsFunctional.

Definition at line 158 of file DefiniteIntegral.cc.

                                                                          {
 
    const Clockwork::QuadratureRule * rule = c->type==OPEN ? 
      (Clockwork::QuadratureRule *) new Clockwork::XtMidpointQuadratureRule() : 
      (Clockwork::QuadratureRule *) new Clockwork::TrapezoidQuadratureRule();
    double xMult=rule->stepMultiplier();
 
    c->nFunctionCalls=0;
    std::vector<double> s(c->MAXITER+2),h(c->MAXITER+2);
    h[1]=1.0;
    for (unsigned int j=1;j<=c->MAXITER;j++) {
      s[j]=rule->integrate(function, c->a, c->b, j);
      c->nFunctionCalls=rule->numFunctionCalls();
      if (j>=c->K) {
    double ss(0.), dss(0.);
    c->polint(h.begin()+j-c->K,s.begin()+j-c->K,0.0,ss, dss);
    if (fabs(dss) <= c->EPS*fabs(ss)) {
      delete rule;
      return ss;
    }
      }
      s[j+1]=s[j];
      h[j+1]=h[j]/xMult/xMult;
    }
    delete rule;
    throw std::runtime_error("DefiniteIntegral:  too many steps.  No convergence");
    return 0.0;                   // Never get here.
  }

setEpsilon()

void Genfun::DefiniteIntegral::setEpsilon

(

double

eps

)

Definition at line 146 of file DefiniteIntegral.cc.

                                              {
    c->EPS=eps;
  }

setMaxIter()

void Genfun::DefiniteIntegral::setMaxIter

(

unsigned int

maxIter

)

Definition at line 150 of file DefiniteIntegral.cc.

                                                        {
    c->MAXITER=maxIter;
  }

setMinOrder()

void Genfun::DefiniteIntegral::setMinOrder

(

unsigned int

order

)

Definition at line 154 of file DefiniteIntegral.cc.

                                                          {
    c->K=(minOrder+1)/2;
  }

---

The documentation for this class was generated from the following files:

• DefiniteIntegral.hh
• DefiniteIntegral.cc