Computing discrete equivariant harmonic maps

We present effective methods to compute equivariant harmonic maps from the universal cover of a surface into a nonpositively curved space. By discretizing the theory appropriately, we show that the energy functional is strongly convex and derive convergence of the discrete heat flow to the energy minimizer, with explicit convergence rate. We also examine center of mass methods, after showing a generalized mean value property for harmonic maps. We feature a concrete illustration of these methods with Harmony, a computer software that we developed in C++, whose main functionality is to numerically compute and display equivariant harmonic maps.