---
title: Optimal Transport Regularization (OTReg)
url: https://www.emergentmind.com/topics/optimal-transport-regularization-otreg
type: topic
---

# Optimal Transport Regularization (OTReg)

Optimal Transport Regularization (OTReg) denotes a family of constructions in which optimal transport is either used as a regularizer inside a larger variational model or is itself regularized by additional convex, nonconvex, or structural terms. In the cited literature, this label covers, among other uses, the dual-Lipschitz or Wasserstein-1 penalty for image restoration, entropic and non-entropic regularizations of Kantorovich transport, adaptive and stochastic regularization schedules in semi-discrete OT, and transport-based alignment losses for representation learning. Across these formulations, the common role of OTReg is to preserve geometric structure while modifying optimization, statistical, or robustness properties in a controlled way [2503.14947] [1610.06447].

## 1. Scope and canonical formulations

In the cited literature, OTReg is used for several mathematically distinct constructions. One recurrent template regularizes the transport plan itself:
\[
W_R(\mu,\nu)\;=\;\min_{\pi\in\Pi(\mu,\nu)}\;\Bigl\{\langle \pi,c_0\rangle+\epsilon\,R(\pi)\Bigr\},
\]
where \(R\) is a convex penalty on the coupling \(\pi\). Another template uses an OT-derived quantity as a penalty inside a separate inverse problem or learning objective, as in image restoration or speech-text alignment. A third template constrains structure directly, for example by imposing per-row regularity budgets or explicit cardinality constraints on the plan [2002.03967] [2310.02925] [2209.15466].

| Usage | Representative formulation | Role |
|---|---|---|
| OT as a regularizer in inverse problems | \(\displaystyle TV(u)+\frac{\alpha}{2}\|f-u-v\|_2^2+\lambda\|v\|_{Lip*}\) | cartoon-texture-noise decomposition |
| Regularized OT on couplings | \(\displaystyle \min_{\pi\in\Pi(\mu,\nu)} \langle \pi,c_0\rangle+\epsilon R(\pi)\) | smoothness, robustness, tractability |
| OT-based alignment loss | \(\displaystyle L_{total}=L_{CE}+\lambda_{OT}L_{OT}\) | speech-text alignment |

The first line is exemplified by the image-restoration model
\[
R_{OT}(u)=\|u-u_0\|_{Lip*}
          = \inf_m\left\{\int_\Omega |m|\,dx:\operatorname{div}m=u-u_0\right\},
\]
which turns a Wasserstein-1 quantity into a regularizer on the image variable. The second line includes entropic, Orlicz-space, quadratic, Rényi, \(\beta\)-potential, and other plan regularizations. The third line appears in spoken language models, where an optimal transport plan between speech and transcript embeddings is computed and then converted into a differentiable regularization loss [2503.14947] [1909.06082] [2508.08131].

## 2. OTReg in image restoration and inverse imaging

In "Image Restoration Models with Optimal Transport and Total Variation Regularization" [2503.14947], OTReg is formulated from the Kantorovich–Rubinstein dual representation of the Wasserstein-1 distance. For nonnegative densities \(\mu,\nu\) on a bounded domain \(\Omega\subset\mathbb R^2\) with the same total mass,
\[
W_1(\mu,\nu)
 = \sup_{\|\phi\|_{Lip}\le 1}\int_\Omega \phi(x)\,[\mu(x)-\nu(x)]\,dx,
\]
equivalently with \(\|\nabla \phi\|_\infty\le 1\). By convex duality,
\[
W_1(\mu,\nu)
= \inf_{m\in [L^1(\Omega)]^2}\left\{\int_\Omega |m(x)|\,dx:\operatorname{div}m=\mu-\nu\right\}.
\]
Setting \(v=u-u_0\), the OT regularizer becomes
\[
R_{OT}(u)=\|u-u_0\|_{Lip*}
=\inf_m\left\{\int_\Omega |m|\,dx:\operatorname{div}m=u-u_0\right\}.
\]

A central theoretical point is the relation to Meyer’s \(G\)-norm. The paper states
\[
\|v\|_{Lip*}=\inf_{v=\operatorname{div}m}\|\,|m|\,\|_{L^1}
\le \inf_{v=\operatorname{div}m}\|\,|m|\,\|_{L^\infty}=\|v\|_G,
\]
and identifies the dual-Lipschitz norm as a negative-Sobolev norm of order \((-1,1)\), namely the case \(p=1\) in
\[
\|v\|_{-1,p}=\inf_{v=\operatorname{div}m}\|m\|_{L^p},\qquad 1\le p\le\infty,
\]
with \(p=\infty\) recovering Meyer’s \(G\)-norm. The combined model studied in practice is
\[
\min_{u\in BV,\;v\in H}\left\{
TV(u)+\frac{\alpha}{2}\|f-u-v\|_2^2+\lambda\|v\|_{Lip*}
\right\},
\]
where \(H=\{v:\int_\Omega v\,dx=0\}\). In the stated interpretation, \(TV(u)=\int_\Omega |\nabla u|\) enforces piecewise constant cartoon structure, \(\|v\|_{Lip*}\) retains fine oscillations, and the \(L^2\) term forces \(u+v\approx f\) [2503.14947].

The numerical scheme alternates between an OT subproblem in \(v\), solved by a Primal–Dual Hybrid Gradient method, and a TV subproblem in \(u\), solved by an Augmented Lagrangian Method. The OT step uses the saddle-point Lagrangian
\[
L(m,v,\phi)=\int |m|+\langle \phi,\operatorname{div}m-v\rangle+\frac{\alpha}{2}\|f_1-v\|^2,
\]
with updates \(m\leftarrow \operatorname{shrink}(m+\mu\nabla\phi;\mu\lambda)\), \(v\leftarrow (v/\nu+\alpha f_1+\phi)/(\alpha+1/\nu)\), and \(\phi\leftarrow \phi+\tau(\operatorname{div}m-v)\). The TV step rewrites \(p=\nabla u\) and applies ALM or ADMM with a closed-form shrinkage update for \(p\). Replacing TV by the nonconvex "Log-TV" of Zhu et al. yields "KR–Log-TV," whose reported effect is to reduce staircasing and keep sharper edges. On "House," "Butterfly," "Parrot," "Barbara," and related tests with additive white Gaussian noise \(\sigma\in\{15/255,25/255\}\), the reported observations are that KR–TV preserves more meaningful contrast than ROF, KR–Log-TV further suppresses staircasing, and PSNR improved \(0.2\)–\(0.5\,\mathrm{dB}\) over ROF in the reported tests [2503.14947].

A related inverse-imaging formulation appears in "Parallel Unbalanced Optimal Transport Regularization for Large Scale Imaging Problems" [1909.00149]. There the unbalanced OT regularizer is written in Beckmann form as
\[
\operatorname{OTReg}_\mu(p,q)
=\min_{M,r}\;\sum_{k=1}^N\|M_{k,:}\|_2+\mu\|r\|_p^p
\quad\text{s.t.}\quad
\operatorname{divergence}(M)-q+p=r.
\]
The stated advantages are linear optimization-variable complexity, a fully parallelizable proximal solver with \(O(N)\) per iteration, and superior empirical performance on dynamical tracking applications in synthetic and real video [1909.00149].

## 3. Regularizing the transport problem itself

A large part of the OTReg literature concerns penalties or constraints placed directly on the transport plan. In "Regularized Optimal Transport and the Rot Mover's Distance" [1610.06447], the primal regularized problem is
\[
d_{C,\lambda,\phi}(p,q)=\min_{X\in\Pi(p,q)}\langle C,X\rangle+\lambda\,\phi(X),
\]
and is shown to be equivalent to a matrix-nearness problem with respect to a Bregman divergence:
\[
\min_{X\in\Pi(p,q)}\langle C,X\rangle+\lambda\phi(X)
=
\lambda\min_{X\in\Pi(p,q)}D_\phi(X\Vert M)+\mathrm{const},
\qquad
M=\nabla\psi(-C/\lambda).
\]
This yields the "rot mover’s distance" and a unified algorithmic framework encompassing Sinkhorn-type and non-Sinkhorn-type regularizers [1610.06447].

Other works broaden the class of admissible penalties. "Orlicz-space regularization for optimal transport and algorithms for quadratic regularization" [1909.06082] replaces Radon plans by densities \(\pi\in L^\Phi(X\times Y)\) and solves
\[
\min_{\substack{\pi\in L^\Phi(X\times Y),\;\pi\ge0\\ \pi_1=\mu,\;\pi_2=\nu}}
\int c\,\pi+\gamma\int \Phi(\pi).
\]
The paper derives the predual
\[
\max_{\alpha\in C(X),\beta\in C(Y)}
\left\{
\int_X\alpha\,d\mu+\int_Y\beta\,d\nu
-\gamma\int \tilde\Phi^*\!\Bigl(\frac{\alpha(x)+\beta(y)-c(x,y)}{\gamma}\Bigr)\,dx\,dy
\right\},
\]
and states strong duality and existence of primal solutions. For \(L^p\) regularization with \(p\ge2\), it also gives cyclic projection, dual gradient descent, simple fixed point, and Nesterov-accelerated methods [1909.06082].

Quadratic regularization is emphasized in graph settings. "Quadratically-Regularized Optimal Transport on Graphs" [1704.08200] studies
\[
\min_{f\ge0,\;Af=b} c^T f+\frac{\alpha}{2}\|f\|_2^2,
\]
states that quadratic regularization preserves sparsity and ensures a unique solution, and derives the dual
\[
\max_\phi\; \phi^T b-\frac{1}{2\alpha}\sum_{e=1}^m\bigl[(A^T\phi)_e-c_e\bigr]_+^2,
\]
whose Hessian is a weighted graph Laplacian. A damped Newton method then reduces each step to a Laplacian linear solve [1704.08200]. "A regularized Interior Point Method for sparse Optimal Transport on Graphs" [2307.05186] introduces coupled primal–dual proximal regularization,
\[
\min_{x\ge0,y}\; c^T x+\frac12\rho\|x-x_k\|^2+\frac12\delta\|y-y_k\|^2
\quad\text{s.t.}\quad
Ax+\delta(y-y_k)=b,
\]
adds a log-barrier, sparsifies the shifted Laplacian normal equations, and proves inner-IPM complexity \(O(n^2\log(1/\varepsilon))\) [2307.05186].

Robustness-motivated OTReg replaces entropy by alternative divergences. "Robust computation of optimal transport by \(\beta\)-potential regularization" [2212.13251] uses
\[
\Phi_\beta(P)=\sum_{i,j}\phi_\beta(P_{ij})
=\frac{1}{\beta(\beta-1)}\sum_{i,j}\bigl(P_{ij}^\beta-\beta P_{ij}+(\beta-1)\bigr),
\]
and solves
\[
\min_{P\in U(a,b)}\;\langle C,P\rangle+\lambda\Phi_\beta(P).
\]
The paper states that, under a cost-separation condition on outliers, the intermediate solution satisfies \(P^{(T)}_{ij}=0\) for all outlier columns \(j\in J\), so no mass is transported onto true outliers [2212.13251]. "Interpolating between Optimal Transport and KL regularized Optimal Transport using Rényi Divergences" [2404.18834] proposes
\[
OT_{1/\lambda,\alpha}(\mu,\nu)
=
\min_{\pi\in\Pi(\mu,\nu)}
\left\{
\langle c,\pi\rangle+\frac1\lambda D_\alpha(\pi\Vert \mu\otimes\nu)
\right\},
\]
with \(\alpha\in(0,1)\). Its interpolation theorem states that, for fixed \(\lambda>0\), \(\alpha\uparrow1\) recovers KL-regularized OT and \(\alpha\downarrow0\) recovers unregularized OT, while for fixed \(\alpha\), \(\lambda\to\infty\) also recovers unregularized OT [2404.18834].

Structural regularization can be imposed more explicitly. "Sparsity-Constrained Optimal Transport" [2209.15466] solves
\[
\min_{P\in\mathbb R_+^{m\times n}}
\langle P,C\rangle+\frac{\gamma}{2}\|P\|_F^2
\quad\text{s.t.}\quad
P\mathbf1_n=a,\;P^T\mathbf1_m=b,\;p_j\in B_k\ \forall j,
\]
where \(B_k=\{t\in\mathbb R^m:\|t\|_0\le k\}\). The paper states that \(k=1\) recovers unregularized OT plus a constant shift, while large \(k\) recovers \(\ell_2\)-regularized OT, thus making the method a middle ground between the two [2209.15466]. "Optimal Transport with Adaptive Regularisation" [2310.02925] replaces a global budget by per-row or per-column budgets. In the row-wise version,
\[
\min_{\pi\in U(\alpha,\beta)}\langle C,\pi\rangle
\quad\text{s.t.}\quad
\pi\in\mathcal B_\psi(\xi),
\]
with
\[
\mathcal B_\psi(\xi)=\{\pi\ge0:\forall i,\ \psi(\pi_{i,\cdot})\le\psi(e_\xi)\}.
\]
Its stated purpose is to enforce a minimum of smoothing per point, in contrast to a single global regularization parameter [2310.02925].

A final conceptual reformulation is given in "Regularized Optimal Transport is Ground Cost Adversarial" [2002.03967]. There any convex regularization is shown, via Fenchel duality, to admit an adversarial ground-cost interpretation:
\[
W_R(\mu,\nu)=
\sup_{u\in C(\mathcal X^2)}
\left\{
\operatorname{OTcost}_{c_0+\epsilon u}(\mu,\nu)-\epsilon R^*(u)
\right\}.
\]
For entropy, the corresponding adversarial cost satisfies
\[
c_R(x,y)=c_0(x,y)+\epsilon\log\frac{d\pi^*}{d(\mu\otimes\nu)}(x,y).
\]
This places OTReg not only in a computational but also in a robust-optimization framework [2002.03967].

## 4. Semi-discrete OTReg and adaptive regularization schedules

Semi-discrete OT provides a distinct OTReg regime in which the source measure is continuous and the target is discrete. In "Semi-Discrete Optimal Transport: Nearly Minimax Estimation With Stochastic Gradient Descent and Adaptive Entropic Regularization" [2405.14459], the entropic semi-dual objective is
\[
H_\epsilon(g)
=
-\int g^{c,\epsilon}(x)\,d\mu(x)-\sum_{j=1}^M g_j w_j,
\]
with
\[
g^{c,\epsilon}(x)
=
-\epsilon\log\sum_{j=1}^M
\exp\!\left(\frac{g_j-\frac12\|x-y_j\|^2}{\epsilon}\right)w_j.
\]
The stochastic gradient is
\[
\nabla_g h_\epsilon(x,g)_j=-w_j+\chi_j^\epsilon(x,g),
\]
and the proposed DRAG algorithm uses a decreasing schedule \(\epsilon_t=t^{-a}\), projected SGD, and Polyak averaging. The paper states that this yields nearly minimax rates, including \(O(t^{-b})\) mean-square error to the true unregularized potential and \(O(t^{-b})\) error for the OT map, with lower bounds \(t^{-2}\) for the potential and \(t^{-1}\) for the map [2405.14459].

"Decreasing Entropic Regularization Averaged Gradient for Semi-Discrete Optimal Transport" [2510.27340] sharpens this perspective. With \(\gamma_k=\gamma_1k^{-b}\), \(\epsilon_k=k^{-a}\), a projection onto a bounded convex set \(\mathcal C\), and Polyak averaging, DRAG is stated to achieve an unbiased \(\mathcal O(1/t)\) sample and iteration complexity for both the OT cost and the potential estimation, and a \(\mathcal O}(1/\sqrt t)\) rate for the OT map. The stated explanation is that early iterations use large \(\epsilon_t\) for strong convexity and low variance, while later iterations use small \(\epsilon_t\) so that the entropic bias vanishes [2510.27340].

A complementary route appears in "Characterizing and computing solutions to regularized semi-discrete optimal transport via an ordinary differential equation" [2504.03030]. Setting \(t=1-\epsilon\), the dual potential \(\psi(t)\) is characterized by
\[
\nabla^2_{\psi\psi}\Phi(\psi(t),t)\,\psi'(t)
+\partial_t\nabla_\psi\Phi(\psi(t),t)=0,
\]
with explicit initial condition at \(t=0\). The paper states global strong convexity of \(\psi\mapsto \Phi(\psi,t)\) on the affine subspace \(\{\psi:\sum\psi_k=0\}\), invertibility of the Hessian on that subspace, and a third-order Runge–Kutta discretization. Its reported numerical finding is that the ODE method is competitive with Newton’s method for squared Euclidean cost and more robust for higher powers of the Euclidean distance [2504.03030].

These semi-discrete developments build on the earlier semi-dual regularization framework of "Semi-dual Regularized Optimal Transport" [1811.05527], which shows that entropic regularization replaces nonsmooth \(c\)-transforms by soft \(c\)-transforms and turns several OT-based variational problems, including barycenters and gradient flows, into smooth convex problems with closed-form gradients and Hessians [1811.05527].

## 5. Algorithms and computational regimes

Algorithmically, OTReg spans matrix scaling, Bregman projection, primal–dual splitting, Newton-type methods, mirror descent, and interior-point methods. For entropic regularization, the baseline is Sinkhorn’s algorithm, which alternates diagonal rescalings of the Gibbs kernel \(K=\exp(-C/\varepsilon)\). "Numerical Methods for Large-Scale Optimal Transport" [2210.11368] states \(O(n^2)\) cost per Sinkhorn iteration and emphasizes that entropy yields a strongly convex surrogate and smooth dual variables. It also presents accelerated primal–dual gradient methods that can use any strongly-convex regularization [2210.11368].

Acceleration of Sinkhorn-type solvers is itself an OTReg topic. "Overrelaxed Sinkhorn-Knopp Algorithm for Regularized Optimal Transport" [1711.01851] introduces overrelaxed Bregman projections with parameter \(\omega\in[1,2)\), proves global convergence using a Lyapunov function \(F(\gamma)=\mathrm{KL}(\gamma^*\Vert\gamma)\), and reports a gain in convergence speed by an order of magnitude in certain low-\(\varepsilon\) regimes [1711.01851].

For general smooth convex regularizers, "Regularized Optimal Transport and the Rot Mover’s Distance" [1610.06447] develops the alternate scaling algorithm and the non-negative alternate scaling algorithm, both based on Dykstra’s algorithm with alternate Bregman projections. In the separable case, row and column projections can be solved by Newton–Raphson updates, and a sparse extension is given for high-dimensional settings [1610.06447]. "Interpolating between Optimal Transport and KL regularized Optimal Transport using Rényi Divergences" [2404.18834] instead uses a nested mirror-descent scheme,
\[
\widetilde P = P^{(k-1)}\circ \exp(-\eta_k\nabla f(P^{(k-1)})),
\qquad
P^{(k)} = SK(\widetilde P;r,c),
\]
combining a mirror step with a Sinkhorn projection onto the transport polytope [2404.18834].

Graph-structured problems admit more specialized second-order methods. In quadratically-regularized graph OT, the dual Hessian is a weighted graph Laplacian, so a damped Newton method requires only Laplacian linear solves [1704.08200]. In sparse graph OT with primal–dual proximal stabilization, the normal equations involve the shifted Laplacian-plus-mass matrix
\[
S_{\rho,\delta}=A(\Theta^{-1}+\rho I)^{-1}A^T+\delta I,
\]
and \(\mu\)-sparsification drops inactive columns of \(A\) while preserving full rank through the \(\delta I\) term [2307.05186].

Outside pure transport problems, algorithm design follows the induced variational structure. The image-restoration KR–TV model alternates PDHG for the OT subproblem and ALM or ADMM for the TV subproblem, with FFT-based linear solves when the imaging operator is convolutional [2503.14947]. The speech alignment OTReg computes an entropic OT plan by Sinkhorn, with each Sinkhorn iteration costing \(O(n_a n_g)\), and reports that with \(20\)–\(50\) Sinkhorn steps the overhead is manageable [2508.08131].

## 6. Applications and empirical behavior

Image restoration is a canonical OT-as-regularizer application. For KR–TV and KR–Log-TV, the reported empirical pattern is contrast preservation, residuals that keep mostly true noise rather than structure, suppression of staircasing, and sharper edges than ROF or pure Log-TV. The reported PSNR gain is \(0.2\)–\(0.5\,\mathrm{dB}\) over ROF on the stated denoising tests [2503.14947].

In spoken language models, "Optimal Transport Regularization for Speech Text Alignment in Spoken Language Models" [2508.08131] formulates speech-text alignment as entropic OT between transformed speech embeddings \(S=\{f_i\}\) and unique transcript embeddings \(T=\{g_j\}\), with cost
\[
C_{i,j}=1-\operatorname{cosine\_similarity}(f_i/\|f_i\|_2,g_j/\|g_j\|_2)\in[0,2].
\]
The OTReg loss is
\[
L_{OT}=L_{cost}(\gamma^*)+\lambda_{spr}L_{spr}(\gamma^*),
\]
and the total objective is \(L_{total}=L_{CE}+\lambda_{OT}L_{OT}\). In the reported English ASR results, Base + OTReg with \(\lambda_{OT}=0.3\) gives \(13.12\) WER on CoVoST-2 and \(6.90\) on FLEURS, compared with \(16.06\) and \(7.63\) for the CE-only base model; the paper states that OTReg enhances speech-text alignment, mitigates the modality gap, and improves cross-domain generalization [2508.08131].

Regularity-based OTReg has also been used for map estimation and domain adaptation. "Regularity as Regularization: Smooth and Strongly Convex Brenier Potentials in Optimal Transport" [1905.10812] imposes \(\alpha I\preceq \nabla^2\varphi \preceq \beta I\), obtains a transport map with two-sided distortion bounds, and reports that the SSNB estimator attains much smaller error constants in high dimension, domain-adaptation accuracy comparable to the best regularized OT methods with mean accuracy \(\approx 86\)–\(87\,\%\), and smooth, contrast-controlled color mappings in a color-transfer task [1905.10812]. "Recovery Bounds on Class-Based Optimal Transport: A Sum-of-Norms Regularization Framework" [1903.03850] adds row and column sum-of-norms penalties to encourage block-diagonal couplings and reports exact class-structure recovery under the stated geometric assumptions; on an MNIST→USPS imbalanced-class task, the best OT-SON score in the excerpt is \(0.550\), compared with \(0.421\) for OT-l1l2 and \(0.414\) for Sinkhorn [1903.03850]. "Optimal Transport with Adaptive Regularisation" [2310.02925] reports that EOTARI-d reaches \(67.4(2.9)\) on MNIST→USPS versus \(64.2(2.8)\) for EOT, while QOTARI-t reaches \(69.3(4.7)\) versus \(68.3(3.9)\) for QOT [2310.02925].

Robustness-oriented OTReg has been evaluated on contaminated data. For \(\beta\)-potential regularization, the reported toy example shows that entropic-OT jumps from approximately \(50\) to approximately \(90\) under a few outliers, whereas \(\beta\)-OT stays at approximately \(50\), and outlier detection via zero column mass achieves \(>95\%\) accuracy on MNIST vs Fashion-MNIST contamination [2212.13251]. On sparse graphs, the proximal-stabilized regularized IPM is reported to outperform Lemon once \(|E|\gtrsim 10^6\), with experiments on graphs up to \(50\)M edges and empirical scaling \(O(|E|^{1.3\div1.4})\) versus Lemon’s \(\sim O(|E|^{2.0})\) [2307.05186]. In large-scale imaging, unbalanced Beckmann OTReg is reported to improve support recovery and rMSE in dynamical tracking, and RPCA+UOT-DF is reported to recover a compressed real video foreground where RPCA, RPCA+L1-DF, and RPCA+BOT-DF fail [1909.00149].

## 7. Interpretive issues, recurring trade-offs, and misconceptions

A persistent misconception is that OTReg is a single method. The cited literature instead uses the term for several non-equivalent operations: penalizing couplings, constraining couplings, using transport-based norms as regularizers in inverse problems, imposing regularity on Brenier potentials, and building alignment losses from optimal plans. This plurality is visible already from the contrast between KR–TV image restoration, entropic and Orlicz-regularized Kantorovich problems, adaptive regularization on rows and columns, and OT-based losses in spoken language models [2503.14947] [1909.06082] [2508.08131].

A second misconception is that regularization is only a computational device. The cited works assign it broader roles. In the Brenier-potential framework, regularity is itself the regularizer, with smoothness and strong convexity used to mitigate the curse of dimensionality and to enable out-of-sample evaluation [1905.10812]. In class-based OT, a sum-of-norms penalty is used to recover block structure under strong cyclical monotonicity and related geometric assumptions [1903.03850]. In the ground-cost adversarial formulation, any convex plan regularization can be reinterpreted as an adversarial perturbation of the ground cost [2002.03967].

A third recurrent issue is the bias–tractability trade-off. Entropic regularization yields smooth objectives and efficient Sinkhorn solvers, but introduces bias and fully dense plans. Quadratic regularization preserves sparsity on graphs; \(\beta\)-potential regularization is designed for robustness to outliers; explicit cardinality constraints give direct control over the number of nonzeros; and adaptive or decreasing regularization schedules are introduced precisely to reduce or remove entropic bias while retaining computational advantages [1704.08200] [2212.13251] [2209.15466] [2510.27340]. This suggests that OTReg is best understood not as a single regularizer, but as a design space for trading off geometry preservation, sparsity, smoothness, robustness, and solver structure.

Source: https://www.emergentmind.com/topics/optimal-transport-regularization-otreg