Current students

Speaker: David Evangelista
Supervisor(s): Joseph Cheriyan and Sophie Spirkl
Committee: Jane Gao, Eric Blais
Location: MC 5417

Abstract: 

A tournament $\T=(V,A)$ on $n$ vertices is an orientation of the complete graph $K_n$. The backedge graph of $T$ with respect to an ordering of $V$ is the undirected graph on vertex set $V$ whose edge set corresponds to the arcs directed from a later vertex to an earlier vertex in the ordering. Backedge graphs provide concise representations of the tournament. The algorithmic problem of determining whether a tournament admits a backedge graph in a given class of undirected graphs varies in complexity, and is often equivalent to computing parameters of tournaments, such as degreewidth when the backedge graph has bounded maximum degree \cite{Davot et al., 2023}. We extend the notion of degreewidth by introducing directional degreewidth, which separately bounds the left-degrees and right-degrees of vertices in addition to bounding the total degrees. We obtain an algorithm for verifying bounds on the directional degreewidth of the tournament, whose runtime is polynomial time when the total degree is unbounded, or fixed-parameter tractable time with respect to the total degree bound otherwise. We also provide a polynomial-time algorithm for computing a $P_3$-free backedge graph of a tournament, if it exists. Together with existing results, the latter result settles the complexity of determining whether a tournament admits an $H$-free backedge graph when $H$ is any graph on three vertices.

Speaker:

Leo Jung

Location: MC 5029

Abstract: 

Difficulties in solving large-scale optimization problems often arise from structural pathologies such as ill-conditioning, Hadamard ill-posedness, and degeneracy, particularly due to the failure of constraint qualifications. While standard algorithms often struggle to address these issues, preprocessing based on structural analysis offers an effective strategy for overcoming such challenges. This thesis investigates several preprocessing methods targeting various sources of these difficulties.
In Part I, we study a nonclassical, average condition number of linear systems, the $\omega$-condition number. Our results demonstrate several advantages of the $\omega$-condition number over the classical $\kappa$-condition number. First, $\omega$ provides a more accurate measure of the conditioning of linear systems by more faithfully capturing the effects of perturbations observed in practice. Second, $\omega$ exhibits superior numerical stability compared to $\kappa$. Third, when used in preconditioner design, $\omega$ more effectively promotes eigenvalue clustering, which is crucial for the efficiency of iterative solvers. Finally, the analytical simplicity of $\omega$ enables the derivation of explicit optimality conditions, allowing for closed-form expressions of optimal preconditioners under various frameworks, including low rank updates of the generalized Jacobian for semismooth Newton methods and diagonal or block-diagonal scaling. For diagonal preconditioning, we further include a comparison between two distinct notions of conditioning.
In Part II, we first answer in the affirmative a long-standing open question of whether the smooth stress function admits local nonglobal minimizers. This quartic nonconvex objective function arises in the exact recovery of a Euclidean distance matrix (EDM) of a given embedding dimension. By eliminating the Hadamard ill-posedness caused by translation and rotation invariance, we stabilize Newton's method and avoid singular Hessians. \\
We then consider the single-element error correction problem as a case study. We first show that the standard nearest EDM formulation based on minimizing the smooth stress function fails to recover the correct EDM in this setting. We then introduce divide-and-conquer strategies based on facial reduction. Our approach efficiently recovers the correct EDM with high accuracy, and we further provide criteria characterizing the existence of multiple solutions.
In Part III, we relate FR to the analysis of the convergence behaviour of a semismooth Newton method for projection onto a spectrahedron, i.e., the intersection of a linear manifold and the semidefinite cone. In this process, we derive an explicit formula for the projection onto a face of the semidefinite cone obtained via regularization and analyze pathologies that arise in the absence of strict feasibility. We further show that ill-conditioning of the Jacobian near optimality characterizes the degeneracy of the projection point. \\
As an application, we consider a simplified Wasserstein barycenter problem, a well-known NP-hard problem. We compute the Wasserstein barycenter by exploiting the structure of the linear constraints to obtain a facially reduced doubly nonnegative (DNN) relaxation. This reduction provides a natural splitting for applying the symmetric alternating direction method of multipliers (sADMM). The resulting algorithm exploits structure in the subproblems to compute strong upper and lower bounds. In most of the instances, we achieve the small gap between these bounds, which means that the original problem is solved.

Speaker:

David Aleman
Affiliation: University of Waterloo
Location: MC 6029

Abstract: 

The multicommodity flow problem involves routing multiple distinct commodities through a shared network. An instance is given by an undirected graph G=(V, E(G) ) with edge capacities, and a collection of source-sink pairs (s_i,t_i) in V with associated nonnegative demands d(s_i, t_i). It will be convenient to think of the source-sink pairs as forming the edges of a demand graph H=( V, E(H) ). A flow is feasible if it routes all demands without exceeding the edge capacities, and it is unsplittable if it routes each demand along a single path. Let C be the smallest value such that the existence of a feasible flow implies the existence of an unsplittable flow that exceeds the edge capacities by at most an additivie amount of C times the maximum demand value. 
We show that if G+H = (V, E(G) U E(H) ) is planar, then  1.5<= C <= 2.
Joint work with Kumar, Poremba, and Shepherd. 
 
Friday, July 31, 2026 11:30 am - 12:30 pm EDT (GMT -04:00)

CombOpt ReadingGroup - Sina Kalantarzadeh-Minimum Bounded Degree Spanning Trees

Speaker:

Sina Kalantarzadeh
Affiliation: University of Waterloo
Location: MC 6029

Abstract: I will present Goemans’s beautiful, though no longer state-of-the-art, 2006 result on the Minimum Bounded-Degree Spanning Tree problem. Given a weighted graph (G=(V,E)) and degree bounds (B\to\mathbb{N}), the goal is to find a minimum-cost spanning tree (T) satisfying (d_T(v)\le B(v)) for every (v\in V). This problem is NP-hard. Fürer and Raghavachari (1992) gave a polynomial-time algorithm for the unweighted setting that produces a spanning tree satisfying (d_T(v)\le B(v)+1). For the weighted problem, Goemans designed an elegant LP-rounding algorithm that returns a tree of cost at most that of the optimal degree-bounded solution while satisfying (d_T(v)\le B(v)+2). I will explain this result and the simple yet beautiful combinatorial optimization ideas underlying it. In 2007, Singh and Lau improved the violation to (B(v)+1) using iterative relaxation, which I might give a talk about in later sessions.

Wednesday, August 5, 2026 2:00 pm - 3:00 pm EDT (GMT -04:00)

Master's Thesis Presentation - Martin Liu

Speaker: Martin Li
Supervisor(s): Vijay Bhattiprolu
Comittee: Jonathan Leake, Levent Tuncel
Location: MC 6483

Abstract:

The $d$-dimensional Grothendieck constant is the smallest constant $K$ such that \begin{align*} \sup\left\{\sum_{i,j=1}^n A_{ij}\langle u_i,v_j\rangle:u_i,v_j\in S^{d-1}\right\}\le K\cdot\sup\left\{\sum_{i,j=1}^n A_{ij}x_iy_j:x_i,y_j\in\{-1,1\}\right\} \end{align*}for any $n\in\mathbb{N}$ and any real $n\times n$ matrix $A$. The inequality above, called the Grothendieck inequality, has made a deep impact in a variety of areas such as functional analysis, quantum information theory, and optimization. Determining the $d$-dimensional Grothendieck constant for any $d\ge 3$ is a long-standing open problem.


In this paper, we propose a worst operator in dimension 3, whose $\infty\to 1$ norm is conjectured to be $1/K_G(3)$. We study a related class of operators with nice geometric interpretations, and we prove the function $f:S^{d-1}\to\{-1,1\}$ corresponding to a hyperplane is uniquely optimal for this class, with the isoperimetric inequality lying at the heart of our proof.

Speaker: David Torregrosa Belén
Affiliation: University of Alicante
Location: MC 5501

Abstract:This talk presents a proximal stochastic subgradient method for minimizing the sum of an expected cost and a lower semicontinuous, prox-bounded function. We target a broad class of nonconvex integrands obeying a nonsmooth, localized variant of the descent lemma in the decision variable, which in particular covers smooth losses with Lipschitz gradient. At each iteration, the expected cost is replaced by a sample average that is progressively refined, and the proximal stepsize is selected by an Armijo-type line search enforcing a
sufficient decrease property up to stochastic errors induced by the sample-based approximation. This framework accommodates more general problem formulations than existing methods and our analysis yields convergence guarantees that, to the best of our knowledge, are new even in the smooth setting. Specifically, we establish almost sure convergence of the sequence of function values and stationarity of every accumulation point of the trajectories under the relaxed requirement that the sample-size sequence be merely nondecreasing and unbounded. Leveraging the Kurdyka-Lojasiewicz property, we further proof convergence of the whole trajectory to a single stationary point. Finally, for exponential-type desingularizing functions and polynomially growing sample sizes, we derive explicit polynomial convergence rates, up to logarithmic factor, for both the function values and the iterates. This is a joint work with Felipe Atenas, Pedro Pérez-Aros and Alejandro
Jofré, from the University of Chile.

Speaker: Ting Kei Pong
Affiliation: The Hong Kong Polytechnic University
Location: MC 5501

Abstract:We consider the problem of minimizing the sum of a Lipschitz differentiable convex function and a proper closed convex function that admits efficient linear minimization oracles, subject to multiple smooth convex inequality constraints. We adapt the classical augmented Lagrangian (AL) method for these problems: in each iteration, our algorithm consists of one step of conditional gradient (CG) method applied to the AL function, followed by an update of the dual variable as in classical AL methods with a diminishing dual stepsize. We study the convergence rate of our algorithm under two standard stepsize rules for the CG method, namely, an open-loop stepsize and the short stepsize, and obtain a rate that matches the best-known complexity for this class of problems. We also establish accelerated rates when the aforementioned proper closed convex function is the indicator function of a uniformly convex set. This is a joint work with Xiaozhou Wang and Zev Woodstock.

Speaker: Maryam Yekta
Affiliation: University of Waterloo
Location: MC 5479

AbstractThe theory of log concave polynomials has recently been developed to study objects and problems in combinatorics and other subfields in mathematics. Particular classes of log concave polynomials called Lorentzian polynomials and denormalized and dually Lorentzian polynomials have been used to prove log concavity statements for various combinatorial sequences. This includes the strongest form of Mason's log concavity conjecture on the independent sets of matroids and the log concavity of sequences of Kostka numbers.
In this talk, we develop an analogous class of power series called denormalized Lorentzian (DL) Laurent series. This class is the natural generalization of DL polynomials to homogeneous power series with the benefit of capturing a number of combinatorial generating series including the Kostant partition function for integer flows of directed graphs. We then analyze specific DL Laurent series to obtain new bounds for integral flows on general directed acyclic graphs and new bounds for the dimensions of weight spaces of parabolic 𝔰𝔩ₙ₊₁(ℂ) Verma modules.

There will be a pre-seminar presenting relevant background at beginning graduate level starting at 1:30pm in MC 5417.

Speaker: Kanstantsin Pashkovich
Affiliation: University of Waterloo
Location: MC 5501

Abstract: Prophet inequalities are a fundamental model for online decision-making under uncertainty. For matroid constraints, Kleinberg and Weinberg gave a tight 1/2-approximation, but their algorithm uses adaptive thresholds that depend on the previously accepted elements. Specifically, the mechanism of Kleinberg and Weinberg accepts an arriving element if and only if it is feasible to add with respect to the matroid constraint and the value of the arriving element passes its threshold; but this threshold depends on the elements accepted so far and on the arriving element itself. Later, Feldman, Svensson, and Zenklusen showed that one can give a 1/4-approximation for general matroids. Their algorithm is almost non-adaptive, i.e., it uses non-adaptive thresholds but changes the underlying matroid to another ``stricter'' matroid. Feldman, Svensson, and Zenklusen also showed that no constant approximation is possible in general matroids using non-adaptive thresholds if one does not change the underlying matroid.

We give a new almost non-adaptive algorithm for matroid prophet inequalities that achieves a 1/2 guarantee. We change the underlying matroid to a ``stricter'' new matroid that is a direct sum of several matroids. For each part of the new matroid, we precompute a single non-adaptive threshold. Once the elements start to arrive, we accept an arriving element as long as it is feasible with respect to the new ``stricter'' matroid and its value passes the precomputed threshold. In addition, we guarantee that our algorithm achieves the 1/2 guarantee not simply with respect to the prophet's expected gain, but with respect to the stronger ex-ante relaxation value.
Thus, we provide the first almost non-adaptive algorithm for the matroid prophet inequality that achieves the best-possible approximation guarantee of 1/2. If time permits, we will also discuss how we used LLM in this project.
Friday, July 24, 2026 11:30 am - 12:30 pm EDT (GMT -04:00)

CombOpt ReadingGroup -Mahtab Alghasi-The Local Dyadic Conjecture

Speaker:

Mahtab Alghasi
Affiliation: University of Waterloo
Location: MC 6029

Abstract: A family of sets $\mathcal{C}$ over a finite ground set $E(\mathcal{C})$ is a clutter if no member of $\mathcal{C}$ properly contains another. A clutter is ideal if its covering polyhedron is integral. A rational number whose denominator is a power of two is called \emph{dyadic}.

A longstanding conjecture of Paul Seymour, known as the Dyadic Conjecture, predicts that, for every ideal clutter, the dual of the set covering linear program admits an optimal solution in which all variables take dyadic values. We present a local version of this conjecture and provide evidence for it by proving the proposed local statement for binary clutters under certain assumptions.
This is joint work with Bertrand Guenin and Levent Tuncel.