PDOT: a Practical Primal-Dual Algorithm and a GPU-Based Solver for Optimal Transport

In this paper, we propose a practical primal-dual algorithm with theoretical guarantees and develop a GPU-based solver, which we dub PDOT, for solving large-scale optimal transport problems. Compared to Sinkhorn algorithm or classic LP algorithms, PDOT can achieve high-accuracy solution while efficiently taking advantage of modern computing architecture, i.e., GPUs. On the theoretical side, we show that PDOT has a data-independent widetilde O(mn(m+n)1.5log(frac1ε)) local flop complexity where ε is the desired accuracy, and m and n are the dimension of the original and target distribution, respectively. We further present a data-dependent widetilde O(mn(m+n)3.5Delta + mn(m+n)1.5log(frac1ε)) global flop complexity of PDOT, where Delta is the precision of the data. On the numerical side, we present PDOT, an open-source GPU solver based on the proposed algorithm. Our extensive numerical experiments consistently demonstrate the well balance of PDOT in computing efficiency and accuracy of the solution, compared to Gurobi and Sinkhorn algorithms.