Parallel Level set algorithm with MPI and accelerated on GPU

Abstract

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.