It involves the evolution of the boundary of the region to the interior.
Its a non-iterative algorithm. I calculate the values only once. ``flag``
helps me determine which pixel value should be computed next.

You can check out the docstrings for ``_init_fmm`` and
``fast_marching_method``.