  A parallel 2D+1 split-step Fourier method with Crank-Nicholson scheme running
on multi-core shared memory architectures is developed to study the propagation
of ultra-short high-intensity laser pulses in air. The parallel method achieves
a near linear speed-up with results for the efficiency of more than 95% on a
24-core machine. This method is of great potential application in studying the
long-distance propagation of the ultra-short high intensity laser pulses.
