Mathematical foundations
Objectives, Stationary Points, and Local Optima
Learn to formulate optimization problems, classify stationary behavior, and recognize when a surrogate objective or parameterization weakens an optimization claim.
By the end you can
- Identify variables, objectives, constraints, and surrogates in an ML optimization problem
- Distinguish local minima, saddles, flat regions, and boundary optima
- Explain why first-order stationarity is necessary but not sufficient
- Audit whether an optimized objective represents the external decision
Optimization starts by declaring what is allowed to improve
An optimization problem contains variables, an objective, a feasible set, and sometimes several competing criteria. The optimizer can only improve the mathematical problem it receives. A stationary point satisfies a local first-order condition. It may be a minimum, maximum, saddle, or flat region. The objective's global structure determines which conclusions are justified.
Before selecting an algorithm, the modeling team should ask whether the objective represents the real decision and whether its variables include every controllable choice. What that question is worth is on the record. A commercial risk-prediction algorithm affecting millions of patients was trained on health-care costs as a surrogate for health need. The decision it fed was who receives extra help. Obermeyer and three colleagues dissected it in Science in 2019, and correcting the substitution changes who is selected: “Remedying this disparity would increase the percentage of Black patients receiving additional help from 17.7 to 46.5%.” No optimizer in that pipeline was failing. The declared objective was.
A precisely optimized surrogate can still be the wrong problem.
Case
The obstacle in high dimensions turned out to be saddle points
The shape of that global structure is an empirical question, not a definitional one. Dauphin and colleagues argued at NIPS 2014 that “a deeper and more profound difficulty originates from the proliferation of saddle points, not local minima, especially in high dimensional problems of practical interest”. A year later, at AISTATS 2015, Choromanska and colleagues reported that “for large-size decoupled networks the lowest critical values of the random loss function form a layered structure and they are located in a well-defined band lower-bounded by the global minimum”, and that “The number of local minima outside that band diminishes exponentially with the size of the network.”
The structural version arrived in 2016, proved on a family where the algebra closes. For the squared loss of deep linear networks — at any depth and any layer widths, given training data for which XX' and XY' have full rank, d_y <= d_x and the relevant covariance has d_y distinct eigenvalues — Kenji Kawaguchi established that “2) every local minimum is a global minimum, 3) every critical point that is not a global minimum is a saddle point, and 4) there exist "bad" saddle points (where the Hessian has no negative eigenvalue) for the deeper networks (with more than three layers), whereas there is no bad saddle point for the shallow networks (with three layers).”
Read the fourth clause slowly. It is this lesson's thesis in its hardest form. Beyond three layers there are points where the gradient vanishes, where the Hessian has no negative eigenvalue — the usual second-order reassurance — and which are still not minima. The first two clauses are good news, and that news survives independent proof. Laurent and von Brecht gave a short, separate argument in 2018, for deep linear networks with an arbitrary convex differentiable loss: “all local minima are global minima if the hidden layers are either 1) at least as wide as the input layer, or 2) at least as wide as the output layer”. Both results describe linear networks. Neither licenses the same claim about the nonlinear model you are actually training.
Example
When a surrogate objective changes the product behavior
Optum's Impact Pro, the risk-prediction system named in the New York regulators' letter of 25 October 2019, separates optimization from utility line by line.
- Business outcome: Which of millions of patients most need extra care — health need itself, the quantity the decision turns on.
- Available label: Health-care costs already recorded for each patient, dense, cheap and complete.
- Trainable surrogate: A risk prediction fitted to those recorded costs. Obermeyer and colleagues identify that substitution as the heart of the failure.
- Hidden bias: Cost and need are not the same quantity, and the gap is not uniform across patients: “Remedying this disparity would increase the percentage of Black patients receiving additional help from 17.7 to 46.5%.”
- Optimization risk: A better fit to logged cost is not a better fit to need, so ordinary model improvement moves the system further into the error rather than out of it.
- Repair: External evaluation against the real outcome — and, failing that, someone else's. Two New York regulators supplied it on the day the study appeared. Linda A. Lacewell, Superintendent, New York State Department of Financial Services, and Howard A. Zucker, Commissioner, New York State Department of Health, wrote jointly to David S. Wichmann, CEO, UnitedHealth Group Incorporated: “These discriminatory results, whether intentional or not, are unacceptable and are unlawful in New York.” Their demand was binary — “We call on you to immediately investigate these reports and demonstrate that this algorithm is not racially discriminatory or to cease using Impact Pro (or any other data analytics program) if you cannot demonstrate that it does not rely on racial biases or perpetuate racially disparate impacts.”
Comparison
Minimum, maximum, saddle, plateau, and boundary solution
A zero gradient is not a complete classification. Supplying the missing classification is not merely tedious. It is provably intractable, for the simplest nonconvex objective there is. Deciding whether a given point of a nonconvex quadratic program is a local minimum is NP-hard: Murty and Kabadi showed it in 1987, and Pardalos and Schnitger showed the same for constrained quadratic programs in 1988. Whether a local minimizer could nonetheless be found efficiently was left open, and it sat on Pardalos and Vavasis's 1992 list of seven open problems in complexity theory for numerical optimization, phrased this way: “Murty and Kabadi (1987, [11]) and Pardalos and Schnitger (1988, [12]) have shown that it is NP-hard to test whether a given point for such a problem is a local minimizer, but that does not rule out the possibility that another point can be found that is easily verified as a local minimizer.” The open problem stood for three decades. Ahmadi and Zhang closed it in 2022.
The flat column of the comparison is not a metaphor either. Sagun and four colleagues measured the Hessian spectrum of trained networks in 2017, and found it comes in two pieces: “we empirically show that the spectrum of the Hessian is composed of two parts: (1) the bulk centered near zero, (2) and outliers away from the bulk”. Most directions in parameter space carry essentially no curvature, which is what weak identification looks like numerically. The same abstract draws the consequence — “the flatness of such landscapes (which can be measured by the singularity of the Hessian) implies that classical notions of basins of attraction may be quite misleading”. Papyan reports the same geometry independently across deep-learning spectra in 2020: “These include spectral outliers, "spikes", and small but distinct continuous distributions, "bumps", often seen beyond the edge of a "main bulk".”
Local minimum
Nearby feasible points have no lower objective.
- May not be global
- Often supported by positive curvature
- Can lie on a boundary
- Depends on the feasible set
Saddle point
Some directions rise while others fall.
- Gradient may be zero
- Hessian can be indefinite
- Common in high-dimensional nonconvex problems
- Can slow first-order methods
Flat region
Many nearby points have nearly equal objective.
- Gradient can be tiny
- Parameters may be weakly identified
- Noise can dominate movement
- Predictions may remain stable
Boundary optimum
Constraint blocks the unconstrained descent direction.
- Gradient need not vanish
- Requires feasible-direction reasoning
- Often described by multipliers
- Common with simplex or norm limits
Visual
The anatomy of an optimization problem
Keeping these pieces separate prevents algorithm choices from hiding modeling decisions. Impact Pro separates cleanly along exactly these lines: the variables and the objective were about recorded cost, the decision the system fed was about need, and only an evaluation rule aimed at the second quantity could ever have shown that the first was standing in for it.
- 1
Variables
The quantities the procedure may change.
- 2
Objective
The scalar or vector criterion to improve.
- 3
Constraints
The allowed region for variables and outputs.
- 4
Regularizers
Additional preferences or penalties.
- 5
Evaluation rule
Evidence used to judge whether optimization helped the real task.
Stationarity depends on the geometry and the feasible set
For an unconstrained differentiable problem, a local interior optimum must have zero gradient. This is a necessary condition, not a sufficient one. With constraints, only feasible directions matter. A nonzero gradient can point outside the allowed region while no feasible descent direction remains. For nonsmooth objectives, gradients may not exist. Subgradients, proximal operators, or other generalized conditions replace the ordinary stationary equation.
The multiplier conditions that describe optimality at a constraint boundary have a disputed date. Kuhn and Tucker published them in 1951. William Karush had already obtained them in 1939, in an unpublished University of Chicago master's thesis, and never advertised the fact. Richard Cottle documented the priority in 2012; his prologue says of Karush's work: “It tells the story of fundamental optimization results that he obtained in his master’s thesis: results that he neither published nor advertised and that were later independently rediscovered and published by Harold W. Kuhn and Albert W. Tucker.” That is why the conditions carry three initials today rather than two.
Analogy
A flawless tuning process can still choose the wrong car
A race car optimized on a short, dry circuit meets a long, wet mountain route in the real competition. The engineering process can be flawless and still select the wrong design. The objective defines the test track. Constraints define which modifications are legal, and the optimizer searches within those rules.
A track keeps its shape no matter how many laps are run on it, while an ML objective also depends on data and feedback generated by previous system behavior. No amount of algorithmic skill rescues a poor objective.
Optimization quality is conditional on problem formulation.
Steps
Write an optimization problem before choosing an optimizer
A clear formulation prevents algorithm tuning from masking a misspecified objective. Step 5 is where a cost column standing in for health need gets written down as a surrogate rather than passed off as the target. Step 6 is where someone measures the system against need instead of against cost. Do that yourself, or a state Superintendent and a state Commissioner will do it for you.
1. Name the decision
State which real action the model will influence.
2. Declare variables
Separate learned parameters, thresholds, and operational controls.
3. Write the objective
Connect each term to a desired behavior or cost.
4. State constraints
Include feasibility, fairness, safety, and resource limits.
5. Identify surrogates
Mark where labels or differentiable losses stand in for real outcomes.
6. Define external evaluation
Specify evidence that can reveal objective misspecification.
Key idea
Stationary points can change appearance under parameterization
Two parameterizations can represent the same predictions with different curvature, gradient scale, and numerical behavior. Deep networks often contain symmetries that create whole families of equivalent parameter values.
In 2017 that became constructive rather than cautionary. The symmetries of rectifier networks let you build models that are exactly equivalent as functions but sit in arbitrarily sharper minima, and Dinh and three colleagues built them. Their abstract states it plainly — “Furthermore, if we allow to reparametrize a function, the geometry of its parameters can change drastically without affecting its generalization properties.” The objective value may remain unchanged while the Hessian spectrum changes under rescaling. Statements about “flatness” therefore require a declared parameterization and neighborhood.
This is not a solved detail that the field has since routed around. Four years later, at NeurIPS 2021, Petzka and four colleagues still recorded that “it is still an open theoretical problem why and under which circumstances flatness is connected to generalization, in particular in light of reparameterizations that change certain flatness measures but leave generalization unchanged”. Optimization geometry belongs to both the function and the coordinates used to express it.
Parameter-space shape is not automatically an invariant property of the represented model.
Key takeaways
- An optimizer improves the declared mathematical objective, not an unstated product goal. Impact Pro was fitted to recorded health-care costs while the decision was who needs care, and repairing the surrogate moves the share of Black patients receiving additional help from 17.7 to 46.5%.
- A zero gradient can describe a minimum, maximum, saddle, or flat region, and telling which is NP-hard for nonconvex quadratic programs (Murty and Kabadi, 1987; Pardalos and Schnitger, 1988). Kawaguchi's deep linear networks even contain "bad" saddle points whose Hessian has no negative eigenvalue.
- Boundary optimality depends on feasible directions rather than an unconstrained zero gradient. The multiplier conditions that express it are Kuhn and Tucker's, published in 1951 and anticipated by Karush's unpublished 1939 thesis.
- Surrogate losses require external evaluation because they can be improved for the wrong reasons — and if you do not supply it, a Superintendent and a Commissioner may, as they did on 25 October 2019.
- Curvature and flatness depend on parameterization and scale. Dinh and three colleagues built exactly equivalent functions sitting in arbitrarily sharper minima in 2017, and measured spectra show a bulk of curvature directions centred near zero (Sagun and colleagues).
- Problem formulation should be reviewed before computational effort is spent on optimizer tuning, because no tuning budget converts a well-fitted cost model into a model of need.