A Polynomial Rate of Convergence for the Dirichlet Problem on Orthodiagonal Maps

We extend recent work of Gurel-Gurevich--Jerison--Nachmias (2020) and Bou-Rabee--Gwynne (2024) by showing that as the mesh of our lattice tends to $0$, we have a polynomial rate of convergence for the Dirichlet problem on orthodiagonal maps with H\"older boundary data to its continuous counterpart. The key idea is that the convolution of a discrete harmonic function on an orthodiagonal map with a smooth mollifier has small Laplacian and so is ``almost harmonic." This also allows us to show that discrete harmonic functions on orthodiagonal maps are Lipschitz in the bulk on a mesoscopic scale.