The propagation of electromagnetic fields between non-parallel planes based on a spectrum-of-plane-wave analysis is discussed and formulations for an efficient numerical implementation are presented in detail. It is shown that with the help of interpolation techniques, the numerical implementation can be done with the standard uniform fast Fourier transform (FFT) of easy access. Different interpolation techniques are numerically examined, and it turns out that the use of cubic interpolation, together with the uniform FFT, brings both significantly increased computational efficiency and high simulation accuracy. Apart from the aspect of computational efficiency, all formulations in this work are generalized in a fully vectorial manner in comparison to previous works.