Papers › Fast Log-Domain Sinkhorn Optimal Transport with Warp-Level GPU Reductions

Fast Log-Domain Sinkhorn Optimal Transport with Warp-Level GPU Reductions

4 Apr 2026arXiv:2605.00837added by Syntology

Hao Xiao

Title, abstract, authors and date from arXiv's metadata (CC0); this paper is not in the Papers with Code archive (frozen 2025-07-28).

Entropic regularized optimal transport (OT) via the Sinkhorn algorithm has become a fundamental tool in machine learning, yet existing implementations either suffer from numerical instability for small regularization parameters or incur significant overhead from deep learning frameworks. We present FastSinkhorn, a lightweight, native CUDA implementation of the log-domain Sinkhorn algorithm that combines warp-level shuffle reductions with shared-memory tiling to achieve high GPU utilization without sacrificing numerical stability. Our solver operates entirely in the log-domain, enabling robust computation for regularization parameters as small as epsilon = 10^{-4} where standard-domain methods fail. On dense OT problems with n = m = 8192, our implementation achieves 12x speedup over the widely-used POT library and 5.9x speedup over GPU-accelerated PyTorch baselines, while consuming only 256 MB of GPU memory. We validate our solver on image color transfer, 3D point cloud matching, and convergence analysis, demonstrating that native CUDA kernels with careful numerical treatment provide a practical and efficient foundation for large-scale optimal transport computation.

PaperPDFCode Syntology ran

In Syntology View this paper on Syntology: its repositories, every harvested function with whether it ran, its licence and the call to fetch it.

Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2605.00837")

Code

Syntology Ran 12 of 12 code samples harvested from 1 repository linked to this paper; 0 have no recorded run. Of those that ran: 12 ran with no contract checked.

By repository: found in paper text by Syntology: 12 samples from 1 repository, 12 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

xiao98/Fast-Sinkhorn-CUDA found in paper text by Syntology report

Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

12 samples harvested; 12 ran; 0 honoured the contract we drafted; 0 have no recorded run. Read from Syntology's graph 2026-09-25; that is when this build read the record, not when the samples ran.

12ran

Licence: 0 of the 12 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from xiao98/Fast-Sinkhorn-CUDA. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

compute_ot_matching xiao98/Fast-Sinkhorn-CUDA/experiments/applications/point_cloud_matching.py found in paper text by Syntology ran fingerprinted MIT (permissive) · f785f92416c75d1d · report
compute_transport_cost xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/bench_pytorch_sinkhorn.py found in paper text by Syntology ran MIT (permissive) · 59b5efe0d457c5fa · report
generate_bunny_like_points xiao98/Fast-Sinkhorn-CUDA/experiments/applications/point_cloud_matching.py found in paper text by Syntology ran fingerprinted MIT (permissive) · f9768a4f1c357138 · report
generate_synthetic_image xiao98/Fast-Sinkhorn-CUDA/experiments/applications/color_transfer.py found in paper text by Syntology ran MIT (permissive) · 95368d0df733a6f4 · report
load_data xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/plot_baselines.py found in paper text by Syntology ran fingerprinted MIT (permissive) · 6829ee1542b9660c · report
make_distributions xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/bench_pytorch_sinkhorn.py found in paper text by Syntology ran MIT (permissive) · ac42230eccda1f28 · report
make_distributions xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/bench_geomloss.py found in paper text by Syntology ran MIT (permissive) · c1fab6a120849fc3 · report
ot_color_transfer xiao98/Fast-Sinkhorn-CUDA/experiments/applications/color_transfer.py found in paper text by Syntology ran MIT (permissive) · cc03a4060e69a92d · report
rotation_matrix xiao98/Fast-Sinkhorn-CUDA/experiments/applications/point_cloud_matching.py found in paper text by Syntology ran MIT (permissive) · 9b6a364baf249a7b · report
run_sinkhorn xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/bench_pot.py found in paper text by Syntology ran MIT (permissive) · ee2f8afaec88b236 · report
sample_pixels xiao98/Fast-Sinkhorn-CUDA/experiments/applications/color_transfer.py found in paper text by Syntology ran MIT (permissive) · ef6edc9926cd2d12 · report
sinkhorn_log_pytorch xiao98/Fast-Sinkhorn-CUDA/experiments/baselines/bench_pytorch_sinkhorn.py found in paper text by Syntology ran MIT (permissive) · f44585a13a67239f · report

Results from the paper

The Papers with Code archive ends with its 2025-07-28 snapshot. This paper's arXiv identifier, 2605.00837, was issued in May 2026, after that date, so the archive has no leaderboard rows for it.

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections