DeepAI AI Chat
Log In Sign Up

A fast semi-discrete optimal transport algorithm for a unique reconstruction of the early Universe

by   Bruno Lévy, et al.

We leverage powerful mathematical tools stemming from optimal transport theory and transform them into an efficient algorithm to reconstruct the fluctuations of the primordial density field, built on solving the Monge-Ampère-Kantorovich equation. Our algorithm computes the optimal transport between an initial uniform continuous density field, partitioned into Laguerre cells, and a final input set of discrete point masses, linking the early to the late Universe. While existing early universe reconstruction algorithms based on fully discrete combinatorial methods are limited to a few hundred thousand points, our algorithm scales up well beyond this limit, since it takes the form of a well-posed smooth convex optimization problem, solved using a Newton method. We run our algorithm on cosmological N-body simulations, from the AbacusCosmos suite, and reconstruct the initial positions of 𝒪(10^7) particles within a few hours with an off-the-shelf personal computer. We show that our method allows a unique, fast and precise recovery of subtle features of the initial power spectrum, such as the baryonic acoustic oscillations.


page 7

page 9

page 11

page 17

page 19

page 20


Accurate Baryon Acoustic Oscillations reconstruction via semi-discrete optimal transport

Optimal transport theory has recently reemerged as a vastly resourceful ...

Partial Optimal Transport for a Constant-Volume Lagrangian Mesh with Free Boundaries

This article introduces a representation of dynamic meshes, adapted to s...

Towards a mathematical theory of trajectory inference

We devise a theoretical framework and a numerical method to infer trajec...

Higher-dimensional power diagrams for semi-discrete optimal transport

Efficient algorithms for solving optimal transport problems are importan...

Dynamic Cell Imaging in PET with Optimal Transport Regularization

We propose a novel dynamic image reconstruction method from PET listmode...

Pictures of Combinatorial Cubes

We prove that the 8-point algorithm always fails to reconstruct a unique...

Computational Semi-Discrete Optimal Transport with General Storage Fees

We propose and analyze a modified damped Newton algorithm to solve the s...