An algorithm for efficient computation of spatial impulse response on the GPU with application in ultrasound simulation