Variational Geometry-aware Neural Network based Method for Solving High-dimensional Diffeomorphic Mapping Problems