REVIEW 3 major objections 3 minor 28 references
Noisy Pairwise-Comparison Random Search for Smooth Nonconvex Optimization
T0 review · 3 major / 3 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read A comparison-only random search can reach approximate first-order stationarity in smooth nonconvex ridge objectives with a number of noisy pairwise comparisons that scales with the intrinsic dimension k, not the ambient dimension d.
desk verdict The central O(k/ε²) squared-gradient claim is not supported by the paper's own lemmas; the paper is promising but needs a major fix before it can be cited. 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 orthogonal projector P=A^T(AA^T)^{-1}A onto the active subspace range(A^T). Its role is to show that the accept/reject decision for f(x)=g(Ax) depends only on Ps_t; since E||Ps_t||² = tr(P) = k, all d-dimensional Gaussian moments in the analysis collapse to k. The proof is carried by a one-step descent inequality for the improve-or-stay update, summed with step size α_t=α0/√(kT); in the confidence model, Bernstein's inequality on the aggregated vote supplies the exponential error control.
What would settle it
On a synthetic ridge family f(x)=||Ax−b||² with fixed k and growing d, run NCRS with the paper's true-k schedule and measure (1/T)Σ||∇f(θ_t)||². If the average squared gradient does not decay as O(√(k/T)) with a prefactor independent of d, Theorem 2.2 fails. A second, more subtle check: run the same algorithm with a k-agnostic schedule (e.g., α_t=1/√t); if the rate then degrades with d, the claim that NCRS adapts without knowing k is refuted.
Extended reading notes
Core claim
On its own terms, the paper's central claim is that accept-or-stay random search driven by a single noisy comparison inherits the intrinsic dimension of a ridge objective. Because f(x+αs) = g(Ax+αAPs) = f(x+αPs) when f(x)=g(Ax), the comparison and the accept/reject decision depend only on the projected direction Ps; the algorithm therefore behaves as a random line search on the k-dimensional active subspace. The one-step descent inequality p α sqrt(2/π) E||∇f(θ_t)||² ≤ E[f(θ_t)-f(θ_{t+1})] + L_f k α_t²/2 turns this into a stationarity bound whose dominant term is O(√(k/T)). For the confidence oracle, the same inequality holds with a majority vote, and the wrong-decision probability decays li
Load-bearing premise
The load-bearing premise is that the optimizer knows the intrinsic dimension k (and, in the uniform-margin rate, the horizon T) so it can set the step sizes α_t=α0/√(kT) and the vote size N=O(k/ε²); without that knowledge, the paper does not establish the stated O(k/ε²) and O(k²/ε⁴) rates.
Editorial extensions
If this is right
- For a smooth nonconvex ridge objective with a k-dimensional active subspace, ε-first-order stationarity is achieved with O(k/(p²ε²)) comparisons; in the deterministic case p=1/2 this is O(k/ε²), independent of d.
- Under the tie-aware confidence oracle, the same guarantee holds with O(k²/ε⁴) total comparisons once the vote size is at least O(L k/ε²), even though the oracle is allowed to become uninformative at near-ties.
- A bounded nuisance gradient outside the active subspace (size τ, living in an m-dimensional subspace) adds only a τ√m/p term, so the stationarity guarantee degrades gracefully under misspecification rather than collapsing.
- The empirical results on masked language modeling fine-tuning show both NCRS and the two-point baseline RSGF improving as intrinsic dimension drops, with NCRS the more sensitive of the two; the paper interprets this as evidence for the k-dependent convergence theory.
Reading between the lines
- The paper's proved schedules require k (and, in Theorem 2.2, T) as inputs, even though the algorithm itself does not observe k; a parameter-free version that estimates k online is the natural next step and is not supplied here.
- Because the descent analysis only uses E||Ps_t||²=tr(P)=k and the Gaussian identity E|⟨∇f,s_t⟩|=sqrt(2/π)||∇f||, the same rate should hold for any direction distribution with those moments; this suggests compressed or adaptive sampling directions could reproduce the guarantee with fewer random bits.
- A controlled synthetic experiment on f(x)=||Ax−b||² with fixed k and varying d would separate the d-independence claim from the real-model confounds in the paper's experiments; the theory predicts flat iteration counts as d grows.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes Noisy-Comparison Random Search (NCRS), a direct-search method that uses one noisy pairwise comparison per iteration, and a variant with a confidence-weighted majority vote for a gap-dependent comparison oracle. For smooth nonconvex objectives with a k-dimensional active subspace, it claims that NCRS attains averaged squared-gradient stationarity with O(k/(p^2 ε^2)) comparisons under a uniform-margin oracle, and that the vote variant attains O(k^2/ε^4) total comparisons under a confidence model. The paper also includes a nearly-ridge robustness extension, a comparison with existing methods, and experiments on language-model fine-tuning and preference-based RL.
Significance. If the stated rates were correct, this would be a notable advance: it would give the first intrinsic-dimension dependence for noisy comparison-based nonconvex optimization and would justify a simple, gradient-free method on low-rank objectives. The algorithmic idea is clean, and the nearly-ridge model is a reasonable way to capture misspecification. However, the manuscript currently proves only a first-order stationarity bound in the gradient norm, not the squared norm used in the theorems and in Table 1. The mismatch is systematic and affects the abstract, the main theorems, the corollaries, and the empirical framing; the stated O(k/(p^2 ε^2)) and O(k^2/ε^4) rates for E||∇f||^2 are not supported by the analysis.
major comments (3)
- [Section 2, Lemma 2.1 / Theorem 2.2; Appendix C] The proof of Lemma 2.1 establishes p α sqrt(2/π) E||∇f(θ_t)|| ≤ E[f(θ_t)-f(θ_{t+1})] + L_f k α_t^2/2. The key identity is E|⟨∇f(θ_t), s_t⟩| = sqrt(2/π)||∇f(θ_t)||, not a squared norm. Averaging with α_t = α_0/√(kT) gives (1/T)Σ E||∇f(θ_t)|| = O((1/p)√(k/T)), i.e., stationarity in the norm sense. Theorem 2.2 states (1/T)Σ E||∇f(θ_t)||^2 ≤ ε at T=O(k/(p^2 ε^2)); this does not follow. The abstract and Table 1 inherit the same unsupported claim. The statements must be corrected to E||∇f||≤ε or a separate argument for the squared-norm measure must be supplied.
- [Section 3, Lemma 3.1, Lemma 3.3, Proposition 3.4] There is an algebraic error in the combination of Lemmas 3.1 and 3.3. Lemma 3.1 contributes -α_t/√(2π) E||∇f(θ_t)||, while Lemma 3.3 contributes +γ_{N,r} α_t √(2/π) E||∇f(θ_t)|| = 2γ_{N,r} α_t/√(2π) E||∇f(θ_t)||. The combined coefficient is therefore (1-2γ_{N,r}), not (1-γ_{N,r}). Proposition 3.4, Theorem 3.5, and Remark 3.6 use the incorrect coefficient; in particular, choosing γ_{N,r} ≤ 1/2 does not guarantee a positive descent coefficient, and the claimed O(k^2/ε^4) complexity in Remark 3.6 is not established.
- [Section 2.1 and Remark 3.6; adaptation claim] Theorem 2.2 requires the step size α_t = α_0/√(kT), and Remark 3.6 chooses α and N using k (and L_f, Δ_f). Algorithm 1 and Algorithm 2 do not receive k as an input, so the statement that NCRS adapts 'despite never observing A or the intrinsic dimension k' is not supported by the theorems. If k is unknown, the stated k-dependent complexity is not achieved by the algorithm as written. The authors should either state that k and T are known to the method, provide a parameter-free or adaptive schedule that preserves the k-dependence, or substantially soften the adaptation claim.
minor comments (3)
- [Throughout] Cross-references are inconsistent: Lemma 2.1 is called 'Theorem 2.1' in the text following it and in Appendix C; Proposition 2.4 is called 'Theorem 2.4' in Appendix D; Section 3 has similar mismatches. Please unify the numbering.
- [Table 1 and abstract] The metric in Table 1 is E||∇f||^2, but the analysis supports E||∇f||. If the authors correct the stationarity measure, the table and abstract should be updated accordingly; the current presentation is internally inconsistent.
- [Remark 2.3] The remark states T=O(k/(p^2 ε^2)) for E||∇f||^2 ≤ ε, but the displayed inequality in Theorem 2.2 has the form O((1/p)√(k/T)). Please make the hidden constants and the norm/squared-norm distinction explicit.
Circularity Check
No circular derivation: NCRS's convergence analysis is self-contained; the only author-overlapping citation (GLD) is background, and the k-dependent rates come from assumptions, not fitted inputs.
full rationale
The derivation chain is self-contained. Lemma 2.1 (Appendix C) proves a one-step descent inequality by conditioning on theta_t, lower-bounding the improvement via Assumption 1.1, and using the Gaussian identity E|⟨∇f(θ_t),s_t⟩| = sqrt(2/π)‖∇f(θ_t)‖; the k-dependence enters through E‖Ps_t‖² = tr(P) = k, an honest computation with the active-subspace projector. Theorem 2.2 then averages this inequality with α_t = α0/√(kT); no coefficient is fitted to the final epsilon and no oracle constant is tuned to produce the stated rate. Section 3 follows the same structure: Lemmas 3.1-3.3 bound the true-improvement and ranking-error terms under Assumption 1.2, and Theorem 3.5/Remark 3.6 aggregate them; the NT = O(k²/epsilon⁴) bound follows from the assumed local growth rho(t) ≥ ct, the second-moment bound, and Bernstein's inequality, not from an ansatz imported by citation. The only author-overlapping citation is GLD [5] (Qiuyi Zhang is a co-author), used in the introduction and Table 1 as conceptual background: 'Methods like Gradientless Descent (GLD) [5] use line search to automatically find the right step size to exploit the inherent optimization dimension.' None of NCRS's convergence theorems invokes GLD's convex result; the self-citation is therefore not load-bearing. Two caveats belong to correctness risk, not circularity: (i) the stated step sizes α_t = α0/√(kT) and vote size N require k and T as inputs, so the paper's claim that NCRS adapts despite never observing A or the intrinsic dimension k is stronger than the implementable algorithm supports; (ii) Lemma 2.1's proof bounds the gradient norm E‖∇f‖, whereas Theorems 2.2/3.5 and Table 1 are written with E‖∇f‖²; a norm bound does not imply the squared-norm bound by Jensen or Cauchy-Schwarz, so the squared-stationarity statements as printed are not established by the displayed inequalities. These are proof and implementation gaps, not reductions of the result to its own inputs.
Assumptions & free parameters
free parameters (1)
- Step-size scale α0
assumptions (7)
- domain assumption Assumption 1.1: uniform-margin comparison oracle with advantage p≤1/2 for all pairs
- domain assumption Assumption 1.2: confidence oracle with E[sign(Δ)R̃]≥ρ(|Δ|), E[R̃²]≤Cρ(|Δ|), and ρ(t)≥ct near 0
- domain assumption Exact ridge structure f(x)=g(Ax) with rank(A)=k
- standard math f is L_f-smooth and bounded below
- standard math Gaussian directional sampling and Isserlis/Wick moment identities
- standard math Bernstein inequality for bounded random variables
- ad hoc to paper Step-size schedule α_t=α0/√(kT) can be set with knowledge of k and T
Cite this review
Pith. "Pith review of Noisy Pairwise-Comparison Random Search for Smooth Nonconvex Optimization." pith.science (2026). https://pith.science/paper/5ATCG3ST
@misc{pith2026260121166,
author = {Pith},
title = {Pith review of: Noisy Pairwise-Comparison Random Search for Smooth Nonconvex Optimization},
year = {2026},
howpublished = {\url{https://pith.science/paper/5ATCG3ST}},
note = {Machine review of arXiv:2601.21166}
}
abstract
We study smooth nonconvex optimization using only noisy pairwise comparisons, without access to gradients or function values. We propose Noisy-Comparison Random Search (NCRS), a simple direct-search method that samples random directions and performs accept/reject updates from comparison feedback. Under a low-dimensional active-subspace structure, NCRS adapts to the intrinsic dimension $k\le d$ rather than the ambient dimension $d$. For a uniform-margin comparison oracle with advantage $p$, NCRS achieves $\epsilon$-first-order stationarity with comparison complexity $\mathcal{O}(k/(p^2\epsilon^2))$. We also introduce a gap-dependent confidence model, where comparison reliability decreases as the objective-value gap between the two candidates becomes small, and analyze a confidence-weighted voting variant of NCRS. For this oracle, the method achieves $\epsilon$-first-order stationarity with total comparison complexity $\mathcal{O}(k^2/\epsilon^4)$. These results provide intrinsic-dimension convergence guarantees for noisy comparison-based random search in smooth nonconvex optimization.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Intrinsic dimensionality explains the effectiveness of language model fine-tuning
Armen Aghajanyan, Sonal Gupta, and Luke Zettlemoyer. Intrinsic dimensionality explains the effectiveness of language model fine-tuning. InProceedings of the 59th annual meeting of the association for computational linguistics and the 11th international joint conference on natural language processing, pages 7319–7328, 2021
2021
-
[2]
Dueling optimization with a monotone adversary
Avrim Blum, Meghal Gupta, Gene Li, Naren Sarayu Manoj, Aadirupa Saha, and Yuanyuan Yang. Dueling optimization with a monotone adversary. InInternational Conference on Algorithmic Learning Theory, pages 221–243. PMLR, 2024
2024
-
[3]
Deep reinforcement learning from human preferences.Advances in neural information processing systems, 30, 2017
Paul F Christiano, Jan Leike, Tom Brown, Miljan Martic, Shane Legg, and Dario Amodei. Deep reinforcement learning from human preferences.Advances in neural information processing systems, 30, 2017
2017
-
[4]
Stochastic first-and zeroth-order methods for nonconvex stochastic programming.SIAM journal on optimization, 23(4):2341–2368, 2013
Saeed Ghadimi and Guanghui Lan. Stochastic first-and zeroth-order methods for nonconvex stochastic programming.SIAM journal on optimization, 23(4):2341–2368, 2013
2013
-
[5]
Gradientless descent: High-dimensional zeroth-order optimization
Daniel Golovin, John Karro, Greg Kochanski, Chansoo Lee, Xingyou Song, and Qiuyi Zhang. Gradientless descent: High-dimensional zeroth-order optimization. InInternational Conference on Learning Representations, 2020
2020
-
[6]
Direct language model alignment from online ai feedback.arXiv preprint arXiv:2402.04792, 2024
Shangmin Guo, Biao Zhang, Tianlin Liu, Tianqi Liu, Misha Khalman, Felipe Llinares, Alexandre Rame, Thomas Mesnard, Yao Zhao, Bilal Piot, et al. Direct language model alignment from online ai feedback.arXiv preprint arXiv:2402.04792, 2024. 13
arXiv 2024
-
[7]
Query complexity of derivative-free optimization.Advances in Neural Information Processing Systems, 25, 2012
Kevin G Jamieson, Robert Nowak, and Ben Recht. Query complexity of derivative-free optimization.Advances in Neural Information Processing Systems, 25, 2012
2012
-
[8]
Kyuyoung Kim, Ah Jeong Seo, Hao Liu, Jinwoo Shin, and Kimin Lee. Margin matching preference optimization: Enhanced model alignment with granular feedback.arXiv preprint arXiv:2410.03145, 2024
arXiv 2024
Show all 28 references
-
[9]
Kingma and Jimmy Ba
Diederik P. Kingma and Jimmy Ba. Adam: A method for stochastic optimization. In Yoshua Bengio and Yann LeCun, editors,ICLR (Poster), 2015
2015
-
[10]
B-pref: Benchmarking preference- based reinforcement learning.arXiv preprint arXiv:2111.03026, 2021
Kimin Lee, Laura Smith, Anca Dragan, and Pieter Abbeel. B-pref: Benchmarking preference- based reinforcement learning.arXiv preprint arXiv:2111.03026, 2021
2021 arXiv
-
[11]
Fine-tuning language models with just forward passes.Advances in Neural Information Processing Systems, 36:53038–53075, 2023
Sadhika Malladi, Tianyu Gao, Eshaan Nichani, Alex Damian, Jason D Lee, Danqi Chen, and Sanjeev Arora. Fine-tuning language models with just forward passes.Advances in Neural Information Processing Systems, 36:53038–53075, 2023
2023
-
[12]
Direct preference optimization: Your language model is secretly a reward model
Rafael Rafailov, Archit Sharma, Eric Mitchell, Christopher D Manning, Stefano Ermon, and Chelsea Finn. Direct preference optimization: Your language model is secretly a reward model. Advances in neural information processing systems, 36:53728–53741, 2023
2023
-
[13]
Dueling convex optimization
Aadirupa Saha, Tomer Koren, and Yishay Mansour. Dueling convex optimization. In Marina Meila and Tong Zhang, editors,Proceedings of the 38th International Conference on Machine Learning, volume 139 ofProceedings of Machine Learning Research, pages 9245–9254. PMLR, 18–24 Jul 2021
2021
-
[14]
Faster convergence with multiway preferences
Aadirupa Saha, Vitaly Feldman, Yishay Mansour, and Tomer Koren. Faster convergence with multiway preferences. InInternational Conference on Artificial Intelligence and Statistics, pages 433–441. PMLR, 2024
2024
-
[15]
Dueling convex optimization with general preferences
Aadirupa Saha, Tomer Koren, and Yishay Mansour. Dueling convex optimization with general preferences. InForty-second International Conference on Machine Learning, 2025
2025
-
[16]
Proximal policy optimization algorithms.arXiv preprint arXiv:1707.06347, 2017
John Schulman, Filip Wolski, Prafulla Dhariwal, Alec Radford, and Oleg Klimov. Proximal policy optimization algorithms.arXiv preprint arXiv:1707.06347, 2017
2017 arXiv
-
[17]
Zeroth-order optimization meets human feedback: Provable learning via ranking oracles
Zhiwei Tang, Dmitry Rybin, and Tsung-Hui Chang. Zeroth-order optimization meets human feedback: Provable learning via ranking oracles. InThe Twelfth International Conference on Learning Representations, 2024
2024
-
[18]
Deepmind control suite
Yuval Tassa, Yotam Doron, Alistair Muldal, Tom Erez, Yazhe Li, Diego de Las Casas, David Budden, Abbas Abdolmaleki, Josh Merel, Andrew Lefrancq, et al. Deepmind control suite. arXiv preprint arXiv:1801.00690, 2018
2018 arXiv
-
[19]
Llama 2: Open foundation and fine-tuned chat models.arXiv preprint arXiv:2307.09288, 2023
Hugo Touvron, Louis Martin, Kevin Stone, Peter Albert, Amjad Almahairi, Yasmine Babaei, Nikolay Bashlykov, Soumya Batra, Prajjwal Bhargava, Shruti Bhosale, et al. Llama 2: Open foundation and fine-tuned chat models.arXiv preprint arXiv:2307.09288, 2023
2023 arXiv
-
[20]
dm control: Software and tasks for continuous control.Software Impacts, 6:100022, 2020
Saran Tunyasuvunakool, Alistair Muldal, Yotam Doron, Siqi Liu, Steven Bohez, Josh Merel, Tom Erez, Timothy Lillicrap, Nicolas Heess, and Yuval Tassa. dm control: Software and tasks for continuous control.Software Impacts, 6:100022, 2020. 14
2020
-
[21]
GLUE: A multi-task benchmark and analysis platform for natural language understanding
Alex Wang, Amanpreet Singh, Julian Michael, Felix Hill, Omer Levy, and Samuel Bowman. GLUE: A multi-task benchmark and analysis platform for natural language understanding. In Tal Linzen, Grzegorz Chrupa la, and Afra Alishahi, editors,Proceedings of the 2018 EMNLP Workshop Bla...
2018 doi
-
[22]
Preference-based optimization from noisy pairwise comparisons.arXiv preprint arXiv:2512.18511, 2025
Siyi Wang, Zifan Wang, and Karl Henrik Johanssson. Preference-based optimization from noisy pairwise comparisons.arXiv preprint arXiv:2512.18511, 2025
2025
-
[23]
Explicit and non-asymptotic query complexities of rank-based zeroth-order algorithms on smooth functions.arXiv preprint arXiv:2512.16200, 2025
Haishan Ye. Explicit and non-asymptotic query complexities of rank-based zeroth-order algorithms on smooth functions.arXiv preprint arXiv:2512.16200, 2025
2025
-
[24]
Zeroth-order policy gradient for reinforcement learning from human feedback without reward inference
Qining Zhang and Lei Ying. Zeroth-order policy gradient for reinforcement learning from human feedback without reward inference. InThe Thirteenth International Conference on Learning Representations, 2025. 15 A Preliminary Lemma A.1.Let f : Rd →R be differentiable and let k∈ {...
2025
-
[25]
(Gradient lives in a fixed k-subspace)There exists a linear subspace V⊂R d with dim(V) =ksuch that: ∀x∈R d,∇f(x)∈V
-
[26]
(2) ⇒ (1).Let x∈R d
(Ridge representation)There exist a full row-rank matrix A∈R k×d and a differentiable functiong:R k →Rsuch that: ∀x∈R d, f(x) =g(Ax), Proof. (2) ⇒ (1).Let x∈R d. We have ⟨∇f (x), v⟩= 0 for all v∈Ker (A). This implies that ∇f(x)∈range(A ⊤) =Ker(A) ⊥,and we have dim(range(A ⊤)) ...
-
[27]
(Gradient leakage outside a k-subspace)There exist a k-dimensional subspace V⊂R d and its orthogonal projectorPsatisfying: ∀x∈R d,∥(I−P)∇f(x)∥ 2 ≤τ
-
[28]
prefer y over x
(Nearly-ridge representation with controlled leakage)There exist a full row-rank matrix A∈R k×d, a differentiable function g : Rk →R , and a differentiable residual η : Rd →R such that, letting P denote the orthogonal projector onto V := range(A⊤), we have: ( ∀x∈R d, f(x) =g(A...
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.