Speeding up code which integrates a 2D gpuArray matrix using Simpson's Rule

Question

I have a (real) 2D gpuArray, which I am using as part of a larger code, and now am trying to also integrate the array using the Composite Simpson Rule inside my main loop (sever