combines Gauss-Seidel iterations with alternating sweeping orderings. In this paper several
parallel implementations of the fast sweeping method are presented. These parallel
algorithms are simple and efficient due to the causality of the underlying partial different
equations. Numerical examples are used to verify our algorithms.