The dynamic formulation of optimal transport has attracted growing interests in scientific computing and machine learning, and its computation requires to solve a PDE-constrained optimization problem. The classical Eulerian discretization based approaches suffer from the curse of dimensionality, which arises from the approximation of high-dimensional velocity field. In this work, we propose a deep learning based method to solve the dynamic optimal transport in high dimensional space. Our method contains three main ingredients: a carefully designed representation of the velocity field, the discretization of the PDE constraint along the characteristics, and the computation of high dimensional integral by Monte Carlo method in each time step. Specifically, in the representation of the velocity field, we apply the classical nodal basis function in time and the deep neural networks in space domain with the H1-norm regularization. This technique promotes the regularity of the velocity field in both time and space such that the discretization along the characteristic remains to be stable during the training process. Extensive numerical examples have been conducted to test the proposed method. Compared to other solvers of optimal transport, our method could give more accurate results in high dimensional cases and has very good scalability with respect to dimension. Finally, we extend our method to more complicated cases such as crowd motion problem.

We undertake a systematic study of the infinitesimal geometry of the Thurston metric, showing that the topology, convex geometry and metric geometry of the tangent and cotangent spheres based at any marked hyperbolic surface representing a point in Teichmüller space can recover the marking and geometry of this marked surface. We then translate the results concerning the infinitesimal structures to global geometric statements for the Thurston metric, most notably deriving rigidity statements for the Thurston metric analogous to the celebrated Royden theorem.

We study Thurston’s Lipschitz and curve metrics, as well as the arc metric on the Teichmüller space of one-hold tori equipped with complete hyperbolic metrics with boundary holonomy of fixed length. We construct natural Lipschitz maps between two surfaces equipped with such hyperbolic metrics that generalize Thurston’s stretch maps and prove the following: (1) On the Teichmüller space of the torus with one boundary component, the Lipschitz and the curve metrics coincide and define a geodesic metric on this space. (2) On the same space, the arc and the curve metrics coincide when the length of the boundary component is ≤4arcsinh(1), but differ when the boundary length is large. We further apply our stretch map generalization to construct novel Thurston geodesics on the Teichmüller spaces of closed hyperbolic surfaces, and use these geodesics to show that the sum-symmetrization of the Thurston metric fails to exhibit Gromov hyperbolicity.

The action of the mapping class group of the thrice-punctured projective plane on its GL(2,C) character variety produces an algorithm for generating the simple length spectra of quasi-Fuchsian thrice-punctured projective planes. We apply this algorithm to quasi-Fuchsian representations of the corresponding fundamental group to prove: a sharp upper-bound for the length its shortest geodesic, a McShane identity and the surprising result of non-polynomial growth for the number of simple closed geodesic lengths.

We prove a McShane-type identity: a series, expressed in terms of geodesic lengths, that sums to 2π for any closed hyperbolic surface with one distinguished point. To do so, we prove a generalized Birman-Series theorem showing that the set of complete geodesics on a hyperbolic surface with large cone angles is sparse.