Classes
class	QuadratureRule

class	TrapezoidQuadratureRule

class	XtMidpointQuadratureRule

Public Member Functions
void	polint (std::vector< double >::iterator xArray, std::vector< double >::iterator yArray, double x, double &y, double &deltay) const

Public Attributes
double	a

double	b

DefiniteIntegral::Type	type

unsigned int	nFunctionCalls

unsigned int	MAXITER

double	EPS

unsigned int	K

Detailed Description

Definition at line 15 of file DefiniteIntegral.cc.

Member Function Documentation

◆ polint()

void Genfun::DefiniteIntegral::Clockwork::polint	(	std::vector< double >::iterator	xArray,
		std::vector< double >::iterator	yArray,
		double	x,
		double &	y,
		double &	deltay
	)		const

Definition at line 187 of file DefiniteIntegral.cc.

                                                                                                                                                        {
    double dif = fabs(x-xArray[1]),dift;
    std::vector<double> cc(K+1),d(K+1);
    unsigned int ns=1;
    for (unsigned int i=1;i<=K;i++) {
      dift=fabs(x-xArray[i]);
      if (dift<dif) {
    ns=i;
    dif=dift;
      }
      cc[i]=d[i]=yArray[i];
    }
    y = yArray[ns--];
    for (unsigned int m=1;m<K;m++) {
      for (unsigned int i=1;i<=K-m;i++) {
    double ho = xArray[i]-x;
    double hp=  xArray[i+m]-x;
    double w=cc[i+1]-d[i];
    double den=ho-hp;
    if (den==0)
      std::cerr
        << "Error in polynomial extrapolation"
        << std::endl;
    den=w/den;
    d[i]=hp*den;
    cc[i]=ho*den;
      }
      deltay = 2*ns < (K-m) ? cc[ns+1]: d[ns--];
      y += deltay;
    }
  }