Computing Teichm\"{u}ller Maps between Polygons

TL;DR

This paper introduces the first constructive iterative method to compute extremal quasiconformal (Teichmüller) maps between polygons, ensuring convergence and minimal angle distortion in both continuous and discrete settings.

Contribution

It presents a novel iterative algorithm for computing extremal quasiconformal maps, with proven convergence and efficiency, extending the problem to punctured sphere mappings.

Findings

The method converges quickly to the extremal map within $O(1/\epsilon^4)$ iterations.

Each iteration involves convex optimization and differential equations, reducing dilatation.

The discrete procedure closely approximates the continuous construction, promising practical applicability.

Abstract

By the Riemann-mapping theorem, one can bijectively map the interior of an $n$-gon $P$ to that of another $n$-gon $Q$ conformally. However, (the boundary extension of) this mapping need not necessarily map the vertices of $P$ to those $Q$. In this case, one wants to find the ``best" mapping between these polygons, i.e., one that minimizes the maximum angle distortion (the dilatation) over \textit{all} points in $P$. From complex analysis such maps are known to exist and are unique. They are called extremal quasiconformal maps, or Teichm\"{u}ller maps. Although there are many efficient ways to compute or approximate conformal maps, there is currently no such algorithm for extremal quasiconformal maps. This paper studies the problem of computing extremal quasiconformal maps both in the continuous and discrete settings. We provide the first constructive method to obtain the extremal quasiconformal map in the continuous setting. Our construction is via an iterative procedure that is proven to converge quickly to the unique extremal map. To get to within $\epsilon$ of the dilatation of the extremal map, our method uses $O(1/\epsilon^{4})$ iterations. Every step of the iteration involves convex optimization and solving differential equations, and guarantees a decrease in the dilatation. Our method uses a reduction of the polygon mapping problem to that of the punctured sphere problem, thus solving a more general problem. We also discretize our procedure. We provide evidence for the fact that the discrete procedure closely follows the continuous construction and is therefore expected to converge quickly to a good approximation of the extremal quasiconformal map.