 
 	assert(rx >= 0 && rx < sourceInfo.size_x);
 	assert(ry >= 0 && ry < sourceInfo.size_y);
-
+	if (d_source[Index(rx,ry,sourceInfo.size_x)] == kFloatMissing) printf("missing value returned\n");
 	return d_source[Index(rx,ry,sourceInfo.size_x)];
 
 }

			// if interpolated grid points are negative, it means that we are outside
			// of the source area
				
			gp.x >= 0 && gp.y >= 0 &&
				
			// if interpolated grid points are larger than source grid in x or y
			// direction, it means again that we are outside of the area
			
			gp.x < (sourceInfo.size_x-1) && gp.y < (sourceInfo.size_y-1)
		)
		{
			

#include "cuda_helper.h"
#include <thrust/sort.h>

const double kEpsilon = 1e-6;

struct point
{
