Genfun::InterpolatingPolynomial Class Reference

#include <InterpolatingPolynomial.hh>

Inheritance diagram for Genfun::InterpolatingPolynomial:

Public Member Functions
	InterpolatingPolynomial ()
	InterpolatingPolynomial (const InterpolatingPolynomial &right)
virtual	~InterpolatingPolynomial ()
virtual double	operator() (double argument) const override
virtual double	operator() (const Argument &a) const override
void	addPoint (double x, double y)
void	getRange (double &min, double &max) const
Public Member Functions inherited from Genfun::AbsFunction
	AbsFunction ()
	AbsFunction (const AbsFunction &right)
virtual	~AbsFunction ()
virtual unsigned int	dimensionality () const
virtual double	operator() (double argument) const =0
virtual double	operator() (const Argument &argument) const =0
virtual AbsFunction *	clone () const =0
virtual FunctionComposition	operator() (const AbsFunction &f) const
virtual ParameterComposition	operator() (const AbsParameter &p) const
Derivative	derivative (const Variable &v) const
Derivative	prime () const
virtual bool	hasAnalyticDerivative () const
virtual Derivative	partial (unsigned int) const

Detailed Description

Author

Definition at line 21 of file InterpolatingPolynomial.hh.

Constructor & Destructor Documentation

InterpolatingPolynomial() [1/2]

Genfun::InterpolatingPolynomial::InterpolatingPolynomial

(

)

Definition at line 13 of file InterpolatingPolynomial.cc.

  :Genfun::AbsFunction()
  {}

InterpolatingPolynomial() [2/2]

Genfun::InterpolatingPolynomial::InterpolatingPolynomial

(

const InterpolatingPolynomial &

right

)

Definition at line 17 of file InterpolatingPolynomial.cc.

    :Genfun::AbsFunction(),xPoints(right.xPoints)
  {}

~InterpolatingPolynomial()

Genfun::InterpolatingPolynomial::~InterpolatingPolynomial ( )

virtual

Definition at line 21 of file InterpolatingPolynomial.cc.

                                                    {
  }

Member Function Documentation

addPoint()

void Genfun::InterpolatingPolynomial::addPoint	(	double	x,
		double	y
	)

Definition at line 60 of file InterpolatingPolynomial.cc.

                                                            {
    xPoints.push_back(std::make_pair(x,y));
  }

getRange()

void Genfun::InterpolatingPolynomial::getRange	(	double &	min,
		double &	max
	)		const

Definition at line 64 of file InterpolatingPolynomial.cc.

                                                                          {
    min=DBL_MAX, max=-DBL_MAX;
    for (unsigned int i=0;i<xPoints.size();i++) {
      min = std::min(min,xPoints[i].first);
      max = std::max(max,xPoints[i].first);
    }
  }

operator()() [1/2]

virtual double Genfun::InterpolatingPolynomial::operator() ( const Argument &  a ) const

inlineoverridevirtual

Implements Genfun::AbsFunction.

Definition at line 38 of file InterpolatingPolynomial.hh.

{return operator() (a[0]);}

operator()() [2/2]

double Genfun::InterpolatingPolynomial::operator() ( double  argument ) const

overridevirtual

Implements Genfun::AbsFunction.

Definition at line 24 of file InterpolatingPolynomial.cc.

                                                            {
    double y=0.0;
    double deltay=0;  // also gives error;
    double dif = fabs(x-xPoints[0].first),dift;
    const unsigned int _K = (unsigned int)xPoints.size(),_KP=_K+1;
    std::vector<double>c(_KP),d(_KP);
    int ns=0;
    for (unsigned int i=0;i<_K;i++) {
      dift=fabs(x-xPoints[i].first);
      if (dift<dif) {
    ns=i;
    dif=dift;
      }
      c[i]=d[i]=xPoints[i].second;
    }
    y = xPoints[ns--].second;
    for (unsigned int m=0;m<_K-1;m++) {
      for (unsigned int i=0;i<_K-m-1;i++) {
    double ho = xPoints[i].first-x;
    double hp=  xPoints[i+m+1].first-x;
    double w=c[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;
    c[i]=ho*den;
      }
      deltay = 2*(ns+1) < (int)(_K-m-1) ? c[ns+1]: d[ns--];
      y += deltay;
    }
    return y;
  }

Referenced by operator()().

---

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

• InterpolatingPolynomial.hh
• InterpolatingPolynomial.cc