Turnstile Streaming Algorithms Might As Well Be Linear Sketches

Total Page:16

File Type:pdf, Size:1020Kb

Turnstile Streaming Algorithms Might As Well Be Linear Sketches Turnstile Streaming Algorithms Might as Well Be Linear Sketches Yi Li Huy L. Nguyên˜ David P. Woodruff Max-Planck Institute for Princeton University IBM Research, Almaden Informatics [email protected] [email protected] [email protected] ABSTRACT Keywords In the turnstile model of data streams, an underlying vector x 2 streaming algorithms, linear sketches, lower bounds, communica- {−m; −m+1; : : : ; m−1; mgn is presented as a long sequence of tion complexity positive and negative integer updates to its coordinates. A random- ized algorithm seeks to approximate a function f(x) with constant probability while only making a single pass over this sequence of 1. INTRODUCTION updates and using a small amount of space. All known algorithms In the turnstile model of data streams [28, 34], there is an under- in this model are linear sketches: they sample a matrix A from lying n-dimensional vector x which is initialized to ~0 and evolves a distribution on integer matrices in the preprocessing phase, and through an arbitrary finite sequence of additive updates to its coor- maintain the linear sketch A · x while processing the stream. At the dinates. These updates are fed into a streaming algorithm, and have end of the stream, they output an arbitrary function of A · x. One the form x x + ei or x x − ei, where ei is the i-th standard cannot help but ask: are linear sketches universal? unit vector. This changes the i-th coordinate of x by an additive In this work we answer this question by showing that any 1- 1 or −1. At the end of the stream, x is guaranteed to satisfy the pass constant probability streaming algorithm for approximating promise that x 2 {−m; −m + 1; : : : ; mgn, where throughout we an arbitrary function f of x in the turnstile model can also be assume that m ≥ 2n. The goal of the streaming algorithm is to implemented by sampling a matrix A from the uniform distribu- make one pass over the data and to use limited memory to compute tion on O(n log m) integer matrices, with entries of magnitude functions of x, such as the frequency moments [1], the number of poly(n), and maintaining the linear sketch Ax. Furthermore, the distinct elements [19], the empirical entropy [12], the heavy hitters logarithm of the number of possible states of Ax, as x ranges over [14, 17], and treating x as a matrix, various quantities in numerical {−m; −m + 1; : : : ; mgn, plus the amount of randomness needed linear algebra such as a low rank approximation [15]. Since com- to store A, is at most a logarithmic factor larger than the space re- puting these quantities exactly or deterministically often requires a quired of the space-optimal algorithm. Our result shows that to prohibitive amount of space, these algorithms are usually random- prove space lower bounds for 1-pass streaming algorithms, it suf- ized and approximate. fices to prove lower bounds in the simultaneous model of commu- Curiously, all known algorithms in the turnstile model have the nication complexity, rather than the stronger 1-way model. More- following form: they sample an r × n matrix A from a distribution over, the fact that we can assume we have a linear sketch with on integer matrices in a preprocessing phase, and then maintain the polynomially-bounded entries further simplifies existing lower bounds, “linear sketch” A · x during the stream1. For known solutions to e.g., for frequency moments we present a simpler proof of the the above problems, the sketch A·x is maintained over the integers Ω(~ n1−2=k) bit complexity lower bound without using communi- (rather, than say, a finite field). From A · x, the algorithm then cation complexity. computes an arbitrary function of A·x, which should with constant probability, be a good approximation to the desired function f(x). Categories and Subject Descriptors Linear sketches are well-suited for turnstile streaming algorithms since given the state A · x and an additive update x x + ei, the F.2.2 [Analysis of Algorithms and Problem Complexity]: Non- new state A · (x + ei) can be computed as A · x + Ai, where Ai numerical Algorithms and Problems is the i-th column of A. This raises the natural question of whether or not all turnstile streaming algorithms need to in fact be linear General Terms sketches. Several works already consider the question of lower bounds Theory on the dimension of linear sketches themselves, rather than try- Permission to make digital or hard copies of all or part of this work for ing to obtain bit complexity lower bounds [4, 27, 32, 36]. For in- personal or classroom use is granted without fee provided that copies are not stance, Andoni et al. [4] prove a lower bound on the dimension of made or distributed for profit or commercial advantage and that copies bear a randomized linear sketch for estimating frequency moments and this notice and the full citation on the first page. Copyrights for components state “We stress that essentially all known algorithms in the general of this work owned by others than ACM must be honored. Abstracting with streaming model are in fact linear estimators.” The relationship be- credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request tween the space complexity of the optimal turnstile algorithm and permissions from [email protected]. 1 STOC ’14, May 31 - June 03 2014, New York, NY, USA There are algorithms in the insertion-only model, where all the Copyright 2014 ACM 978-1-4503-2710-7/14/05 ...$15.00. updates are required to be positive, which are not linear sketches, http://dx.doi.org/10.1145/2591796.2591812 e.g., the algorithm for frequency moments given by Alon et al. [1]. that of a linear sketch is thus important for understanding the appli- and lower bounds the communication required among the players 2 cability of lower bounds for linear sketches . to compute a function f(X1;:::;Xk). For 1-pass lower bounds, it Some progress has been made on the above question. Ganguly suffices to consider 1-way communication, in which P1 speaks to [22] shows that for approximating the `1-heavy hitters, any deter- P2 who speaks to P3, etc., and Pk outputs the answer. To obtain ministic streaming algorithm might as well be a linear sketch over lower bounds in the turnstile model, each player Pi creates a data the integers. Ganguly also generalized this reduction to determin- stream σi from his/her input Xi. Then P1 runs a data stream al- istic streaming algorithms for approximating the `p-heavy hitters, gorithm A on σ1, and transmits the memory contents of A to P2, for every 1 ≤ p ≤ 2 [23]. who continues the computation of A on σ2, etc. At the end, A will In a related work, Feldman et al. [18] introduced the MUD have been executed on the stream σ1 ◦ σ2 ◦ · · · ◦ σk. If the output model for computational tasks, which stands for massive, unordered, of A determines f with constant probability, then the space of the distributed algorithms, which contains linear sketches as a spe- streaming algorithm is at least the randomized 1-way communica- cial case. The authors show that any deterministic streaming al- tion complexity of f divided by k. gorithm for computing a symmetric (order-invariant) total single- While for some problems randomized 1-way communication lower valued function exactly can be simulated by a MUD algorithm. bounds are easy via a reduction from the Indexing problem [31], They partially extend this to approximation algorithms, but require for other problems such as k-player Set Disjointness and Gap-`1 that the approximating function is also a symmetric function of its [7, 13, 26, 30], 1-way lower bounds are not much easier to prove input, rather than only requiring that the function being approxi- than 2-way communication lower bounds with current techniques. mated is symmetric. For randomized algorithms they require that A weaker model than the 1-way model is the simultaneous commu- for each fixed random string, the streaming algorithm computes a nication model. In this model there are k players P1;:::;Pk, each symmetric function of its input. with an input Xi, but the communication is even more restricted. Thus, existing results on this problem were either for determin- There is an additional player, called the referee, and each of the istic algorithms and specific problems, or did not work for arbitrary players Pi can only send a single message, and only to the referee, approximation algorithms. though they may share a common random string. The referee an- nounces the output, which should equal f(X1;:::;Xk) with con- Our Contributions: We show that any 1-pass constant probability stant probability. Before our work, one could not use randomized streaming algorithm for approximating an arbitrary function f of x communication lower bounds for f to obtain space lower bounds in the turnstile model can be implemented by sampling a matrix A for streaming algorithms since there is no simulation in which the from the uniform distribution with support on O(n log m) integer state of a streaming algorithm can be passed sequentially among the matrices, each with entries of magnitude poly(n log m), and main- players. However, since our result shows that the optimal stream- taining the linear sketch A · x over the integers. Furthermore, the ing algorithm can be implemented using a sketching matrix A, up logarithm of the number of possible states of A·x, as x ranges over to a logarithmic factor, each player Pi can compute A · xi, where n {−m; −m + 1; : : : ; mg , plus the number of random bits needed xi is the underlying vector associated with the stream σi that Pi to sample A, is at most an O(log m) factor larger than the space creates from his/her input Xi to the communication game.
Recommended publications
  • Data Streams and Applications in Computer Science
    Data Streams and Applications in Computer Science David P. Woodruff IBM Research Almaden [email protected] Abstract This is a short survey of my work in data streams and related appli- cations, such as communication complexity, numerical linear algebra, and sparse recovery. The goal is give a non-technical overview of results in data streams, and highlight connections between these different areas. It is based on my Presburger lecture award given at ICALP, 2014. 1 The Data Stream Model Informally speaking, a data stream is a sequence of data that is too large to be stored in available memory. The data may be a sequence of numbers, points, edges in a graph, and so on. There are many examples of data streams, such as internet search logs, network traffic, sensor networks, and scientific data streams (such as in astronomics, genomics, physical simulations, etc.). The abundance of data streams has led to new algorithmic paradigms for processing them, which often impose very stringent requirements on the algorithm’s resources. Formally, in the streaming model, there is a sequence of elements a1;:::; am presented to an algorithm, where each element is drawn from a universe [n] = f1;:::; ng. The algorithm is allowed a single or a small number of passes over the stream. In network applications, the algorithm is typically only given a single pass, since if data on a network is not physically stored somewhere, it may be impossible to make a second pass over it. In other applications, such as when data resides on external memory, it may be streamed through main memory a small number of times, each time constituting a pass of the algorithm.
    [Show full text]
  • Efficient Semi-Streaming Algorithms for Local Triangle Counting In
    Efficient Semi-streaming Algorithms for Local Triangle Counting in Massive Graphs ∗ Luca Becchetti Paolo Boldi Carlos Castillo “Sapienza” Università di Roma Università degli Studi di Milano Yahoo! Research Rome, Italy Milan, Italy Barcelona, Spain [email protected] [email protected] [email protected] Aristides Gionis Yahoo! Research Barcelona, Spain [email protected] ABSTRACT Categories and Subject Descriptors In this paper we study the problem of local triangle count- H.3.3 [Information Systems]: Information Search and Re- ing in large graphs. Namely, given a large graph G = (V; E) trieval we want to estimate as accurately as possible the number of triangles incident to every node v 2 V in the graph. The problem of computing the global number of triangles in a General Terms graph has been considered before, but to our knowledge this Algorithms, Measurements is the first paper that addresses the problem of local tri- angle counting with a focus on the efficiency issues arising Keywords in massive graphs. The distribution of the local number of triangles and the related local clustering coefficient can Graph Mining, Semi-Streaming, Probabilistic Algorithms be used in many interesting applications. For example, we show that the measures we compute can help to detect the 1. INTRODUCTION presence of spamming activity in large-scale Web graphs, as Graphs are a ubiquitous data representation that is used well as to provide useful features to assess content quality to model complex relations in a wide variety of applica- in social networks. tions, including biochemistry, neurobiology, ecology, social For computing the local number of triangles we propose sciences, and information systems.
    [Show full text]
  • Streaming Quantiles Algorithms with Small Space and Update Time
    Streaming Quantiles Algorithms with Small Space and Update Time Nikita Ivkin Edo Liberty Kevin Lang [email protected] [email protected] [email protected] Amazon Amazon Yahoo Research Zohar Karnin Vladimir Braverman [email protected] [email protected] Amazon Johns Hopkins University ABSTRACT An approximate version of the problem relaxes this requirement. Approximating quantiles and distributions over streaming data has It is allowed to output an approximate rank which is off by at most been studied for roughly two decades now. Recently, Karnin, Lang, εn from the exact rank. In a randomized setting, the algorithm and Liberty proposed the first asymptotically optimal algorithm is allowed to fail with probability at most δ. Note that, for the for doing so. This manuscript complements their theoretical result randomized version to provide a correct answer to all possible by providing a practical variants of their algorithm with improved queries, it suffices to amplify the success probability by running constants. For a given sketch size, our techniques provably reduce the algorithm with a failure probability of δε, and applying the 1 the upper bound on the sketch error by a factor of two. These union bound over O¹ ε º quantiles. Uniform random sampling of ¹ 1 1 º improvements are verified experimentally. Our modified quantile O ε 2 log δ ε solves this problem. sketch improves the latency as well by reducing the worst-case In network monitoring [16] and other applications it is critical 1 1 to maintain statistics while making only a single pass over the update time from O¹ ε º down to O¹log ε º.
    [Show full text]
  • Finding Graph Matchings in Data Streams
    Finding Graph Matchings in Data Streams Andrew McGregor⋆ Department of Computer and Information Science, University of Pennsylvania, Philadelphia, PA 19104, USA [email protected] Abstract. We present algorithms for finding large graph matchings in the streaming model. In this model, applicable when dealing with mas- sive graphs, edges are streamed-in in some arbitrary order rather than 1 residing in randomly accessible memory. For ǫ > 0, we achieve a 1+ǫ 1 approximation for maximum cardinality matching and a 2+ǫ approxi- mation to maximum weighted matching. Both algorithms use a constant number of passes and O˜( V ) space. | | 1 Introduction Given a graph G = (V, E), the Maximum Cardinality Matching (MCM) problem is to find the largest set of edges such that no two adjacent edges are selected. More generally, for an edge-weighted graph, the Maximum Weighted Matching (MWM) problem is to find the set of edges whose total weight is maximized subject to the condition that no two adjacent edges are selected. Both problems are well studied and exact polynomial solutions are known [1–4]. The fastest of these algorithms solves MWM with running time O(nm + n2 log n) where n = V and m = E . | | | | However, for massive graphs in real world applications, the above algorithms can still be prohibitively computationally expensive. Examples include the vir- tual screening of protein databases. (See [5] for other examples.) Consequently there has been much interest in faster algorithms, typically of O(m) complexity, that find good approximate solutions to the above problems. For MCM, a linear time approximation scheme was given by Kalantari and Shokoufandeh [6].
    [Show full text]
  • Lecture 1 1 Graph Streaming
    CS 671: Graph Streaming Algorithms and Lower Bounds Rutgers: Fall 2020 Lecture 1 September 1, 2020 Instructor: Sepehr Assadi Scribe: Sepehr Assadi Disclaimer: These notes have not been subjected to the usual scrutiny reserved for formal publications. They may be distributed outside this class only with the permission of the Instructor. 1 Graph Streaming Motivation Massive graphs appear in most application domains nowadays: web-pages and hyperlinks, neurons and synapses, papers and citations, or social networks and friendship links are just a few examples. Many of the computing tasks in processing these graphs, at their core, correspond to classical problems such as reachability, shortest path, or maximum matching|problems that have been studied extensively for almost a century now. However, the traditional algorithmic approaches to these problems do not accurately capture the challenges of processing massive graphs. For instance, we can no longer assume a random access to the entire input on these graphs, as their sheer size prevents us from loading them into the main memory. To handle these challenges, there is a rapidly growing interest in developing algorithms that explicitly account for the restrictions of processing massive graphs. Graph streaming algorithms have been particularly successful in this regard; these are algorithms that process input graphs by making one or few sequential passes over their edges while using a small memory. Not only streaming algorithms can be directly used to address various challenges in processing massive graphs such as I/O-efficiency or monitoring evolving graphs in realtime, they are also closely connected to other models of computation for massive graphs and often lead to efficient algorithms in those models as well.
    [Show full text]
  • Streaming Algorithms
    3 Streaming Algorithms Great Ideas in Theoretical Computer Science Saarland University, Summer 2014 Some Admin: • Deadline of Problem Set 1 is 23:59, May 14 (today)! • Students are divided into two groups for tutorials. Visit our course web- site and use Doodle to to choose your free slots. • If you find some typos or like to improve the lecture notes (e.g. add more details, intuitions), send me an email to get .tex file :-) We live in an era of “Big Data”: science, engineering, and technology are producing in- creasingly large data streams every second. Looking at social networks, electronic commerce, astronomical and biological data, the big data phenomenon appears everywhere nowadays, and presents opportunities and challenges for computer scientists and mathematicians. Com- paring with traditional algorithms, several issues need to be considered: A massive data set is too big to be stored; even an O(n2)-time algorithm is too slow; data may change over time, and algorithms need to cope with dynamic changes of the data. Hence streaming, dynamic and distributed algorithms are needed for analyzing big data. To illuminate the study of streaming algorithms, let us look at a few examples. Our first example comes from network analysis. Any big website nowadays receives millions of visits every minute, and very frequent visit per second from the same IP address are probably gen- erated automatically by computers and due to hackers. In this scenario, a webmaster needs to efficiently detect and block those IPs which generate very frequent visits per second. Here, instead of storing all IP record history which are expensive, we need an efficient algorithm that detect most such IPs.
    [Show full text]
  • Anytime Algorithms for Stream Data Mining
    Anytime Algorithms for Stream Data Mining Von der Fakultat¨ fur¨ Mathematik, Informatik und Naturwissenschaften der RWTH Aachen University zur Erlangung des akademischen Grades eines Doktors der Naturwissenschaften genehmigte Dissertation vorgelegt von Diplom-Informatiker Philipp Kranen aus Willich, Deutschland Berichter: Universitatsprofessor¨ Dr. rer. nat. Thomas Seidl Visiting Professor Michael E. Houle, PhD Tag der mundlichen¨ Prufung:¨ 14.09.2011 Diese Dissertation ist auf den Internetseiten der Hochschulbibliothek online verfugbar.¨ Contents Abstract / Zusammenfassung1 I Introduction5 1 The Need for Anytime Algorithms7 1.1 Thesis structure......................... 16 2 Knowledge Discovery from Data 17 2.1 The KDD process and data mining tasks ........... 17 2.2 Classification .......................... 25 2.3 Clustering............................ 36 3 Stream Data Mining 43 3.1 General Tools and Techniques................. 43 3.2 Stream Classification...................... 52 3.3 Stream Clustering........................ 59 II Anytime Stream Classification 69 4 The Bayes Tree 71 4.1 Introduction and Preliminaries................. 72 4.2 Indexing density models.................... 76 4.3 Experiments........................... 87 4.4 Conclusion............................ 98 i ii CONTENTS 5 The MC-Tree 99 5.1 Combining Multiple Classes.................. 100 5.2 Experiments........................... 111 5.3 Conclusion............................ 116 6 Bulk Loading the Bayes Tree 117 6.1 Bulk loading mixture densities . 117 6.2 Experiments..........................
    [Show full text]
  • Regular Path Query Evaluation on Streaming Graphs
    Regular Path ery Evaluation on Streaming Graphs∗ Anil Pacaci Angela Bonifati M. Tamer Özsu University of Waterloo Lyon 1 University University of Waterloo [email protected] [email protected] [email protected] ABSTRACT on the entire stream. Most existing work focus on the snap- We study persistent query evaluation over streaming graphs, shot model, which assumes that graphs are static and fully which is becoming increasingly important. We focus on available, and adhoc queries reect the current state of the navigational queries that determine if there exists a path database (e.g., [24, 45, 62, 64, 68–70]). The dynamic graph between two entities that satises a user-specied constraint. model addresses the evolving nature of these graphs; how- We adopt the Regular Path Query (RPQ) model that species ever, algorithms in this model assume that the entire graph navigational patterns with labeled constraints. We propose is fully available and they compute how the output changes deterministic algorithms to eciently evaluate persistent as the graph is updated [15, 42, 47, 60]. RPQs under both arbitrary and simple path semantics in a In this paper, we study the problem of persistent query uniform manner. Experimental analysis on real and synthetic processing over streaming graphs, addressing the limitations streaming graphs shows that the proposed algorithms can of existing approaches. We adopt the Regular Path Query process up to tens of thousands of edges per second and (RPQ) model that focuses on path navigation, e.g., nding eciently answer RPQs that are commonly used in real- pairs of users in a network connected by a path whose label world workloads.
    [Show full text]
  • Stochastic Streams: Sample Complexity Vs. Space Complexity
    Stochastic Streams: Sample Complexity vs. Space Complexity Michael Crouch1, Andrew McGregor2, Gregory Valiant3, and David P. Woodruff4 1 Bell Labs Ireland 2 University of Massachusetts Amherst 3 Stanford University 4 IBM Research Almaden Abstract We address the trade-off between the computational resources needed to process a large data set and the number of samples available from the data set. Specifically, we consider the following abstraction: we receive a potentially infinite stream of IID samples from some unknown distribu- tion D, and are tasked with computing some function f(D). If the stream is observed for time t, how much memory, s, is required to estimate f(D)? We refer to t as the sample complexity and s as the space complexity. The main focus of this paper is investigating the trade-offs between the space and sample complexity. We study these trade-offs for several canonical problems stud- ied in the data stream model: estimating the collision probability, i.e., the second moment of a distribution, deciding if a graph is connected, and approximating the dimension of an unknown subspace. Our results are based on techniques for simulating different classical sampling proced- ures in this model, emulating random walks given a sequence of IID samples, as well as leveraging a characterization between communication bounded protocols and statistical query algorithms. 1998 ACM Subject Classification F.2.3 Tradeoffs between Complexity Measures Keywords and phrases data streams, sample complexity, frequency moments, graph connectivity, subspace approximation Digital Object Identifier 10.4230/LIPIcs.ESA.2016.258 1 Introduction Big data systems must process data at multiple levels of a hierarchy, starting at local per-process or per-container monitoring and extending up to aggregate statistics for entire databases or data centers.
    [Show full text]
  • Finding Matchings in the Streaming Model
    University of Pennsylvania ScholarlyCommons Technical Reports (CIS) Department of Computer & Information Science January 2004 Finding Matchings in the Streaming Model Andrew McGregor University of Pennsylvania Follow this and additional works at: https://repository.upenn.edu/cis_reports Recommended Citation Andrew McGregor, "Finding Matchings in the Streaming Model", . January 2004. University of Pennsylvania Department of Computer and Information Science Technical Report No. MS-CIS-04-19. This paper is posted at ScholarlyCommons. https://repository.upenn.edu/cis_reports/21 For more information, please contact [email protected]. Finding Matchings in the Streaming Model Abstract This report presents algorithms for finding large matchings in the streaming model. In this model, applicable when dealing with massive graphs, edges are streamed-in in some arbitrary order rather than residing in randomly accessible memory. For ε > 0, we achieve a 1/(1+ε) approximation for maximum cardinality matching and a 1/(2+ε) approximation to maximum weighted matching. Both algorithms use a constant number of passes. Comments University of Pennsylvania Department of Computer and Information Science Technical Report No. MS- CIS-04-19. This technical report is available at ScholarlyCommons: https://repository.upenn.edu/cis_reports/21 Finding Matchings in the Streaming Model Andrew McGregor∗ Abstract This report presents algorithms for finding large matchings in the streaming model. In this model, applicable when dealing with massive graphs, edges are streamed-in in some arbitrary order rather than 1 residing in randomly accessible memory. For > 0, we achieve a 1+ approximation for maximum 1 cardinality matching and a 2+ approximation to maximum weighted matching. Both algorithms use a constant number of passes.
    [Show full text]
  • Suffix Tree, Minwise Hashing and Streaming Algorithms for Big Data Analysis in Bioinformatics
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Computer Science and Engineering: Theses, Computer Science and Engineering, Department Dissertations, and Student Research of Fall 12-4-2020 SUFFIX TREE, MINWISE HASHING AND STREAMING ALGORITHMS FOR BIG DATA ANALYSIS IN BIOINFORMATICS Sairam Behera University of Nebraska-Lincoln, [email protected] Follow this and additional works at: https://digitalcommons.unl.edu/computerscidiss Part of the Computer Engineering Commons, and the Computer Sciences Commons Behera, Sairam, "SUFFIX TREE, MINWISE HASHING AND STREAMING ALGORITHMS FOR BIG DATA ANALYSIS IN BIOINFORMATICS" (2020). Computer Science and Engineering: Theses, Dissertations, and Student Research. 201. https://digitalcommons.unl.edu/computerscidiss/201 This Article is brought to you for free and open access by the Computer Science and Engineering, Department of at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Computer Science and Engineering: Theses, Dissertations, and Student Research by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. SUFFIX TREE, MINWISE HASHING AND STREAMING ALGORITHMS FOR BIG DATA ANALYSIS IN BIOINFORMATICS by Sairam Behera A DISSERTATION Presented to the Faculty of The Graduate College at the University of Nebraska In Partial Fulfilment of Requirements For the Degree of Doctor of Philosophy Major: Computer Science Under the Supervision of Professor Jitender S. Deogun Lincoln, Nebraska December, 2020 SUFFIX TREE, MINWISE HASHING AND STREAMING ALGORITHMS FOR BIG DATA ANALYSIS IN BIOINFORMATICS Sairam Behera, Ph.D. University of Nebraska, 2020 Adviser: Jitender S. Deogun In this dissertation, we worked on several algorithmic problems in bioinformatics using mainly three approaches: (a) a streaming model, (b) suffix-tree based indexing, and (c) minwise-hashing (minhash) and locality-sensitive hashing (LSH).
    [Show full text]
  • Streaming Algorithms for Halo Finders
    Streaming Algorithms for Halo Finders Zaoxing Liu∗ 1, Nikita Ivkin∗ 1, Lin F. Yangy 2, Mark Neyrincky 3, Gerard Lemsony 3, Alexander S. Szalayy 3, Vladimir Braverman∗ 4, Tamas Budavariy, Randal Burns∗, Xin Wangy 3 ∗Department of Computer Science yDepartment of Physics & Astronomy Johns Hopkins University Baltimore, MD 21218, USA Abstract—Cosmological N-body simulations are essential for I. INTRODUCTION studies of the large-scale distribution of matter and galaxies in the Universe. This analysis often involves finding clusters of The goal of astrophysics is to explain the observed particles and retrieving their properties. Detecting such “halos” properties of the universe we live in. In cosmology in among a very large set of particles is a computationally inten- sive problem, usually executed on the same super-computers particular, one tries to understand how matter is distributed that produced the simulations, requiring huge amounts of on the largest scales we can observe. In this effort, advanced memory. computer simulations play an ever more important role. Sim- Recently, a new area of computer science emerged. This area, ulations are currently the only way to accurately understand called streaming algorithms, provides new theoretical methods the nonlinear processes that produce cosmic structures such to compute data analytics in a scalable way using only a single as galaxies and patterns of galaxies. Hence a large amount of pass over a data sets and logarithmic memory. effort is spent on running simulations modelling representa- The main contribution of this paper is a novel connection tive parts of the universe in ever greater detail. A necessary between the N-body simulations and the streaming algorithms.
    [Show full text]