“…Nevertheless, the impedance matrix filling time remains O((N a N e ) 2 ), the same as for the original MoM, which is prohibitive for large arrays. To overcome this obstacle, different techniques have been proposed, such as the multipole expansion algorithm [20], [21], the interpolatory technique [22], the perturbation technique based on the domain Green's function [23], or the combination of multipole and Contour-FFT approaches for printed structures [24]. The multipole expansion, however, suffers from the instability of near-field interaction, which appears when antennas are too close to each other, i.e.…”