Abstract: We present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In a nutshell, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers and is suitable for PDE-constrained optimization and other large-scale applications.
We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary, i.e., in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary.
In the context of optimal control problems, e.g., in PDE-constrained optimization, ALESQP provides a unified framework to efficiently handle general constraints on both the state variables and the control variables. A key algorithmic feature is a constraint decomposition strategy that allows ALESQP to exploit problem-specific variable scalings and inner products. We present several examples with state and control inequality constraints where ALESQP shows remarkable mesh-independent performance, requiring only a handful of outer (AL) iterations to meet constraint tolerances at the level of machine precision. At the same time, ALESQP uses the inner (SQP) loop economically, requiring only a few dozen SQP iterations in total.
Abstract: Multiscale substitution tilings are a new family of tilings of Euclidean space that are generated by multiscale substitution rules. Unlike the standard setup of substitution tilings, which is a basic object of study within the aperiodic order community and includes examples such as the Penrose and the pinwheel tilings, multiple distinct scaling constants are allowed, and the defining process of inflation and subdivision is a continuous one. Under a certain irrationality assumption on the scaling constants, this construction gives rise to a new class of tilings, tiling spaces and tiling dynamical systems, which are intrinsically different from those that arise in the standard setup. In the talk I will describe these new objects and discuss various structural, geometrical, statistical and dynamical results. Based on joint work with Yaar Solomon.