The package OptimalTransport.jl implements the sinkhorn method.
-
Implement comparisons to Sinkhorn. See the Pluto notebook test/optimal-transport.jl for a start. This implements a small test problem shown in the documentation of OptimalTransport.jl.
-
Note from the screenshot below that the @btime results in the Pluto notebook: KLLS has significant overhead relative to sinkhorn. (See screenshot.) Do you think you can spot where are the majority of the allocations? Is it KLLS or Trunk or ...?
-
Additional test problems. See Section 6 of Blondel et al 2018.

The package
OptimalTransport.jlimplements thesinkhornmethod.Implement comparisons to Sinkhorn. See the Pluto notebook
test/optimal-transport.jlfor a start. This implements a small test problem shown in the documentation ofOptimalTransport.jl.Note from the screenshot below that the
@btimeresults in the Pluto notebook: KLLS has significant overhead relative to sinkhorn. (See screenshot.) Do you think you can spot where are the majority of the allocations? Is it KLLS or Trunk or ...?Additional test problems. See Section 6 of Blondel et al 2018.