Large scale numerical experiments are commonplace today in theoretical physics. The high performance algorithms described herein are the most compact, efficient methods known for representing and analyzing systems modeled well by sets or graphs. After studying how these implementations maximize instruction throughput and optimize memory access patterns, we apply them to causal set quantum gravity, in which spacetime is represented by a partially ordered set. We build upon the low-level set and graph algorithms to optimize the calculation of the causal set action, and then discuss how to measure boundaries of a discrete spacetime. We then examine the broader applicability of these algorithms to greedy information routing in random geometric graphs embedded in Lorentzian manifolds, which requires us to find new closed-form solutions to the geodesic differential equations in Friedmann-Lemaitre-Robertson-Walker spacetimes. Finally, we consider the vacuum selection problem in string theory, where we show a network-centered approach yields a dynamical mechanism for vacuum selection in the context of multiverse cosmology. These algorithms have broad applicability to many physical systems, and they improve existing methods by reducing simulation runtimes by orders of magnitude.

Geodesics are used in a wide array of applications in cosmology and astrophysics. However, it is not a trivial task to efficiently calculate exact geodesic distances in an arbitrary spacetime. We show that in spatially flat (3+1)-dimensional Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetimes, it is possible to integrate the second-order geodesic differential equations, and derive a general method for finding both timelike and spacelike distances given initial-value or boundary-value constraints. In flat spacetimes with either dark energy or matter, whether dust, radiation, or a stiff fluid, we find an exact closed-form solution for geodesic distances. In spacetimes with a mixture of dark energy and matter, including spacetimes used to model our physical universe, there exists no closed-form solution, but we provide a fast numerical method to compute geodesics. A general method is also described for determining the geodesic connectedness of an FLRW manifold, provided only its scale factor.