Merging discrete Morse vector fields: a case of stubborn geometric parallelization
Abstract
We address the basic question in discrete Morse theory of combining discrete gradient fields that are partially defined on subsets of the given complex. This is a well-posed question when the discrete gradient field is generated using a fixed algorithm which has a local nature. One example is ProcessLowerStars, a widely used algorithm for computing persistent homology associated to a grey-scale image in 2D or 3D. While the algorithm for may be inherently local, being computed within stars of vertices and so embarrassingly parallelizable, in practical use it is natural to want to distribute the computation over patches , apply the chosen algorithm to compute the fields associated to each patch, and then assemble the ambient field from these. Simply merging the fields from the patches, even when that makes sense, gives a wrong answer. We develop both very general merging procedures and leaner versions designed for specific, easy to arrange covering patterns.
Cite
@article{arxiv.2111.09295,
title = {Merging discrete Morse vector fields: a case of stubborn geometric parallelization},
author = {Douglas Lenseth and Boris Goldfarb},
journal= {arXiv preprint arXiv:2111.09295},
year = {2022}
}
Comments
18 pages, 3 figures