@article{530, abstract = {Inclusion–exclusion is an effective method for computing the volume of a union of measurable sets. We extend it to multiple coverings, proving short inclusion–exclusion formulas for the subset of Rn covered by at least k balls in a finite set. We implement two of the formulas in dimension n=3 and report on results obtained with our software.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, journal = {Computational Geometry: Theory and Applications}, pages = {119 -- 133}, publisher = {Elsevier}, title = {{Multiple covers with balls I: Inclusion–exclusion}}, doi = {10.1016/j.comgeo.2017.06.014}, volume = {68}, year = {2018}, } @article{1072, abstract = {Given a finite set of points in Rn and a radius parameter, we study the Čech, Delaunay–Čech, Delaunay (or alpha), and Wrap complexes in the light of generalized discrete Morse theory. Establishing the Čech and Delaunay complexes as sublevel sets of generalized discrete Morse functions, we prove that the four complexes are simple-homotopy equivalent by a sequence of simplicial collapses, which are explicitly described by a single discrete gradient field.}, author = {Bauer, Ulrich and Edelsbrunner, Herbert}, journal = {Transactions of the American Mathematical Society}, number = {5}, pages = {3741 -- 3762}, publisher = {American Mathematical Society}, title = {{The Morse theory of Čech and delaunay complexes}}, doi = {10.1090/tran/6991}, volume = {369}, year = {2017}, } @article{1173, abstract = {We introduce the Voronoi functional of a triangulation of a finite set of points in the Euclidean plane and prove that among all geometric triangulations of the point set, the Delaunay triangulation maximizes the functional. This result neither extends to topological triangulations in the plane nor to geometric triangulations in three and higher dimensions.}, author = {Edelsbrunner, Herbert and Glazyrin, Alexey and Musin, Oleg and Nikitenko, Anton}, issn = {02099683}, journal = {Combinatorica}, number = {5}, pages = {887 -- 910}, publisher = {Springer}, title = {{The Voronoi functional is maximized by the Delaunay triangulation in the plane}}, doi = {10.1007/s00493-016-3308-y}, volume = {37}, year = {2017}, } @inproceedings{836, abstract = {Recent research has examined how to study the topological features of a continuous self-map by means of the persistence of the eigenspaces, for given eigenvalues, of the endomorphism induced in homology over a field. This raised the question of how to select dynamically significant eigenvalues. The present paper aims to answer this question, giving an algorithm that computes the persistence of eigenspaces for every eigenvalue simultaneously, also expressing said eigenspaces as direct sums of “finite” and “singular” subspaces.}, author = {Ethier, Marc and Jablonski, Grzegorz and Mrozek, Marian}, booktitle = {Special Sessions in Applications of Computer Algebra}, isbn = {978-331956930-7}, location = {Kalamata, Greece}, pages = {119 -- 136}, publisher = {Springer}, title = {{Finding eigenvalues of self-maps with the Kronecker canonical form}}, doi = {10.1007/978-3-319-56932-1_8}, volume = {198}, year = {2017}, } @article{718, abstract = {Mapping every simplex in the Delaunay mosaic of a discrete point set to the radius of the smallest empty circumsphere gives a generalized discrete Morse function. Choosing the points from a Poisson point process in ℝ n , we study the expected number of simplices in the Delaunay mosaic as well as the expected number of critical simplices and nonsingular intervals in the corresponding generalized discrete gradient. Observing connections with other probabilistic models, we obtain precise expressions for the expected numbers in low dimensions. In particular, we obtain the expected numbers of simplices in the Poisson–Delaunay mosaic in dimensions n ≤ 4.}, author = {Edelsbrunner, Herbert and Nikitenko, Anton and Reitzner, Matthias}, issn = {00018678}, journal = {Advances in Applied Probability}, number = {3}, pages = {745 -- 767}, publisher = {Cambridge University Press}, title = {{Expected sizes of poisson Delaunay mosaics and their discrete Morse functions}}, doi = {10.1017/apr.2017.20}, volume = {49}, year = {2017}, } @article{1433, abstract = {Phat is an open-source C. ++ library for the computation of persistent homology by matrix reduction, targeted towards developers of software for topological data analysis. We aim for a simple generic design that decouples algorithms from data structures without sacrificing efficiency or user-friendliness. We provide numerous different reduction strategies as well as data types to store and manipulate the boundary matrix. We compare the different combinations through extensive experimental evaluation and identify optimization techniques that work well in practical situations. We also compare our software with various other publicly available libraries for persistent homology.}, author = {Bauer, Ulrich and Kerber, Michael and Reininghaus, Jan and Wagner, Hubert}, issn = { 07477171}, journal = {Journal of Symbolic Computation}, pages = {76 -- 90}, publisher = {Academic Press}, title = {{Phat - Persistent homology algorithms toolbox}}, doi = {10.1016/j.jsc.2016.03.008}, volume = {78}, year = {2017}, } @article{1662, abstract = {We introduce a modification of the classic notion of intrinsic volume using persistence moments of height functions. Evaluating the modified first intrinsic volume on digital approximations of a compact body with smoothly embedded boundary in Rn, we prove convergence to the first intrinsic volume of the body as the resolution of the approximation improves. We have weaker results for the other modified intrinsic volumes, proving they converge to the corresponding intrinsic volumes of the n-dimensional unit ball.}, author = {Edelsbrunner, Herbert and Pausinger, Florian}, journal = {Advances in Mathematics}, pages = {674 -- 703}, publisher = {Academic Press}, title = {{Approximation and convergence of the intrinsic volume}}, doi = {10.1016/j.aim.2015.10.004}, volume = {287}, year = {2016}, } @article{1295, abstract = {Voronoi diagrams and Delaunay triangulations have been extensively used to represent and compute geometric features of point configurations. We introduce a generalization to poset diagrams and poset complexes, which contain order-k and degree-k Voronoi diagrams and their duals as special cases. Extending a result of Aurenhammer from 1990, we show how to construct poset diagrams as weighted Voronoi diagrams of average balls.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, journal = {Electronic Notes in Discrete Mathematics}, pages = {169 -- 174}, publisher = {Elsevier}, title = {{Multiple covers with balls II: Weighted averages}}, doi = {10.1016/j.endm.2016.09.030}, volume = {54}, year = {2016}, } @article{1805, abstract = {We consider the problem of deciding whether the persistent homology group of a simplicial pair (K,L) can be realized as the homology H∗(X) of some complex X with L ⊂ X ⊂ K. We show that this problem is NP-complete even if K is embedded in double-struck R3. As a consequence, we show that it is NP-hard to simplify level and sublevel sets of scalar functions on double-struck S3 within a given tolerance constraint. This problem has relevance to the visualization of medical images by isosurfaces. We also show an implication to the theory of well groups of scalar functions: not every well group can be realized by some level set, and deciding whether a well group can be realized is NP-hard.}, author = {Attali, Dominique and Bauer, Ulrich and Devillers, Olivier and Glisse, Marc and Lieutier, André}, journal = {Computational Geometry: Theory and Applications}, number = {8}, pages = {606 -- 621}, publisher = {Elsevier}, title = {{Homological reconstruction and simplification in R3}}, doi = {10.1016/j.comgeo.2014.08.010}, volume = {48}, year = {2015}, } @article{2035, abstract = {Considering a continuous self-map and the induced endomorphism on homology, we study the eigenvalues and eigenspaces of the latter. Taking a filtration of representations, we define the persistence of the eigenspaces, effectively introducing a hierarchical organization of the map. The algorithm that computes this information for a finite sample is proved to be stable, and to give the correct answer for a sufficiently dense sample. Results computed with an implementation of the algorithm provide evidence of its practical utility. }, author = {Edelsbrunner, Herbert and Jablonski, Grzegorz and Mrozek, Marian}, journal = {Foundations of Computational Mathematics}, number = {5}, pages = {1213 -- 1244}, publisher = {Springer}, title = {{The persistent homology of a self-map}}, doi = {10.1007/s10208-014-9223-y}, volume = {15}, year = {2015}, } @inproceedings{1495, abstract = {Motivated by biological questions, we study configurations of equal-sized disks in the Euclidean plane that neither pack nor cover. Measuring the quality by the probability that a random point lies in exactly one disk, we show that the regular hexagonal grid gives the maximum among lattice configurations. }, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel and Kurlin, Vitaliy}, booktitle = {Proceedings of the 27th Canadian Conference on Computational Geometry}, location = {Ontario, Canada}, pages = {128--135}, publisher = {Queen's University}, title = {{Relaxed disk packing}}, volume = {2015-August}, year = {2015}, } @inbook{10817, abstract = {The Morse-Smale complex can be either explicitly or implicitly represented. Depending on the type of representation, the simplification of the Morse-Smale complex works differently. In the explicit representation, the Morse-Smale complex is directly simplified by explicitly reconnecting the critical points during the simplification. In the implicit representation, on the other hand, the Morse-Smale complex is given by a combinatorial gradient field. In this setting, the simplification changes the combinatorial flow, which yields an indirect simplification of the Morse-Smale complex. The topological complexity of the Morse-Smale complex is reduced in both representations. However, the simplifications generally yield different results. In this chapter, we emphasize properties of the two representations that cause these differences. We also provide a complexity analysis of the two schemes with respect to running time and memory consumption.}, author = {Günther, David and Reininghaus, Jan and Seidel, Hans-Peter and Weinkauf, Tino}, booktitle = {Topological Methods in Data Analysis and Visualization III.}, editor = {Bremer, Peer-Timo and Hotz, Ingrid and Pascucci, Valerio and Peikert, Ronald}, isbn = {9783319040981}, issn = {2197-666X}, pages = {135--150}, publisher = {Springer Nature}, title = {{Notes on the simplification of the Morse-Smale complex}}, doi = {10.1007/978-3-319-04099-8_9}, year = {2014}, } @inbook{10893, abstract = {Saddle periodic orbits are an essential and stable part of the topological skeleton of a 3D vector field. Nevertheless, there is currently no efficient algorithm to robustly extract these features. In this chapter, we present a novel technique to extract saddle periodic orbits. Exploiting the analytic properties of such an orbit, we propose a scalar measure based on the finite-time Lyapunov exponent (FTLE) that indicates its presence. Using persistent homology, we can then extract the robust cycles of this field. These cycles thereby represent the saddle periodic orbits of the given vector field. We discuss the different existing FTLE approximation schemes regarding their applicability to this specific problem and propose an adapted version of FTLE called Normalized Velocity Separation. Finally, we evaluate our method using simple analytic vector field data.}, author = {Kasten, Jens and Reininghaus, Jan and Reich, Wieland and Scheuermann, Gerik}, booktitle = {Topological Methods in Data Analysis and Visualization III }, editor = {Bremer, Peer-Timo and Hotz, Ingrid and Pascucci, Valerio and Peikert, Ronald}, isbn = {9783319040981}, issn = {2197-666X}, pages = {55--69}, publisher = {Springer}, title = {{Toward the extraction of saddle periodic orbits}}, doi = {10.1007/978-3-319-04099-8_4}, volume = {1}, year = {2014}, } @inproceedings{2043, abstract = {Persistent homology is a popular and powerful tool for capturing topological features of data. Advances in algorithms for computing persistent homology have reduced the computation time drastically – as long as the algorithm does not exhaust the available memory. Following up on a recently presented parallel method for persistence computation on shared memory systems [1], we demonstrate that a simple adaption of the standard reduction algorithm leads to a variant for distributed systems. Our algorithmic design ensures that the data is distributed over the nodes without redundancy; this permits the computation of much larger instances than on a single machine. Moreover, we observe that the parallelism at least compensates for the overhead caused by communication between nodes, and often even speeds up the computation compared to sequential and even parallel shared memory algorithms. In our experiments, we were able to compute the persistent homology of filtrations with more than a billion (109) elements within seconds on a cluster with 32 nodes using less than 6GB of memory per node.}, author = {Bauer, Ulrich and Kerber, Michael and Reininghaus, Jan}, booktitle = {Proceedings of the Workshop on Algorithm Engineering and Experiments}, editor = { McGeoch, Catherine and Meyer, Ulrich}, location = {Portland, USA}, pages = {31 -- 38}, publisher = {Society of Industrial and Applied Mathematics}, title = {{Distributed computation of persistent homology}}, doi = {10.1137/1.9781611973198.4}, year = {2014}, } @inbook{2044, abstract = {We present a parallel algorithm for computing the persistent homology of a filtered chain complex. Our approach differs from the commonly used reduction algorithm by first computing persistence pairs within local chunks, then simplifying the unpaired columns, and finally applying standard reduction on the simplified matrix. The approach generalizes a technique by Günther et al., which uses discrete Morse Theory to compute persistence; we derive the same worst-case complexity bound in a more general context. The algorithm employs several practical optimization techniques, which are of independent interest. Our sequential implementation of the algorithm is competitive with state-of-the-art methods, and we further improve the performance through parallel computation.}, author = {Bauer, Ulrich and Kerber, Michael and Reininghaus, Jan}, booktitle = {Topological Methods in Data Analysis and Visualization III}, editor = {Bremer, Peer-Timo and Hotz, Ingrid and Pascucci, Valerio and Peikert, Ronald}, pages = {103 -- 117}, publisher = {Springer}, title = {{Clear and Compress: Computing Persistent Homology in Chunks}}, doi = {10.1007/978-3-319-04099-8_7}, year = {2014}, } @inproceedings{2153, abstract = {We define a simple, explicit map sending a morphism f : M → N of pointwise finite dimensional persistence modules to a matching between the barcodes of M and N. Our main result is that, in a precise sense, the quality of this matching is tightly controlled by the lengths of the longest intervals in the barcodes of ker f and coker f . As an immediate corollary, we obtain a new proof of the algebraic stability theorem for persistence barcodes [5, 9], a fundamental result in the theory of persistent homology. In contrast to previous proofs, ours shows explicitly how a δ-interleaving morphism between two persistence modules induces a δ-matching between the barcodes of the two modules. Our main result also specializes to a structure theorem for submodules and quotients of persistence modules. Copyright is held by the owner/author(s).}, author = {Bauer, Ulrich and Lesnick, Michael}, booktitle = {Proceedings of the Annual Symposium on Computational Geometry}, location = {Kyoto, Japan}, pages = {355 -- 364}, publisher = {ACM}, title = {{Induced matchings of barcodes and the algebraic stability of persistence}}, doi = {10.1145/2582112.2582168}, year = {2014}, } @inproceedings{2155, abstract = {Given a finite set of points in Rn and a positive radius, we study the Čech, Delaunay-Čech, alpha, and wrap complexes as instances of a generalized discrete Morse theory. We prove that the latter three complexes are simple-homotopy equivalent. Our results have applications in topological data analysis and in the reconstruction of shapes from sampled data. Copyright is held by the owner/author(s).}, author = {Bauer, Ulrich and Edelsbrunner, Herbert}, booktitle = {Proceedings of the Annual Symposium on Computational Geometry}, location = {Kyoto, Japan}, pages = {484 -- 490}, publisher = {ACM}, title = {{The morse theory of Čech and Delaunay filtrations}}, doi = {10.1145/2582112.2582167}, year = {2014}, } @inproceedings{2156, abstract = {We propose a metric for Reeb graphs, called the functional distortion distance. Under this distance, the Reeb graph is stable against small changes of input functions. At the same time, it remains discriminative at differentiating input functions. In particular, the main result is that the functional distortion distance between two Reeb graphs is bounded from below by the bottleneck distance between both the ordinary and extended persistence diagrams for appropriate dimensions. As an application of our results, we analyze a natural simplification scheme for Reeb graphs, and show that persistent features in Reeb graph remains persistent under simplification. Understanding the stability of important features of the Reeb graph under simplification is an interesting problem on its own right, and critical to the practical usage of Reeb graphs. Copyright is held by the owner/author(s).}, author = {Bauer, Ulrich and Ge, Xiaoyin and Wang, Yusu}, booktitle = {Proceedings of the Annual Symposium on Computational Geometry}, location = {Kyoto, Japan}, pages = {464 -- 473}, publisher = {ACM}, title = {{Measuring distance between Reeb graphs}}, doi = {10.1145/2582112.2582169}, year = {2014}, } @article{2255, abstract = {Motivated by applications in biology, we present an algorithm for estimating the length of tube-like shapes in 3-dimensional Euclidean space. In a first step, we combine the tube formula of Weyl with integral geometric methods to obtain an integral representation of the length, which we approximate using a variant of the Koksma-Hlawka Theorem. In a second step, we use tools from computational topology to decrease the dependence on small perturbations of the shape. We present computational experiments that shed light on the stability and the convergence rate of our algorithm.}, author = {Edelsbrunner, Herbert and Pausinger, Florian}, issn = {09249907}, journal = {Journal of Mathematical Imaging and Vision}, number = {1}, pages = {164 -- 177}, publisher = {Springer}, title = {{Stable length estimates of tube-like shapes}}, doi = {10.1007/s10851-013-0468-x}, volume = {50}, year = {2014}, } @inproceedings{10897, abstract = {Taking images is an efficient way to collect data about the physical world. It can be done fast and in exquisite detail. By definition, image processing is the field that concerns itself with the computation aimed at harnessing the information contained in images [10]. This talk is concerned with topological information. Our main thesis is that persistent homology [5] is a useful method to quantify and summarize topological information, building a bridge that connects algebraic topology with applications. We provide supporting evidence for this thesis by touching upon four technical developments in the overlap between persistent homology and image processing.}, author = {Edelsbrunner, Herbert}, booktitle = {Graph-Based Representations in Pattern Recognition}, isbn = {9783642382208}, issn = {1611-3349}, location = {Vienna, Austria}, pages = {182--183}, publisher = {Springer Nature}, title = {{Persistent homology in image processing}}, doi = {10.1007/978-3-642-38221-5_19}, volume = {7877}, year = {2013}, }