The main reason behind this behaviour is that you are trying ot use a general purpose integrator on a function which is not smooth enough to benefit from the high order integration methods used. The lack of smoothness between each piece forces the integrator to go to great expense. Because the newer integration functions (like "integral2") are much more conservative and reliable than the old ones, the expense may be great indeed.
The best (and by far more efficient) way to integrate over an interpolated function is to split between the various pieces, in fact the integral is extremely difficult to integrate at the borders of each piece , while the piece itself is very easy to integrate. You can then iterate through the pieces and add up all the contributions.
As an example:
A2 = 0;
for ix = 1:N-1
xleft = xData(ix);
xright = xData(ix+1);
for iy = 1:N-1
yleft = yData(iy);
yright = yData(iy+1);
A2 = A2 + integral2(G,xleft,xright,yleft,yright,'RelTol',1e-10,'AbsTol',1e-10,'method', 'auto');
end
end
