 The underlying method constructs a signed distance field that extends from the boundary, and can be used to solve the motion of the boundary in this field.The powerful level-set method is due to Osher and Sethian 1988.However, the straightforward implementation via a dense d-dimensional array of values, results in both time and storage complexity of formula_1, where formula_2 is the cross sectional resolution of the spatial extents of the domain and formula_3 is the number of spatial dimensions of the domain.