 #include <NFmiStereographicArea.h>
 #include <NFmiGrid.h>
 #include "cuda_helper.h"
+#include <thrust/sort.h>
+
+const double kEpsilon = 1e-5;
 
 struct point
 {

}

__device__
double Linear(double dx, double left, double right)
{
	// return (1 - dx) * left + dx * right;

		)
		{
			
			//targetInfo.interpolation = himan::kNearestPointValue;

			switch (targetInfo.interpolation)
			{
				case himan::kBiLinear:
					interp = BiLinearInterpolation(d_source, sourceInfo, gp);
					break;
					
				case himan::kNearestPoint:
					interp = NearestPointInterpolation(d_source, sourceInfo, gp);
					break;
					
				case himan::kNearestPointValue:
					interp = NearestPointValueInterpolation(d_source, sourceInfo, gp);
					break;
			}
		}

