The basic of the BPM technique in the frequency domain relies on treating the slowly varying envelope of the monochromatic electromagnetic field under paraxial propagation, thus allowing efficient numerical computation in terms of speed and allocated memory. In addition, the BPM based on finite differences is an easy way to implement robust and efficient computer codes.