In this paper, we propose a parallel and scalable approach for geodesic distance computation on triangle meshes. Our key observation is that the recovery of with heat method [1] can be reformulated as optimization its gradients subject to integrability, which solved using an efficient first-order requires no linear system solving converges quickly. Afterward, efficiently recovered by integratio...