The performances of different methods used for the computation of the gradient of the free space Green's function for doubly periodic infinite arrays are investigated, for spacings smaller than the wavelength. An optimal transition between plane wave and cylindrical wave decompositions is proposed and an extension of the method is proposed for the large wavelength case. The joint use of different wave decompositions, combined with the Levin accelerator, allows to obtain exponential convergence for all observation points.
Guerin, N., Craeye, C., & Dardenne, X. (2009). Accelerated Computation of the Free Space Green’s Function Gradient of Infinite Phased Arrays of Dipoles. IEEE Transactions on Antennas and Propagation, 57(10), 3430-3434. https://doi.org/10.1109/TAP.2009.2028684 (Original work published 2009)