Georgia Institute of Technology College of Sciences

Accelerated gradient methods play a central role in optimization, achieving the optimal convergence rates in many settings. While many extensions of Nesterov's original acceleration method have been proposed, it is not yet clear what is the natural scope of the acceleration concept. In this work, we study accelerated methods from a continuous-time perspective. We show there is a Bregman Lagrangian functional that generates a large class of accelerated methods in continuous time, including (but not limited to) accelerated gradient descent, its non-Euclidean extension, and accelerated higher-order gradient methods. We show that in continuous time, these accelerated methods correspond to traveling the same curve in spacetime at different speeds. This is in contrast to the family of rescaled gradient flows, which correspond to changing the distance in space. We show how to implement both the rescaled and accelerated gradient methods as algorithms in discrete time with matching convergence rates. These algorithms achieve faster convergence rates for convex optimization under higher-order smoothness assumptions. We will also discuss lower bounds and some open questions. Joint work with Ashia Wilson and Michael Jordan.

The Virtual Element Method (VEM), is a very recent technology introduced in [Beirao da Veiga, Brezzi, Cangiani, Manzini, Marini, Russo, 2013, M3AS] for the discretization of partial differential equations, that has shared a good success in recent years. The VEM can be interpreted as a generalization of the Finite Element Method that allows to use general polygonal and polyhedral meshes, still keeping the same coding complexity and allowing for arbitrary degree of accuracy. The Virtual Element Method makes use of local functions that are not necessarily polynomials and are defined in an implicit way. Nevertheless, by a wise choice of the degrees of freedom and introducing a novel construction of the associated stiffness matrixes, the VEM avoids the explicit integration of such shape functions. In addition to the possibility to handle general polytopal meshes, the flexibility of the above construction yields other interesting properties with respect to more standard Galerkin methods. For instance, the VEM easily allows to build discrete spaces of arbitrary C^k regularity, or to satisfy exactly the divergence-free constraint for incompressible fluids. The present talk is an introduction to the VEM, aiming at showing the main ideas of the method. We consider for simplicity a simple elliptic model problem (that is pure diffusion with variable coefficients) but set ourselves in the more involved 3D setting. In the first part we introduce the adopted Virtual Element space and the associated degrees of freedom, first by addressing the faces of the polyhedrons (i.e. polygons) and then building the space in the full volumes. We then describe the construction of the discrete bilinear form and the ensuing discretization of the problem. Furthermore, we show a set of theoretical and numerical results. In the very final part, we will give a glance at more involved problems, such as magnetostatics (mixed problem with more complex spaces interacting) and large deformation elasticity (nonlinear problem).

We assume the stock is modeled by a Markov regime-switching diffusion process and that, the benchmark depends on the economic factor. Then, we solve a risk-sensitive benchmarked asset management problem of a firm. Our method consists of finding the portfolio strategy that minimizes the risk sensitivity of an investor in such environment, using the general maximum principle.After the above presentation, the speaker will discuss some of her ongoing research.

Large-scale geophysical flows, i.e. the ocean and atmosphere, evolve on spatial scales ranging from meters to thousands of kilometers, and on temporal scales ranging from seconds to decades. These scales interact in a highly nonlinear fashion, making it extremely challenging to reliably and accurately capture the long-term dynamics of these flows on numerical models. In fact, this problem is closely associated with the grand challenges of long-term weather and climate predictions. Unstructured meshes have been gaining popularity in recent years on geophysical models, thanks to its being almost free of polar singularities, and remaining highly scalable even at eddy resolving resolutions. However, to unleash the full potential of these meshes, new schemes are needed. This talk starts with a brief introduction to large-scale geophysical flows. Then it goes over the main considerations, i.e. various numerical and algorithmic choices, that one needs to make in deisgning numerical schemes for these flows. Finally, a new vorticity-divergence based finite volume scheme will be introduced. Its strength and challenges, together with some numerical results, will be presented and discussed.