__ Summary and Contributions__:
This paper proposes a framework for online computation of entropy-regularized optimal transport. Compared to current methods that operate on full or mini-batches, this method allows for streaming computation. On the other hand, in contrast with most current methods that operate exclusively on discrete distributions, this one can be used on continuous distributions, for which it relies on a non-parametric representation of the OT dual potentials. The paper proposes various versions of this “online Sinkhorn” scheme, provides convergence results under some assumptions, and presents empirical results on two settings: continuous potential estimation and warm-starting batch Sinkhorn. The results show this is a promising method for online OT computation.
--- POST-REBUTTAL --
I thank the authors for their response to my questions/concerns. I do hope they take into account to suggestion to include high-dimensional experiments (even if via embedding), to make the results even more compelling.

__ Strengths__: - Novelty/Originality. Streaming / online computation of optimal transport distances has received much less attention that its batch-computation counterpart, so work which is makes advances in this front is likely to have a strong impact in the community.
- Soundness. The proposed method is sound, grounded on solid theory and skillfully leverages ideas from Stochastic Approximation for handling potentials as continuous functions.
- Broad/thorough experimental evaluation. Except for the caveat on dimensionality (see weaknesses) The paper does a very good job of validating and comparing the various versions of the proposed approach on synthetic tasks.

__ Weaknesses__: - Motivation. While I think online OT estimation is an interesting challenge from a theoretical standpoint, I feel like the paper is missing a practical/concrete motivation for why/when online OT computation is needed or preferred over batch computation, e.g., perhaps emphasizing the “progressive enrichment” aspect of this framework
- No experiments on high-dimensional settings, arguably those for which a streaming setting would be most compelling. Related to the point above, a “real” application that truly benefits from streaming computation of OT would definitely strengthen the paper.

__ Correctness__: The main claims and methods presented in the paper are overall correct. A few smaller issues:
- I’m not entirely convinced that using “number of computations” in the x-axis makes sense for Figs 1, 3 etc. Batch-methods that involve some initial pre-computation (e.g., distance matrix computation for batch Sinkhorn), will of course make little to no progress early on. I’m not sure what the the “right” comparison would be, but I don’t think it’s this one - It’s comparing apples to oranges. At any rate, it would be useful to higlight in the plots what part of the computation corresponds to these “no-progress” steps.
- In L62 it is claimed that the memory complexity increases lineraly on n_t. Should this be O(n_t^2)?
- I could not find a discussion or details on how the learning rate \eta_n is chosen in practice, not sure if I missed it. This is particularly relevant considering Assumption 2.

__ Clarity__:
The paper is mostly well written, is relatively easy to follow and understand. Some of the proofs in the Appendix are a bit dry and could use more hand-holding.

__ Relation to Prior Work__:
Some related work is missing and/or could be discussed in more detail:
* [39] is one of the few recent OT methods for the streaming setting (albeit focused on barycenters). While cited here, it is not described in detailed nor compared against. It would be useful for the reader to have such a comparison in order to better set this work in context.
* There is a rich literature on streaming algorithms for the (classic, non-regularized) EMD problem (e.g. [1,2]), which is quite relevant but completely missing from the related work section here.
[1] Andoni et al., “Efficient Sketches for Earth-Mover Distance, with Applications“, 2009
[2] Indyk, “Efficient Sketches for Earth-Mover Distance, with Applications”, 2004

__ Reproducibility__: Yes

__ Additional Feedback__: Typos/suggestions:
- The thickness of the lines in the the figures makes it hard to compare them. Consider using thinner lines or transparency.
- Inconsistent big-O notation in L51.

__ Summary and Contributions__: The paper introduces an algorithm for computing optimal transport (OT) distances between two distributions, given access only to streams of samples from the distributions. (this is in contrast to the usual setup, where you take fixed samples from the distributions, then compute OT between those empirical distributions).
The authors prove appropriate convergence, sample complexity, and computational complexity results for the algorithm (the latter of which are comparable to Sinkhorn). They also provide compelling experiments showing that the algorithm yields better estimates using the same number of samples, and that it can even effectively warm-start Sinkhorn.

__ Strengths__: The algorithm in this paper seems convincingly better than all alternatives I am aware of for OT between general (e.g. continuous) distributions.
The authors back the algorithm with theory that is sensible and reasonably well-explained. The experiments have reasonable breadth and are presented in a clear way (I particularly like the plots in Figure 2 comparing the potentials to those from the baseline from [19]). The authors are also reasonably up-front about the limitations of the algorithm (e.g. that you need to either increase the batch size, or use the fully-corrective scheme in section 3.2).
Faster/better OT solvers are highly relevant to NeurIPS, as OT is a tool with wide applicability in ML. As for the significant of this specific algorithm, though it is not likely to displace Sinkhorn e.g. in OT loss functions in deep learning, it is a nice treat that Sinkhorn can be cheaply warm-started via this paper (I imagine this could be used widely).

__ Weaknesses__: There are a few ways in which the experiments could be further improved:
Generally, I think more could be done to empirically understand how convergence is affected by various properties of C and the distributions involved.
The experiments seem to exclusively show convergence plots in terms of the potentials -- many readers will ask: how does the algorithm here compare in terms of the estimated OT distance itself?
I am always a bit suspicious of optimization papers that label the x-axis by iterations (or in this case samples) and not by time, though it may be unfair to compare the algorithm here to highly optimized (e.g. GPU-accelerated) Sinkhorn iterations. Could the authors comment on actual runtime?

__ Correctness__: The theoretical claims seem correct, or at least correct enough that they could be fixed should a bug be found. The experiments also seem correct.

__ Clarity__: For the most part I found the paper a joy to read, though the notation got a bit heavy on page 4.
I also got a bit confused at lines 288-289 -- I had trouble figuring out whether it was a two step process (f and g are fit until C is formed as a consequence, then run Sinkhorn) or whether the two steps happened in parallel (f and g are fit, while C is being built on another thread, then when C is done, start running Sinkhorn with whatever f and g are now).

__ Relation to Prior Work__: The paper is well contextualized; I found it clearly demarcated how it is different than prior work.

__ Reproducibility__: Yes

__ Additional Feedback__: The constants in the bounds, e.g. A in proposition 2, look a bit scary -- "It behaves like exp(1/eps)." Any intuition as to whether this can be improved?

__ Summary and Contributions__: This paper considers the computation of the entropy regularized optimal transport distance. In particular, unlike the classic Sinkhorn algorithm which assumes an emprical approximation of the underlying measures is available beforehand, this paper proposes an online Sinkhorn method which is able to compute the said distance using a growing sample from the underlying measures. Under the premise that the batch size per round is increasing over time, the proposed method is proved to find the true Sinkhorn potential with a nearly-O(1/n) asymptotic sample complexity.

__ Strengths__: The proposed method has a lower per-iteration computational complexity than the classic Sinkhorn algorithm as it handles the empirical approximation of the underlying measures like streaming data. Concretely, the per-iteration computational complexity is $O(N n_t d)$ where $N$ is the total number of available samples, $d$ is the dimension of the ground space, and $n_t$ is the (growing) batch size in round t.
In contrast, under an $O(Nd)$ memory constraint, the classic Sinkhorn algorithm would require $O(N^2d)$ computations per round.
Note that the per-iteration cost of the classic Sinkhorn algorithm can be reduced to $O(N^2)$ if the pairwise cost matrix is pre-computed and stored (which of course would entail a $O(N^2 d)$ computation expense and a $O(N^2)$ memory expense).
Besides, the author proves that if the size of the sample batch available per iteration grows (but can be arbitrarily slow), the proposed online Sinkhorn algorithm is able to compute the Sinkhorn potential at the rate of $O(1/N)$ where N is the number of the total available samples.

__ Weaknesses__: I find the Proposition 4 surprising. In the limiting case $\iota$ -> 0, P4 states that the sample complexity is $O(1/N)$, but P2 means that for $\iota=0$, the proposed online Sinkhorn algorithm only converges to a neighbor of the true Sinkhorn potential. So there is a bit jump here.
======update after author response=======
In Prop 4., if we take \iota to be arbitrarily close to 0, the proposed algorithm has the best performance O(1/N). However, if we take \iota to be exactly zero, according to prop 2, there is no convergence guarantee.
So what should we expect in practice? When taking a sufficiently small \iota, all finite iterations have the same batch size, should it converge or not? Either way, one of the proposition is violated.
The authors should address this in their revision.

__ Correctness__: I am not sure about that. Please see the discussion in the Weaknesses section.

__ Clarity__: I do not see how Proposition 2 agrees with the limiting case of Proposition 4. Maybe the author can elaborate on that.

__ Relation to Prior Work__: Yes

__ Reproducibility__: Yes

__ Additional Feedback__:

__ Summary and Contributions__: The authors propose a novel scheme for online computation of the sinkhorn divergence, of great interest for the ML community. Claims (different modes of convergence of the algorithm) are substantiated by thorough theoretical analysis and experiments on simple datasets of up to 10d.
Post author feedback
Authors have addressed my points but I won't update score since it is already very good.

__ Strengths__: The paper has many strengths and is bound to generate impact in the machine learning community. Authors thoroughly elaborate the underlying theory and make the distinction (through a number of convergence results) of random and online sinkhorn. Specifically, authors appeal to stochastic optimization to 'solve' the lack of almost sure convergence of the 'naive' randomized algorithm.
There is an is in depth discussion about the trade-offs entailed by the fact that the batch size has to grow to obtain convergence.
Authors also show an interesting connection with mirror descent.
This is a timely contribution, this paper addresses a relevant computational issue with Sinkhorn divergences.
Comparisons with existing methods are sensible and experiments insightful (but see below). A number of computational 'tricks' (as warm starts) are also presented and the impactful is measured. This adds value.
The supplement is comprehensive and contains very nice additional experiments.

__ Weaknesses__: Experiments are performed in only simple cases. I understand this is mostly a theory paper but for the sake of impact the paper would really benefit from including experiments with more complex dataset, where the theoretical gains are observed in practice.

__ Correctness__: As far as I can judge there are no mistakes

__ Clarity__: Yes. Perhaps the notation is a bit hard to understand with n_t and \eta_t
Line 170 is a bit confusing. Could you clarify what you mean by "Soft C transform"

__ Relation to Prior Work__: The paper cites most of exciting work but it misses a citation to relevant similar work, though from a more statistic perspective https://arxiv.org/pdf/1812.09150.pdf

__ Reproducibility__: Yes

__ Additional Feedback__: Authors may discuss on the relation between their work and the one from Bercu and Bigot (see above). Also, does your intuition about stochastic EM relates somehow to the recently develop Sinkhorn EM algorithm? https://arxiv.org/abs/2006.16548
Authors may also cite https://arxiv.org/abs/1905.11882 as this works is stated for general spaces. The above paper contains a more general sample complexiy statement.
In Figure 1, authors may provide entire traces for all methods.