Parallel Methods and Higher Dimensional NLS Equations

Alternating direction implicit (ADI) schemes are proposed for the solution of the two-dimensional coupled nonlinear Schrödingerequation. These schemes are of second- and fourth-order accuracy in spaceand second order in time. The resulting schemes in each ADI computation step correspond to a block tridiagonal system which can be solvedby using one-dimensional block tridiagonal algorithm with a considerablesaving in computational time. These schemes are very well suited for parallel implementation on a high performance system with many processorsdue to the nature of the computation that involves solving the same blocktridiagonal systems with many right hand sides. Numerical experimentson one processor system are conducted to demonstrate the efficiency andaccuracy of these schemes by comparing them with the analytic solutions. The results show that the proposed schemes give highly accurate results