function [ g ] = push_forward( f, a, resolutions )

[n,m]=size(f);
x=0.5-m/2:m/2-0.5;
y=0.5-n/2:n/2-0.5;
[x,y]=meshgrid(x*resolutions(1),y*resolutions(2));

xt=x-(a(1));
yt=y-(a(2));

g=interp2(x,y,f,xt,yt,'*cubic',0);

end

