In this article, we present an implementation of a flexible framework for non-parametric variational image registration, realized as part of ITK's finite difference solver hierarchy. In a variational registration setting, the transformation is found by [...]