Source author record

Primoz Skraba

Primoz Skraba appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Central Limit Theorem for Euclidean Minimal Spanning Acycles

We investigate asymptotics for the minimal spanning acycles of the (Alpha)-Delaunay complex on a stationary Poisson process on $\mathbb{R}^d, d \geq 2$. Minimal spanning acycles are topological (or higher-dimensional) generalization of minimal spanning trees. We establish a central limit theorem for total weight of the minimal spanning acycle on a Poisson-Delaunay complex. Our approach also allows us to establish central limit theorems for sum of birth times and lifetimes in the persistent diagram of the Delaunay complex. The key to our proof is in showing the so-called weak stabilization of minimal spanning acycles which proceeds by establishing suitable chain maps and uses matroidal properties of minimal spanning acycles. In contrast to the proof of weak-stabilization for Euclidean minimal spanning trees via percolation-theoretic estimates, our weak-stabilization proof is algebraic in nature and provides an alternative proof even in the case of minimal spanning trees.

preprint2022arXiv

On the Universality of Random Persistence Diagrams

One of the most elusive challenges within the area of topological data analysis is understanding the distribution of persistence diagrams. Despite much effort, this is still largely an open problem. In this paper, we present a series of novel conjectures regarding the behavior of persistence diagrams arising from random point-clouds. We claim that these diagrams obey a universal probability law, and include an explicit expression as a candidate for what this law is. We back these conjectures with an exhaustive set of experiments, including both simulated and real data. We demonstrate the power of these conjectures by proposing a new hypothesis testing framework for individual features within persistence diagrams.

preprint2021arXiv

Euler Characteristic Surfaces

We study the use of the Euler characteristic for multiparameter topological data analysis. Euler characteristic is a classical, well-understood topological invariant that has appeared in numerous applications, including in the context of random fields. The goal of this paper is to present the extension of using the Euler characteristic in higher-dimensional parameter spaces. While topological data analysis of higher-dimensional parameter spaces using stronger invariants such as homology continues to be the subject of intense research, Euler characteristic is more manageable theoretically and computationally, and this analysis can be seen as an important intermediary step in multi-parameter topological data analysis. We show the usefulness of the techniques using artificially generated examples, and a real-world application of detecting diabetic retinopathy in retinal images.

preprint2020arXiv

A Topology Layer for Machine Learning

Topology applied to real world data using persistent homology has started to find applications within machine learning, including deep learning. We present a differentiable topology layer that computes persistent homology based on level set filtrations and edge-based filtrations. We present three novel applications: the topological layer can (i) regularize data reconstruction or the weights of machine learning models, (ii) construct a loss on the output of a deep generative network to incorporate topological priors, and (iii) perform topological adversarial attacks on deep networks trained with persistence features. The code (www.github.com/bruel-gabrielsson/TopologyLayer) is publicly available and we hope its availability will facilitate the use of persistent homology in deep learning and other gradient based applications.

preprint2020arXiv

Homological Percolation: The Formation of Giant k-Cycles

In this paper we introduce and study a higher-dimensional analogue of the giant component in continuum percolation. Using the language of algebraic topology, we define the notion of giant k-dimensional cycles (with 0-cycles being connected components). Considering a continuum percolation model in the flat d-dimensional torus, we show that all the giant k-cycles (k=1,...,d-1) appear in the regime known as the thermodynamic limit. We also prove that the thresholds for the emergence of the giant k-cycles are increasing in k and are tightly related to the critical values in continuum percolation. Finally, we provide bounds for the exponential decay of the probabilities of giant cycles appearing.

preprint2020arXiv

Randomly Weighted $d-$complexes: Minimal Spanning Acycles and Persistence Diagrams

A weighted $d-$complex is a simplicial complex of dimension $d$ in which each face is assigned a real-valued weight. We derive three key results here concerning persistence diagrams and minimal spanning acycles (MSAs) of such complexes. First, we establish an equivalence between the MSA face-weights and \emph{death times} in the persistence diagram. Next, we show a novel stability result for the MSA face-weights which, due to our first result, also holds true for the death and birth times, separately. Our final result concerns a perturbation of a mean-field model of randomly weighted $d-$complexes. The $d-$face weights here are perturbation of some i.i.d. distribution while all the lower-dimensional faces have a weight of $0$. If the perturbations decay sufficiently quickly, we show that suitably scaled extremal nearest face-weights, face-weights of the $d-$MSA, and the associated death times converge to an inhomogeneous Poisson point process. This result completely characterizes the extremal points of persistence diagrams and MSAs. The point process convergence and the asymptotic equivalence of three point processes are new for any weighted random complex model, including even the non-perturbed case. Lastly, as a consequence of our stability result, we show that Frieze's $ζ(3)$ limit for random minimal spanning trees and the recent extension to random MSAs by Hino and Kanazawa also hold in suitable noisy settings.

preprint2019arXiv

Homological Percolation and the Euler Characteristic

In this paper we study the connection between the phenomenon of homological percolation (the formation of "giant" cycles in persistent homology), and the zeros of the expected Euler characteristic curve. We perform an experimental study that covers four different models: site-percolation on the cubical and permutahedral lattices, the Poisson-Boolean model, and Gaussian random fields. All the models are generated on the flat torus $T^d$, for $d=2,3,4$. The simulation results strongly indicate that the zeros of the expected Euler characteristic curve approximate the critical values for homological-percolation. Our results also provide some insight about the approximation error. Further study of this connection could have powerful implications both in the study of percolation theory, and in the field of Topological Data Analysis.

preprint2016arXiv

Maximally Persistent Cycles in Random Geometric Complexes

We initiate the study of persistent homology of random geometric simplicial complexes. Our main interest is in maximally persistent cycles of degree-$k$ in persistent homology, for a either the \cech or the Vietoris--Rips filtration built on a uniform Poisson process of intensity $n$ in the unit cube $[0,1]^d$. This is a natural way of measuring the largest "$k$-dimensional hole" in a random point set. This problem is in the intersection of geometric probability and algebraic topology, and is naturally motivated by a probabilistic view of topological inference. We show that for all $d \ge 2$ and $1 \le k \le d-1$ the maximally persistent cycle has (multiplicative) persistence of order $$ Θ\left(\left(\frac{\log n}{\log \log n} \right)^{1/k} \right),$$ with high probability, characterizing its rate of growth as $n \to \infty$. The implied constants depend on $k$, $d$, and on whether we consider the Vietoris--Rips or \cech filtration.

preprint2013arXiv

A Comparison of Relaxations of Multiset Cannonical Correlation Analysis and Applications

Canonical correlation analysis is a statistical technique that is used to find relations between two sets of variables. An important extension in pattern analysis is to consider more than two sets of variables. This problem can be expressed as a quadratically constrained quadratic program (QCQP), commonly referred to Multi-set Canonical Correlation Analysis (MCCA). This is a non-convex problem and so greedy algorithms converge to local optima without any guarantees on global optimality. In this paper, we show that despite being highly structured, finding the optimal solution is NP-Hard. This motivates our relaxation of the QCQP to a semidefinite program (SDP). The SDP is convex, can be solved reasonably efficiently and comes with both absolute and output-sensitive approximation quality. In addition to theoretical guarantees, we do an extensive comparison of the QCQP method and the SDP relaxation on a variety of synthetic and real world data. Finally, we present two useful extensions: we incorporate kernel methods and computing multiple sets of canonical vectors.

preprint2013arXiv

Approximating Local Homology from Samples

Recently, multi-scale notions of local homology (a variant of persistent homology) have been used to study the local structure of spaces around a given point from a point cloud sample. Current reconstruction guarantees rely on constructing embedded complexes which become difficult in high dimensions. We show that the persistence diagrams used for estimating local homology, can be approximated using families of Vietoris-Rips complexes, whose simple constructions are robust in any dimension. To the best of our knowledge, our results, for the first time, make applications based on local homology, such as stratification learning, feasible in high dimensions.

preprint2013arXiv

Persistence modules: Algebra and algorithms

Persistent homology was shown by Carlsson and Zomorodian to be homology of graded chain complexes with coefficients in the graded ring $\kk[t]$. As such, the behavior of persistence modules -- graded modules over $\kk[t]$ is an important part in the analysis and computation of persistent homology. In this paper we present a number of facts about persistence modules; ranging from the well-known but under-utilized to the reconstruction of techniques to work in a purely algebraic approach to persistent homology. In particular, the results we present give concrete algorithms to compute the persistent homology of a simplicial complex with torsion in the chain complex.

preprint2011arXiv

A spectral sequence for parallelized persistence

We approach the problem of the computation of persistent homology for large datasets by a divide-and-conquer strategy. Dividing the total space into separate but overlapping components, we are able to limit the total memory residency for any part of the computation, while not degrading the overall complexity much. Locally computed persistence information is then merged from the components and their intersections using a spectral sequence generalizing the Mayer-Vietoris long exact sequence. We describe the Mayer-Vietoris spectral sequence and give details on how to compute with it. This allows us to merge local homological data into the global persistent homology. Furthermore, we detail how the classical topology constructions inherent in the spectral sequence adapt to a persistence perspective, as well as describe the techniques from computational commutative algebra necessary for this extension. The resulting computational scheme suggests a parallelization scheme, and we discuss the communication steps involved in this scheme. Furthermore, the computational scheme can also serve as a guideline for which parts of the boundary matrix manipulation need to co-exist in primary memory at any given time allowing for stratified memory access in single-core computation. The spectral sequence viewpoint also provides easy proofs of a homology nerve lemma as well as a persistent homology nerve lemma. In addition, the algebraic tools we develop to approch persistent homology provide a purely algebraic formulation of kernel, image and cokernel persistence (D. Cohen-Steiner, H. Edelsbrunner, J. Harer, and D. Morozov. Persistent homology for kernels, images, and cokernels. In Proceedings of the twentieth Annual ACM-SIAM Symposium on Discrete Algorithms, pages 1011-1020. Society for Industrial and Applied Mathematics, 2009.)