阿克玛插值(2D)

#ifndef INTERPOLATION_HEADER

#define INTERPOLATION_HEADER

 

#include<cmath>

#include<set>

#include<vector>

#include<algorithm>

 

using namespace std;

 

#define DELTA 1e-6

#define VERY_CLOSE(x,y) fabs((x)-(y))<DELTA

 

template<typename T>

class Point2D

{

public:

     Point2D(T _X=0,T _Y=0)

         :X(_X),Y(_Y)

     {

     }

     T X,Y;

};

 

bool MyFunc(const Point2D<double>& P1, const Point2D<double>& P2)

{

     return P1.X+DELTA<P2.X;

}

class Point2D_X_Less_Delta

{