We present a numerical algorithm for an accurate and efficient computation of the convolution frequency domain layered media Green's function with given density function. Instead compressing matrix directly as in classical fast multipole, direct solver, H-matrix algorithms, this new considers translated form original so existing blocks from highly optimized free-space multipole method can be ea...