massive interpolation...

[Nico]

Hello,

I have a 2D grid 'vxf' with 'nx+1' by 'ny+1' points, representing a
discrete representation of a scalar field which size is 'xm' by 'ym'
The grid is regular and the distance between each point is supposed to
be given by 'dx'

example :

vxf = dblarr(nx+1,ny+1)
vxf = get_values() ; personal function..

x = dx*indgen(nx+1)
y = dy*indgen(ny+1)

then I have about 6 million positions (xp,yp) within the range [0:xm,
0:ym] where I would like to interpolate the value of 'vxf'

Would it be an easy an quick way to do that, without the loop over all
the 6 million positions ?

Thanks,
Nico