Abstract
A rigorous and efficient algorithm which combines the fast multipole method and conjugate gradient fast Fourier transform is presented for the full-wave analysis of large-scale finite-sized periodic structures. The memory requirement and computational complexity of the new algorithm are only of the order of O(N), where N denotes the number of unknowns. Hence the new method can handle very large electromagnetic problem of periodic structures. For example, a problem consisting of 12 millions unknowns can be solved in a personal computer. Numerical results prove the accuracy and efficiency of the new algorithm.