  Level set method has been used to capture interface motion. Narrow band
algorithm is applied to localize the solving of level-set PDE on global domain
to a tube around interface. Due to the unknown evolving interface, narrow band
algorithm brings load balance problem for parallelizing computing. This work
presents a tool for evenly distributing work loads on CPUs. On the other hand,
numerically solving level-set PDE only needs simple operations but on large
grid points. This work also presents a GPU acceleration for solving level-set
PDE using finite difference method.
