Added bilinear interpolation function for many points and a simple test