No description available.

This dissertation seeks to describe advancements made in computational methods for multi-messenger astrophysics (MMA) using gravitational waves GW and neutrinos during Advanced LIGO (aLIGO)’s first through third observing runs (O1-O3) and, looking forward, to describe novel computational techniques suited to the challenges of both the burgeoning MMA field and high-performance computing as a whole. The first two chapters provide an overview of MMA as it pertains to gravitational wave/high energy neutrino (GWHEN) searches, including a summary of expected astrophysical sources as well as GW, neutrino, and gamma-ray detectors used in their detection. These are followed in the third chapter by an in-depth discussion of LIGO’s timing system, particularly the diagnostic subsystem, describing both its role in MMA searches and the author’s contributions to the system itself. The fourth chapter provides a detailed description of the Low-Latency Algorithm for Multi-messenger Astrophysics (LLAMA), the GWHEN pipeline developed by the author and used in O2 and O3. Relevant past multi-messenger searches are described first, followed by the O2 and O3 analysis methods, the pipeline’s performance, scientific results, and finally, an in-depth account of the library’s structure and functionality. In particular, the author’s high-performance multi-order coordinates (MOC) HEALPix image analysis library, HPMOC, is described. HPMOC increases performance of HEALPix image manipulations by several orders of magnitude vs. naive single-resolution approaches while presenting a simple high-level interface and should prove useful for diverse future MMA searches. The performance improvements it provides for LLAMA are also covered. The final chapter of this dissertation builds on the approaches taken in developing HPMOC, presenting several novel methods for efficiently storing and analyzing large data sets, with applications to MMA and other data-intensive fields. A family of depth-first multi-resolution ordering of HEALPix images — DEPTH9, DEPTH19, and DEPTH40 — is defined, along with algorithms and use cases where it can improve on current approaches, including high-speed streaming calculations suitable for serverless compute or FPGAs. For performance-constrained analyses on HEALPix data (e.g. image analysis in multi-messenger search pipelines) using SIMD processors, breadth-first data structures can provide short-circuiting calculations in a data-parallel way on compressed data; a simple compression method is described with application to further improving LLAMA performance. A new storage scheme and associated algorithms for efficiently compressing and contracting tensors of varying sparsity is presented; these demuxed tensors (D-Tensors) have equivalent asymptotic time and space complexity to optimal representations of both dense and sparse matrices, and could be used as a universal drop-in replacement to reduce code complexity and developer effort while improving performance of existing non-optimized numerical code. Finally, the big bucket hash table (B-Table), a novel type of hash table making guarantees on data layout (vs. load factor), is described, along with optimizations it allows for (like hardware acceleration, online rebuilds, and hard realtime applications) that are not possible with existing hash table approaches. These innovations are presented in the hope that some will prove useful for improving future MMA searches and other data-intensive applications.

Includes bibliographical references (p. 259-261). / This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections. / Thesis (Ph. D.)--Joint Program in Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Mechanical Engineering; and the Woods Hole Oceanographic Institution), 2007. / In this thesis, a numerically effcient three-dimensional propagation and scattering model is developed based on the three-dimensional coupled mode theory for axisymmetric bathymetry. The three-dimensional coupled mode approach applied in this thesis is fundamentally identical to the one applied in earlier models, such as the one presented by Taroudakis [20]. Thus, it is based on a Fourier expansion of the acoustic field around a seamount, with each azimuthal expansion coefficient being represented by a two-way coupled mode formulation. However, earlier formulations were severely limited in terms of frequency, size and geometry of the seamount, the seabed composition, and the distance between the source and the seamount, and are totally inadequate for modeling high-frequency, large-scale seamount problems. By introducing a number of changes in the numerical formulation and using a standard normal mode model (C-SNAP) for determining the fundamental modal solutions and coupling coefficients, orders of magnitude improvement in efficiency and fidelity has been achieved, allowing for realistic propagation and scattering scenarios to be modeled, including effects of seamount roughness and realistic sedimentary structure. / (cont.) Also, by the simple superposition principle, the computational requirements are made independent of the distance between the seamount and the source and receivers, and dependent only on the geometry of the seamount and the frequency of the source. First, this thesis investigates the scattering from a cylindrical island, which is the simplest case of a conical seamount problem. Our model, using the superposition method, can solve the cylindrical problem in Athanassoulis and Prospathopoulos's paper [3] with the same accuracy while saving about 4/5 computational effort. Second, this thesis demonstrates the spectral coupled mode approach, which includes a two-way coupled mode model and a superposition representation of the field. Third, this thesis applies the three-dimensional model to investigate some physics issues of three-dimensional seamount scattering. As a result of the investigation, we learn that the Nx2D model is a poor approximation of the true three-dimensional model when the three-dimensional effects are significant, though it is a good approximation of the three-dimensional model otherwise. The convergence of the model in terms of the seamount discretization is also discussed and demonstrated. / (cont.) Finally, our three-dimensional spectral coupled mode model is tested by the application of the Kermit Seamount problem with realistic ocean environmental data from the 2004 BASSEX experiment. / by Wenyu Luo. / Ph.D.