网格测地线算法（Geodesics in Heat）附源码

　　测地线又称为大地线，可以定义为空间曲面上两点的局部最短路径。测地线具有广泛的应用，例如在工业上测地线最短的性质就意味着最优最省，在航海和航空中，轮船和飞机的运行路线就是测地线。[Crane et al. 2013]提出了利用热运动方程来计算网格测地线的方法，可以想象一下，当一根烫的针尖接触到曲面上的一点时，热量会随着时间的推移而扩散，测地距离因此可以和热运动相联系。具体算法过程如下图所示：

　　第一步：热运动方程用来描述热的传播状态：

　　将热运动方程离散化并整理后得到：

其中：id为单位矩阵，t为时间间隔，Δ为离散Laplacian算子，u_t为t时刻的热状态，u₀为初始时刻的热状态。

　　第二步：第一步计算得到的热梯度方向与测地距离的梯度方向相同，由Eikonal方程知道测地距离的梯度为单位向量，于是通过归一化热梯度我们得到测地距离的梯度：

<