In our previous paper (Gaster et al., 2018, arXiv:1810.11932), we showed that the theory of harmonic maps between Riemannian manifolds, especially hyperbolic surfaces, may be discretized by introducing a triangulation of the domain manifold with independent vertex and edge weights. In the present paper, we study convergence of the discrete theory back to the smooth theory when taking finer and finer triangulations, in the general Riemannian setting. We present suitable conditions on the weighted triangulations that ensure convergence of discrete harmonic maps to smooth harmonic maps, introducing the notion of (almost) asymptotically Laplacian weights, and we offer a systematic method to construct such weighted triangulations in the two-dimensional case. Our computer software Harmony successfully implements these methods to compute equivariant harmonic maps in the hyperbolic plane.
