Optimization over the Diffeomorphism Group Using Riemannian BFGS with Application

Applied and Computational Mathematics Seminar
Wednesday, March 4, 2020 - 1:00pm for 1 hour (actually 50 minutes)
Skiles 006
Dr. Darshan Bryner – Naval Surface Warfare Center, Panama City Division
Haomin Zhou

Please Note: This is a part of IEEE Signal Processing Society Lecture Series, organized by Dr. Alessio Medda (alessiomedda@ieee.org). PLEASE RSVP to https://events.vtools.ieee.org/m/222947

The set of diffeomorphisms of the unit interval, or “warping functions,” plays an important role in many in functional data analysis applications. Most prominently, the problem of registering, or aligning, pairs of functions depends on solving for an element of the diffeomorphism group that, when applied to one function, optimally aligns it to the other.
The registration problem is posed as the unconstrained minimization of a cost function over the infinite dimensional diffeomorphism function space. We make use of its well-known Riemannian geometry to implement an efficient, limited memory Riemannian BFGS optimization scheme. We compare performance and results to the benchmark algorithm, Dynamic Programming, on several functional datasets. Additionally, we apply our methodology to the problem of non-parametric density estimation and compare to the benchmark performance of MATLAB’s built-in kernel density estimator ‘ksdensity’.