Parallel Level set algorithm with MPI and accelerated on GPU

TL;DR

This paper introduces a parallel level set algorithm that balances load distribution on CPUs and accelerates PDE solving on GPUs, enhancing efficiency in interface motion simulations.

Contribution

It presents a load balancing tool for parallel CPU computations and implements GPU acceleration for solving level-set PDEs using finite difference methods.

Findings

Effective load balancing on CPUs for narrow band algorithms.

Significant speed-up in PDE solving on GPUs.

Improved efficiency in interface motion simulations.

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.