Flux-based level set method on rectangular grids and computation of first arrival time functions