24 double dif = fabs(
x-xPoints[0].first),dift;
25 const unsigned int _K=xPoints.size(),_KP=_K+1;
26 std::vector<double>c(_KP),d(_KP);
28 for (
unsigned int i=0;
i<_K;
i++) {
29 dift=fabs(
x-xPoints[
i].first);
34 c[
i]=d[
i]=xPoints[
i].second;
36 y = xPoints[ns--].second;
37 for (
unsigned int m=0;m<_K-1;m++) {
38 for (
unsigned int i=0;
i<_K-m-1;
i++) {
39 double ho = xPoints[
i].first-
x;
40 double hp= xPoints[
i+m+1].first-
x;
45 <<
"Error in polynomial extrapolation"
51 deltay = 2*(ns+1) < (
int)(_K-m-1) ? c[ns+1]: d[ns--];
58 xPoints.push_back(std::make_pair(
x,y));
62 min=DBL_MAX, max=-DBL_MAX;
63 for (
unsigned int i=0;
i<xPoints.size();
i++) {
64 min = std::min(min,xPoints[
i].first);
65 max = std::max(max,xPoints[
i].first);