Source author record

Arun Jambulapati

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

6works
3topics
4close collaborators

Published work

6 publications

preprint2026arXiv

Convex optimization with pp-norm oracles

In recent years, there have been significant advances in efficiently solving ℓs\ell_s-regression using linear system solvers and ℓ2\ell_2-regression [Adil-Kyng-Peng-Sachdeva, J. ACM'24]. Would efficient smoothed ℓp\ell_p-norm solvers lead to even faster rates for solving ℓs\ell_s-regression when 2≤p<s2 \leq p < s? In this paper, we give an affirmative answer to this question and show how to solve ℓs\ell_s-regression using O~(nν1+ν)\tilde{O}(n^{\fracν{1+ν}}) iterations of solving smoothed ℓp\ell_p regression problems, where ν:=1p−1sν:= \frac{1}{p} - \frac{1}{s}. To obtain this result, we provide improved accelerated rates for convex optimization problems when given access to an ℓps(λ)\ell_p^s(λ)-proximal oracle, which, for a point cc, returns the solution of the regularized problem min⁡xf(x)+λ∣∣x−c∣∣ps\min_{x} f(x) + λ||x-c||_p^s. Additionally, we show that these rates for the ℓps(λ)\ell_p^s(λ)-proximal oracle are optimal for algorithms that query in the span of the outputs of the oracle, and we further apply our techniques to settings of high-order and quasi-self-concordant optimization.

preprint2022arXiv

Optimal Methods for Higher-Order Smooth Monotone Variational Inequalities

In this work, we present new simple and optimal algorithms for solving the variational inequality (VI) problem for pthp^{th}-order smooth, monotone operators -- a problem that generalizes convex optimization and saddle-point problems. Recent works (Bullins and Lai (2020), Lin and Jordan (2021), Jiang and Mokhtari (2022)) present methods that achieve a rate of O~(ε−2/(p+1))\tilde{O}(ε^{-2/(p+1)}) for p≥1p\geq 1, extending results by (Nemirovski (2004)) and (Monteiro and Svaiter (2012)) for p=1,2p=1,2. A drawback to these approaches, however, is their reliance on a line search scheme. We provide the first pthp^{\textrm{th}}-order method that achieves a rate of O(ε−2/(p+1)).O(ε^{-2/(p+1)}). Our method does not rely on a line search routine, thereby improving upon previous rates by a logarithmic factor. Building on the Mirror Prox method of Nemirovski (2004), our algorithm works even in the constrained, non-Euclidean setting. Furthermore, we prove the optimality of our algorithm by constructing matching lower bounds. These are the first lower bounds for smooth MVIs beyond convex optimization for p>1p > 1. This establishes a separation between solving smooth MVIs and smooth convex optimization, and settles the oracle complexity of solving pthp^{\textrm{th}}-order smooth MVIs.

preprint2022arXiv

RECAPP: Crafting a More Efficient Catalyst for Convex Optimization

The accelerated proximal point algorithm (APPA), also known as "Catalyst", is a well-established reduction from convex optimization to approximate proximal point computation (i.e., regularized minimization). This reduction is conceptually elegant and yields strong convergence rate guarantees. However, these rates feature an extraneous logarithmic term arising from the need to compute each proximal point to high accuracy. In this work, we propose a novel Relaxed Error Criterion for Accelerated Proximal Point (RECAPP) that eliminates the need for high accuracy subproblem solutions. We apply RECAPP to two canonical problems: finite-sum and max-structured minimization. For finite-sum problems, we match the best known complexity, previously obtained by carefully-designed problem-specific algorithms. For minimizing max⁡yf(x,y)\max_y f(x,y) where ff is convex in xx and strongly-concave in yy, we improve on the best known (Catalyst-based) bound by a logarithmic factor.

preprint2020arXiv

Acceleration with a Ball Optimization Oracle

Consider an oracle which takes a point xx and returns the minimizer of a convex function ff in an ℓ2\ell_2 ball of radius rr around xx. It is straightforward to show that roughly r−1log⁡1εr^{-1}\log\frac{1}ε calls to the oracle suffice to find an εε-approximate minimizer of ff in an ℓ2\ell_2 unit ball. Perhaps surprisingly, this is not optimal: we design an accelerated algorithm which attains an εε-approximate minimizer with roughly r−2/3log⁡1εr^{-2/3} \log \frac{1}ε oracle queries, and give a matching lower bound. Further, we implement ball optimization oracles for functions with locally stable Hessians using a variant of Newton's method. The resulting algorithm applies to a number of problems of practical and theoretical import, improving upon previous results for logistic and ℓ∞\ell_\infty regression and achieving guarantees comparable to the state-of-the-art for ℓp\ell_p regression.

preprint2020arXiv

Robust Sub-Gaussian Principal Component Analysis and Width-Independent Schatten Packing

We develop two methods for the following fundamental statistical task: given an εε-corrupted set of nn samples from a dd-dimensional sub-Gaussian distribution, return an approximate top eigenvector of the covariance matrix. Our first robust PCA algorithm runs in polynomial time, returns a 1−O(εlog⁡ε−1)1 - O(ε\logε^{-1})-approximate top eigenvector, and is based on a simple iterative filtering approach. Our second, which attains a slightly worse approximation factor, runs in nearly-linear time and sample complexity under a mild spectral gap assumption. These are the first polynomial-time algorithms yielding non-trivial information about the covariance of a corrupted sub-Gaussian distribution without requiring additional algebraic structure of moments. As a key technical tool, we develop the first width-independent solvers for Schatten-pp norm packing semidefinite programs, giving a (1+ε)(1 + ε)-approximate solution in O(plog⁡(ndε)ε−1)O(p\log(\tfrac{nd}ε)ε^{-1}) input-sparsity time iterations (where nn, dd are problem dimensions).

preprint2020arXiv

Towards Optimal Running Times for Optimal Transport

In this work, we provide faster algorithms for approximating the optimal transport distance, e.g. earth mover's distance, between two discrete probability distributions μ,ν∈Δnμ, ν\in Δ^n. Given a cost function C:[n]×[n]→R≥0C : [n] \times [n] \to \mathbb{R}_{\geq 0} where C(i,j)≤1C(i,j) \leq 1 quantifies the penalty of transporting a unit of mass from ii to jj, we show how to compute a coupling XX between rr and cc in time O~(n2/ε)\widetilde{O}\left(n^2 /ε\right) whose expected transportation cost is within an additive εε of optimal. This improves upon the previously best known running time for this problem of O~(min{n9/4/ε,n2/ε2})\widetilde{O}\left(\text{min}\left\{ n^{9/4}/ε, n^2/ε^2 \right\}\right). We achieve our results by providing reductions from optimal transport to canonical optimization problems for which recent algorithmic efforts have provided nearly-linear time algorithms. Leveraging nearly linear time algorithms for solving packing linear programs and for solving the matrix balancing problem, we obtain two separate proofs of our stated running time. Further, one of our algorithms is easily parallelized and can be implemented with depth O~(1/ε)\widetilde{O}(1/ε). Moreover, we show that further algorithmic improvements to our result would be surprising in the sense that any improvement would yield an o(n2.5)o(n^{2.5}) algorithm for \textit{maximum cardinality bipartite matching}, for which currently the only known algorithms for achieving such a result are based on fast-matrix multiplication.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.