Research
Theory and algorithms for optimization in function spaces
My work combines functional and variational analysis, PDE theory and probability with the design and analysis of numerical methods. Four threads run through it.
Optimization under uncertainty and risk
Design and control problems for physical systems are typically posed with uncertain data. Optimizing an expected value is often not enough: one needs to control the tail of the distribution of outcomes. I work on the analytical foundations of risk-averse PDE-constrained optimization, including existence, optimality conditions and stability with respect to the underlying probability law. I also develop algorithms that remain efficient in infinite dimensions.
Examples include the epi-regularization of risk measures, primal–dual and interior-point methods for risk minimization, and the asymptotic analysis of Monte Carlo (sample average) approximations, much of it in long-standing collaboration with Drew P. Kouri (Sandia National Laboratories).
Selected papers
- Risk-averse PDE-constrained optimization using the conditional value-at-risk SIAM Journal on Optimization, Vol. 26, No. 1, pp. 365–396 (2016) DOI
- Existence and optimality conditions for risk-averse PDE-constrained optimization SIAM/ASA Journal on Uncertainty Quantification, Vol. 6, No. 2, pp. 787–815 (2018) DOI
- Epi-regularization of risk measures Mathematics of Operations Research, Vol. 45, No. 2, pp. 774–795 (2020) DOI
- A primal–dual algorithm for risk minimization Mathematical Programming, Vol. 193, No. 1, pp. 337–363 (2022) DOI
- Asymptotic consistency for nonconvex risk-averse stochastic optimization with infinite-dimensional decision spaces Mathematics of Operations Research, Vol. 49, No. 3, pp. 1403–1418 (2024) DOIarXiv
- Asymptotic properties of Monte Carlo methods in elliptic PDE-constrained optimization under uncertainty Numerische Mathematik, Vol. 156, No. 5, pp. 1887–1914 (2024) DOIarXiv
- Risk-averse optimal control of random elliptic variational inequalities ESAIM: Control, Optimisation and Calculus of Variations, Vol. 31, Art. 71 (2025) DOIarXiv
Structure-preserving numerical methods
Pointwise bounds, obstacles and other inequality constraints appear throughout mechanics, physics and design. Proximal Galerkin and the latent variable proximal point (LVPP) method, developed with Brendan Keith and coauthors, enforce such constraints exactly at the discrete level by working with a latent variable inside a Bregman proximal point iteration. The result is a family of finite element solvers whose iteration counts are largely insensitive to mesh refinement.
The same ideas underlie SiMPL, a simple first-order method for density-based topology optimization, and recent extensions to phase-field fracture and multi-material design. Proximal Galerkin received the 2025 Langtangen Prize. The methods are implemented in FEniCS, Firedrake and MFEM.
Selected papers
- Proximal Galerkin: A structure-preserving finite element method for pointwise bound constraints Foundations of Computational Mathematics, Vol. 26, No. 1, pp. 385–481 (2026) Langtangen Prize 2025DOIarXivCode
- The latent variable proximal point algorithm for variational problems with inequality constraints Computer Methods in Applied Mechanics and Engineering, Vol. 445, Art. 118181 (2025) DOIarXivCode
- Analysis of the SiMPL method for density-based topology optimization SIAM Journal on Optimization, Vol. 35, No. 2, pp. 1134–1164 (2025) DOIarXiv
- A simple introduction to the SiMPL method for density-based topology optimization Structural and Multidisciplinary Optimization, Vol. 68, No. 4, Art. 74 (2025) DOIarXiv
- Proximal Galerkin for phase field fracture PreprintarXiv
- The SiMPL method for multi-material topology optimization PreprintarXiv
Variational analysis and optimal control
Many control and design problems have a nonsmooth structure: the state is defined by a variational inequality, a complementarity system or the solution of a lower-level optimization problem. I use tools from variational analysis, such as generalized derivatives, coderivatives and directional differentiability, to derive stationarity conditions for such problems in function spaces and to design solution algorithms.
A related strand concerns generalized Nash equilibrium problems with PDE constraints, including risk-neutral and mean-field formulations.
Selected papers
- First-order optimality conditions for elliptic mathematical programs with equilibrium constraints via variational analysis SIAM Journal on Optimization, Vol. 21, No. 4, pp. 1561–1593 (2011) DOI
- Several approaches for the derivation of stationarity conditions for elliptic MPECs with upper-level control constraints Mathematical Programming, Vol. 146, No. 1–2, pp. 555–582 (2014) DOI
- Generalized Nash equilibrium problems in Banach spaces: Theory, Nikaido–Isoda-based path-following methods, and applications SIAM Journal on Optimization, Vol. 25, No. 3, pp. 1826–1856 (2015) DOI
- A bundle-free implicit programming approach for a class of elliptic MPECs in function space Mathematical Programming, Vol. 160, No. 1–2, pp. 271–305 (2016) DOI
- On the directional differentiability of the solution mapping for a class of variational inequalities of the second kind Set-Valued and Variational Analysis, Vol. 26, No. 3, pp. 631–642 (2018) DOI
- Optimal control of the Landau–de Gennes model of nematic liquid crystals SIAM Journal on Control and Optimization, Vol. 61, No. 4, pp. 2546–2570 (2023) DOIarXiv
- On a differential generalized Nash equilibrium problem with mean field interaction SIAM Journal on Optimization, Vol. 34, No. 3, pp. 2005–2030 (2024) DOI
Optimization for trustworthy AI
At SURE-AI, the Norwegian Centre for Sustainable, Risk-Averse and Ethical AI, I lead the work package on optimization in AI. The aim is to bring the mathematics of risk-averse optimization to bear on learning systems: risk-aware and explainable inference, the stability of statistical estimators, and energy-efficient algorithms with guarantees.
Selected papers
- Inference as risk identification Computational Management Science, Vol. 23, No. 2, Art. 34 (2026) DOIarXiv
- Exponential map models as an interpretable framework for generating neural spatial representations Advances in Neural Information Processing Systems (NeurIPS 2026), to appear (2026) ConferenceDOI Link is to the bioRxiv preprint