Georgia Institute of Technology College of Sciences

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’.