← Back to Research

Structured convex bodies + applications

Overview

My research develops geometric frameworks for understanding structured convexity with applications in optimization and statistics. This includes the study of classical objects like zonotopes and polytopes in high-dimensions and generalizations to operator-theoretic settings.

Structured Convex Bodies

In joint work with Venkat Chandarsekaran, we introduce operatopes, which are Minkowski sums of affine images of operator norm balls—a natural generalization of zonotopes (Minkowski sums of line segments). Taking limits yields operanoids, defined as expectations of random affine images of operator norm balls. These convex bodies generalize limits of zonotopes, called zonoids, which have applications in statistics, stochastic processes, engineering, and the approximation capabilities of neural networks. This work also opens a new asymptotic regimes when the dimension of the matrices become large, leading to the study of noncommutative zonoids using free probability theory.

Operatope structure Operatope structure Operatope structure
Figure 2: Operatopes as generalizations of zonotopes

In joint work with Gilles Bonnet, we study facets of convex hulls of random vectors sampled uniformly from the unit sphere in high dimensions. We establish asymptotic formulas for the height of facets and expected number of facets as dimension grows, identifying different regimes with distinct asymptotic behavior.

Spectrahedral Regression

Convex regression—fitting a convex function to input-output data—is fundamental in statistics and machine learning. We introduce spectrahedral regression, which fits spectrahedral functions (maximum eigenvalue of an affine matrix expression) to data.

This significantly generalizes polyhedral (max-affine) regression by allowing the fitted function to capture both smooth and singular features. We prove approximation bounds showing how well spectrahedral functions approximate arbitrary convex functions via statistical risk analysis. We also develop an alternating minimization algorithm that converges geometrically with high probability to near-optimal solutions given good initialization.

Spectrahedral function
Figure 3: Spectrahedral functions generalize piecewise linear convex functions

Related Publications

Operatopes, operanoids, and noncommutative zonoids
Eliza O'Reilly, Venkat Chandrasekaran
Preprint, 2026
Spectrahedral regression
Eliza O'Reilly, Venkat Chandrasekaran
SIAM Journal on Optimization, 33(2), 2023
Facets of spherical random polytopes
Gilles Bonnet, Eliza O'Reilly
Mathematische Nachrichten, 295(10):1901–1933, 2022