FFT-based convolution algorithm for fast and precise numerical evaluating diffracted field by photon sieve

FFT-based convolution is proposed to numerical solve Fresnel–Kirchhoff integral in Fresnel regime carefully and in a very shorter time in comparison to direct solving convolution. To show its capability, the algorithm was implemented to evaluate amplitude of a diffracted plane wave at the focal plane of photon sieves with different focal lengths. The calculated amplitudes are completely the same calculated via convolving operation but has advantageous of taking very very shorter time. The calculation was also repeated using single-FFT algorithm that produce same result for all ranges either below or upper the sampling criteria and different results in comparison to the other two methods.