Rewrite multiprocessing in `ldgm/reduction.py`

Currently, multiprocessing in `ldgm/reduction.py` is implemented with `pool.imap_unordered` with zipped parameters being passed to `find_reach_set()`. This works, but is pretty inelegant and should be rewritten using a different multiprocessing function`.