← Search

Robert M. Gower

16 accepted papers

2026

The Polar Express: Optimal Matrix Sign Methods and their Application to the Muon Algorithm

ICLR 2026oral

Computing the polar decomposition and the related matrix sign function has been a well-studied problem in numerical analysis for decades. Recently, it has emerged as an important subroutine within the Muon algorithm for training deep neural networks. However, the requirements of this application di…

Cited by 0SourcecodeScholar
2025

Fisher meets Feynman: score-based variational inference with a product of experts

NeurIPS 2025spotlight

We introduce a highly expressive yet distinctly tractable family for black-box variational inference (BBVI). Each member of this family is a weighted product of experts (PoE), and each weighted expert in the product is proportional to a multivariate $t$-distribution. These products of experts can…

Cited by 0SourceScholar
2024

Batch and match: black-box variational inference with a score-based divergence

ICML 2024spotlight

Most leading implementations of black-box variational inference (BBVI) are based on optimizing a stochastic evidence lower bound (ELBO). But such approaches to BBVI often converge slowly due to the high variance of their gradient estimates and their sensitivity to hyperparameters. In this work, we p…

Cited by 7SourcePDFScholar
2024

Directional Smoothness and Gradient Methods: Convergence and Adaptivity

NeurIPS 2024poster

We develop new sub-optimality bounds for gradient descent (GD) that depend on the conditioning of the objective along the path of optimization, rather than on global, worst-case constants. Key to our proofs is directional smoothness, a measure of gradient variation that we use to develop upper-boun…

Cited by 5SourcePDFScholar
2024

EigenVI: score-based variational inference with orthogonal function expansions

NeurIPS 2024spotlight

We develop EigenVI, an eigenvalue-based approach for black-box variational inference (BBVI). EigenVI constructs its variational approximations from orthogonal function expansions. For distributions over $\mathbb{R}^D$, the lowest order term in these expansions provides a Gaussian variational approxi…

Cited by 3SourcePDFScholar
2024

Improving Convergence and Generalization Using Parameter Symmetries

ICLR 2024oral

In many neural networks, different values of the parameters may result in the same loss value. Parameter space symmetries are loss-invariant transformations that change the model parameters. Teleportation applies such transformations to accelerate optimization. However, the exact mechanism behind th…

2024

MoMo: Momentum Models for Adaptive Learning Rates

ICML 2024poster

Training a modern machine learning architecture on a new task requires extensive learning-rate tuning, which comes at a high computational cost. Here we develop new Polyak-type adaptive learning rates that can be used on top of any momentum method, and require less tuning to perform well. We first d…

2023

Linear Convergence of Natural Policy Gradient Methods with Log-Linear Policies

ICLR 2023poster

We consider infinite-horizon discounted Markov decision processes and study the convergence rates of the natural policy gradient (NPG) and the Q-NPG methods with the log-linear policy class. Using the compatible function approximation framework, both methods with log-linear policies can be written a…

Cited by 51SourcePDFScholar
2023

Provable convergence guarantees for black-box variational inference

NeurIPS 2023poster

Black-box variational inference is widely used in situations where there is no proof that its stochastic optimization succeeds. We suggest this is due to a theoretical gap in existing stochastic optimization proofs—namely the challenge of gradient estimators with unusual noise bounds, and a composit…

Cited by 26SourcePDFScholar
2023

SP2 : A Second Order Stochastic Polyak Method

ICLR 2023poster

Recently the SP (Stochastic Polyak step size) method has emerged as a competitive adaptive method for setting the step sizes of SGD. SP can be interpreted as a method specialized to interpolated models, since it solves the interpolation equations. SP solves these equation by using local linearizati…

Cited by 13SourcePDFScholar
2023

Variational Inference with Gaussian Score Matching

NeurIPS 2023poster

Variational inference (VI) is a method to approximate the computationally intractable posterior distributions that arise in Bayesian statistics. Typically, VI fits a simple parametric distribution to be close to the target posterior, optimizing an appropriate objective such as the evidence lower b…

Cited by 13SourcePDFScholar
2022

A general sample complexity analysis of vanilla policy gradient

AISTATS 2022poster

We adapt recent tools developed for the analysis of Stochastic Gradient Descent (SGD) in non-convex optimization to obtain convergence and sample complexity guarantees for the vanilla policy gradient (PG). Our only assumptions are that the expected return is smooth w.r.t. the policy parameters, that…

Cited by 83SourcePDFScholar
2022

SAN: Stochastic Average Newton Algorithm for Minimizing Finite Sums

AISTATS 2022poster

We present a principled approach for designing stochastic Newton methods for solving finite sum optimization problems. Our approach has two steps. First, we re-write the stationarity conditions as a system of nonlinear equations that associates each data point to a new row. Second, we apply a Subsam…