Tianyi Lin
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Tianyi Lin has a total of 24 co-authors. Here are the top co-authors by number of shared publications:
- Michael. I. Jordan
University of California Berkeley USA
7 shared publications - Shuzhong Zhang
Department of Industrial and Systems Engineering, University of Minnesota, Minneapolis, MN 55455 USA
2 shared publications - Zhengyuan Zhou
Stern School of Business, New York University, New York, New York 10012
2 shared publications Aldo Pacchiano
Department of Electrical Engineering and Computer Sciences Department of Statistics University of California Berkeley Microsoft Research NYC
1 shared publicationAsuman E. Ozdaglar
1 shared publication- Bo Jiang
Research Institute for Interdisciplinary Sciences, School of Information Management and Engineering, Shanghai University of Finance and Economics, Shanghai 200433, China; and Key Laboratory of Interdisciplinary Research of Computation and Economics (Shanghai University of Finance and Economics), Ministry of Education, Shanghai 200433, China; and Dishui Lake Advanced Finance Institute, Shanghai University of Finance and Economics, Shanghai 200120, China
1 shared publication
Latest Publications
- Open Access
- 30 May 2026
- Published by Elsevier in SSRN Electronic Journal
Robust Learning in Creator Markets: Algorithm and Optimal Regret Guarantee
We study robust learning in creator markets under the ε-contamination model. Here, a platform sequentially allocates advertiser demand to creators, collects rewards, and incentivizes them through revenue sharing. For each advertiser, rewards follow an unknown fixed distribution on at least T − εT rounds of a horizon of length T , while on the remaining at most εT rounds they can take arbitrary values in [0, 1]. Our first contribution is to formalize this robust learning problem in creator markets and demonstrate its algorithmic challenges. We prove that any learning algorithm suffers a regret of Ω(εT + T 2/3), and that the vanilla UCB update rule underlying Zhu et al. (2023) admits a linear-regret failure mode in an oracle-aided recommendation-only sub- problem even when εT < T 2/3. We then develop a new algorithm based on adaptive elimination, for which we establish a high-probability regret bound of ˜O(εT + T 2/3), optimal up to logarithmic factors. The algorithm is transparent and clarifies the key explore-exploit trade-off under contamination. Numerical experiments provide a focused stress test under front-loaded adversarial contamination, showing that the proposed elim- ination rule is more stable than vanilla UCB and Thompson sampling in this setting. These experiments complement the theoretical robustness guarantee rather than serve as an exhaustive empirical comparison across contamination models and baseline families. Using the outcomes induced by our algorithm, we also analyze how the distribution of creator capabilities and contamination shape long-term market concentration and revenue sharing, uncovering new insights that complement prior economic models (Bhargava 2022).
- 30 May 2026
- Published by Elsevier in SSRN Electronic Journal
Open AccessRobust Learning in Creator Markets: Algorithm and Optimal Regret Guarantee
We study robust learning in creator markets under the ε-contamination model. Here, a platform sequentially allocates advertiser demand to creators, collects rewards, and incentivizes them through revenue sharing. For each advertiser, rewards follow an unknown fixed distribution on at least T − εT rounds of a horizon of length T , while on the remaining at most εT rounds they can take arbitrary values in [0, 1]. Our first contribution is to formalize this robust learning problem in creator markets and demonstrate its algorithmic challenges. We prove that any learning algorithm suffers a regret of Ω(εT + T 2/3), and that the vanilla UCB update rule underlying Zhu et al. (2023) admits a linear-regret failure mode in an oracle-aided recommendation-only sub- problem even when εT < T 2/3. We then develop a new algorithm based on adaptive elimination, for which we establish a high-probability regret bound of ˜O(εT + T 2/3), optimal up to logarithmic factors. The algorithm is transparent and clarifies the key explore-exploit trade-off under contamination. Numerical experiments provide a focused stress test under front-loaded adversarial contamination, showing that the proposed elim- ination rule is more stable than vanilla UCB and Thompson sampling in this setting. These experiments complement the theoretical robustness guarantee rather than serve as an exhaustive empirical comparison across contamination models and baseline families. Using the outcomes induced by our algorithm, we also analyze how the distribution of creator capabilities and contamination shape long-term market concentration and revenue sharing, uncovering new insights that complement prior economic models (Bhargava 2022).
- Open Access
- 1 January 2026
- Published by Elsevier in SSRN Electronic Journal
How AI Aggregation Affects Knowledge
Artificial intelligence (AI) changes social learning when aggregated outputs become training data for future predictions. To study this, we extend the DeGroot model by introducing an AI aggregator that trains on population beliefs and feeds synthesized signals back to agents. We define the learning gap as the deviation of long-run beliefs from the efficient benchmark, allowing us to capture how AI aggregation affects learning. Our main result identifies a threshold in the speed of updating: when the aggregator updates too quickly, there is no positive-measure set of training weights that robustly improves learning across a broad class of environments, whereas such weights exist when updating is sufficiently slow. We then compare global and local architectures. Local aggregators trained on proximate or topic-specific data robustly improve learning in all environments. Consequently, replacing specialized local aggregators with a single global aggregator worsens learning in at least one dimension of the state.Institutional subscribers to the NBER working paper series, and residents of developing countries may download this paper without additional charge at www.nber.org.
- 1 January 2026
- Published by Elsevier in SSRN Electronic Journal
Open AccessHow AI Aggregation Affects Knowledge
Artificial intelligence (AI) changes social learning when aggregated outputs become training data for future predictions. To study this, we extend the DeGroot model by introducing an AI aggregator that trains on population beliefs and feeds synthesized signals back to agents. We define the learning gap as the deviation of long-run beliefs from the efficient benchmark, allowing us to capture how AI aggregation affects learning. Our main result identifies a threshold in the speed of updating: when the aggregator updates too quickly, there is no positive-measure set of training weights that robustly improves learning across a broad class of environments, whereas such weights exist when updating is sufficiently slow. We then compare global and local architectures. Local aggregators trained on proximate or topic-specific data robustly improve learning in all environments. Consequently, replacing specialized local aggregators with a single global aggregator worsens learning in at least one dimension of the state.Institutional subscribers to the NBER working paper series, and residents of developing countries may download this paper without additional charge at www.nber.org.
- Top Cited Papers
- 1 November 2025
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Operations Research
- Vol. 73 (6) , 3219-3244
Doubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback
Doubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback Curious about how players can learn and adapt in unknown games without knowing the game’s dynamics? In “Doubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback,” Ba, Lin, Zhang, and Zhou present a novel bandit learning algorithm for no-regret learning in games where each player only observes its reward determined by all players’ current joint action, not its gradient. Focusing on smooth and strongly monotone games, they introduce a bandit learning algorithm using self-concordant barrier functions. This algorithm achieves optimal single-agent regret and optimal last-iterate convergence rate in multiagent learning to the Nash equilibrium. Their work significantly improves previous methods and demonstrates the algorithm’s effectiveness through numerical results in various applications.
- 1 November 2025
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Operations Research
- Vol. 73 (6) , 3219-3244
Top Cited PapersDoubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback
Doubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback Curious about how players can learn and adapt in unknown games without knowing the game’s dynamics? In “Doubly Optimal No-Regret Online Learning in Strongly Monotone Games with Bandit Feedback,” Ba, Lin, Zhang, and Zhou present a novel bandit learning algorithm for no-regret learning in games where each player only observes its reward determined by all players’ current joint action, not its gradient. Focusing on smooth and strongly monotone games, they introduce a bandit learning algorithm using self-concordant barrier functions. This algorithm achieves optimal single-agent regret and optimal last-iterate convergence rate in multiagent learning to the Nash equilibrium. Their work significantly improves previous methods and demonstrates the algorithm’s effectiveness through numerical results in various applications.
- 1 May 2025
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Operations Research
- Vol. 73 (3) , 1675-1702
Adaptive, Doubly Optimal No-Regret Learning in Strongly Monotone and Exp-Concave Games with Gradient Feedback
Feasible Online Learning with Gradient Feedback Online gradient descent (OGD) is well-known to be doubly optimal under strong convexity or monotonicity assumptions: (1) in the single-agent setting, it achieves an optimal regret of [ ] for strongly convex cost functions, and (2) in the multiagent setting of strongly monotone games with each agent employing OGD we obtain last-iterate convergence of the joint action to a unique Nash equilibrium at an optimal rate of [ ]. Whereas these finite-time guarantees highlight its merits, OGD has the drawback that it requires knowing the strong convexity/monotonicity parameters. In “Adaptive, Doubly Optimal No-Regret Learning in Strongly Monotone and Exp-Concave Games with Gradient Feedback,” M. Jordan, T. Lin, and Z. Zhou design a fully adaptive OGD algorithm, AdaOGD, that does not require a priori knowledge of these parameters. In the single-agent setting, the algorithm achieves [ ] regret under strong convexity, which is optimal up to a log factor. Further, if each agent employs AdaOGD in strongly monotone games, the joint action converges in a last-iterate sense to a unique Nash equilibrium at a rate of [ ], again optimal up to log factors. The algorithms are illustrated in a learning version of the classic newsvendor problem, in which, because of lost sales, only (noisy) gradient feedback can be observed. The results immediately yield the first feasible and near-optimal algorithm for both the single-retailer and multiretailer settings.
- 1 May 2025
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Operations Research
- Vol. 73 (3) , 1675-1702
Adaptive, Doubly Optimal No-Regret Learning in Strongly Monotone and Exp-Concave Games with Gradient Feedback
Feasible Online Learning with Gradient Feedback Online gradient descent (OGD) is well-known to be doubly optimal under strong convexity or monotonicity assumptions: (1) in the single-agent setting, it achieves an optimal regret of [ ] for strongly convex cost functions, and (2) in the multiagent setting of strongly monotone games with each agent employing OGD we obtain last-iterate convergence of the joint action to a unique Nash equilibrium at an optimal rate of [ ]. Whereas these finite-time guarantees highlight its merits, OGD has the drawback that it requires knowing the strong convexity/monotonicity parameters. In “Adaptive, Doubly Optimal No-Regret Learning in Strongly Monotone and Exp-Concave Games with Gradient Feedback,” M. Jordan, T. Lin, and Z. Zhou design a fully adaptive OGD algorithm, AdaOGD, that does not require a priori knowledge of these parameters. In the single-agent setting, the algorithm achieves [ ] regret under strong convexity, which is optimal up to a log factor. Further, if each agent employs AdaOGD in strongly monotone games, the joint action converges in a last-iterate sense to a unique Nash equilibrium at a rate of [ ], again optimal up to log factors. The algorithms are illustrated in a learning version of the classic newsvendor problem, in which, because of lost sales, only (noisy) gradient feedback can be observed. The results immediately yield the first feasible and near-optimal algorithm for both the single-retailer and multiretailer settings.
- Open Access
Explicit Second-Order Min-Max Optimization Methods with Optimal Convergence Guarantee
We propose and analyze several inexact regularized Newton-type methods for finding a global saddle point of \emph{convex-concave} unconstrained min-max optimization problems. Compared to first-order methods, our understanding of second-order methods for min-max optimization is relatively limited, as obtaining global rates of convergence with second-order information is much more involved. In this paper, we examine how second-order information can be used to speed up extra-gradient methods, even under inexactness. Specifically, we show that the proposed methods generate iterates that remain within a bounded set and that the averaged iterates converge to an $\epsilon$-saddle point within $O(\epsilon^{-2/3})$ iterations in terms of a restricted gap function. This matched the theoretically established lower bound in this context. We also provide a simple routine for solving the subproblem at each iteration, requiring a single Schur decomposition and $O(\log\log(1/\epsilon))$ calls to a linear system solver in a quasi-upper-triangular system. Thus, our method improves the existing line-search-based second-order min-max optimization methods by shaving off an $O(\log\log(1/\epsilon))$ factor in the required number of Schur decompositions. Finally, we present numerical experiments on synthetic and real data that demonstrate the efficiency of the proposed methods.Open AccessExplicit Second-Order Min-Max Optimization Methods with Optimal Convergence Guarantee
We propose and analyze several inexact regularized Newton-type methods for finding a global saddle point of \emph{convex-concave} unconstrained min-max optimization problems. Compared to first-order methods, our understanding of second-order methods for min-max optimization is relatively limited, as obtaining global rates of convergence with second-order information is much more involved. In this paper, we examine how second-order information can be used to speed up extra-gradient methods, even under inexactness. Specifically, we show that the proposed methods generate iterates that remain within a bounded set and that the averaged iterates converge to an $\epsilon$-saddle point within $O(\epsilon^{-2/3})$ iterations in terms of a restricted gap function. This matched the theoretically established lower bound in this context. We also provide a simple routine for solving the subproblem at each iteration, requiring a single Schur decomposition and $O(\log\log(1/\epsilon))$ calls to a linear system solver in a quasi-upper-triangular system. Thus, our method improves the existing line-search-based second-order min-max optimization methods by shaving off an $O(\log\log(1/\epsilon))$ factor in the required number of Schur decompositions. Finally, we present numerical experiments on synthetic and real data that demonstrate the efficiency of the proposed methods. - 13 March 2024
- Published by Springer Nature in Mathematical Programming
- Vol. 209 (1-2) , 609-650
Perseus: a simple and optimal high-order method for variational inequalities
- 13 March 2024
- Published by Springer Nature in Mathematical Programming
- Vol. 209 (1-2) , 609-650
Perseus: a simple and optimal high-order method for variational inequalities
- Open Access
Perseus: A Simple and Optimal High-Order Method for Variational Inequalities
This paper settles an open and challenging question pertaining to the design of simple and optimal high-order methods for solving smooth and monotone variational inequalities (VIs). A VI involves finding $x^\star \in \mathcal{X}$ such that $\langle F(x), x - x^\star\rangle \geq 0$ for all $x \in \mathcal{X}$. We consider the setting in which $F$ is smooth with up to $(p-1)^{th}$-order derivatives. For $p = 2$, the cubic regularized Newton method was extended to VIs with a global rate of $O(\epsilon^{-1})$. An improved rate of $O(\epsilon^{-2/3}\log\log(1/\epsilon))$ can be obtained via an alternative second-order method, but this method requires a nontrivial line-search procedure as an inner loop. Similarly, high-order methods based on line-search procedures have been shown to achieve a rate of $O(\epsilon^{-2/(p+1)}\log\log(1/\epsilon))$. As emphasized by Nesterov, however, such procedures do not necessarily imply practical applicability in large-scale applications, and it would be desirable to complement these results with a simple high-order VI method that retains the optimality of the more complex methods. We propose a $p^{th}$-order method that does \textit{not} require any line search procedure and provably converges to a weak solution at a rate of $O(\epsilon^{-2/(p+1)})$. We prove that our $p^{th}$-order method is optimal in the monotone setting by establishing a matching lower bound under a generalized linear span assumption. Our method with restarting attains a linear rate for smooth and uniformly monotone VIs and a local superlinear rate for smooth and strongly monotone VIs. Our method also achieves a global rate of $O(\epsilon^{-2/p})$ for solving smooth and nonmonotone VIs satisfying the Minty condition and when augmented with restarting it attains a global linear and local superlinear rate for smooth and nonmonotone VIs satisfying the uniform/strong Minty condition.Open AccessPerseus: A Simple and Optimal High-Order Method for Variational Inequalities
This paper settles an open and challenging question pertaining to the design of simple and optimal high-order methods for solving smooth and monotone variational inequalities (VIs). A VI involves finding $x^\star \in \mathcal{X}$ such that $\langle F(x), x - x^\star\rangle \geq 0$ for all $x \in \mathcal{X}$. We consider the setting in which $F$ is smooth with up to $(p-1)^{th}$-order derivatives. For $p = 2$, the cubic regularized Newton method was extended to VIs with a global rate of $O(\epsilon^{-1})$. An improved rate of $O(\epsilon^{-2/3}\log\log(1/\epsilon))$ can be obtained via an alternative second-order method, but this method requires a nontrivial line-search procedure as an inner loop. Similarly, high-order methods based on line-search procedures have been shown to achieve a rate of $O(\epsilon^{-2/(p+1)}\log\log(1/\epsilon))$. As emphasized by Nesterov, however, such procedures do not necessarily imply practical applicability in large-scale applications, and it would be desirable to complement these results with a simple high-order VI method that retains the optimality of the more complex methods. We propose a $p^{th}$-order method that does \textit{not} require any line search procedure and provably converges to a weak solution at a rate of $O(\epsilon^{-2/(p+1)})$. We prove that our $p^{th}$-order method is optimal in the monotone setting by establishing a matching lower bound under a generalized linear span assumption. Our method with restarting attains a linear rate for smooth and uniformly monotone VIs and a local superlinear rate for smooth and strongly monotone VIs. Our method also achieves a global rate of $O(\epsilon^{-2/p})$ for solving smooth and nonmonotone VIs satisfying the Minty condition and when augmented with restarting it attains a global linear and local superlinear rate for smooth and nonmonotone VIs satisfying the uniform/strong Minty condition. - 4 January 2023
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Mathematics of Operations Research
- Vol. 48 (4) , 2353-2382
Monotone Inclusions, Acceleration, and Closed-Loop Control
We propose and analyze a new dynamical system with a closed-loop control law in a Hilbert space [ ], aiming to shed light on the acceleration phenomenon for monotone inclusion problems, which unifies a broad class of optimization, saddle point, and variational inequality (VI) problems under a single framework. Given an operator [ ] that is maximal monotone, we propose a closed-loop control system that is governed by the operator [ ], where a feedback law [ ] is tuned by the resolution of the algebraic equation [ ] for some [ ]. Our first contribution is to prove the existence and uniqueness of a global solution via the Cauchy–Lipschitz theorem. We present a simple Lyapunov function for establishing the weak convergence of trajectories via the Opial lemma and strong convergence results under additional conditions. We then prove a global ergodic convergence rate of [ ] in terms of a gap function and a global pointwise convergence rate of [ ] in terms of a residue function. Local linear convergence is established in terms of a distance function under an error bound condition. Further, we provide an algorithmic framework based on the implicit discretization of our system in a Euclidean setting, generalizing the large-step hybrid proximal extragradient framework. Even though the discrete-time analysis is a simplification and generalization of existing analyses for a bounded domain, it is largely motivated by the aforementioned continuous-time analysis, illustrating the fundamental role that the closed-loop control plays in acceleration in monotone inclusion. A highlight of our analysis is a new result concerning [ ]-order tensor algorithms for monotone inclusion problems, complementing the recent analysis for saddle point and VI problems. Funding: This work was supported in part by the Mathematical Data Science Program of the Office of Naval Research [Grant N00014-18-1-2764] and by the Vannevar Bush Faculty Fellowship Program [Grant N00014-21-1-2941].
- 4 January 2023
- Published by Institute for Operations Research and the Management Sciences (INFORMS) in Mathematics of Operations Research
- Vol. 48 (4) , 2353-2382
Monotone Inclusions, Acceleration, and Closed-Loop Control
We propose and analyze a new dynamical system with a closed-loop control law in a Hilbert space [ ], aiming to shed light on the acceleration phenomenon for monotone inclusion problems, which unifies a broad class of optimization, saddle point, and variational inequality (VI) problems under a single framework. Given an operator [ ] that is maximal monotone, we propose a closed-loop control system that is governed by the operator [ ], where a feedback law [ ] is tuned by the resolution of the algebraic equation [ ] for some [ ]. Our first contribution is to prove the existence and uniqueness of a global solution via the Cauchy–Lipschitz theorem. We present a simple Lyapunov function for establishing the weak convergence of trajectories via the Opial lemma and strong convergence results under additional conditions. We then prove a global ergodic convergence rate of [ ] in terms of a gap function and a global pointwise convergence rate of [ ] in terms of a residue function. Local linear convergence is established in terms of a distance function under an error bound condition. Further, we provide an algorithmic framework based on the implicit discretization of our system in a Euclidean setting, generalizing the large-step hybrid proximal extragradient framework. Even though the discrete-time analysis is a simplification and generalization of existing analyses for a bounded domain, it is largely motivated by the aforementioned continuous-time analysis, illustrating the fundamental role that the closed-loop control plays in acceleration in monotone inclusion. A highlight of our analysis is a new result concerning [ ]-order tensor algorithms for monotone inclusion problems, complementing the recent analysis for saddle point and VI problems. Funding: This work was supported in part by the Mathematical Data Science Program of the Office of Naval Research [Grant N00014-18-1-2764] and by the Vannevar Bush Faculty Fellowship Program [Grant N00014-21-1-2941].
- Open Access
Online Nonsubmodular Minimization with Delayed Costs: From Full Information to Bandit Feedback
Motivated by applications to online learning in sparse estimation and Bayesian optimization, we consider the problem of online unconstrained nonsubmodular minimization with delayed costs in both full information and bandit feedback settings. In contrast to previous works on online unconstrained submodular minimization, we focus on a class of nonsubmodular functions with special structure, and prove regret guarantees for several variants of the online and approximate online bandit gradient descent algorithms in static and delayed scenarios. We derive bounds for the agent's regret in the full information and bandit feedback setting, even if the delay between choosing a decision and receiving the incurred cost is unbounded. Key to our approach is the notion of $(\alpha, \beta)$-regret and the extension of the generic convex relaxation model from~\citet{El-2020-Optimal}, the analysis of which is of independent interest. We conduct and showcase several simulation studies to demonstrate the efficacy of our algorithms.
Open AccessOnline Nonsubmodular Minimization with Delayed Costs: From Full Information to Bandit Feedback
Motivated by applications to online learning in sparse estimation and Bayesian optimization, we consider the problem of online unconstrained nonsubmodular minimization with delayed costs in both full information and bandit feedback settings. In contrast to previous works on online unconstrained submodular minimization, we focus on a class of nonsubmodular functions with special structure, and prove regret guarantees for several variants of the online and approximate online bandit gradient descent algorithms in static and delayed scenarios. We derive bounds for the agent's regret in the full information and bandit feedback setting, even if the delay between choosing a decision and receiving the incurred cost is unbounded. Key to our approach is the notion of $(\alpha, \beta)$-regret and the extension of the generic convex relaxation model from~\citet{El-2020-Optimal}, the analysis of which is of independent interest. We conduct and showcase several simulation studies to demonstrate the efficacy of our algorithms.
- Open Access
- 22 October 2021
- Published by Springer Nature in Mathematical Programming
- Vol. 195 (1-2) , 929-975
A control-theoretic perspective on optimal high-order optimization
We provide a control-theoretic perspective on optimal tensor algorithms for minimizing a convex function in a finite-dimensional Euclidean space. Given a function $$\varPhi : {\mathbb {R}}^d \rightarrow {\mathbb {R}}$$ that is convex and twice continuously differentiable, we study a closed-loop control system that is governed by the operators $$\nabla \varPhi $$ and $$\nabla ^2 \varPhi $$ together with a feedback control law $$\lambda (\cdot )$$ satisfying the algebraic equation $$(\lambda (t))^p\Vert \nabla \varPhi (x(t))\Vert ^{p-1} = \theta $$ for some $$\theta \in (0, 1)$$ . Our first contribution is to prove the existence and uniqueness of a local solution to this system via the Banach fixed-point theorem. We present a simple yet nontrivial Lyapunov function that allows us to establish the existence and uniqueness of a global solution under certain regularity conditions and analyze the convergence properties of trajectories. The rate of convergence is $$O(1/t^{(3p+1)/2})$$ in terms of objective function gap and $$O(1/t^{3p})$$ in terms of squared gradient norm. Our second contribution is to provide two algorithmic frameworks obtained from discretization of our continuous-time system, one of which generalizes the large-step A-HPE framework of Monteiro and Svaiter (SIAM J Optim 23(2):1092–1125, 2013) and the other of which leads to a new optimalp-th order tensor algorithm. While our discrete-time analysis can be seen as a simplification and generalization of Monteiro and Svaiter (2013), it is largely motivated by the aforementioned continuous-time analysis, demonstrating the fundamental role that the feedback control plays in optimal acceleration and the clear advantage that the continuous-time perspective brings to algorithmic design. A highlight of our analysis is that we show that all of thep-th order optimal tensor algorithms that we discuss minimize the squared gradient norm at a rate of $$O(k^{-3p})$$ , which complements the recent analysis in Gasnikov et al. (in: COLT, PMLR, pp 1374–1391, 2019), Jiang et al. (in: COLT, PMLR, pp 1799–1801, 2019) and Bubeck et al. (in: COLT, PMLR, pp 492–507, 2019).- 22 October 2021
- Published by Springer Nature in Mathematical Programming
- Vol. 195 (1-2) , 929-975
Open AccessA control-theoretic perspective on optimal high-order optimization
We provide a control-theoretic perspective on optimal tensor algorithms for minimizing a convex function in a finite-dimensional Euclidean space. Given a function $$\varPhi : {\mathbb {R}}^d \rightarrow {\mathbb {R}}$$ that is convex and twice continuously differentiable, we study a closed-loop control system that is governed by the operators $$\nabla \varPhi $$ and $$\nabla ^2 \varPhi $$ together with a feedback control law $$\lambda (\cdot )$$ satisfying the algebraic equation $$(\lambda (t))^p\Vert \nabla \varPhi (x(t))\Vert ^{p-1} = \theta $$ for some $$\theta \in (0, 1)$$ . Our first contribution is to prove the existence and uniqueness of a local solution to this system via the Banach fixed-point theorem. We present a simple yet nontrivial Lyapunov function that allows us to establish the existence and uniqueness of a global solution under certain regularity conditions and analyze the convergence properties of trajectories. The rate of convergence is $$O(1/t^{(3p+1)/2})$$ in terms of objective function gap and $$O(1/t^{3p})$$ in terms of squared gradient norm. Our second contribution is to provide two algorithmic frameworks obtained from discretization of our continuous-time system, one of which generalizes the large-step A-HPE framework of Monteiro and Svaiter (SIAM J Optim 23(2):1092–1125, 2013) and the other of which leads to a new optimalp-th order tensor algorithm. While our discrete-time analysis can be seen as a simplification and generalization of Monteiro and Svaiter (2013), it is largely motivated by the aforementioned continuous-time analysis, demonstrating the fundamental role that the feedback control plays in optimal acceleration and the clear advantage that the continuous-time perspective brings to algorithmic design. A highlight of our analysis is that we show that all of thep-th order optimal tensor algorithms that we discuss minimize the squared gradient norm at a rate of $$O(k^{-3p})$$ , which complements the recent analysis in Gasnikov et al. (in: COLT, PMLR, pp 1374–1391, 2019), Jiang et al. (in: COLT, PMLR, pp 1799–1801, 2019) and Bubeck et al. (in: COLT, PMLR, pp 492–507, 2019).