The plane wave expansion method for calculating magneto-photonic crystals is modified via complex Fourier factorization whose basic principles are briefly explained. The new approach is used to simulate the photonic band structure of a two-dimensional magneto-photonic crystal made as Bi:YIG rods arranged with square periodicity.
The considerable improvement of the convergence properties is demonstrated for selected points of the bands.