REVIEW 2 major objections 6 minor 33 references
Accelerated Bregman gradient methods for relatively smooth and relatively Lipschitz continuous minimization problems
T0 review · 2 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Accelerated Bregman proximal gradient methods keep their convergence rate when the oracle is inexact, up to an additive noise term, provided the Bregman divergence obeys a triangular scaling inequality with exponent between 1 and 2.
desk verdict Useful new accelerated Bregman methods with inexact oracles, but the advertised Poisson inverse problem application is not covered by the theory because its divergence has the wrong scaling exponent. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the triangular scaling property of the Bregman divergence $V(x,y)=d(x)-d(y)-\langle\nabla d(y),x-y\rangle$: the inequality $V((1-\theta)x+\theta z,(1-\theta)x+\theta \tilde z)\le \theta^\gamma V(z,\tilde z)$ for every $\theta\in[0,1]$. This inequality is what lets the proofs bound the Bregman distance between the convex combination points $y_k=(1-\theta_k)x_k+\theta_k z_k$ and $x_{k+1}=(1-\theta_k)x_k+\theta_k z_{k+1}$ by $\theta_k^\gamma V(z_{k+1},z_k)$, converting the algorithm's updates into a telescoping Lyapunov inequality. The second ingredient is the $(\delta,L)$-oracle model $0\le f(x)-\left(f_\delta(y)+\langle\nabla f_\delta(y),x-y\rangle\right)\le LV(x,y)+\delta$, and the third is relative smoothness, $f(y)\le f(x)+\langle\nabla f(x),y-x\rangle+LV(y,x)$, which replaces the usual Lipschitz-gradient upper model.
What would settle it
Evaluate the triangular scaling inequality numerically for the prox-function used in an application: draw random triples $(x,z,\tilde z)$ and compare $V((1-\theta)x+\theta z,(1-\theta)x+\theta\tilde z)$ with $\theta^\gamma V(z,\tilde z)$. If the smallest exponent that fits is $\le1$, then the proof step behind Theorem 4.5 is invalid for that problem and the accelerated guarantee cannot be certified; for the Poisson experiments with the log-barrier prox, the paper's own Example 2.9 indicates exactly this situation.
Extended reading notes
Core claim
The central claim is that the inexact-oracle accelerated Bregman proximal gradient method, under the triangular scaling property with exponent $\gamma\in(1,2]$, satisfies $f(x_N)-f(x_*) \le 2L\left(\frac{\gamma}{\gamma+N-1}\right)^\gamma V(x_*,x_0)+(2(N-1)L+1)\delta$ after $N$ iterations; in the noiseless limit $\delta=0$ this recovers the $O(1/N^\gamma)$ rate of the exact-oracle method. The paper also proves a similar-triangles fast gradient method with rate $8LR^2/(N+1)^2 + \frac{2}{A_N}\sum_{k=0}^{N-1}A_{k+1}\delta_k$, and an adaptive intermediate Bregman method whose rate $(k+2)^{-(p-1)(\gamma-1)-1}$ interpolates between non-accelerated ($p=1$) and accelerated ($p=2$) behaviour. The universality statements extend both accelerated families to relatively Lipschitz continuous objectives, at the cost of a slightly degraded oracle in the line-search test.
Load-bearing premise
The entire argument rests on the chosen Bregman divergence satisfying the triangular scaling inequality with an exponent strictly between 1 and 2; the paper's own examples show the Kullback-Leibler divergence has exponent 1 and the Itakura-Saito divergence has exponent below 1, so the accelerated guarantees do not cover those natural measures.
Editorial extensions
If this is right
- For any relatively smooth problem whose prox-function has a triangular scaling factor $\gamma\in(1,2]$, inexact oracles preserve the accelerated rate; noise only adds a term proportional to $\delta$.
- The same iterations solve relatively Lipschitz continuous problems, so the methods are universal in the sense the paper states.
- The adaptive intermediate method gives a user a dial $p\in[1,2]$ between non-accelerated and accelerated convergence, without switching algorithms.
- In the Euclidean case $\gamma=2$, the intermediate method's rate is $(k+2)^{-p}$ plus $(k+2p)^{p-1}\delta$, making the interpolation explicit.
- With fixed noise $\delta$, the bound in Theorem 4.5 grows linearly in $N$, so the method is best viewed as achieving accuracy of order $\delta$ rather than converging to zero.
Reading between the lines
- A design lesson the paper leaves implicit is that the triangular scaling exponent should be treated as a selection criterion for the prox-function, since the Kullback-Leibler and Itakura-Saito divergences lie outside the accelerated regime.
- One testable extension would be to estimate $\gamma$ numerically for a candidate divergence and use it as a pre-screening test before committing to an accelerated Bregman method.
- The growth of the noise term in Theorem 4.5 suggests that when $\delta$ is fixed and not small, a user should switch to the intermediate method; this is a practical reading of the interpolation the authors construct.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes accelerated Bregman gradient methods with an inexact oracle for minimizing convex functions that are relatively smooth or relatively Lipschitz continuous with respect to a reference prox-function. It presents four algorithms: an adaptive fast gradient method based on similar triangles (Algorithm 1), a non-adaptive accelerated Bregman proximal gradient method (Algorithm 2), an adaptive version of it (Algorithm 3), and an adaptive intermediate Bregman method (Algorithm 4) that interpolates between accelerated and non-accelerated behavior. Convergence-rate theorems are stated for the relatively smooth case (Theorems 3.1, 4.1, 4.5, and 5.4) in the presence of an inexact oracle, and the universality of the methods for relatively Lipschitz continuous problems is argued informally. The paper concludes with numerical experiments on a Poisson inverse problem comparing the proposed algorithms with an adaptive Bregman proximal gradient method.
Significance. If the results are correct, they would constitute a useful extension of accelerated Bregman proximal methods to the inexact-oracle setting and introduce an interesting noise-robust intermediate method. The proofs are self-contained, and the non-adaptive Theorem 4.1 appears coherent. However, the actual scope is considerably narrower than advertised: the accelerated guarantees require the triangular scaling property with exponent γ strictly in (1,2], while the motivating Poisson inverse problem uses a divergence with γ<1, and numerical experiments treat γ as a freely tunable parameter. In addition, the adaptive line-search condition in Algorithm 3 is scale-inconsistent and breaks the proof of Theorem 4.5. These are load-bearing issues that prevent acceptance in the present form.
major comments (2)
- [§4.2, Algorithm 3, line 8] The line-search condition L_{k+1} ≥ (θ_{k+1}/θ_k)^γ in Algorithm 3 is scale-inconsistent: the right-hand side is dimensionless, while L_{k+1} scales with the relative smoothness parameter. If the true parameter L is less than 1, then for large k this condition forces L_{k+1} to be at least about 1, and the assumption L_1 < 2L need not imply L_N < 2L. The proof of Theorem 4.5 explicitly uses this condition in the chain following Eq. (47) to replace θ_k^γ L_{k+1} by θ_{k+1}^γ, so the final bound (46) is not justified. The condition should be replaced by a scale-corrected version, for example L_{k+1} ≥ (θ_{k+1}/θ_k)^γ L_k or another form that respects the scaling of L.
- [§6 and Examples 2.8–2.9] The accelerated convergence results in Theorems 3.1, 4.1, and 4.5 assume the triangular scaling property with γ ∈ (1,2]. The paper's own Examples 2.8 and 2.9 state that the KL divergence has γ=1 and the Itakura–Saito divergence has γ<1, while Example 2.3, which motivates the paper, uses the log-barrier prox-function d(x)=−Σ ln x_i, whose Bregman divergence is VIS. Hence γ<1 for the Poisson inverse problem, and the convergence theorems do not apply to the problem advertised in the abstract. The numerical experiments in Section 6 run the same Poisson problem with γ=2, 1.4, and 1.1, but γ is a structural property of the divergence rather than an algorithmic input; no prox-function satisfying triangular scaling with those exponent values is specified for the experiments. The claims and experiments should be restricted to divergences with γ∈(1,2], or the experimental setup must be justified by explicit prox-functions for each γ.
minor comments (6)
- [Example 2.9] The name 'Itakuro–Saito divergence' is a typo; it should be 'Itakura–Saito divergence'.
- [Algorithm 1, line 4] The update L_{k+1}=2^{i_k−1}L_k gives L_{k+1}=L_k/2 when i_k=0; clarify whether the exponent is intended to be i_k or explain why halving is allowed, and state that the line search terminates in finite time.
- [Algorithm 4, lines 3–5] The while condition at line 3 refers to y0 before y0 is computed in line 5; the pseudocode should be reordered so that y0 is initialized before the condition is evaluated, or a repeat-until form should be used.
- [Theorem 5.4 and Corollary 5.3] The statement uses V(x0,x*)=d(x*)≤R0, but this equality holds only when x0 minimizes d and d(x0)=0; the assumption should be stated as V(x0,x*)≤R0 or d(x*)≤R0.
- [Eq. (52) and Theorem 5.4 proof] The assumption 0<α_k≤B_k≤A_k is stated for all k but is not verified in the paper; please add a short proof or a reference showing that the updates of Algorithm 4 maintain these inequalities.
- [General] The manuscript contains numerous typographical and grammatical errors (e.g., 'we concluded its convergence rate', 'Itakuro–Saito', 'with constant γ >0'); a thorough language revision is recommended.
Circularity Check
No significant circularity: the convergence proofs are self-contained and reduce to the stated assumptions, not to fitted constants or self-citation chains.
full rationale
The derivation chain is self-contained. Theorems 3.1, 4.1, 4.5, and 5.4 each prove their rates by explicit telescoping arguments (e.g., summing inequalities (24), (35), and (42)) from the standing assumptions: relative smoothness (10), the inexact-oracle inequality (16), and the triangular scaling property (17). No fitted parameter is renamed as a prediction: the parameters L, delta, and gamma are inputs or structural assumptions, not quantities tuned to the target bound. The cited works ([6] for the triangular scaling property and the APBG analogy, [12] for the similar-triangles framework, [28] for the relatively-Lipschitz inequality) are external lemmas with stated assumptions and are used as tools rather than as substitutes for the target convergence proofs. The numerical experiments' use of gamma = 2, 1.4, and 1.1 for a problem whose Bregman divergence is the Itakura-Saito divergence (gamma < 1 by Example 2.9) is an application-consistency issue, not a circular one. Thus no step reduces the paper's conclusions to its own inputs by construction.
Assumptions & free parameters
free parameters (1)
- gamma (triangular scaling factor) in numerical experiments =
2, 1.4, 1.1
assumptions (4)
- domain assumption Bregman divergence satisfies the triangular scaling property with exponent gamma in (1,2].
- domain assumption The pair (f_delta, grad f_delta) is a (delta,L)-oracle satisfying Definition 2.5.
- ad hoc to paper The relatively Lipschitz inequality from [28] (f(x) <= f(y)+<grad f(y),x-y>+L V(x,y)+L V(y,x)+delta) supports the inexact-oracle line search.
- ad hoc to paper The adaptive condition L_{k+1} >= (theta_{k+1}/theta_k)^gamma in Algorithm 3 is scale-invariant.
Cite this review
Pith. "Pith review of Accelerated Bregman gradient methods for relatively smooth and relatively Lipschitz continuous minimization problems." pith.science (2026). https://pith.science/paper/V3FS5UXV
@misc{pith2026241116743,
author = {Pith},
title = {Pith review of: Accelerated Bregman gradient methods for relatively smooth and relatively Lipschitz continuous minimization problems},
year = {2026},
howpublished = {\url{https://pith.science/paper/V3FS5UXV}},
note = {Machine review of arXiv:2411.16743}
}
read the original abstract
In this paper, we propose some accelerated methods for solving optimization problems under the condition of relatively smooth and relatively Lipschitz continuous functions with an inexact oracle. We consider the problem of minimizing the convex differentiable and relatively smooth function concerning a reference convex function. The first proposed method is based on a similar triangles method with an inexact oracle, which uses a special triangular scaling property for the used Bregman divergence. The other proposed methods are non-adaptive and adaptive (tuning to the relative smoothness parameter) accelerated Bregman proximal gradient methods with an inexact oracle. These methods are universal in the sense that they are applicable not only to relatively smooth but also to relatively Lipschitz continuous optimization problems. We also introduced an adaptive intermediate Bregman method which interpolates between slower but more robust algorithms non-accelerated and faster, but less robust accelerated algorithms. We conclude the paper with the results of numerical experiments demonstrating the advantages of the proposed algorithms for the Poisson inverse problem.
Figures
Reference graph
Works this paper leans on
-
[1]
Lectures on convex optimization
Nesterov Y. Lectures on convex optimization. Vol. 137, Springer International Publishing, 2018. 30
work page 2018
-
[2]
First-Order Methods in Optimization
Beck A. First-Order Methods in Optimization. Philadelphia: Society for Industrial and Applied Mathematics; 2017
work page 2017
-
[3]
Lu H. relative continuity for non-lipschitz nonsmooth convex optimization using stochastic (or deterministic) mirror descent. IJOO. 2019;1(4):288-303
work page 2019
-
[4]
A descent Lemma beyond Lipschitz gradient continuity: first-order method revisited and applications
Bauschke H, Bolte J, Teboulle M. A descent Lemma beyond Lipschitz gradient continuity: first-order method revisited and applications. Math. Oper. Res. 2017;42(2):330-348. doi: https://doi.org/10.1287/moor.2016.0817
arXiv 2017
-
[5]
Relatively smooth convex optimization by first-order methods, and applications
Lu H, Freund R, Nesterov Y. Relatively smooth convex optimization by first-order methods, and applications. SIOPT. 2018;28(1):333-354. doi: https://doi.org/10.1137/16M1099546
-
[6]
Accelerated Bregman proximal gradient methods for relatively smooth convex optimization
Hanzely F, Richtarik P, Xiao L. Accelerated Bregman proximal gradient methods for relatively smooth convex optimization. Comput. Optim. Appl. 2021;79:405–440. doi: https://doi.org/10.1007/s10589-021-00273-8
-
[7]
Implementable tensor methods in unconstrained convex optimization
Nesterov Y. Implementable tensor methods in unconstrained convex optimization. Math. Program. 2021;186:157–183. doi: https://doi.org/10.1007/s10107-019-01449-1
-
[8]
Inexact accelerated high-order proximal-point methods
Nesterov Y. Inexact accelerated high-order proximal-point methods. Math. Program. 2023;197:1-26. doi: https://doi.org/10.1007/s10107-021-01727-x
Show all 33 references
-
[9]
Optimization methods for large-scale machine learning
L´ eon B, Curtis F, Nocedal J. Optimization methods for large-scale machine learning. SIAM review. 2018;60(2):223-311. doi: https://doi.org/10.1137/16M1080173
2018 doi
-
[10]
Random search for hyper-parameter optimization
Bergstra J, Bengio Y. Random search for hyper-parameter optimization. J. Mach. Learn. Res. 2012;13(10):281-305
2012
-
[11]
First-order methods of smooth convex optimization with inexact oracle
Devolder O, Glineur A, Nesterov Y. First-order methods of smooth convex optimization with inexact oracle. Math. Program. 2014;146:37–75. doi: https://doi.org/10.1007/s10107- 013-0677-5
2014 doi
-
[12]
Inexact model: A framework for optimiza- tion and variational inequalities
Stonyakin F, Tyurin A, Gasnikov A, et al. Inexact model: A framework for optimiza- tion and variational inequalities. Optim. Methods Softw. 2021;36(6):1155–1201. doi: https://doi.org/10.1080/10556788.2021.1924714
2021
-
[13]
Universal intermediate gradient method for convex problems with inexact oracle
Kamzolov D, Dvurechensky P, Gasnikov A. Universal intermediate gradient method for convex problems with inexact oracle. Optim. Methods Softw. 2021;36(6):1289-1316. doi: https://doi.org/10.1080/10556788.2019.1711079
2021
-
[14]
Intermediate gradient methods for smooth convex problems with inexact oracle
Devolder O, Glineur F, Nesterov Y. Intermediate gradient methods for smooth convex problems with inexact oracle. (No. UCL-Universite Catholique de Louvain). Technical report, CORE-2013017, 2013
2013
-
[15]
Stochastic intermediate gradient method for convex prob- lems with stochastic inexact oracle
Dvurechensky P, Gasnikov A. Stochastic intermediate gradient method for convex prob- lems with stochastic inexact oracle. J. Optim. Theory Appl. 2016;171:121-145. doi: https://doi.org/10.1007/s10957-016-0999-6
2016 doi
-
[16]
A fast iterative shrinkage-thresholding algorithm for linear inverse problems
Beck A, Teboulle M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems. SIIMS. 2009;2(1):183-202. doi: https://doi.org/10.1137/080716542
2009 doi
-
[17]
Gradient methods for minimizing composite functions
Nesterov Y. Gradient methods for minimizing composite functions. Math. Program. 2013;140:125-161. doi: https://doi.org/10.1007/s10107-012-0629-5
2013 doi
-
[18]
Interior gradient and proximal methods for convex and conic op- timization
Auslender A, Teboulle M. Interior gradient and proximal methods for convex and conic op- timization. SIOPT. 2006;16(3):697-725. doi: https://doi.org/10.1137/S105262340342782
2006 doi
-
[19]
Problem Complexity and Method Efficiency in Optimization
Nemirovsky A, Yudin D. Problem Complexity and Method Efficiency in Optimization. New York: J. Wiley and Sons; 1983
1983
-
[20]
Image deblurring with Poisson data: from cells to galaxies
Bertero M, Boccacci P, Desidera G, et al. Image deblurring with Poisson data: from cells to galaxies. Inverse Problems. 2009;25(12). doi: 10.1088/0266-5611/25/12/123006
2009 doi
-
[21]
Why least squares and maximum entropy an axiomatic approach to inference for linear iverse problems
Csiszar I. Why least squares and maximum entropy an axiomatic approach to inference for linear iverse problems. Ann. Stat. 1991;19(4):2032-2066. doi: 10.1214/aos/1176348385
1991
-
[22]
Optimum design in regression problems
Kiefer J, Wolfowitz J. Optimum design in regression problems. Ann. Math. Stat. 1959;30(2):271-294. doi: 10.1214/aoms/1177706252
1959
-
[23]
Optimal and efficient designs of experiments
Atwood C. Optimal and efficient designs of experiments. Ann. Math. Stat. 1969;40(5):1570-1602. doi: https://doi.org/10.1214/aoms/1177697374
1969
-
[24]
A simple convergence analysis of Bregman proximal gradient algorithm
Zhou Y, Liang Y, Shen L. A simple convergence analysis of Bregman proximal gradient algorithm. Comput. Optim. Appl. 2019;93:903-912. doi: https://doi.org/10.1007/s10589- 31 019-00092-y
2019 doi
-
[25]
Convergence analysis of a proximal-like minimization algorithm us- ing Bregman functions
Chen G, Teboulle M. Convergence analysis of a proximal-like minimization algorithm us- ing Bregman functions. SIOPT. 1993;3(3):538–543. doi: https://doi.org/10.1137/0803026
1993 doi
-
[26]
Learning supervised pagerank with gradient-based and gradient-free optimization methods
Bogolubsky L, Dvurechensky P, Gasnikov A, et al. Learning supervised pagerank with gradient-based and gradient-free optimization methods. In Advances in Neural Informa- tion Processing Systems 29, D.D. Lee, M. Sugiyama, U.V. Luxburg, I. Guyon, and R. Garnett, eds. Curran Assoc...
2016
-
[27]
Universal gradient methods for convex optimization problems
Nesterov Y. Universal gradient methods for convex optimization problems. Math. Pro- gram. 2015;152:381-404. doi: https://doi.org/10.1007/s10107-014-0790-0
2015 doi
-
[28]
Adaptive Algorithms for Relatively Lips- chitz Continuous Convex Optimization Problems
Stonyakin F, Alkousa M, Titov T, et al. Adaptive Algorithms for Relatively Lips- chitz Continuous Convex Optimization Problems. Pure and Applied Functional Analysis. 2023;8(5):1505-1526
2023
-
[29]
Adaptive Methods or Variational Inequalities with Relatively Smooth and Reletively Strongly Monotone Operators
Ablaev S, Stonyakin F, Alkousa M, et al. Adaptive Methods or Variational Inequalities with Relatively Smooth and Reletively Strongly Monotone Operators. Program Comput Soft. 2023;(49):485–492. doi: https://doi.org/10.1134/S0361768823060026
2023 doi
-
[30]
Lectures on modern convex optimization: analysis, algorithms, and engineering applications
Ben-Tal A, Nemirovski A. Lectures on modern convex optimization: analysis, algorithms, and engineering applications. SIAM; 2001
2001
-
[31]
A simplified view of first order methods for optimization
Teboulle M. A simplified view of first order methods for optimization. Math. Program. 2018;170:67-96. doi: https://doi.org/10.1007/s10107-018-1284-2
2018 doi
-
[32]
Pegasos: primal estimated sub-gradient solver for SVM
Shwartz S, Singer Y, Srebro N, et al. Pegasos: primal estimated sub-gradient solver for SVM. Math. Program. 2011;127:3–30. doi: https://doi.org/10.1007/s10107-010-0420-4
2011 doi
-
[33]
Bregman gradient methods for relatively-smooth optimization
Dragomir A. Bregman gradient methods for relatively-smooth optimization. PhD thesis. UT1 Capitole, 2021. 32
2021
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.