Parallel simulations of femtosecond laser pulses propagation in air based on OpenMP

The propagation of axis-symmetric ultra-short high-intensity laser pulse in air can be described by the 2D+1 non-linear Schrodinger equation, which is often solved by the FCN method (split-step Fourier method with Crank-Nicholson scheme in the transverse direction). However, since the FCN method is very time-consuming, the series FCN code can only deal with the simulations of the propagation of laser pulse up to hundreds meters. We develop a parallel algorithm for the FCN method basing on OpenMP. Numerical simulations demonstrate the correctness and high parallel-efficiency of the parallel FCN code. The speed-up ratio is about 12 for 15 threads. The parallel FCN code can be used to simulate the long-distance propagation of the ultra-short high-intensity laser pulses.