Fast Topological Signal Identification and Persistent Cohomological Cycle Matching

Within the context of topological data analysis, the problems of identifying topological significance and matching signals across datasets are important and useful inferential tasks in many applications. The limitation of existing solutions to these problems, however, is computational speed. In this paper, we harness the state-of-the-art for persistent homology computation by studying the problem of determining topological prevalence and cycle matching using a cohomological approach, which increases their feasibility and applicability to a wider variety of applications and contexts. We demonstrate this on a wide range of real-life, large-scale, and complex datasets. We extend existing notions of topological prevalence and cycle matching to include general non-Morse filtrations. This provides the most general and flexible state-of-the-art adaptation of topological signal identification and persistent cycle matching, which performs comparisons of orders of ten for thousands of sampled points in a matter of minutes on


Introduction
Persistent homology, one of the cornerstones of topological data analysis, studies the lifespan of the topological features in a nested sequence of topological spaces by tracking the changes in its homology groups.It provides a robust statistical summary of data, capturing its "shape" and "size" and has been been applied to many scientific disciplines in recent years with great success.The diagram consisting of the homology groups of the filtration connected by the maps induced by the inclusions is called the persistence module.From this, the persistence barcode (or simply, barcode) is derived-a canonical summary of the aforementioned lifespans as a set of half-open intervals.
A natural question is whether it is possible to compare the barcodes obtained from different filtrations, which would, for instance, provide a correspondence between some of their intervals.Several solutions have been proposed.Gonzalez-Diaz & Soriano-Trigueros (2020) derive a basis-independent partial matching for ladder modules and zigzag persistence.A different method of persistent extension to find analogous barsespecially interesting if there is no known mapping between the persistence modules-was very recently introduced by Yoon et al. (2022).Bauer & Lesnick (2015) match intervals of barcodes using a known mapping between the persistence modules.This notion was recently reinterpreted in statistical terms by Reani & Bobrowski (2021b), who propose a similar interval matching using image-persistence, which was first introduced by Cohen-Steiner et al. (2009).This matching is applied to define a prevalence score-a measure of the significance of a given interval in a barcode.Typically, persistence (i.e., the length of the interval) is interpreted as topological significance or signal: longer intervals correspond to "true" features, while short ones are attributed to topological noise.However, this practice can be misleading since persistence is highly affected by the distance of sampling points and usually has a higher value for cycles created at a larger scale.The prevalence score proposed by Reani & Bobrowski (2021b) bypasses this shortcoming by taking into account the statistical heuristics of the problem: it is obtained by matching persistence intervals across diagrams of several resamplings of the data.
A limitation common to all previously proposed barcode comparison techniques is that they are all computationally very expensive, which significantly limits their practicality in many applications and to many real datasets.In this paper, we address this specific issue by leveraging the current state-of-the-art in persistent homology computation, Ripser (Bauer, 2021), which studies the dual perspective and computes persistent cohomology, thus taking advantage of its equivalence to persistent homology (de Silva et al., 2011a).Furthermore, recently, Ripser was adapted to the setting of image-persistence via Ripser-image (Bauer & Schmahl, 2022).We apply this technology to the interval matching approach proposed by Reani & Bobrowski (2021b), as well as specialize and extend their definitions to allow for greater flexibility and applicability.The final result of our contributions is state-of-the-art software for interval matching, executable in a matter of minutes using only standard institutional high performance computing facilities, which we showcase on a wide variety of complex and large-scale datasets, such as static and time-lapse imaging and video data from biomedical and astrophysical applications.

Contributions
• We specialize the definition of interval matching proposed by Reani & Bobrowski (2021b) to simplexwise filtrations (see Definition 15), making it compatible with the output of Ripser-image (Bauer & Schmahl, 2022) and more widely applicable.
• We present a comprehensive case study of different definitions for a matching affinity score that extends the original score proposed by Reani & Bobrowski (2021b).
• We comprehensively showcase and demonstrate representative applications of our specialized definitions and code to complex and large-scale datasets.

Outline
We begin by introducing the fundamentals of persistent homology and set relevant notations in section 1.We also review image-persistence and use it to present interval matching as proposed in Reani & Bobrowski (2021b).In section 2, we adapt the definition of image-persistence to the various homology settings and study how these frameworks are related.Here, we propose our specialized definition of interval matching and revisit the notion of matching affinity by Reani & Bobrowski (2021b) in a case study of alternative formulations.In section 3, we present applications of the notion of cycle matching to a variety of data sets diverse in nature and aimed at different objectives.We close with a discussion of our contributions and proposals for future work in section 4.

Preliminaries
In this section we introduce the fundamental concepts underlying our work and establish some relevant notations that we will use throughout the rest of the paper.

The Four Standard Persistence Modules
A filtration is a family of nested subspaces {X t : t ∈ T } of some space X X t ⊂ X s ⊂ X, for t ≤ s, where T ⊂ R is a totally ordered indexing set.In this paper, we work with filtered complexes, specifically, we further assume that X is a finite simplicial complex and the spaces X t are simplicial subcomplexes of X.
Filtered complexes can also be interpreted as diagrams X • : T → Simp of simplicial complexes indexed over some finite totally ordered set T , such that all maps in the diagram are inclusions.
A re-indexing of a filtration changes the indexing set T to another totally ordered set I using some monotonic map r : I → T so that X r(i) = X t .For instance, if {X t : t ∈ T } is a filtered complex, T = {t 1 , . . ., t n } is finite and we have the re-indexing r(i) = t i that allows for a reparameterization of the filtration over the natural numbers Applying the corresponding homology functor to the simplicial complexes in a filtered complex and the inclusions X i ⊂ X i+1 between consecutive spaces gives us the following diagrams (1.4) Following de Silva et al. (2011a), we will call these the four standard persistence modules.The first persistence module (1.1) corresponds to absolute homology and is the one most often used.The expressions following are the persistence modules for absolute cohomology (1.2), relative homology (1.3), and relative cohomology (1.4).Unless otherwise stated, homology and cohomology will have field coefficients, so that these persistence modules are made up of vector spaces and linear maps.
The assumption of field coefficients allows us to invoke the structure theorem (Zomorodian & Carlsson, 2005).This is one of the foundational results in persistent homology and ensures that, up to isomorphism, any persistence module, such as the ones above, can be decomposed in a direct sum of interval modules.An interval module consists of copies of the field of coefficients over an interval range of indices; these copies are connected by the identity map and the trivial vector space outside of that interval.This allows for the interpretation that some (co)cycle is born at the beginning of the interval and dies at the end of it.For instance, for the absolute homology module, where the sub-index denotes the range of indices over which the interval module is nontrivial.
The collection of intervals that appears in the decomposition of the structure theorem is an invariant of the isomorphism type of the persistence module.This collection is the persistence barcode of the filtration . The intervals from the barcode are called persistence intervals and the start and end points of the intervals are the persistence pairs.
The persistence barcode provides a summary of the lifespans of the topological features of the filtration.Persistence pairs are often interpreted with real indices t bm and t dm associated to the natural indices b m and d m via re-indexing.In this case, the convention dictates that the barcode be represented by half-open intervals.These intervals exclude the real-valued death time of the subsequent step of the filtration, in relation to when the feature was last designated "alive" using natural indices: This convention also involves setting t 0 = −∞-notice that the index i = 0 might appear in the barcodes of relative (co)homology-and t M +1 = ∞.It is also customary to discard intervals where t bm = t dm+1 .
Figure 1: Filtration of a triangle indexed with natural numbers.
The barcodes for such a filtration for the 0-and 1-dimensional homologies are Considering the re-indexing given by the diameter of the larger simplex in the complex t i = max{diam (σ) : σ ∈ X i }, the previous barcodes would then become so that we may discard the 1-dimensional persistent homology.
The four standard persistence modules carry the same information: the barcode of the setting of absolute (or relative) homology setting is the same as the barcode of the setting of absolute (or relative) cohomology.
We can also find a bijection between the bars in the barcodes of the relative setting and the bars in the corresponding absolute setting.For further the details on this equivalence, see de Silva et al. (2011a).

Image-Persistence and Interval Matching
The idea underlying image-persistence is to study the persistent homology of some filtered complex inside another larger filtered complex.Let X and Z be finite simplicial complexes and f : X → Z an injective map between them.Let {X i : 1 ≤ i ≤ n} and {Z i : 1 ≤ i ≤ n} be filtrations associated to the previous complexes and denote the restrictions to the steps in the filtrations by Note that these are also injective maps, which gives rise to the following commutative diagram for all 1 ≤ i ≤ n − 1: where ι X i and ι Z i are the corresponding inclusion maps between consecutive steps in the corresponding filtration.
Applying the homology functor to the previous diagram gives rise to another commutative diagram: which now involves the homology groups and the induced linear maps.The commutativity of this diagram allows for the following definition.
Definition 2 (Image-persistent homology).The persistence module given by the subspaces Im(H * (f i )) ⊂ H * (Z i ) and the restrictions of the maps H * (ι Z i ) is called image-persistent homology.
The elements in Im (H * (f i )) can be seen as cycles in X i up to boundaries in Z i , which we gain by studying one filtration inside another.
Remark 3. Since Im H * (f i ) is a subspace of H * (Z i ), a death in the image-persistence module implies a death in the persistent homology of the space Z.At a topological level, this implication can be linked to the fact that every cycle from an image-persistence module is in fact a cycle in the codomain of the function used to define the image-persistence module.This cycle may have been born before as a cycle of the codomain, however, it gets bounded at the same time for both the image-persistence and the persistence module of the codomain.On the other hand, a birth in the image-persistence module implies a birth in the persistent homology of the space X.This is a result from the fact that every cycle in the image persistence corresponds to a cycle of the domain, which may get bounded in the image-persistence module at an earlier time as it is studied within the "larger" codomain.See Cohen-Steiner et al. (2009) for further details on the relations between these three modules.
The definition of interval matching introduced by Reani & Bobrowski (2021b) is based on image-persistence to compare the persistence bars of two diagrams and is restricted to Morse filtrations.
Definition 4 (Morse filtration).A filtration {X t : t ∈ R} is a Morse filtration if there exists a finite set T = {t 1 , ..., t n } ⊂ R such that the following are satisfied: 1.For all t / ∈ T , there exists ϵ > 0 small enough such that for every 0 < ϵ ′ < ϵ the map induced by inclusion is an isomorphism for every homology group.Equivalently, the homology does not change at t.
2. For all t ∈ T , there exists ϵ > 0 small enough so that for any 0 < ϵ ′ < ϵ either ) is injective and the dimension of the vector space increases by one, or ) is surjective and the dimension decreases by one.
Equivalently, the homology changes allowed are either the creation of a single new cycle or the termination of a single existing cycle.
We now review how to match the persistence intervals of two filtered complexes inside a third comparison space.Let X, Y, Z be finite simplicial complexes with Morse filtrations and {Z i : 1 ≤ i ≤ n}.Assume we have injective maps for every 1 ≤ i ≤ n such that f j | Xi = f i and g j | Yi = g i for every i ≤ j.With these assumptions, Reani & Bobrowski (2021b) match persistence intervals as follows.
Definition 5 (Matching intervals, Reani & Bobrowski (2021b) From the definition of image-persistence, this also holds in the respective image-persistence modules.Recall that a birth in the image-persistence module means a birth in the corresponding persistent homology module of X or Y .This allows each bar from the image-persistence module to have the same birth time as exactly one bar in the associated persistent homology.
From the death perspective, recall that a death happening in any of the image-persistence modules means a death in H * (Z • ).Thus, there is also at most one bar in each image-persistence diagram dying at any given time.Consequently, each bar from an image-persistence module can share the same death time with at most one bar from the other image-persistence module.These notions of uniqueness induced by assuming Morse filtrations guarantee that there are no ambiguous matchings.

Matching Affinity and Prevalence Score
The prevalence score was proposed by Reani & Bobrowski (2021b) as an alternative measure to persistencein the sense of interval length-as an indicator for topological significance in noisy data.It takes inspiration from bootstrapping techniques, which is a well-known and powerful subsampling with replacement method, originally proposed in the statistical literature by Efron (1982).
The formulation of the prevalence score takes into account the inherent tendency that as the sample size grows, noisy generators tend to reappear frequently.Due to this tendency, the affinity of a match must first be discussed before prevalence may be considered.Affinity is a score assigned to every match that considers the lifetimes of the persistent cycles and image-persistent cycles involved in the definition of interval matching.
Recall that the Jaccard index of two intervals I and J is given by where ∥•∥ denotes the length of an interval.
Definition 7 (Matching affinity, Reani & Bobrowski (2021b)).The matching affinity of two bars α, β matching through their image-bars α, β is defined as the product With this definition, the prevalence score may now be formally introduced.
Definition 8 (Prevalence score, Reani & Bobrowski (2021b)).Given some reference space X = X ref and resampling spaces X (1) , . . ., X (K) , any α ∈ Pers(H * (X • )) has a prevalence score defined by where k) as a comparison space and the inclusions into the union as the connecting maps f : X ref → Z and g : X (k) → Z (if a matching does not exist, set ρ = 0).

Clearing Algorithm and Cohomology
The basic algorithm to compute the persistent homology of a filtered complex is based on reducing each column of the matrix of the boundary operator on the complex by adding columns on its left, from left to right, to obtain a reduced matrix.From this reduced matrix, the barcode can be readily attained.
Chen & Kerber (2011) proposed an optimization of this process called the clearing algorithm, based on the observation that since the matrix to reduce comes from a boundary operator, some columns in the reduced matrix must be null after the reduction and do not play a role in the reduction process.The clearing algorithm then reduces the boundary matrix in blocks from right to left so that it becomes possible to detect these null columns beforehand and set them directly to 0. In that way, it is possible to avoid reducing these columns and thereby accelerate the computation.However, the increase in speed is burdened by the large number of columns in the first block that must be reduced in the boundary matrix.
In an application to compute circular coordinates, de Silva et al. (2011b) observed that computing persistent cohomology was generally faster than computing persistent homology.This phenomenon was later confirmed by Bauer et al. (2017), who also realized that this increase in speed was coming from an implicit use of the clearing algorithm to compute persistent cohomology by de Silva et al. (2011b).One of the contributions of Ripser (Bauer, 2021), which also implements this optimization, is to provide a formal argument for this increase in speed: the advantage of using the clearing algorithm in persistent cohomology stems from the fact that in the relative coboundary matrix, the first block to reduce is significantly smaller.From the reduced coboundary matrix, the barcode for the relative cohomology setting (1.4) is read off, which is equivalent to the barcode for the homology setting (1.1) as established by de Silva et al. (2011a).
We note in particular that Ripser only considers Vietoris-Rips persistent homology, which is based on the Vietoris-Rips filtration, and does not compute other filtrations.The Vietoris-Rips filtration is a standard filtration often considered in computational applications and settings.Recall that the Vietoris-Rips filtration of a finite metric space (P, d) is VR • (P, d) where the simplicial complex at filtration value ϵ is By applying the homology functor, we obtain the Vietoris-Rips persistent homology H * (VR • (P, d)).

Cycle Matching in the Setting of Cohomology
In this section, we present our theoretical contributions.Specifically, we address the current gaps in the literature by providing a comprehensive account of the extension of image-persistence to the outstanding settings of the four standard persistence modules and the relations between such modules.Subsequently, we specialize the definition of interval matching to simplex-wise filtrations and outline how to implement our specialization using Ripser-image.We finish the section with a case study of alternative definitions of the matching affinity.

The Four Image-Persistence Modules
We only need functoriality applied to the following commutative diagram for 1 ≤ i ≤ n − 1 to define image-persistence, which can then be easily extended to obtain four imagepersistence modules as parallels to the four standard persistence modules presented previously in Section 1.1.
In the setting of absolute cohomology, applying the corresponding homology functor gives us the following commutative diagram for every 1 ≤ i ≤ n − 1.Since we are working with field coefficients, the objects with superscripts are dual to the objects with subscripts from the diagram for homology.The commutativity of the diagram allows for the following definition.
Definition 9 (Image-persistent cohomology).The image-persistent cohomology is defined as the persistence module Im H * (f given by the subspaces Im(H * (f i )) ⊂ H * (X i ) and the restrictions of the maps H * (ι X i ).We now consider the relative settings.Recall that we have a map f : for every 1 ≤ i ≤ n.Since we also have X i ⊂ X j ⊂ X, for i ≤ j, we can write for the identity in X.The same can be written for the identity ι Z This gives us the following commutative diagram Using functoriality in the relative setting we obtain the subsequent commutative diagrams where again, the subscripts and superscripts mean duality.Observe the notation of the homology functor applied to f : (X, X i ) → (Z, Z i ).Commutativity again allows for the following definitions.
Definition 10 (Image-persistent relative homology).The image-persistent relative homology is the persistence module given by the vector spaces Im(H * (f, f i )) ⊂ H * (Z, Z i ) and the restrictions of the linear maps H * (ι Z ).
Definition 11 (Image-persistent relative cohomology).The image-persistent relative cohomology is the persistence module given by the vector spaces Im(H * (f, f i )) ⊂ H * (X, X i ) and the restrictions of the linear maps H * (ι X ).
Remark 12.The four image-persistence modules above currently appear in Bauer & Schmahl (2021) as an example of application of a broader theory of lifespan functors, i.e., endofunctors in the category of persistence modules and matching diagrams-equivalent to barcodes-which are related to boundedness properties of the intervals.Their work is precisely motivated by the extension of the duality results by de Silva et al.
(2011a) to the setting of image-persistence.A version of a result from Bauer & Schmahl ( 2021) is then revisited in Bauer & Schmahl (2022) in Proposition 3.12 in order to implement Ripser-image, which plays a central role in our work.This result will be further explained in Section 2.2 and is fully referenced in Proposition 14.

Equivalence Among the Four Image-Persistence Settings
A natural question to ask after introducing the four image-persistence modules is whether we can expect equivalences among them akin to the ones proved by de Silva et al. (2011a) for the standard persistence modules.In search of a first immediate answer to this question, we check directly whether the persistence modules of homology and cohomology provide the same information.
Proposition 13.The following equalities hold: Proof.It is sufficient to prove that the maps naturally induced between the images, , for all i ≤ j have the same rank.This is true since where the second equality (2.1) is given by duality and the third equality (2.2) is given by the commutativity of the diagram in absolute cohomology.Since persistence barcodes are uniquely determined by dimensions and ranks, we have shown that the image-persistent absolute homology and image-persistent cohomology barcodes are the same.The same proof applies to prove equality of the barcodes in the relative setting.
As for equivalence between the absolute and relative settings, the arguments used in de Silva et al. (2011a) are not directly applicable for image-persistence.However, if Pers 0 denotes the finite intervals and Pers ∞ the infinite intervals of a given persistence barcode, we have the following correspondence.
Proposition 14 (Bauer & Schmahl (2022), Proposition 3.12).We have Additionally, the map I → T \ I defines the bijections This result means that in order to determine the barcode of Im H * (f • ), it suffices to compute Pers ∞ (H * (X, X • )) and Pers 0 (Im H * (f, f • )).Both of these persistence diagrams may be computed applying a matrix reduction algorithm to appropriate boundary matrices.Bauer & Schmahl (2022) also show that the clearing algorithm implemented in Ripser (see Section 1.4) can be applied to compute image-persistence.In this way, the code for Ripser can be fully adapted to this setting to achieve state-of-the-art computations for image-persistence.

Matching Intervals in Non-Morse Filtrations
Ripser-image provides the barcode of the image-persistent homology for Vietoris-Rips filtrations.As noted previously in Remark 6, this presents a significant obstacle to implement efficient interval matching using Ripser-image: Vietoris-Rips filtrations are not Morse filtrations.To overcome this limitation, we introduce a specialization of Definition 5 that resolves the matches between bars with shared birth or death time.We first recall the definition of simplex-wise filtration.
Additionally, it is a simplex-wise filtration if for every i ∈ I such that X i ̸ = ∅ there is some simplex σ i and some index j < i, such that X i ∖ X j = {σ i }.
Observe that in simplex-wise filtrations, there is a bijection between the indices of the filtration and the simplices of the complex X.Consequently, the persistence pairs in the intervals of the barcode can be associated with particular simplices.We call the simplices corresponding to birth times positive simplices and those corresponding to death times negative simplices; we say that an interval is created by its positive simplex and destroyed by the negative simplex.
In addition, simplex-wise filtrations are Morse filtrations, which allows us to directly apply the definition of interval matching proposed by Reani & Bobrowski (2021b) (Definition 5): given the correspondence between birth (resp.death) times and positive (resp.negative) simplices, we can rephrase Definition 5 in the following manner.Let X, Y, Z be finite simplicial complexes with simplex-wise filtrations {X i : i ∈ I}, {Y i : i ∈ I}, and {Z i : i ∈ I}.Assume we have injective maps f : X → Z and g : Y → Z with the usual notation for the restrictions for every i ∈ I. such that the following conditions are satisfied: • α and α are created by the same simplex (seen in X and in f (X), respectively); • β and β are created by the same simplex (seen in Y and in g(Y ), respectively); • α and β are destroyed by the same simplex in Z.
In Definition 16, the use of the phrase "by the same simplex" in relation to intervals α and α means that if σ is the positive simplex associated to α, then f (σ) is the positive simplex associated to α; this meaning also applies to the bars β and β with their positive simplices connected by the function g.Notice that this association is well defined since both f and g are assumed to be injective.Similarly, the notion of α and β being destroyed by the same simplex means that the negative simplices associated to these intervals are precisely represented by a single simplex within the larger complex Z.
Again, many constructions of filtered complexes do not comply with Definition 15, in particular, Vietoris-Rips filtrations are not simplex-wise in general.However, for any filtered complex we can always find a re-indexing that refines the filtration and turns it into an essential simplex-wise filtration.This can be done by finding a partial ordering of the simplices in each of the steps of the original filtration that extends to a total ordering in the whole complex.For instance, Ripser uses a lexicographic refinement of the Vietoris-Rips filtration, which orders the simplices by dimension, diameter and a combinatorial system (see Bauer (2021) for further detail).
Consequently, restricting the matching to simplex-wise filtrations does not introduce any additional difficulties; rather, it enables the application of cycle matching to general filtrations.The process described above is a standard procedure in persistent homology solvers (i.e., algorithms that compute persistent homology barcodes), which compute the barcode for the refined simplex-wised filtration and then retrieve the barcode of the original filtration by relabelling the endpoints of the intervals as explained in Section 1.1.This process can be extended to interval matching in a similar manner: given a general filtration, we refine it to a simplex-wise filtration, compute the interval matching following Definition 16, and then recover the interval matching for the original filtration.This approach addresses a gap in the work of Reani & Bobrowski (2021b) and thus makes their ideas more general and widely applicable, increasing their practical utility.
Remark 17.Notice that to implement the interval matching in Definition 16 for general filtrations, we need to further assume that underlying simplex-wise refinements are compatible in the three filtrations.This means that the simplices are added to the persistence modules H * (X • ) and H * (Y • ) and to their image-persistence modules Im H * (f • ) and Im H * (g • ) in the same order.This can always be achieved by first setting the orders in X and Y and then in Z accordingly.

Implementing Cycle Matching with Ripser-Image
The input to Ripser-image consists of two Vietoris-Rips filtrations where the two metrics d, d ′ on a finite set P satisfy d(p, q) ≥ d ′ (p, q) for all p, q ∈ P.However, the setting for interval matching is slightly different.
From section 1.2, to implement interval matching on finite point clouds X and Y living in the same ambient space, we can consider their union P = X ∪ Y and any metric d ′ on P induced from that ambient space.This will induce metrics X is obtained by setting a very large distance between any point in X and any point in Y, and any pair of points in Y, all seen in the union.Then, up to a threshold corresponding to that large distance and up to the points in P \ X = Y, we have X ) which puts us in the setting of Ripser-image.Notice that in the matrix representation of d X there are rows and columns corresponding to points in Y.This construction can be also applied to (Y, d Y ), preserving the same order as before in rows and columns to ensure compatibility.In this manner, we obtain the three Vietoris-Rips filtrations and we consider the inclusions as the connecting functions for the matching.
The code for Ripser and Ripser-image assigns a unique index to any simplex of the input filtered complex using a lexicographic refinement.However, it does not provide the indices associated to the positive and negative simplices of the persistence intervals of the barcode in its output.These values can be readily retrieved with a slight change the original code.By arranging the matrices representing the finite metric spaces as explained above, these indices allow for the implementation of Definition 16 without affecting the computational runtime of both programs.
It is crucial to observe at this point that a change in the order of the columns and rows within the distance matrices used as input for Ripser-image, which are in correspondence with the points of X and Y in the setting described above, will alter the indices assigned to simplices through the lexicographic refinement.This could, in principle, alter the outcome of the interval matching using Ripser-image as described thus far.Nonetheless, as long as we are consistent with the ordering in all of the three matrices involved in the computations, this observation poses no problem for the implementation of the interval matching.In our code, the points in X occupy the foremost positions, preserving the same order in all three matrices, and the points in Y assume the terminal positions, also preserving their order across matrices.
A Note on Terminology: Cycle Matching.Reani & Bobrowski (2021b) refer to this framework of matching intervals in persistent homology as cycle registration.In this paper, whenever we use the term "cycle matching," we assume that there is a way of finding cycles in the final simplicial complex of the filtration that correspond to the intervals in its barcode.Note that, in general, these representative cycles are not unique-in fact, the persistence intervals are associated to homology classes, i.e., equivalence classes of cycles.However, there are methods to uniquely determine the representative cycles.One such method considers the columns of the reduced boundary matrix corresponding to the killing simplices, which provides the simplices that make up a cycle killed by that simplex.
Further on in section 3, we implement cycle matching by using the version ripser-tight-representative-cycles of Ripser.This feature provides the representative cycles corresponding to the intervals in the barcode computed by Ripser.Note that Ripser does not reduce the boundary matrix but the relative coboundary matrix (see the previous discussion from section 1.4), and thus we cannot implement the aforementioned method directly.However, Čufar & Virk (2021) develop an adaptation of this idea to obtain state-of-the-art computations of barcodes and representatives.The method proposed in Čufar & Virk (2021) uses persistent cohomology to obtain the persistence pairs and reduces the boundary matrix only using the columns corresponding to death indices.This is the technique that ripser-tight-representative-cycles implements.
Processing Multiple Jobs in Parallel.A computational advantage of the bootstrapping approach proposed by Reani & Bobrowski (2021b) is that this technique is parallelizable, despite the inherently nonparallelizable nature of persistent homology computations.Recall that once the barcode of the reference sample is computed, to obtain the prevalence score of its intervals, these intervals are matched with the intervals in the barcodes of K-many resamplings.These matchings can be processed in parallel jobs using a high performance computer cluster (HPC) with a workload manager and job scheduling system, such as SLURM or OpenPBS.In each of these jobs, first, the barcode of the corresponding resampling and the barcodes of the image-persistence modules involved are computed, then the generalized interval matching is implemented.This allows for a dramatic increase in efficiency in the computational runtime, with respect to a sequential execution of the code: the total computational time corresponds to the one of the slowest job, instead of the sum of the computational times of the individual jobs.

Revisiting Matching Affinity
The matching affinity introduced in Definition 7 relies on a particular choice of pairs of intervals to compare through their Jaccard index.In principle, other selections are also valid to obtain different definitions of the matching affinity.We now study the behavior of four such affinities in the example of two circles with same radius but diverging centers.This will allow us to conclude that only one of these definitions exhibits a significant difference with respect to the others.
The following concrete example provides an intuition on how the different affinities behave.Consider two circles of radius 1 and centers shifted by a distance s.We expect that the matching affinity decreases as the center-to-center distance s increases, until reaching 0 (i.e., no match) beyond a certain value.The result of this experiment is displayed in Figure 2, where we see that all matching affinities decrease with respect to s and that the cutoff value is 1.Affinities A, B, and C follow very similar decreasing behaviors in a linear fashion, whereas affinity D has a distinct plateau-like behavior.We now further investigate this phenomenon.
Assume that we have two persistence bars α and β matched via their image-bars α and β.We know that the birth times of α and α, and β and β coincide, respectively, and that the bars α and β also have the same death time.Thus, having high affinity A, for instance, depends on the two following phenomena.Firstly, the bars α and β must have similar birth and death times.This means that the cycles in X and Y that generated them should be similar in size.Secondly, the death time of the image-bars should also be similar to the death time of the original bars.Geometrically, this means that the cycles that generated the bars should have high overlapping surface when considered in the union of the point clouds.These ideas are illustrated in Figure 3, where the affinity A of a match of two circles decreases significantly when the circles have different centers or radii.
Returning to Figure 2 with these ideas in mind, we see that there are few noticeable but not fundamental differences between the affinities A, B, and C. Indeed, the Jaccard indices Jac(α, β) and Jac(α, β) have similar magnitude-if a pair of cycles are similar in size in the spaces X and Y , the cycles corresponding to their image-bars in the union will also have similar sizes.This implies that the affinities A and B are very similar.Affinity C is slightly lower than affinities A and B only because it has an extra multiplicative factor with magnitude less than one.
It also makes sense that the matching affinities A, B, and C drop when the spatial overlap between cycles decreases.This happens because these affinities include the Jaccard indices between bars and image-bars, which are influenced by this overlap.The matching affinity D does not consider such a comparison, and thus, remains at a higher value until there is no match anymore, when it drops abruptly.Such a behavior could be useful in certain situations.In real-life applications, it might happen that we are interested in matching topological features that shrink or enlarge significantly, or which appear misplaced in the samples.Matching affinity D would then be more sensitive to these matches, by assigning a higher prevalence score to them.However, this feature could be undesirable in other contexts, as we discuss in the next observation.
We now need to check whether the four affinities are consistent with the original motivation, which is the condition proposed in Reani & Bobrowski (2021b).Namely, random cycles that appear in resamplings and get matched several times should be assigned a low prevalence score.Similar to Reani & Bobrowski (2021b), we consider a uniform sampling of the unit square with N ref = 1000 points and compute the prevalence scores of its bars by finding matches with K = 20 different resamplings of N = 1000 points from that same distribution.The results of this experiment are given in Figure 4.For affinity D, some random cycles are assigned quite high prevalence scores in the range 0.6-0.7.This must be taken into account when interpreting the affinity scores for applications using affinity D.

Applications
In this section, we demonstrate with numerous examples that cycle matching can be applied to real-life, large-scale, complex datasets from biology and astrophysics.Several applications motivate the usage of cycle matching and prevalence.First, we can identify common topological features shared by two spaces as a direct application of cycle matching.We can use this to track features both spatially and over time, on consecutive slices of an object or on consecutive time frames.We demonstrate both of these applications in this section.
Second, the most prevalent features in data can be identifed by applying cycle matching repeatedly after resampling from the same distribution.By doing this, we can detect prevalent cycles in large-scale and complex data.Here prevalence is computed using Affinity A exclusively.We demonstrate this on cosmic web data and cell actin network data.Prevalence gives rise to an enriched visualization via the prevalenceaugmented barcode, where length corresponds to persistence while thickness and color both correspond to prevalence.
Comparison with Classical Computer Vision Tools for Feature Tracking.It is important to note at this point that the application we are suggesting here differs from existing algorithms for tracking features in computer vision, such as those implemented in OpenCV (Bradski, 2000).These algorithms take images of an object of interest as input and use different methods to track the object through subsequent frames, identifying an area in the frame in which the object is located.The technique we are proposing fundamentally differs from the procedure above in four aspects.First, we track topological features, which are structures that are inherently different from objects and object locations in images.Moreover, we do not need any prior knowledge of what to track: persistent homology directly detects the topological features of the image, and our only input is the video or stack of images where the topological features are present.In addition, our output is the chain of features from the persistent homology of the images that are matched in subsequent frames, which we can represent through cycle representatives as explained above, but which is not in principle related to any specific area in the image, as is the case of the output of the feature tracking algorithms.Lastly, notice that we can track several features concurrently, without needing to re-run the whole algorithm, which is what usually happens in feature tracking algorithms in computer vision.

Tunneling: Tracking Intervals Over Slices
As a first application of cycle matching, we tracked intervals over two-dimensional slices of three-dimensional objects.In biomedical imaging, for instance, it is common that data are made up of spherical or tubular elements, such as in vessels and other biological organs with channeling functions.Using cycle matching, we can match the closed contours delimiting the spherical or tubular elements across slices.
Data: Lateral Line in Zebrafish.To demonstrate this application, we used a biological imaging dataset, in particular the dataset with image ID 9836972 provided in Hartmann et al. (2020).This is a stack of twodimensional confocal images from the zebrafish posterior lateral line primordium (pLLP).The pLLP is a primitive expression of the lateral line-an organ in fish that allows them to detect the pattern of water flow over their body surface.It appears at the embryonic stage in the form of a rosette-shaped cluster of cells.
The circular contours that we can see in the images of the dataset (see Figure 5) are precisely these cells; we will track these along the height of the stack.
We considered a stack of 15 images with a 0.66 µm gap and size of 300 × 300 pixels in each image, with a resolution of 0.1 µm per pixel.We thresholded the images with the Otsu method and uniformly sampled them with N = 1000 points.We matched persistent intervals between pairs of Vietoris-Rips filtrations on consecutive slices.Some features were matched on consecutive frames and formed a tunnel: we were able to detect 32 such tunnels, each stained with a different color.We computed the geometric generators drawn here to represent the persistence intervals that were matched using the ripser-tight-representative-cycles module of Ripser (Bauer, 2021).The results are shown in Figure 5.As a byproduct of this approach, we can identify slices on which a cell appears then disappears.

Video Data: Tracking Features Over Time
Cycle matching can be used to track topological features over time, by matching the barcodes of consecutive frames in a video, or, in a biological context, at different stages of disease development.This method detects common topological patterns surviving across consecutive time points and quantifies the quality of the match through the affinity scores.
Data: Heart Valves in Zebrafish.To illustrate this application, we analyzed a video of the atrioventricular valve (AVV) of a wild-type AB zebrafish, from Scherz et al. (2008).This video is taken 76 hours post fecundation and at a rate of 50 milliseconds.The specimen studied comes from a transgenic line that allows the monitoring of the two chambers that make up the primitive heart of the zebrafish embryo, and how its contraction over time generates embryonic heartbeats.The contraction is especially pronounced for the right chamber, as can be seen from Figure 6.
We selected 10 frames capturing one contraction and matched cycles on consecutive frames (Figure 6).We sampled N = 500 points on each of the images after applying a thresholding technique based on the mean of gray-scale values.We successfully detected the persistent intervals delimiting the two chambers, and tracked them on all consecutive frames.We also tracked their size variation.Note that the matching affinities are much more variable for the cycle on the right (in red), which is expected since the right chamber changes abruptly in shape.As before, we show on each frame of Figure 6 the generator associated to each persistence interval, obtained with the ripser-tight-representative-cycles module of Ripser, while using the same color to stain matched cycles.
Data: Time-Lapse Images of Human Embryos.As another example to demonstrate feature tracking over time, we tracked intervals on 10 consecutive frames of time-lapse embryo data from Gomez et al. (2022) (Figure 7).A time-lapse imaging (TLI) system with a special camera captures images of a human embryo every 50 to 100 minutes and features different stages of cellular division.We matched samples of N = 500 points on the images after applying a Sato operator and a threshold using the Otsu method.We were able in particular to detect cell division as the appearance of a new topological feature (in red).

Prevalent Cycles
Our next application is to find prevalent cycles in order to reveal significantly organized topological patterns in data.We will demonstrate this on cosmic web data (Figure 8) and cell actin data.Recall that this consists of comparing multiple resamplings X (1) , . . ., X (K) to a reference space X ref and finding all possible matching pairs of persistence intervals between X ref and any This approach becomes especially useful in situations where the initial (unknown) distribution has high topological complexity, and for which we only have access to point cloud samplings.In other situations, the distribution may be partially or even entirely known already-for instance if we are given an image from which to sample points.Oftentimes, an image may suffer from noise or contrast variations, but we can recover the true cycles of the object.Even in the absence of noise, brighter and weaker signals may still capture interesting information (for instance, depth in 2D images)-prevalence takes this into account.
We can also study the profile of prevalence scores corresponding to an image to characterize its topological structure in comparison to another image.We may visualize this using prevalence-augmented persistence barcodes, where the length of a bar still represents persistence in the usual sense, while its thickness and color represents prevalence.
Data: The Cosmic Web.First, we identified prevalent cycles in the cosmic web, based on the point distribution of galaxies from the BOSS CMASS database (Dawson et al., 2013).Matter in the universe is arranged along an intricate pattern involving filamentary structures (de Lapparent et al., 1986;York et al., 2000).However, the reconstruction of these filaments is still challenging; multiple methods to address this reconstruction problem have been proposed (Malavasi et al., 2020).Instead of detecting filaments with an uncertainty score (e.g., Duque et al., 2022), we propose to detect cycles with a prevalence score.
The final version of the BOSS CMASS dataset used here is included in the current data release SDSS DR17 (Abdurrouf et al., 2022) and was released in SDSS DR12 (Alam et al., 2015).We selected galaxies with right ascension 170 < RA < 190, declination 30 < dec < 50 and redshift range 0.564 < z < 0.57 and projected the points onto the (RA, dec) 2D space.
We sampled the reference space X ref using N ref = 1000 points and resampled the dataset K = 20 times with 300 points in each resampling X (k) , by adding Gaussian noise of magnitude 0.1 to each point.The noise scale, required for the prevalence score bootstrapping technique, is determined based on both the data magnitude and the width of the filaments present in the original sample that we want to detect.We performed 20 comparisons of barcodes to find matching cycles.The results are shown in Figure 8, where cycles with different prevalence scores can be visualized.
Data: Cell Actin.Next, we computed the prevalence of cycles in biological imaging data of cell actin.Actin networks are essential in scaffolding the inner structure of cells, enabling in particular cell motility and reshaping.Figure 9, featuring data from Svitkina & Borisy (1999), shows significant loss of actin filaments in the rear of the lamellipodium network due to the absence of some stabilizing chemicals during extraction.We selected three crops I, II and III, to study from this image, where the actin filaments are sparse, half-sparse and half-dense, and dense, respectively.We then thresholded each cropped image using the Otsu thresholding method to segment filaments, and restricted the original image pixel intensities to these filaments to finally obtain a discrete probability distribution.Points were sampled from this distribution and their spatial coordinates were perturbed with a Gaussian noise of standard deviation 10% of a pixel side.
We show the results of our computations in Figures 10,11,12,13,14,and 15.We found that larger voids in the structure of the actin mesh led to larger cycles of higher prevalence, such as in crop I, whereas small voids in denser parts led to smaller cycles of lower prevalence, such as in crop III.As seen in Figures 10 and  11, we correctly identified large cycle representatives in crop I, small ones in crop III, and a mix of small and large ones in crop II, at the transitional region between rear and front of the lamellipodium.The usual persistence barcodes (Figure 12) show numerous short-lived bars in crop III, but some longer-lived features for crops I and II.This barcode information can be enriched by visualizing the prevalence as the thickness and color of a bar (Figure 13).It is interesting to note that the highest prevalence scores throughout selected data were found in crop II, corresponding to two highly prevalent cycles (see Figure 14).Their scores are higher than those from crop I, due to contrast variations of larger amplitude along filaments of crop II.Indeed, prevalence can be interpreted as a certainty measure of topological features.Finally, a scatter plot of prevalence versus persistence scores (Figure 15) confirmed that prevalence is not monotonous with respect to persistence and longer intervals do not necessarily correspond to more prevalent features.

A Note on Computational Runtime
As mentioned previously in section 2.4, the advantage of the framework proposed by Reani & Bobrowski (2021b) is that the procedure of matching is easily parallelizable.With access to standard institutional high performance computing (HPC) resources requiring simply CPU processing, use of a single node, one CPU per task, and max 30 GB of memory per CPU, the problem of identifying prevalent cycles or matching intervals in spaces with a range of 100-1000 points generally reduces to a matter of minutes, ranging from seconds to a few hours for the applications showcased here.
We present in Table 1 the runtimes corresponding to the real datasets from Section 3.1 and Section 3.2, where we track topological features over a set of frames.We include the number of points in the samples N and the number of subsamples K.In Table 2, we exhibit the runtimes associated to the real datasets shown in Section 3.3, where we compute prevalent features.We include the number of points on the reference space N ref , the number of points in the resamples N , and the number of resamples K. Table 3 collects runtimes for computing prevalent features on synthetic datasets that consist of point clouds uniformly sampled in the unit square of the plane.Note that we took N ref = N in the synthetic examples.
In Table 1, one runtime corresponds to the computation of 1. the barcodes of the two samplings on consecutive frames; 2. the two image-barcodes of those samplings in their union; and 3. the matching itself, which compares the barcodes.
In Table 2 and Table 3, step (i) only covers the computation of the barcode on the resampling X (k) corresponding to that job, and in step (ii), we compute the image-persistence of the reference space X ref and the resampling X (k) in their union.The runtime needed to compute the barcode of the reference space X ref , which is just a few seconds and is computed once before the parallel jobs, is not included in Tables 2 and 3.
The computational bottleneck here is not the number of matchings K but the number of points N and N ref sampled in the respective spaces for each matching instead.The number of matchings could be increased arbitrarily (up to the capacity of the HPC) without affecting the computational runtime while maintaining it on the order of minutes or a few hours.Increasing it here allows for more precise estimations of, for instance, the average runtime needed to find all possible matchings between two spaces.However, median runtime increases following a power law T ∼ cst • N p where cst is a constant and p ≃ 3.266 (see Figure 16), so that increasing the number of sampled points N by a factor of 10 would result in increasing runtime by a factor of 10 p ≃ 1845.16).

Software and Data Availability
The code used to perform all experiments here is freely and publicly available at our project GitHub repository https://github.com/inesgare/interval-matching.It is fully adaptable for individual user customization.
Where possible, the data we have used in this section are also provided on the same GitHub repository so that all experiments and examples in our paper are fully reproducible.Note that some of the data we have used here required an institutional materials transfer agreement, so these data were not made available on our repository.

Discussion
In this paper, we studied the problem of identifying topologically significant features in noisy data where the usual measure of persistence in the sense of the length of an interval in a persistence barcode is unsatisfactory.We also studied the problem of comparing barcodes over different filtrations and identifying correspondences between persistence intervals.To date, the various existing proposals to these problems have faced significant computational limitations.The main contribution of our work is an extension of existing notions of topological significance and cycle matching to provide the most general and flexible definitions, which we then implement using the dual perspective of cohomology to achieve the fastest available identification of prevalent cycles as well as cycle matching.Our implementation now makes these approaches practical and applicable to real-life, large-scale, complex datasets with execution times ranging from a matter of minutes to a few hours for 100-1000 sampled points, using only standard institutional HPC facilities.Our work inspires several directions for future research, which we now discuss.
First, one natural question is to understand the behavior of the prevalence-augmented barcodes as the number of resampling spaces K and the number of sampling points N for each space increase to infinity.It will be important to understand how at the limit augmented barcodes are related to the original distribution, quantify the rate and type of convergence, and study how the choice of filtration (Rips, Čech, coupled-Alpha (Reani & Bobrowski, 2021a)) and of affinity (A, B, C, D) affects the convergence, if at all.This would then allow us to design probabilistically-founded statistical tests, for example, to determine whether a collection of point clouds has been sampled from a specific distribution with known barcode, or not, based on some confidence intervals.Developing such a framework would provide a practical approach to choosing a threshold to identify "true" cycles in the original data as bars whose prevalence are higher than 1 − ϵ.This could also have practical implications on certain applications, for example, by directly extracting "true" topological signal in complex settings such as imaging data, and perhaps even bypassing the need for computationally expensive procedures, such as image segmentation.
Another theoretical question of interest is the study of the new metric d IMp on persistence modules introduced by Reani & Bobrowski (2021b) in section 6.2.This metric involves comparing matched intervals between two modules directly, and not comparing birth-death times between persistence diagrams which discards important information about spatial correspondence, as explained previously.It is reasonable to expect that as N increases, the distance between two persistence modules converges to zero if the point clouds are drawn from the same distribution.Likewise, this could be an alternative approach towards the design of a (non-parametric) statistical test to determine whether two points clouds were sampled from a same distribution.

Figures
Figure 2: Mean and standard deviation of the affinities of the matches between two circles of radius 1 with centers shifted according to the horizontal axis.The circles were sampled with N = 100 points without noise added.We considered 15 equidistant distances between 0 and 1 and took 15 samples at each step.3. Median runtime increases as a power law T ∼ cst N 3.266 with respect to N (a linear regression on the logarithmic scale gave a score of 0.996, coefficient 3.266 and intercept −14.087 with respect to the data points of Table 3).Left: linear-linear scale.Right: log-log scale.The dashed line corresponds to the result of the linear regression on the log-log scale.

Figure 3 :
Figure3: Three matches between samples of 100 points of two circles with Gaussian noise of magnitude 0.1 added.In Figure3a; the circles have the same radii and centres; in Figure3bthe circles have the same radii and centres 0.7 units of length apart; and in Figure3c, the circles have coinciding centres and radii 1 and 1.5.

Figure 4 :
Figure 4: Most matched intervals in resamplings of N = 1000 points of the uniform distribution in the unit square.Top left: Frequency of reappearance of the 15 most frequently matched intervals.Remainder of the figure: Cycles representing the persistence intervals of X ref stained by their prevalence score using the affinity specified above of the image.

Figure 5 :
Figure 5: Cycle matching to track cell contours on images of slices of the posterior Lateral Line Primordium (pLLP) of the zebrafish.We applied an Otsu threshold and took samples of N = 1000 points on the images.Cycles matched across consecutive slices can be grouped into tunnels, each stained with a different color.Data courtesy of Hartmann et al. (2020).

Figure 6 :
Figure6: Cycle matching on 10 frames from a video courtesy ofScherz et al. (2008).We took samples of N = 500 points and applied a threshold based on the mean of the gray-scale values before sampling.The intervals matched are stained in the same color.Below each image we display the affinity of the match between the interval in the image and the corresponding interval in the subsequent image.

Figure 7 :
Figure 7: Cycle matching on the time-lapse embryo dataset(Gomez et al., 2022) between samples of N = 500 points in the images after applying a Sato operator and a threshold using the Otsu method.Cycles that get matched are stained in the same color.Below each image one can find the affinity of the match between the interval in the image and the corresponding interval in the previous image.

Figure 8 :
Figure 8: Prevalent cycles in the cosmic web.Cycle representatives are stained by prevalence score and galaxies (from the original BOSS CMASS data) are shown as blue dots.We formed the reference space by sampling the galaxies with N ref = 1000 perturbed points and performed K = 20 comparisons to spaces formed by sampling with N = 300 perturbed points each.

Figure 9 :
Figure 9: Electron micrograph of actin network in Xenopus keratocyte lamellipodium, whose rear part disassembled in the course of unprotected extraction and front part remained dense as in control cells.Selected crops I, II, III (from left to right) are shown as red rectangles.Original image CIL:24800 is from the Cell Image Library database (Ellisman et al., 2021), available under CC BY-NC-SA 3.0 License, corresponding to Figure 6b of Svitkina & Borisy (1999).

Figure 10 :
Figure 10: Prevalent cycles of reference space X with N ref = 1200 points (shown as blue dots), based on K = 30 resampling spaces of N = 500 points.Cycle representatives are stained by prevalence score.From left to right: crops I, II, III.Colorbar describes prevalence scores.

Figure 11 :
Figure 11: Prevalent cycles overlayed on original image (see Figures 9 and 10).From left to right: crops I, II, III.Colorbar describes prevalence scores.
Figure 11: Prevalent cycles overlayed on original image (see Figures 9 and 10).From left to right: crops I, II, III.Colorbar describes prevalence scores.

Figure 12 :
Figure 12: Persistence barcodes of reference space X.Values on the horizontal axis correspond to birth time and length of a bar to persistence.From left to right: crops I, II, III.

Figure 13 :
Figure 13: Prevalence-augmented persistence barcodes of reference space X.Thickness and color of a bar correspond to prevalence score.Values on the horizontal axis correspond to birth time and length of the bars to persistence.From left to right: crops I, II, III.Colorbar describes prevalence scores.

Figure 14 :
Figure 14: Prevalence scores, sorted by birth time of the persistence interval.One dot stands for one interval.From left to right: crops I, II, III.

Figure 15 :
Figure 15: Persistence vs prevalence scores in the augmented barcode of the reference space X.One dot stands for one persistence interval.From left to right: crops I, II, III.

Figure 16 :
Figure 16: Median runtime v.s.number of sampling points N for the synthetic examples of Table3.Median runtime increases as a power law T ∼ cst N 3.266 with respect to N (a linear regression on the logarithmic scale gave a score of 0.996, coefficient 3.266 and intercept −14.087 with respect to the data points of Table3).Left: linear-linear scale.Right: log-log scale.The dashed line corresponds to the result of the linear regression on the log-log scale.
(f • ) arises from a bar in the barcode of the persistence module H * (X • ) so that both share the same birth time; we have the same respective coincidence for the modules H * (g • ) and H * (Y • ).This justifies our procedure to first match a bar α ∈ Pers(H * (X • )) and a bar β ∈ Pers(H * (Y • )) with the bars α ∈ Pers(H * (f • )) and β ∈ Pers(H * (g • )) when their birth times coincide.Similarly, the death of a bar in an image-persistence module implies that there is a bar in the barcode of H * (Z • ) sharing that same death time.Consequently, once we have identified the bars α and β, the rule to match them is that they both share the same death time, and thus are related to the same persistence interval of the module H * (Z • ).Remark 6.The Morse assumption is crucial in order for the notion in Definition 5 to be well-defined.Having Morse filtrations for X and Y ensures that there is at most one birth at each time in H * (X • ) and H * (Y • ).

Table 1 :
Computational Runtimes of Tracking Experiments on Real Data.

Table 2 :
Computational Runtimes of Prevalence Experiments on Real Data.

Table 3 :
Computational Runtimes of Prevalence Experiments on Synthetic Data.Median runtime increases as a power law T ∼ cst • N 3.266 with respect to N (see Figure